Rare genetic diseases are typically studied in referral populations, resulting in underdiagnosis and biased assessment of penetrance and phenotype. To address this, we develop a generalizable method of genotype inference based on distant relatedness and deploy this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population. We identify 9 LQT5 families referred to a single specialty clinic, each carrying p.Asp76Asn, the most common LQT5 variant. We uncover recent common ancestry and a single shared haplotype among probands. Application to a non-referral population of 69,819 BioVU biobank subjects identifies 22 additional subjects sharing this haplotype, which we confirm to carry p.Asp76Asn. Referral and non-referral carriers have prolonged QT interval corrected for heart rate (QTc) compared to controls, and, among carriers, the QTc polygenic score is independently associated with QTc prolongation. Thus, our innovative analysis of shared chromosomal segments identifies undiagnosed cases of genetic disease and refines the understanding of LQT5 penetrance and phenotype.

Introduction: Transcriptome-wide Association Studies (TWAS) extend genome-wide association studies (GWAS) by integrating genetically-regulated gene expression models. We performed the most powerful AD-TWAS to date, using summary statistics from cis -eQTL meta-analyses and the largest clinically-adjudicated Alzheimer's Disease (AD) GWAS.

Methods: We implemented the OTTERS TWAS pipeline, leveraging cis -eQTL data from cortical brain tissue (MetaBrain; N=2,683) and blood (eQTLGen; N=31,684) to predict gene expression, then applied these models to AD-GWAS data (Cases=21,982; Controls=44,944).

Results: We identified and validated five novel gene associations in cortical brain tissue ( PRKAG1 , C3orf62 , LYSMD4 , ZNF439 , SLC11A2 ) and six genes proximal to known AD-related GWAS loci (Blood: MYBPC3 ; Brain: MTCH2 , CYB561 , MADD , PSMA5 , ANXA11 ). Further, using causal eQTL fine-mapping, we generated sparse models that retained the strength of the AD-TWAS association for MTCH2 , MADD , ZNF439 , CYB561 , and MYBPC3 .

Discussion: Our comprehensive AD-TWAS discovered new gene associations and provided insights into the functional relevance of previously associated variants.

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Background: The X chromosome is often omitted in disease association studies despite containing thousands of genes that may provide insight into well-known sex differences in the risk of Alzheimer's disease (AD).

Objective: To model the expression of X chromosome genes and evaluate their impact on AD risk in a sex-stratified manner.

Methods: Using elastic net, we evaluated multiple modeling strategies in a set of 175 whole blood samples and 126 brain cortex samples, with whole genome sequencing and RNA-seq data. SNPs (MAF > 0.05) within the cis-regulatory window were used to train tissue-specific models of each gene. We apply the best models in both tissues to sex-stratified summary statistics from a meta-analysis of Alzheimer's Disease Genetics Consortium (ADGC) studies to identify AD-related genes on the X chromosome.

Results: Across different model parameters, sample sex, and tissue types, we modeled the expression of 217 genes (95 genes in blood and 135 genes in brain cortex). The average model R2 was 0.12 (range from 0.03 to 0.34). We also compared sex-stratified and sex-combined models on the X chromosome. We further investigated genes that escaped X chromosome inactivation (XCI) to determine if their genetic regulation patterns were distinct. We found ten genes associated with AD at p < 0.05, with only ARMCX6 in female brain cortex (p = 0.008) nearing the significance threshold after adjusting for multiple testing (α = 0.002).

Conclusions: We optimized the expression prediction of X chromosome genes, applied these models to sex-stratified AD GWAS summary statistics, and identified one putative AD risk gene, ARMCX6.

Keywords: Alzheimer’s disease; X chromosome; bioinformatics; elastic net regression; gene expression; gene prediction; sex differences; transcriptome.

Background: Preeclampsia is a hypertensive disorder of pregnancy characterized by widespread vascular inflammation. It occurs frequently in pregnancy, often without known risk factors, and has high rates of maternal and fetal morbidity and mortality. Identification of biomarkers that predict preeclampsia and its cardiovascular sequelae before clinical onset, or even before pregnancy, is a critical unmet need for the prevention of adverse pregnancy outcomes.

Methods: We explored differences in cardiovascular proteomics (Olink Explore 384) in 256 diverse pregnant persons across 2 centers (26% Hispanic, 21% Black).

Results: We identified significant differences in plasma abundance of markers associated with angiogenesis, blood pressure, cell adhesion, inflammation, and metabolism between individuals delivering with preeclampsia and controls, some of which have not been widely described previously and are not represented in the preeclampsia placental transcriptome. While we observed a broadly similar pattern in early (<34 weeks) versus late (≥34 weeks) preeclampsia, several proteins related to hemodynamic stress, hemostasis, and immune response appeared to be more highly dysregulated in early preeclampsia relative to late preeclampsia.

Conclusions: These results demonstrate the value of performing targeted proteomics using a panel of cardiovascular biomarkers to identify biomarkers relevant to preeclampsia pathophysiology and highlight the need for larger multiomic studies to define modifiable pathways of surveillance and intervention upstream to preeclampsia diagnosis.

Keywords: blood pressure; multiomic; placenta; proteomics; risk factors.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Metabolic dysfunction-associated steatotic liver disease (MASLD) prevalence is increasing in parallel with an obesity pandemic, calling for novel strategies for prevention and treatment. We defined a circulating proteome of human MASLD across ≈7000 proteins in ≈5000 individuals from diverse, at-risk populations across the metabolic health spectrum, demonstrating reproducible diagnostic performance and specifying both known and novel metabolic pathways relevant to MASLD (central carbon and amino acid metabolism, hepatocyte regeneration, inflammation, fibrosis, insulin sensitivity). A parsimonious proteomic signature of MASLD was associated with a protection from MASLD and its related multi-system metabolic consequences in >26000 free-living individuals, with an additive effect to polygenic risk. The MASLD proteome was encoded by genes that demonstrated transcriptional enrichment in liver, with spatial transcriptional activity in areas of steatosis in human liver biopsy and dynamicity for select targets in human liver across stages of steatosis. We replicated several top relations from proteomics and spatial tissue transcriptomics in a humanized "liver-on-a-chip" model of MASLD, highlighting the power of a full translational approach to discovery in MASLD. Collectively, these results underscore utility of blood-based proteomics as a dynamic "liquid biopsy" of human liver relevant to clinical biomarker and mechanistic applications.

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet β cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and β cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by β cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the β cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by β cells. RFX6 perturbation in primary human islet cells alters β cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

Rare genetic diseases are typically studied in referral populations, resulting in underdiagnosis and biased assessment of penetrance and phenotype. To address this, we developed a generalizable method of genotype inference based on distant relatedness and deployed this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population. We identified 9 LQT5 families referred to a single specialty clinic, each carrying p.Asp76Asn, the most common LQT5 variant. We uncovered recent common ancestry and a single shared haplotype among probands. Application to a non-referral population of 69,819 BioVU biobank subjects identified 22 additional subjects sharing this haplotype, subsequently confirmed to carry p.Asp76Asn. Referral and non-referral carriers had prolonged QTc compared to controls, and, among carriers, QTc polygenic score additively associated with QTc prolongation. Thus, our novel analysis of shared chromosomal segments identified undiagnosed cases of genetic disease and refined the understanding of LQT5 penetrance and phenotype.

Background: Adipokines are hormones secreted from adipose tissue and are associated with cardiometabolic diseases (CMD). Functional differences between adipokines (leptin, adiponectin, and resistin) are known, but inconsistently reported associations with CMD and lack of studies in Hispanic populations are research gaps. We investigated the relationship between subclinical atherosclerosis and multiple adipokine measures.

Methods: Cross-sectional data from the Cameron County Hispanic Cohort (N = 624; mean age = 50; Female = 70.8%) were utilized to assess associations between adipokines [continuous measures of adiponectin, leptin, resistin, leptin-to-adiponectin ratio (LAR), and adiponectin-resistin index (ARI)] and early atherosclerosis [carotid-intima media thickness (cIMT)]. We adjusted for sex, age, body mass index (BMI), smoking status, cytokines, fasting blood glucose levels, blood pressure, lipid levels, and medication usage in the fully adjusted linear regression model. We conducted sexes-combined and sex-stratified analyses to account for sex-specificity and additionally tested whether stratification of participants by their metabolic status (metabolically elevated risk for CMD as defined by having two or more of the following conditions: hypertension, dyslipidemia, insulin resistance, and inflammation vs. not) influenced the relationship between adipokines and cIMT.

Results: In the fully adjusted analyses, adiponectin, leptin, and LAR displayed significant interaction by sex (p < 0.1). Male-specific associations were between cIMT and LAR [β(SE) = 0.060 (0.016), p = 2.52 × 10-4], and female-specific associations were between cIMT and adiponectin [β(SE) = 0.010 (0.005), p = 0.043] and ARI [β(SE) = - 0.011 (0.005), p = 0.036]. When stratified by metabolic health status, the male-specific positive association between LAR and cIMT was more evident among the metabolically healthy group [β(SE) = 0.127 (0.015), p = 4.70 × 10-10] (p for interaction by metabolic health < 0.1). However, the female-specific associations between adiponectin and cIMT and ARI and cIMT were observed only among the metabolically elevated risk group [β(SE) = 0.014 (0.005), p = 0.012 for adiponectin; β(SE) = - 0.015 (0.006), p = 0.013 for ARI; p for interaction by metabolic health < 0.1].

Conclusion: Associations between adipokines and cIMT were sex-specific, and metabolic health status influenced the relationships between adipokines and cIMT. These heterogeneities by sex and metabolic health affirm the complex relationships between adipokines and atherosclerosis.

Keywords: Adipokines; Cardiometabolic health; Carotid intima-media thickness; Health disparities; Mexican–American; Population-based cohort; Sex-differences.

Minority populations are largely absent from clinical research trials. The neglect of these populations has become increasingly apparent, with escalating cancer burdens and chronic disease. The challenges to recruitment of minorities in the United States are multiple including trust or lack thereof. Keys to successful recruitment are responding to community issues, its history, beliefs, and its social and economic pressures. The strategy we have used in many low-income, sometimes remote, communities is to recruit staff from the same community and train them in the required basic research methods. They are the first line of communication. After our arrival in the Texas Rio Grande Valley in 2001, we applied these principles learned over years of global research, to studies of chronic diseases. Beginning in 2004, we recruited and trained a team of local women who enrolled in a cohort of over five thousand Mexican Americans from randomly selected households. This cohort is being followed, and the team has remained, acquiring not only advanced skills (ultrasound, FibroScan, retinal photos, measures of cognition, etc.) but capacity to derive key health information. Currently, we are participating in multiple funded studies, including an NIH clinical trial, liver disease, obesity, and diabetes using multiomics aimed at developing precision medicine approaches to chronic disease prevention and treatment.

Keywords: Chronic disease; Recruitment; community engagement; health disparities; minority clinical studies.

Introduction: Apical periodontitis (AP) is a common consequence of root canal infection leading to periapical bone resorption. Microbial and host genetic factors and their interactions have been shown to play a role in AP development and progression. Variations in a few genes have been reported in association with AP; however, the lack of genome-wide studies has hindered progress in understanding the molecular mechanisms involved. Here, we report the first genome-wide association study of AP in a large and well-characterized population.

Methods: Male and female adults (n = 932) presenting with deep caries and AP (cases), or deep caries without AP (controls) were included. Genotyping was performed using the Illumina Expanded Multi-Ethnic Genotyping Array (MEGA). Single-variant association testing was performed adjusting for sex and 5 principal components. Subphenotype association testing, analyses of genetically regulated gene expression, polygenic risk score, and phenome-wide association (PheWAS) analyses were also conducted.

Results: Eight loci reached near genome-wide significant association with AP (P < 5 × 10-6); gene-focused analyses replicated 3 previously reported associations (P < 8.9 × 10-5). Sex-specific and subphenotype-specific analyses revealed additional significant associations with variants genome-wide. Functionally oriented gene-based analyses revealed 8 genes significantly associated with AP (P < 5 × 10-5), and PheWAS analysis revealed 33 phecodes associated with AP risk score (P < 3.08 × 10-5).

Conclusions: This study identified novel genes/loci contributing to AP and specific contributions to AP risk in men and women. Importantly, we identified additional systemic conditions significantly associated with AP risk. Our findings provide strong evidence for host-mediated effects on AP susceptibility.

Keywords: Gene association; periapical lesion; pulp necrosis; susceptibility.

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

Background: The X chromosome is often omitted in disease association studies despite containing thousands of genes which may provide insight into well-known sex differences in the risk of Alzheimer's Disease.

Objective: To model the expression of X chromosome genes and evaluate their impact on Alzheimer's Disease risk in a sex-stratified manner.

Methods: Using elastic net, we evaluated multiple modeling strategies in a set of 175 whole blood samples and 126 brain cortex samples, with whole genome sequencing and RNA-seq data. SNPs (MAF>0.05) within the cis-regulatory window were used to train tissue-specific models of each gene. We apply the best models in both tissues to sex-stratified summary statistics from a meta-analysis of Alzheimer's disease Genetics Consortium (ADGC) studies to identify AD-related genes on the X chromosome.

Results: Across different model parameters, sample sex, and tissue types, we modeled the expression of 217 genes (95 genes in blood and 135 genes in brain cortex). The average model R2 was 0.12 (range from 0.03 to 0.34). We also compared sex-stratified and sex-combined models on the X chromosome. We further investigated genes that escaped X chromosome inactivation (XCI) to determine if their genetic regulation patterns were distinct. We found ten genes associated with AD at p 0.05, with only ARMCX6 in female brain cortex (p = 0.008) nearing the significance threshold after adjusting for multiple testing (α = 0.002).

Conclusions: We optimized the expression prediction of X chromosome genes, applied these models to sex-stratified AD GWAS summary statistics, and identified one putative AD risk gene, ARMCX6.

Keywords: Alzheimer’s disease; Bioinformatics; Elastic Net Regression; Gene Expression; Gene Prediction; Sex Differences; Transcriptome; X Chromosome.

Importance: The diagnosis and study of rare genetic disease is often limited to referral populations, leading to underdiagnosis and a biased assessment of penetrance and phenotype.

Objective: To develop a generalizable method of genotype inference based on distant relatedness and to deploy this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population.

Participants: We identified 9 LQT5 probands and 3 first-degree relatives referred to a single Genetic Arrhythmia clinic, each carrying D76N (p.Asp76Asn), the most common variant implicated in LQT5. The non-referral population consisted of 69,879 ancestry-matched subjects in BioVU, a large biobank that links electronic health records to dense array data. Participants were enrolled from 2007-2022. Data analysis was performed in 2022.

Exposures: We developed and applied a novel approach to genotype inference (Distant Relatedness for Identification and Variant Evaluation, or DRIVE) to identify shared, identical-by-descent (IBD) large chromosomal segments in array data.

Main outcomes and measures: We sought to establish genetic relatedness among the probands and to use genomic segments underlying D76N to identify other potential carriers in BioVU. We then further studied the role of D76N in LQT5 pathogenesis.

Results: Genetic reconstruction of pedigrees and distant relatedness detection among clinic probands using DRIVE revealed shared recent common ancestry and identified a single long shared haplotype. Interrogation of the non-referral population in BioVU identified a further 23 subjects sharing this haplotype, and sequencing confirmed D76N carrier status in 22, all previously undiagnosed with LQT5. The QTc was prolonged in D76N carriers compared to BioVU controls, with 40% penetrance of QTc ≥ 480 msec. Among D76N carriers, a QTc polygenic score was additively associated with QTc prolongation.

Conclusions and relevance: Detection of IBD shared chromosomal segments around D76N enabled identification of distantly related and previously undiagnosed rare-variant carriers, demonstrated the contribution of polygenic risk to monogenic disease penetrance, and further established LQT5 as a primary arrhythmia disorder. Analysis of shared chromosomal regions spanning disease-causing mutations can identify undiagnosed cases of genetic diseases.

Aims/hypothesis: The Latino population has been systematically underrepresented in large-scale genetic analyses, and previous studies have relied on the imputation of ungenotyped variants based on the 1000 Genomes (1000G) imputation panel, which results in suboptimal capture of low-frequency or Latino-enriched variants. The National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) released the largest multi-ancestry genotype reference panel representing a unique opportunity to analyse rare genetic variations in the Latino population. We hypothesise that a more comprehensive analysis of low/rare variation using the TOPMed panel would improve our knowledge of the genetics of type 2 diabetes in the Latino population.

Methods: We evaluated the TOPMed imputation performance using genotyping array and whole-exome sequence data in six Latino cohorts. To evaluate the ability of TOPMed imputation to increase the number of identified loci, we performed a Latino type 2 diabetes genome-wide association study (GWAS) meta-analysis in 8150 individuals with type 2 diabetes and 10,735 control individuals and replicated the results in six additional cohorts including whole-genome sequence data from the All of Us cohort.

Results: Compared with imputation with 1000G, the TOPMed panel improved the identification of rare and low-frequency variants. We identified 26 genome-wide significant signals including a novel variant (minor allele frequency 1.7%; OR 1.37, p=3.4 × 10-9). A Latino-tailored polygenic score constructed from our data and GWAS data from East Asian and European populations improved the prediction accuracy in a Latino target dataset, explaining up to 7.6% of the type 2 diabetes risk variance.

Conclusions/interpretation: Our results demonstrate the utility of TOPMed imputation for identifying low-frequency variants in understudied populations, leading to the discovery of novel disease associations and the improvement of polygenic scores.

Data availability: Full summary statistics are available through the Common Metabolic Diseases Knowledge Portal ( https://t2d.hugeamp.org/downloads.html ) and through the GWAS catalog ( https://www.ebi.ac.uk/gwas/ , accession ID: GCST90255648). Polygenic score (PS) weights for each ancestry are available via the PGS catalog ( https://www.pgscatalog.org , publication ID: PGP000445, scores IDs: PGS003443, PGS003444 and PGS003445).

Keywords: GWAS meta-analysis; Latino population; Polygenic score; TOPMed imputation; Type 2 diabetes.

Background: Haptoglobin (HP) is an antioxidant of apolipoprotein E (APOE), and previous reports have shown HP binds with APOE and amyloid beta (Aβ) to aid its clearance. A common structural variant of the HP gene distinguishes it into two alleles: HP1 and HP2.

Methods: HP genotypes were imputed in 29 cohorts from the Alzheimer's Disease Genetics Consortium (N = 20,512). Associations between the HP polymorphism and Alzheimer's disease (AD) risk and age of onset through APOE interactions were investigated using regression models.

Results: The HP polymorphism significantly impacts AD risk in European-descent individuals (and in meta-analysis with African-descent individuals) by modifying both the protective effect of APOE ε2 and the detrimental effect of APOE ε4. The effect is particularly significant among APOE ε4 carriers.

Discussion: The effect modification of APOE by HP suggests adjustment and/or stratification by HP genotype is warranted when APOE risk is considered. Our findings also provided directions for further investigations on potential mechanisms behind this association.

Keywords: Alzheimer's disease; age at onset; apolipoprotein E; haptoglobin polymorphism.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Uncovering the genetic underpinnings of musical ability and engagement is a foundational step for exploring their wide-ranging associations with cognition, health, and neurodevelopment. Prior studies have focused on using twin and family designs, demonstrating moderate heritability of musical phenotypes. The current study used genome-wide complex trait analysis and polygenic score (PGS) approaches utilizing genotype data to examine genetic influences on two musicality traits (rhythmic perception and music engagement) in N = 1792 unrelated adults in the Vanderbilt Online Musicality Study. Meta-analyzed heritability estimates (including a replication sample of Swedish individuals) were 31% for rhythmic perception and 12% for self-reported music engagement. A PGS derived from a recent study on beat synchronization ability predicted both rhythmic perception (β = 0.11) and music engagement (β = 0.19) in our sample, suggesting that genetic influences underlying self-reported beat synchronization ability also influence individuals' rhythmic discrimination aptitude and the degree to which they engage in music. Cross-trait analyses revealed a modest contribution of PGSs from several nonmusical traits (from the cognitive, personality, and circadian chronotype domains) to individual differences in musicality (β = -0.06 to 0.07). This work sheds light on the complex relationship between the genetic architecture of musical rhythm processing, beat synchronization, music engagement, and other nonmusical traits.

Purpose: Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk.

Methods: To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results.

Results: GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022.

Conclusion: Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

Trial registration: ClinicalTrials.gov NCT05277116.

Keywords: Common variants; Family history; Genotyping; Monogenic risks; Polygenic risk scores.

Importance: Developmental language disorder (DLD) is a common (with up to 7% prevalence) yet underdiagnosed childhood disorder whose underlying biological profile and comorbidities are not fully understood, especially at the population level.

Objective: To identify clinically relevant conditions that co-occur with DLD at the population level.

Design, setting, and participants: This case-control study used an electronic health record (EHR)-based population-level approach to compare the prevalence of comorbid health phenotypes between DLD cases and matched controls. These cases were identified using the Automated Phenotyping Tool for Identifying Developmental Language Disorder algorithm of the Vanderbilt University Medical Center EHR, and a phenome enrichment analysis was used to identify comorbidities. An independent sample was selected from the Geisinger Health System EHR to test the replication of the phenome enrichment using the same phenotyping and analysis pipeline. Data from the Vanderbilt EHR were accessed between March 2019 and October 2020, while data from the Geisinger EHR were accessed between January and March 2022.

Main outcomes and measures: Common and rare comorbidities of DLD at the population level were identified using EHRs and a phecode-based enrichment analysis.

Results: Comorbidity analysis was conducted for 5273 DLD cases (mean [SD] age, 16.8 [7.2] years; 3748 males [71.1%]) and 26 353 matched controls (mean [SD] age, 14.6 [5.5] years; 18 729 males [71.1%]). Relevant phenotypes associated with DLD were found, including learning disorder, delayed milestones, disorders of the acoustic nerve, conduct disorders, attention-deficit/hyperactivity disorder, lack of coordination, and other motor deficits. Several other health phenotypes not previously associated with DLD were identified, such as dermatitis, conjunctivitis, and weight and nutrition, representing a new window into the clinical complexity of DLD.

Conclusions and relevance: This study found both rare and common comorbidities of DLD. Comorbidity profiles may be leveraged to identify risk of additional health challenges, beyond language impairment, among children with DLD.

Summary: Genomic data are often processed in batches and analyzed together to save time. However, it is challenging to combine multiple large VCFs and properly handle imputation quality and missing variants due to the limitations of available tools. To address these concerns, we developed IMMerge, a Python-based tool that takes advantage of multiprocessing to reduce running time. For the first time in a publicly available tool, imputation quality scores are correctly combined with Fisher's z transformation.

Availability and implementation: IMMerge is an open-source project under MIT license. Source code and user manual are available at https://github.com/belowlab/IMMerge.

Objective: Mexican Americans are disproportionally affected by non-alcoholic fatty liver disease (NAFLD), liver fibrosis and hepatocellular carcinoma. Noninvasive means to identify those in this population at high risk for these diseases are urgently needed. Approach: The Cameron County Hispanic Cohort (CCHC) is a population-based cohort with high rates of obesity (51%), type 2 diabetes (28%) and NAFLD (49%). In a subgroup of 564 CCHC subjects, we evaluated 339 genetic variants previously reported to be associated with liver injury markers aspartate aminotransferase (AST) and alanine aminotransferase (ALT) in United Kingdom and Japanese cohorts. Results: Association was confirmed for 86 variants. Among them, 27 had higher effect allele frequency in the CCHC than in the United Kingdom and Japanese cohorts, and 16 had stronger associations with AST and ALT than rs738409 (PNPLA3). These included rs17710008 (MYCT1), rs2519093 (ABO), rs1801690 (APOH), rs10409243 (S1PR2), rs1800759 (LOC100507053) and rs2491441 (RGL1), which were also associated with steatosis and/or liver fibrosis measured by vibration-controlled transient elastography. Main contributors to advanced fibrosis risk were rs11240351 (CNTN2), rs1800759 (LOC100507053), rs738409 (PNPLA3) and rs1801690 (APOH), with advanced fibrosis detected in 37.5% of subjects with 3 of these 4 variants [AOR = 11.6 (95% CI) = 3.8-35.3]. AST- and ALT-associated variants implicated distinct pathways (ethanol and galactose degradation versus antigen presentation and B cell development). Finally, 8 variants, including rs62292950 (DNAJC13), were associated with gut microbiome changes. Conclusion: These genotype-phenotype findings may have utility in risk modeling and disease prevention in this high-risk population.

[This corrects the article DOI: 10.1016/j.xhgg.2022.100099.].

Background: Concurrent variation in adiposity and inflammation suggests potential shared functional pathways and pleiotropic disease underpinning. Yet, exploration of pleiotropy in the context of adiposity-inflammation has been scarce, and none has included self-identified Hispanic/Latino populations. Given the high level of ancestral diversity in Hispanic American population, genetic studies may reveal variants that are infrequent/monomorphic in more homogeneous populations.

Methods: Using multi-trait Adaptive Sum of Powered Score (aSPU) method, we examined individual and shared genetic effects underlying inflammatory (CRP) and adiposity-related traits (Body Mass Index [BMI]), and central adiposity (Waist to Hip Ratio [WHR]) in HLA participating in the Population Architecture Using Genomics and Epidemiology (PAGE) cohort (N = 35,871) with replication of effects in the Cameron County Hispanic Cohort (CCHC) which consists of Mexican American individuals.

Results: Of the > 16 million SNPs tested, variants representing 7 independent loci were found to illustrate significant association with multiple traits. Two out of 7 variants were replicated at statistically significant level in multi-trait analyses in CCHC. The lead variant on APOE (rs439401) and rs11208712 were found to harbor multi-trait associations with adiposity and inflammation.

Conclusions: Results from this study demonstrate the importance of considering pleiotropy for improving our understanding of the etiology of the various metabolic pathways that regulate cardiovascular disease development.

Keywords: Genetic pleiotropy; Hispanic Americans; Inflammation; Obesity.

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Erratum for

Identification of genetic effects underlying type 2 diabetes in South Asian and European populations.

Mexican Americans have a high prevalence of diabetes and burden of diabetes-related complications, highlighting the need for novel preventive strategies and noninvasive predictors of diabetes risk tailored to this population. Changes in the gut microbiome have the potential to predict diabetes. Here, we aimed to identify alterations in the gut microbiome associated with diabetes in the high-risk population of Mexican Americans in South Texas. Stool samples were collected from 216 subjects from the population-based Cameron County Hispanic Cohort. Among them, 75 had type 2 diabetes. Taxonomic and functional profiling of the stool samples were assessed by 16S and shotgun metagenomic sequencing, and the influence of genetic factors was explored. The gut microbiome of subjects with diabetes was enriched with proinflammatory Proteobacteria members (Enterobacteriaceae, Escherichia-Shigella) and depleted of butyrate-producing Clostridiales members (Faecalibacterium prausnitzii, Peptostreptococcaceae, and Clostridium sensu stricto 1). The accompanying metagenomic changes in subjects with diabetes suggested dysregulated amino acid metabolism, reduced galacturonate and glucuronate catabolism (correlating with Faecalibacterium prausnitzii abundance), and enriched heme biosynthesis (correlating with Enterobacteriaceae abundance). Polymorphism rs7129790 near MMP27 was strongly associated with high Proteobacteria abundance and was more frequent in this cohort and in individuals of Mexican ancestry than in Europeans. In conclusion, Mexican Americans in South Texas with diabetes display distinct gut microbiome and metagenomic signatures. These signatures may have utility in risk modeling and disease prevention in this high-risk population. IMPORTANCE The gut microbiome composition varies across ethnicities and geographical locations, yet studies on diabetes-associated microbiome changes specific to high-risk Mexican Americans are lacking. Here, we aimed to identify specific alterations associated with diabetes in this population, as well as host genetic factors that may explain increased disease susceptibility in this ethnic group. Using samples from a population-based cohort of Mexican Americans with a high prevalence of obesity and diabetes, we confirmed findings from studies on other ethnicities that suggested promotion of a chronic proinflammatory environment, loss of butyrate production, and compromised intestinal barrier integrity. High abundance of proinflammatory Proteobacteria was associated with a polymorphism that was more frequent in this cohort and in individuals of Mexican ancestry than in Europeans. Validation of microbiome-based risk models for diabetes should be evaluated in prospective cohort studies.

Keywords: Mexican American population; diabetes; gut microbiome; health disparity; metagenome.

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Keywords: Hispanic/Latino; anthropometrics; diversity; fine-mapping; obesity; population stratification; trans-ancestral or trans-ethnic.

South Asians are at high risk of developing type 2 diabetes (T2D). We carried out a genome-wide association meta-analysis with South Asian T2D cases (n = 16,677) and controls (n = 33,856), followed by combined analyses with Europeans (neff = 231,420). We identify 21 novel genetic loci for significant association with T2D (P = 4.7 × 10-8 to 5.2 × 10-12), to the best of our knowledge at the point of analysis. The loci are enriched for regulatory features, including DNA methylation and gene expression in relevant tissues, and highlight CHMP4B, PDHB, LRIG1 and other genes linked to adiposity and glucose metabolism. A polygenic risk score based on South Asian-derived summary statistics shows ~4-fold higher risk for T2D between the top and bottom quartile. Our results provide further insights into the genetic mechanisms underlying T2D, and highlight the opportunities for discovery from joint analysis of data from across ancestral populations.

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele.

Prosody perception is fundamental to spoken language communication as it supports comprehension, pragmatics, morphosyntactic parsing of speech streams, and phonological awareness. A particular aspect of prosody: perceptual sensitivity to speech rhythm patterns in words (i.e., lexical stress sensitivity), is also a robust predictor of reading skills, though it has received much less attention than phonological awareness in the literature. Given the importance of prosody and reading in educational outcomes, reliable and valid tools are needed to conduct large-scale health and genetic investigations of individual differences in prosody, as groundwork for investigating the biological underpinnings of the relationship between prosody and reading. Motivated by this need, we present the Test of Prosody via Syllable Emphasis ("TOPsy") and highlight its merits as a phenotyping tool to measure lexical stress sensitivity in as little as 10 min, in scalable internet-based cohorts. In this 28-item speech rhythm perception test [modeled after the stress identification test from Wade-Woolley (2016)], participants listen to multi-syllabic spoken words and are asked to identify lexical stress patterns. Psychometric analyses in a large internet-based sample shows excellent reliability, and predictive validity for self-reported difficulties with speech-language, reading, and musical beat synchronization. Further, items loaded onto two distinct factors corresponding to initially stressed vs. non-initially stressed words. These results are consistent with previous reports that speech rhythm perception abilities correlate with musical rhythm sensitivity and speech-language/reading skills, and are implicated in reading disorders (e.g., dyslexia). We conclude that TOPsy can serve as a useful tool for studying prosodic perception at large scales in a variety of different settings, and importantly can act as a validated brief phenotype for future investigations of the genetic architecture of prosodic perception, and its relationship to educational outcomes.

Keywords: dyslexia; phenotyping; prosody; rhythm; speech.

Type 2 diabetes is a complex, systemic disease affected by both genetic and environmental factors. Previous research has identified genetic variants associated with type 2 diabetes risk; however, gene regulatory changes underlying progression to metabolic dysfunction are still largely unknown. We investigated RNA expression changes that occur during diabetes progression using a two-stage approach. In our discovery stage, we compared changes in gene expression using two longitudinally collected blood samples from subjects whose fasting blood glucose transitioned to a level consistent with type 2 diabetes diagnosis between the time points against those who did not with a novel analytical network approach. Our network methodology identified 17 networks, one of which was significantly associated with transition status. This 822-gene network harbors many genes novel to the type 2 diabetes literature but is also significantly enriched for genes previously associated with type 2 diabetes. In the validation stage, we queried associations of genetically determined expression with diabetes-related traits in a large biobank with linked electronic health records. We observed a significant enrichment of genes in our identified network whose genetically determined expression is associated with type 2 diabetes and other metabolic traits and validated 31 genes that are not near previously reported type 2 diabetes loci. Finally, we provide additional functional support, which suggests that the genes in this network are regulated by enhancers that operate in human pancreatic islet cells. We present an innovative and systematic approach that identified and validated key gene expression changes associated with type 2 diabetes transition status and demonstrated their translational relevance in a large clinical resource.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common, severe craniofacial malformation that imposes significant medical, psychosocial and financial burdens. NSCL/P is a multifactorial disorder with genetic and environmental factors playing etiologic roles. Currently, only 25% of the genetic variation underlying NSCL/P has been identified by linkage, candidate gene and genome-wide association studies. In this study, whole-genome sequencing and genome-wide genotyping followed by polygenic risk score (PRS) and linkage analyses were used to identify the genetic etiology of NSCL/P in a large three-generation family. We identified a rare missense variant in PDGFRA (c.C2740T; p.R914W) as potentially etiologic in a gene-based association test using pVAAST (P = 1.78 × 10-4) and showed decreased penetrance. PRS analysis suggested that variant penetrance was likely modified by common NSCL/P risk variants, with lower scores found among unaffected carriers. Linkage analysis provided additional support for PRS-modified penetrance, with a 7.4-fold increase in likelihood after conditioning on PRS. Functional characterization experiments showed that the putatively causal variant was null for signaling activity in vitro; further, perturbation of pdgfra in zebrafish embryos resulted in unilateral orofacial clefting. Our findings show that a rare PDGFRA variant, modified by additional common NSCL/P risk variants, have a profound effect on NSCL/P risk. These data provide compelling evidence for multifactorial inheritance long postulated to underlie NSCL/P and may explain some unusual familial patterns.

Despite a lifetime prevalence of at least 5%, developmental stuttering, characterized by prolongations, blocks, and repetitions of speech sounds, remains a largely idiopathic speech disorder. Family, twin, and segregation studies overwhelmingly support a strong genetic influence on stuttering risk; however, its complex mode of inheritance combined with thus-far underpowered genetic studies contribute to the challenge of identifying and reproducing genes implicated in developmental stuttering susceptibility. We conducted a trans-ancestry genome-wide association study (GWAS) and meta-analysis of developmental stuttering in two primary datasets: The International Stuttering Project comprising 1,345 clinically ascertained cases from multiple global sites and 6,759 matched population controls from the biobank at Vanderbilt University Medical Center (VUMC), and 785 self-reported stuttering cases and 7,572 controls ascertained from The National Longitudinal Study of Adolescent to Adult Health (Add Health). Meta-analysis of these genome-wide association studies identified a genome-wide significant (GWS) signal for clinically reported developmental stuttering in the general population: a protective variant in the intronic or genic upstream region of SSUH2 (rs113284510, protective allele frequency = 7.49%, Z = -5.576, p = 2.46 × 10-8) that acts as an expression quantitative trait locus (eQTL) in esophagus-muscularis tissue by reducing its gene expression. In addition, we identified 15 loci reaching suggestive significance (p < 5 × 10-6). This foundational population-based genetic study of a common speech disorder reports the findings of a clinically ascertained study of developmental stuttering and highlights the need for further research.

Keywords: complex trait; genome-wide assocation study; population-study; stuttering.

Aims/hypothesis: Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study.

Methods: We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci.

Results: Four novel associations were identified (p < 5 × 10-9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis.

Conclusions/interpretation: Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations.

Data availability: Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog ( https://www.ebi.ac.uk/gwas/downloads/summary-statistics ).

Keywords: Fine-mapping; Genome-wide association study; Glucose; Glycaemic traits; HbA1c; Insulin; Transethnic population.

Late-onset Alzheimer disease (LOAD) is highly polygenic, with a heritability estimated between 40 and 80%, yet risk variants identified in genome-wide studies explain only ~8% of phenotypic variance. Due to its increased power and interpretability, genetically regulated expression (GReX) analysis is an emerging approach to investigate the genetic mechanisms of complex diseases. Here, we conducted GReX analysis within and across 51 tissues on 39 LOAD GWAS data sets comprising 58,713 cases and controls from the Alzheimer's Disease Genetics Consortium (ADGC) and the International Genomics of Alzheimer's Project (IGAP). Meta-analysis across studies identified 216 unique significant genes, including 72 with no previously reported LOAD GWAS associations. Cross-brain-tissue and cross-GTEx models revealed eight additional genes significantly associated with LOAD. Conditional analysis of previously reported loci using established LOAD-risk variants identified eight genes reaching genome-wide significance independent of known signals. Moreover, the proportion of SNP-based heritability is highly enriched in genes identified by GReX analysis. In summary, GReX-based meta-analysis in LOAD identifies 216 genes (including 72 novel genes), illuminating the role of gene regulatory models in LOAD.

Developmental stuttering is a speech disorder characterized by disruption in the forward movement of speech. This disruption includes part-word and single-syllable repetitions, prolongations, and involuntary tension that blocks syllables and words, and the disorder has a life-time prevalence of 6-12%. Within Vanderbilt's electronic health record (EHR)-linked biorepository (BioVU), only 142 individuals out of 92,762 participants (0.15%) are identified with diagnostic ICD9/10 codes, suggesting a large portion of people who stutter do not have a record of diagnosis within the EHR. To identify individuals affected by stuttering within our EHR, we built a PheCode-driven Gini impurity-based classification and regression tree model, PheML, by using comorbidities enriched in individuals affected by stuttering as predicting features and imputing stuttering status as the outcome variable. Applying PheML in BioVU identified 9,239 genotyped affected individuals (a clinical prevalence of ∼10%) for downstream genetic analysis. Ancestry-stratified GWAS of PheML-imputed affected individuals and matched control individuals identified rs12613255, a variant near CYRIA on chromosome 2 (B = 0.323; p value = 1.31 × 10-8) in European-ancestry analysis and rs7837758 (B = 0.518; p value = 5.07 × 10-8), an intronic variant found within the ZMAT4 gene on chromosome 8, in African-ancestry analysis. Polygenic-risk prediction and concordance analysis in an independent clinically ascertained sample of developmental stuttering cases validate our GWAS findings in PheML-imputed affected and control individuals and demonstrate the clinical relevance of our population-based analysis for stuttering risk.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Background and aims: Hispanics are disproportionately affected by NAFLD, liver fibrosis, cirrhosis, and HCC. Preventive strategies and noninvasive means to identify those in this population at high risk for liver fibrosis, are urgently needed. We aimed to characterize the gut microbiome signatures and related biological functions associated with liver fibrosis in Hispanics and identify environmental and genetic factors affecting them.

Approach and results: Subjects of the population-based Cameron County Hispanic Cohort (CCHC; n = 217) were screened by vibration-controlled transient elastography (FibroScan). Among them, 144 (66.7%) had steatosis and 28 (13.0%) had liver fibrosis. The gut microbiome of subjects with liver fibrosis was enriched with immunogenic commensals (e.g., Prevotella copri, Holdemanella, Clostridiaceae 1) and depleted of Bacteroides caccae, Parabacteroides distasonis, Enterobacter, and Marinifilaceae. The liver fibrosis-associated metagenome was characterized by changes in the urea cycle, L-citrulline biosynthesis and creatinine degradation pathways, and altered synthesis of B vitamins and lipoic acid. These metagenomic changes strongly correlated with the depletion of Parabacteroides distasonis and enrichment of Prevotella and Holdemanella. Liver fibrosis was also associated with depletion of bacterial pathways related to L-fucose biosynthesis. Alcohol consumption, even moderate, was associated with high Prevotella abundance. The single-nucleotide polymorphisms rs3769502 and rs7573751 in the NCK adaptor protein 2 (NCK2) gene positively associated with high Prevotella abundance.

Conclusion: Hispanics with liver fibrosis display microbiome profiles and associated functional changes that may promote oxidative stress and a proinflammatory environment. These microbiome signatures, together with NCK2 polymorphisms, may have utility in risk modeling and disease prevention in this high-risk population.

Genomic discovery and characterization of risk loci for type 2 diabetes (T2D) have been conducted primarily in individuals of European ancestry. We conducted a multiethnic genome-wide association study of T2D among 53,102 cases and 193,679 control subjects from African, Hispanic, Asian, Native Hawaiian, and European population groups in the Population Architecture Genomics and Epidemiology (PAGE) and Diabetes Genetics Replication and Meta-analysis (DIAGRAM) Consortia. In individuals of African ancestry, we discovered a risk variant in the TGFB1 gene (rs11466334, risk allele frequency (RAF) = 6.8%, odds ratio [OR] = 1.27, p = 2.06 × 10-8), which replicated in independent studies of African ancestry (p = 6.26 × 10-23). We identified a multiethnic risk variant in the BACE2 gene (rs13052926, RAF = 14.1%, OR = 1.08, p = 5.75 × 10-9), which also replicated in independent studies (p = 3.45 × 10-4). We also observed a significant difference in the performance of a multiethnic genetic risk score (GRS) across population groups (pheterogeneity = 3.85 × 10-20). Comparing individuals in the top GRS risk category (40%-60%), the OR was highest in Asians (OR = 3.08) and European (OR = 2.94) ancestry populations, followed by Hispanic (OR = 2.39), Native Hawaiian (OR = 2.02), and African ancestry (OR = 1.57) populations. These findings underscore the importance of genetic discovery and risk characterization in diverse populations and the urgent need to further increase representation of non-European ancestry individuals in genetics research to improve genetic-based risk prediction across populations.

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

In the 40 years since Some Statistical Methods Useful in Circulation Research was published, many of the same battery of statistical tests and concepts, such as t tests, ANOVA, P values, effect sizes, and SEs are still abundantly employed in hypothesis-driven research.1 Newer methods, too, have emerged to address the challenges of big data analysis. Some methods now routinely employed to extract insights from data include regression analysis, supervised and unsupervised machine learning for clustering, density estimation, and dimensionality reduction (eg, viSNE), as well as prediction modeling and enrichment analyses. Additionally, in basic science research, it is now common to encounter hypothesis-free analyses, in marked contrast to traditional statistical analyses that begin with an explicit hypothesis. To encourage reproducibility, rigor, interpretability, and transparency, many editorial teams, including those of Circulation Research and the AHA journals, have developed statistical guidelines for authors. Given the rapidly changing data landscape, such guidelines must extend beyond “What statistical test should I use?” (a question that often can be addressed by a decision tree diagram in applied statistical analysis textbooks), to address higher-level challenges that frequently face authors including multiple testing, standards of reporting, robustness to violations of assumptions, and the limitations of conventional measures of significance. To better support authors and readers, we have assembled some topics that warrant particular attention in basic and clinical scientific publications such as those published in Circulation Research. These guidelines are intended to complement those outlined by the American Heart Association’s Statistical Taskforce in their concurrent “Guidelines for Statistical Reporting in Cardiovascular Medicine: A Special Report from the American Heart Association.”

Statistical analyses are a crucial component of the biomedical research process and are necessary to draw inferences from biomedical research data. The application of sound statistical methodology is a prerequisite for publication in the American Heart Association (AHA) journal portfolio. The objective of this document is to summarize key aspects of statistical reporting that might be most relevant to the authors, reviewers, and readership of AHA journals. The AHA Scientific Publication Committee convened a task force to inventory existing statistical standards for publication in biomedical journals and to identify approaches suitable for the AHA journal portfolio. The experts on the task force were selected by the AHA Scientific Publication Committee, who identified 12 key topics that serve as the section headers for this document. For each topic, the members of the writing group identified relevant references and evaluated them as a resource to make the standards summarized herein. Each section was independently reviewed by an expert reviewer who was not part of the task force. Expert reviewers were also permitted to comment on other sections if they chose. Differences of opinion were adjudicated by consensus. The standards presented in this report are intended to serve as a guide for high-quality reporting of statistical analyses methods and results.

Purpose: This study aimed to identify cases of developmental stuttering and associated comorbidities in de-identified electronic health records (EHRs) at Vanderbilt University Medical Center, and, in turn, build and test a stuttering prediction model.

Methods: A multi-step process including a keyword search of medical notes, a text-mining algorithm, and manual review was employed to identify stuttering cases in the EHR. Confirmed cases were compared to matched controls in a phenotype code (phecode) enrichment analysis to reveal conditions associated with stuttering (i.e., comorbidities). These associated phenotypes were used as proxy variables to phenotypically predict stuttering in subjects within the EHR that were not otherwise identifiable using the multi-step identification process described above.

Results: The multi-step process resulted in the manually reviewed identification of 1,143 stuttering cases in the EHR. Highly enriched phecodes included codes related to childhood onset fluency disorder, adult-onset fluency disorder, hearing loss, sleep disorders, atopy, a multitude of codes for infections, neurological deficits, and body weight. These phecodes were used as variables to create a phenome risk classifier (PheRC) prediction model to identify additional high likelihood stuttering cases. The PheRC prediction model resulted in a positive predictive value of 83 %.

Conclusions: This study demonstrates the feasibility of using EHRs in the study of stuttering and found phenotypic associations. The creation of the PheRC has the potential to enable future studies of stuttering using existing EHR data, including investigations into the genetic etiology.

Keywords: Developmental stuttering; Electronic health records; Machine learning; Stuttering comorbidities.

A primary challenge in the analysis of free-ranging animal populations is the accurate estimation of relatedness among individuals. Many aspects of population analysis rely on knowledge of relatedness patterns, including socioecology, demography, heritability and gene mapping analyses, wildlife conservation and the management of breeding colonies. Methods for determining relatedness using genome-wide data have improved our ability to determine kinship and reconstruct pedigrees in humans. However, methods for reconstructing complex pedigree structures and estimating distant relatedness (beyond third-degree) have not been widely applied to other species. We sequenced the genomes of 150 male rhesus macaques from the Tulane National Primate Research Center colony to estimate pairwise relatedness, reconstruct closely related pedigrees, estimate more distant relationships and augment colony records. Methods for determining relatedness developed for human genetic data were applied and evaluated in the analysis of nonhuman primates, including identity-by-descent-based methods for pedigree reconstruction and shared segment-based inference of more distant relatedness. We compared the genotype-based pedigrees and estimated relationships to available colony pedigree records and found high concordance (95.5% agreement) between expected and identified relationships for close relatives. In addition, we detected distant relationships not captured in colony records, including some as distant as twelfth-degree. Furthermore, while deep sequence coverage is preferable, we show that this approach can also provide valuable information when only low-coverage (5×) sequence data is available. Our findings demonstrate the value of these methods for determination of relatedness in various animal populations, with diverse applications to conservation biology, evolutionary and ecological research and biomedical studies.

Keywords: bioinfomatics/phyloinfomatics; kinship; mammals; quantitative genetics.

Given the coronavirus disease 2019 (COVID-19) pandemic, investigations into host susceptibility to infectious diseases and downstream sequelae have never been more relevant. Pneumonia is a lung disease that can cause respiratory failure and hypoxia and is a common complication of infectious diseases, including COVID-19. Few genome-wide association studies (GWASs) of host susceptibility and severity of pneumonia have been conducted. We performed GWASs of pneumonia susceptibility and severity in the Vanderbilt University biobank (BioVU) with linked electronic health records (EHRs), including Illumina Expanded Multi-Ethnic Global Array (MEGAEX)-genotyped European ancestry (EA, n= 69,819) and African ancestry (AA, n = 15,603) individuals. Two regions of large effect were identified: the CFTR locus in EA (rs113827944; OR = 1.84, p value = 1.2 × 10-36) and HBB in AA (rs334 [p.Glu7Val]; OR = 1.63, p value = 3.5 × 10-13). Mutations in these genes cause cystic fibrosis (CF) and sickle cell disease (SCD), respectively. After removing individuals diagnosed with CF and SCD, we assessed heterozygosity effects at our lead variants. Further GWASs after removing individuals with CF uncovered an additional association in R3HCC1L (rs10786398; OR = 1.22, p value = 3.5 × 10-8), which was replicated in two independent datasets: UK Biobank (n = 459,741) and 7,985 non-overlapping BioVU subjects, who are genotyped on arrays other than MEGAEX. This variant was also validated in GWASs of COVID-19 hospitalization and lung function. Our results highlight the importance of the host genome in infectious disease susceptibility and severity and offer crucial insight into genetic effects that could potentially influence severity of COVID-19 sequelae.

Keywords: GWAS; biobank; electronic health record; host genetic effect; pneumonia.

No absract available

Purpose Data mining algorithms using electronic health records (EHRs) are useful in large-scale population-wide studies to classify etiology and comorbidities (Casey et al., 2016). Here, we apply this approach to developmental language disorder (DLD), a prevalent communication disorder whose risk factors and epidemiology remain largely undiscovered. Method We first created a reliable system for manually identifying DLD in EHRs based on speech-language pathologist (SLP) diagnostic expertise. We then developed and validated an automated algorithmic procedure, called, Automated Phenotyping Tool for identifying DLD cases in health systems data (APT-DLD), that classifies a DLD status for patients within EHRs on the basis of ICD (International Statistical Classification of Diseases and Related Health Problems) codes. APT-DLD was validated in a discovery sample (N = 973) using expert SLP manual phenotype coding as a gold-standard comparison and then applied and further validated in a replication sample of N = 13,652 EHRs. Results In the discovery sample, the APT-DLD algorithm correctly classified 98% (concordance) of DLD cases in concordance with manually coded records in the training set, indicating that APT-DLD successfully mimics a comprehensive chart review. The output of APT-DLD was also validated in relation to independently conducted SLP clinician coding in a subset of records, with a positive predictive value of 95% of cases correctly classified as DLD. We also applied APT-DLD to the replication sample, where it achieved a positive predictive value of 90% in relation to SLP clinician classification of DLD. Conclusions APT-DLD is a reliable, valid, and scalable tool for identifying DLD cohorts in EHRs. This new method has promising public health implications for future large-scale epidemiological investigations of DLD and may inform EHR data mining algorithms for other communication disorders. Supplemental Material

One of the forerunners that pioneered the revolution of high-throughput genomic technologies is the genotyping microarray technology, which can genotype millions of single-nucleotide variants simultaneously. Owing to apparent benefits, such as high speed, low cost, and high throughput, the genotyping array has gained lasting applications in genome-wide association studies (GWAS) and thus accumulated an enormous amount of data. Empowered by continuous manufactural upgrades and analytical innovation, unconventional applications of genotyping array data have emerged to address more diverse genetic problems, holding promise of boosting genetic research into human diseases through the re-mining of the rich accumulated data. Here, we review several unconventional genotyping array analysis techniques that have been built on the idea of large-scale multivariant analysis and provide empirical application examples. These unconventional outcomes of genotyping arrays include polygenic score, runs of homozygosity (ROH)/heterozygosity ratio, distant pedigree computation, and mitochondrial DNA (mtDNA) copy number inference.

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Genome-wide association studies have revolutionized our understanding of the genetic underpinnings of cardiometabolic disease. Yet, the inadequate representation of individuals of diverse ancestral backgrounds in these studies may undercut their ultimate potential for both public health and precision medicine. The goal of this review is to describe the imperativeness of studying the populations who are most affected by cardiometabolic disease, to the aim of better understanding the genetic underpinnings of the disease. We support this premise by describing the current variation in the global burden of cardiometabolic disease and emphasize the importance of building a globally and ancestrally representative genetics evidence base for the identification of population-specific variants, fine-mapping, and polygenic risk score estimation. We discuss the important ethical, legal, and social implications of increasing ancestral diversity in genetic studies of cardiometabolic disease and the challenges that arise from the (1) lack of diversity in current reference populations and available analytic samples and the (2) unequal generation of health-associated genomic data and their prediction accuracies. Despite these challenges, we conclude that additional, unprecedented opportunities lie ahead for public health genomics and the realization of precision medicine, provided that the gap in diversity can be systematically addressed. Achieving this goal will require concerted efforts by social, academic, professional and regulatory stakeholders and communities, and these efforts must be based on principles of equity and social justice.

Keywords: cardiovascular diseases; genomics; global burden of disease; metabolic diseases; minority health; precision medicine; social justice.

Lipid levels are important markers for the development of cardio-metabolic diseases. Although hundreds of associated loci have been identified through genetic association studies, the contribution of genetic factors to variation in lipids is not fully understood, particularly in U.S. minority groups. We performed genome-wide association analyses for four lipid traits in over 45,000 ancestrally diverse participants from the Population Architecture using Genomics and Epidemiology (PAGE) Study, followed by a meta-analysis with several European ancestry studies. We identified nine novel lipid loci, five of which showed evidence of replication in independent studies. Furthermore, we discovered one novel gene in a PrediXcan analysis, minority-specific independent signals at eight previously reported loci, and potential functional variants at two known loci through fine-mapping. Systematic examination of known lipid loci revealed smaller effect estimates in African American and Hispanic ancestry populations than those in Europeans, and better performance of polygenic risk scores based on minority-specific effect estimates. Our findings provide new insight into the genetic architecture of lipid traits and highlight the importance of conducting genetic studies in diverse populations in the era of precision medicine.

Genetic ataxias are associated with mutations in hundreds of genes with high phenotypic overlap complicating the clinical diagnosis. Whole-exome sequencing (WES) has increased the overall diagnostic rate considerably. However, the upper limit of this method remains ill-defined, hindering efforts to address the remaining diagnostic gap. To further assess the role of rare coding variation in ataxic disorders, we reanalyzed our previously published exome cohort of 76 predominantly adult and sporadic-onset patients, expanded the total number of cases to 260, and introduced analyses for copy number variation and repeat expansion in a representative subset. For new cases (n = 184), our resulting clinically relevant detection rate remained stable at 47% with 24% classified as pathogenic. Reanalysis of the previously sequenced 76 patients modestly improved the pathogenic rate by 7%. For the combined cohort (n = 260), the total observed clinical detection rate was 52% with 25% classified as pathogenic. Published studies of similar neurological phenotypes report comparable rates. This consistency across multiple cohorts suggests that, despite continued technical and analytical advancements, an approximately 50% diagnostic rate marks a relative ceiling for current WES-based methods and a more comprehensive genome-wide assessment is needed to identify the missing causative genetic etiologies for cerebellar ataxia and related neurodegenerative diseases.

Keywords: ataxia; cerebellar ataxia; cerebellum; diagnostic testing; exome; gait disorders; genetics; genomics; neurogenetics; spastic paraparesis; spastic paraplegia; spinocerebellar ataxia.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Background: The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

Methods: We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

Results: We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

Conclusions: Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

No abstract available

Keywords: atrial fibrillation; atrial fibrosis; family history; genetics.

Sleep disordered breathing (SDB)-related overnight hypoxemia is associated with cardiometabolic disease and other comorbidities. Understanding the genetic bases for variations in nocturnal hypoxemia may help understand mechanisms influencing oxygenation and SDB-related mortality. We conducted genome-wide association tests across 10 cohorts and 4 populations to identify genetic variants associated with three correlated measures of overnight oxyhemoglobin saturation: average and minimum oxyhemoglobin saturation during sleep and the percent of sleep with oxyhemoglobin saturation under 90%. The discovery sample consisted of 8,326 individuals. Variants with p < 1 × 10(-6) were analyzed in a replication group of 14,410 individuals. We identified 3 significantly associated regions, including 2 regions in multi-ethnic analyses (2q12, 10q22). SNPs in the 2q12 region associated with minimum SpO2 (rs78136548 p = 2.70 × 10(-10)). SNPs at 10q22 were associated with all three traits including average SpO2 (rs72805692 p = 4.58 × 10(-8)). SNPs in both regions were associated in over 20,000 individuals and are supported by prior associations or functional evidence. Four additional significant regions were detected in secondary sex-stratified and combined discovery and replication analyses, including a region overlapping Reelin, a known marker of respiratory complex neurons.These are the first genome-wide significant findings reported for oxyhemoglobin saturation during sleep, a phenotype of high clinical interest. Our replicated associations with HK1 and IL18R1 suggest that variants in inflammatory pathways, such as the biologically-plausible NLRP3 inflammasome, may contribute to nocturnal hypoxemia.

Purpose of review: Neurodegenerative diseases, neuropsychiatric disorders, and related traits have highly complex etiologies but are also highly heritable and identifying the causal genes and biological pathways underlying these traits may advance the development of treatments and preventive strategies. While many genome-wide association studies (GWAS) have successfully identified variants contributing to polygenic neurodegenerative and neuropsychiatric phenotypes including Alzheimer's disease (AD), schizophrenia (SCZ), and bipolar disorder (BPD) amongst others, interpreting the biological roles of significantly-associated variants in the genetic architecture of these traits remains a significant challenge. Here we review several 'omics' approaches which attempt to bridge the gap from associated genetic variants to phenotype by helping define the functional roles of GWAS loci in the development of neuropsychiatric disorders and traits.

Recent findings: Several common 'omics' approaches have been applied to examine neuropsychiatric traits, such as nearest-gene mapping, trans-ethnic fine mapping, annotation enrichment analysis, transcriptomic analysis, and pathway analysis, and each of these approaches has strengths and limitations in providing insight into biological mechanisms. One popular emerging method is the examination of tissue-specific genetically-regulated gene expression (GReX), which aggregates the genetic variants' effects at the gene-level. Furthermore, proteomic, metabolomic, and microbiomic studies and phenome-wide association studies will further enhance our understanding of neuropsychiatric traits.

Summary: GWAS has been applied to neuropsychiatric traits for a decade, but our understanding about the biological function of identified variants remains limited. Today, technological advancements have created analytical approaches for integrating transcriptomics, metabolomics, proteomics, pharmacology and toxicology as tools for understanding the functional roles of genetics variants. These data, as well as the broader clinical information provided by electronic health records, can provide additional insight and complement genomic analyses.

Keywords: Functional Annotation; Functional Interpretation; Genetically Regulated Expression; Genome-Wide Association Studies; Omics.

Interpretation of genetic association results is difficult because signals often lack biological context. To generate hypotheses of the functional genetic etiology of complex cardiometabolic traits, we estimated the genetically determined component of gene expression from common variants using PrediXcan (1) and determined genes with differential predicted expression by trait. PrediXcan imputes tissue-specific expression levels from genetic variation using variant-level effect on gene expression in transcriptome data. To explore the value of imputed genetically regulated gene expression (GReX) models across different ancestral populations, we evaluated imputed expression levels for predictive accuracy genome-wide in RNA sequence data in samples drawn from European-ancestry and African-ancestry populations and identified substantial predictive power using European-derived models in a non-European target population. We then tested the association of GReX on 15 cardiometabolic traits including blood lipid levels, body mass index, height, blood pressure, fasting glucose and insulin, RR interval, fibrinogen level, factor VII level and white blood cell and platelet counts in 15 755 individuals across three ancestry groups, resulting in 20 novel gene-phenotype associations reaching experiment-wide significance across ancestries. In addition, we identified 18 significant novel gene-phenotype associations in our ancestry-specific analyses. Top associations were assessed for additional support via query of S-PrediXcan (2) results derived from publicly available genome-wide association studies summary data. Collectively, these findings illustrate the utility of transcriptome-based imputation models for discovery of cardiometabolic effect genes in a diverse dataset.

Obstructive sleep apnea (OSA) is a common disorder associated with increased risk of cardiovascular disease and mortality. Its prevalence and severity vary across ancestral background. Although OSA traits are heritable, few genetic associations have been identified. To identify genetic regions associated with OSA and improve statistical power, we applied admixture mapping on three primary OSA traits [the apnea hypopnea index (AHI), overnight average oxyhemoglobin saturation (SaO2) and percentage time SaO2 < 90%] and a secondary trait (respiratory event duration) in a Hispanic/Latino American population study of 11 575 individuals with significant variation in ancestral background. Linear mixed models were performed using previously inferred African, European and Amerindian local genetic ancestry markers. Global African ancestry was associated with a lower AHI, higher SaO2 and shorter event duration. Admixture mapping analysis of the primary OSA traits identified local African ancestry at the chromosomal region 2q37 as genome-wide significantly associated with AHI (P < 5.7 × 10-5), and European and Amerindian ancestries at 18q21 suggestively associated with both AHI and percentage time SaO2 < 90% (P < 10-3). Follow-up joint ancestry-SNP association analyses identified novel variants in ferrochelatase (FECH), significantly associated with AHI and percentage time SaO2 < 90% after adjusting for multiple tests (P < 8 × 10-6). These signals contributed to the admixture mapping associations and were replicated in independent cohorts. In this first admixture mapping study of OSA, novel associations with variants in the iron/heme metabolism pathway suggest a role for iron in influencing respiratory traits underlying OSA.

Background: Genome-wide association studies have implicated the transcription factor 7-like 2 (TCF7L2) gene in type 2 diabetes risk, and more recently, in decreased body mass index. Given the contrary direction of genetic effects on these two traits, it has been suggested that the observed association with body mass index may reflect either selection bias or a complex underlying biology at TCF7L2.

Methods: Using 9031 Hispanic/Latino adults (21-76 years) with complete weight history and genetic data from the community-based Hispanic Community Health Study/Study of Latinos (HCHS/SOL, Baseline 2008-2011), we estimated the multivariable association between the additive number of type 2 diabetes increasing-alleles at TCF7L2 (rs7903146-T) and body mass index. We then used structural equation models to simultaneously model the genetic association on changes in body mass index across the life course and estimate the odds of type 2 diabetes per TCF7L2 risk allele.

Results: We observed both significant increases in type 2 diabetes prevalence at examination (independent of body mass index) and decreases in mean body mass index and waist circumference across genotypes at rs7903146. We observed a significant multivariable association between the additive number of type 2 diabetes-risk alleles and lower body mass index at examination. In our structured modeling, we observed non-significant inverse direct associations between rs7903146-T and body mass index at ages 21 and 45 years, and a significant positive association between rs7903146-T and type 2 diabetes onset in both middle and late adulthood.

Conclusions: Herein, we replicated the protective effect of rs7930146-T on body mass index at multiple time points in the life course, and observed that these effects were not explained by past type 2 diabetes status in our structured modeling. The robust replication of the negative effects of TCF7L2 on body mass index in multiple samples, including in our diverse Hispanic/Latino community-based sample, supports a growing body of literature on the complex biologic mechanism underlying the functional consequences of TCF7L2 on obesity and type 2 diabetes across the life course.

Keywords: Diabetes; Genetics; Hispanic/Latinos; Obesity; TCF7L2.

We expanded GWAS discovery for type 2 diabetes (T2D) by combining data from 898,130 European-descent individuals (9% cases), after imputation to high-density reference panels. With these data, we (i) extend the inventory of T2D-risk variants (243 loci, 135 newly implicated in T2D predisposition, comprising 403 distinct association signals); (ii) enrich discovery of lower-frequency risk alleles (80 index variants with minor allele frequency <5%, 14 with estimated allelic odds ratio >2); (iii) substantially improve fine-mapping of causal variants (at 51 signals, one variant accounted for >80% posterior probability of association (PPA)); (iv) extend fine-mapping through integration of tissue-specific epigenomic information (islet regulatory annotations extend the number of variants with PPA >80% to 73); (v) highlight validated therapeutic targets (18 genes with associations attributable to coding variants); and (vi) demonstrate enhanced potential for clinical translation (genome-wide chip heritability explains 18% of T2D risk; individuals in the extremes of a T2D polygenic risk score differ more than ninefold in prevalence).

The Alzheimer's Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10-7), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10-7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10-6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

No abstract available

Keywords: clinical trials; sample size; statistics

No abstract available

Background/aims: The Alzheimer's Disease Sequencing Project (ADSP) aims to identify novel genes influencing Alzheimer's disease (AD). Variants within genes known to cause dementias other than AD have previously been associated with AD risk. We describe evidence of co-segregation and associations between variants in dementia genes and clinically diagnosed AD within the ADSP.

Methods: We summarize the properties of known pathogenic variants within dementia genes, describe the co-segregation of variants annotated as "pathogenic" in ClinVar and new candidates observed in ADSP families, and test for associations between rare variants in dementia genes in the ADSP case-control study. The participants were clinically evaluated for AD, and they represent European, Caribbean Hispanic, and isolate Dutch populations.

Results/conclusions: Pathogenic variants in dementia genes were predominantly rare and conserved coding changes. Pathogenic variants within ARSA, CSF1R, and GRN were observed, and candidate variants in GRN and CHMP2B were nominated in ADSP families. An independent case-control study provided evidence of an association between variants in TREM2, APOE, ARSA, CSF1R, PSEN1, and MAPT and risk of AD. Variants in genes which cause dementing disorders may influence the clinical diagnosis of AD in a small proportion of cases within the ADSP.

Keywords: Alzheimer genetics; Alzheimer’s disease; Arylsulfatase A pseudodeficiency; Candidate genes; ClinVar; Frontotemporal dementia; Pathogenicity; Rare variants.

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis-expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants.

Keywords: eQTL; genetics; rare variants; sequencing; type 2 diabetes.

QT interval prolongation is a heritable risk factor for ventricular arrhythmias and can predispose to sudden death. Most genome-wide association studies (GWAS) of QT were performed in European ancestral populations, leaving other groups uncharacterized. Herein we present the first QT GWAS of Hispanic/Latinos using data on 15,997 participants from four studies. Study-specific summary results of the association between 1000 Genomes Project (1000G) imputed SNPs and electrocardiographically measured QT were combined using fixed-effects meta-analysis. We identified 41 genome-wide significant SNPs that mapped to 13 previously identified QT loci. Conditional analyses distinguished six secondary signals at NOS1AP (n = 2), ATP1B1 (n = 2), SCN5A (n = 1), and KCNQ1 (n = 1). Comparison of linkage disequilibrium patterns between the 13 lead SNPs and six secondary signals with previously reported index SNPs in 1000G super populations suggested that the SCN5A and KCNE1 lead SNPs were potentially novel and population-specific. Finally, of the 42 suggestively associated loci, AJAP1 was suggestively associated with QT in a prior East Asian GWAS; in contrast BVES and CAP2 murine knockouts caused cardiac conduction defects. Our results indicate that whereas the same loci influence QT across populations, population-specific variation exists, motivating future trans-ethnic and ancestrally diverse QT GWAS.

Obstructive sleep apnea (OSA) is a common heritable disorder displaying marked sexual dimorphism in disease prevalence and progression. Previous genetic association studies have identified a few genetic loci associated with OSA and related quantitative traits, but they have only focused on single ethnic groups, and a large proportion of the heritability remains unexplained. The apnea-hypopnea index (AHI) is a commonly used quantitative measure characterizing OSA severity. Because OSA differs by sex, and the pathophysiology of obstructive events differ in rapid eye movement (REM) and non-REM (NREM) sleep, we hypothesized that additional genetic association signals would be identified by analyzing the NREM/REM-specific AHI and by conducting sex-specific analyses in multiethnic samples. We performed genome-wide association tests for up to 19,733 participants of African, Asian, European, and Hispanic/Latino American ancestry in 7 studies. We identified rs12936587 on chromosome 17 as a possible quantitative trait locus for NREM AHI in men (N = 6,737; P = 1.7 × 10-8) but not in women (P = 0.77). The association with NREM AHI was replicated in a physiological research study (N = 67; P = 0.047). This locus overlapping the RAI1 gene and encompassing genes PEMT1, SREBF1, and RASD1 was previously reported to be associated with coronary artery disease, lipid metabolism, and implicated in Potocki-Lupski syndrome and Smith-Magenis syndrome, which are characterized by abnormal sleep phenotypes. We also identified gene-by-sex interactions in suggestive association regions, suggesting that genetic variants for AHI appear to vary by sex, consistent with the clinical observations of strong sexual dimorphism.

Keywords: genetics; genome-wide association studies; multiethnic; obstructive sleep apnea; sexual dimorphism.

We carried out an admixture mapping study of lipid traits in two samples from Mexico City. Native American locus ancestry was significantly associated with triglyceride levels in a broad region of chromosome 11 overlapping the BUD13, ZNF259 and APOA5 genes. In our fine-mapping analysis of this region using dense genome-wide data, rs964184 is the only marker included in the 99% credible set of SNPs, providing strong support for rs964184 as the causal variant within this region. The frequency of the allele associated with increased triglyceride concentrations (rs964184-G) is between 30-40% higher in Native American populations from Mexico than in European populations. The evidence currently available for this variant indicates that it may be exerting its effect through three potential mechanisms: 1) modification of enhancer activity, 2) regulation of the expression of several genes in cis and/or trans, or 3) modification of the methylation patterns of the promoter of the APOA5 gene.

Genetic determinants of sleep-disordered breathing (SDB), a common set of disorders that contribute to significant cardiovascular and neuropsychiatric morbidity, are not clear. Overnight nocturnal oxygen saturation (SaO2) is a clinically relevant and easily measured indicator of SDB severity but its genetic contribution has never been studied. Our recent study suggests nocturnal SaO2 is heritable. We performed linkage analysis, association analysis and haplotype analysis of average nocturnal oxyhaemoglobin saturation in participants in the Cleveland Family Study (CFS), followed by gene-based association and additional tests in four independent samples. Linkage analysis identified a peak (LOD = 4.29) on chromosome 8p23. Follow-up association analysis identified two haplotypes in angiopoietin-2 (ANGPT2) that significantly contributed to the variation of SaO2 (P = 8 × 10-5) and accounted for a portion of the linkage evidence. Gene-based association analysis replicated the association of ANGPT2 and nocturnal SaO2. A rare missense SNP rs200291021 in ANGPT2 was associated with serum angiopoietin-2 level (P = 1.29 × 10-4), which was associated with SaO2 (P = 0.002). Our study provides the first evidence for the association of ANGPT2, a gene previously implicated in acute lung injury syndromes, with nocturnal SaO2, suggesting that this gene has a broad range of effects on gas exchange, including influencing oxygenation during sleep.

Background: There is an increasing appreciation for a series of less traditional risk factors that should not be ignored when considering type 2 diabetes, obesity, hypertension, and cardiovascular disease. These include aortic stiffness, cardiac structure, impaired endothelial function and obstructive sleep apnea. They are associated to varying degrees with each disease categorization and with each other. It is not clear whether they represent additional complications, concomitants or antecedents of disease. Starr County, Texas, with its predominantly Mexican American population has been shown previously to bear a disproportionate burden of the major disease categories, but little is known about the distribution of these less traditional factors.

Methods: Type 2 diabetes, obesity and hypertension frequencies were determined through a systematic survey of Starr County conducted from 2002 to 2006. Individuals from this examination and an enriched set with type 2 diabetes were re-examined from 2010 to 2014 including assessment of cardiac structure, sleep apnea, endothelial function and aortic stiffness. Individual and combined frequencies of these inter-related (i.e., axis) conditions were estimated and associations evaluated.

Results: Household screening of 5230 individuals aged 20 years and above followed by direct physical assessment of 1610 identified 23.7 % of men and 26.7 % of women with type 2 diabetes, 46.2 and 49.5 % of men and women, respectively with obesity and 32.1 and 32.4 % with hypertension. Evaluation of pulse wave velocity, left ventricular mass, endothelial function and sleep apnea identified 22.3, 12.7, 48.6 and 45.2 % of men as having "at risk" values for each condition, respectively. Corresponding numbers in women were 16.0, 17.9, 23.6 and 28.8 %. Cumulatively, 88 % of the population has one or more of these while 50 % have three or more.

Conclusions: The full axis of conditions is high among Mexican Americans in Starr County, Texas. Individual and joint patterns suggest a genesis well before overt disease. Whether they are all mediated by common underlying factors or whether there exist multiple mechanisms remains to be seen.

Keywords: Aortic stiffness; Endothelial function; Hispanic; Hypertension; Left ventricular mass; Mexican American; Obesity; Prevalence; Sleep apnea; Type 2 diabetes.

Rationale: Obstructive sleep apnea is a common disorder associated with increased risk for cardiovascular disease, diabetes, and premature mortality. Although there is strong clinical and epidemiologic evidence supporting the importance of genetic factors in influencing obstructive sleep apnea, its genetic basis is still largely unknown. Prior genetic studies focused on traits defined using the apnea-hypopnea index, which contains limited information on potentially important genetically determined physiologic factors, such as propensity for hypoxemia and respiratory arousability.

Objectives: To define novel obstructive sleep apnea genetic risk loci for obstructive sleep apnea, we conducted genome-wide association studies of quantitative traits in Hispanic/Latino Americans from three cohorts.

Methods: Genome-wide data from as many as 12,558 participants in the Hispanic Community Health Study/Study of Latinos, Multi-Ethnic Study of Atherosclerosis, and Starr County Health Studies population-based cohorts were metaanalyzed for association with the apnea-hypopnea index, average oxygen saturation during sleep, and average respiratory event duration.

Measurements and main results: Two novel loci were identified at genome-level significance (rs11691765, GPR83, P = 1.90 × 10-8 for the apnea-hypopnea index, and rs35424364; C6ORF183/CCDC162P, P = 4.88 × 10-8 for respiratory event duration) and seven additional loci were identified with suggestive significance (P < 5 × 10-7). Secondary sex-stratified analyses also identified one significant and several suggestive associations. Multiple loci overlapped genes with biologic plausibility.

Conclusions: These are the first genome-level significant findings reported for obstructive sleep apnea-related physiologic traits in any population. These findings identify novel associations in inflammatory, hypoxia signaling, and sleep pathways.

Keywords: GWAS; Prader–Willi syndrome; genetics; hypoxia; sleep apnea.

To gain insight into potential regulatory mechanisms through which the effects of variants at four established type 2 diabetes (T2D) susceptibility loci (CDKAL1, CDKN2A-B, IGF2BP2 and KCNQ1) are mediated, we undertook transancestral fine-mapping in 22 086 cases and 42 539 controls of East Asian, European, South Asian, African American and Mexican American descent. Through high-density imputation and conditional analyses, we identified seven distinct association signals at these four loci, each with allelic effects on T2D susceptibility that were homogenous across ancestry groups. By leveraging differences in the structure of linkage disequilibrium between diverse populations, and increased sample size, we localised the variants most likely to drive each distinct association signal. We demonstrated that integration of these genetic fine-mapping data with genomic annotation can highlight potential causal regulatory elements in T2D-relevant tissues. These analyses provide insight into the mechanisms through which T2D association signals are mediated, and suggest future routes to understanding the biology of specific disease susceptibility loci.

We performed genome-wide meta-analysis of lipid traits on three samples of Mexican and Mexican American ancestry comprising 4,383 individuals, and followed up significant and highly suggestive associations in three additional Hispanic samples comprising 7,876 individuals. Genome-wide significant signals were observed in or near CELSR2, ZNF259/APOA5, KANK2/DOCK6 and NCAN/MAU2 for total cholesterol, LPL, ABCA1, ZNF259/APOA5, LIPC and CETP for HDL cholesterol, CELSR2, APOB and NCAN/MAU2 for LDL cholesterol, and GCKR, TRIB1, ZNF259/APOA5 and NCAN/MAU2 for triglycerides. Linkage disequilibrium and conditional analyses indicate that signals observed at ABCA1 and LIPC for HDL cholesterol and NCAN/MAU2 for triglycerides are independent of previously reported lead SNP associations. Analyses of lead SNPs from the European Global Lipids Genetics Consortium (GLGC) dataset in our Hispanic samples show remarkable concordance of direction of effects as well as strong correlation in effect sizes. A meta-analysis of the European GLGC and our Hispanic datasets identified five novel regions reaching genome-wide significance: two for total cholesterol (FN1 and SAMM50), two for HDL cholesterol (LOC100996634 and COPB1) and one for LDL cholesterol (LINC00324/CTC1/PFAS). The top meta-analysis signals were found to be enriched for SNPs associated with gene expression in a tissue-specific fashion, suggesting an enrichment of tissue-specific function in lipid-associated loci.