MUSSEL: Enhanced Bayesian polygenic risk prediction leveraging information across multiple ancestry groups

Polygenic risk scores (PRSs) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in summary statistics from genome-wide association studies (GWASs) across multiple ancestry groups via Bayesian hierarchical modeling and ensemble learning. In our simulation studies and data analyses across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. For example, MUSSEL has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, trait architecture, and linkage disequilibrium reference samples; thus, ultimately a combination of methods may be needed to generate the most robust PRSs across diverse populations.

NIH Toolbox performance of persons with Parkinson's disease according to GBA1 and STN-DBS status

OBJECTIVE

Mutations in the glucocerebrosidase (GBA1) gene and subthalamic nucleus deep brain stimulation (STN-DBS) are independently associated with cognitive dysfunction in persons with Parkinson's disease (PwP). We hypothesized that PwP with both GBA1 mutations and STN-DBS are at greater risk of cognitive dysfunction than PwP with only GBA1 mutations or STN-DBS, or neither. In this study, we determined the pattern of cognitive dysfunction in PwP based on GBA1 mutation status and STN-DBS treatment.

METHODS

PwP who are GBA1 mutation carriers with or without DBS (GBA1+DBS+, GBA1+DBS-), and noncarriers with or without DBS (GBA1-DBS+, GBA1-DBS-) were included. Using the NIH Toolbox, cross-sectional differences in response inhibition, processing speed, and episodic memory were compared using analysis of variance with adjustment for relevant covariates.

RESULTS

Data were available for 9 GBA1+DBS+, 14 GBA1+DBS-, 17 GBA1-DBS+, and 26 GBA1-DBS- PwP. In this cross-sectional study, after adjusting for covariates, we found that performance on the Flanker test (measure of response inhibition) was lower in GBA1+DBS+ PwP compared with GBA1-DBS+ PwP (P = 0.030).

INTERPRETATION

PwP who carry GBA1 mutations and have STN-DBS have greater impaired response inhibition compared with PwP with STN-DBS but without GBA1 mutations. Longitudinal data, including preoperative scores, are required to definitively determine whether GBA1 mutation carriers respond differently to STN-DBS, particularly in the domain of response inhibition.

MUSSEL: Enhanced Bayesian Polygenic Risk Prediction Leveraging Information across Multiple Ancestry Groups

Polygenic risk scores (PRS) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across different populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in the summary statistics from genome-wide association studies (GWAS) across multiple ancestry groups. MUSSEL conducts Bayesian hierarchical modeling under a MUltivariate Spike-and-Slab model for effect-size distribution and incorporates an Ensemble Learning step using super learner to combine information across different tuning parameter settings and ancestry groups. In our simulation studies and data analyses of 16 traits across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. The method, for example, has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African Ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, underlying trait architecture, and the choice of reference samples for LD estimation, and thus ultimately, a combination of methods may be needed to generate the most robust PRS across diverse populations.

The Rosetta Phenotype Harmonization Method Facilitates Finding a Relationship Quantitative Trait Locus for a Complex Cognitive Trait

Genetics researchers increasingly combine data across many sources to increase power and to conduct analyses that cross multiple individual studies. However, there is often a lack of alignment on outcome measures when the same constructs are examined across studies. This inhibits comparison across individual studies and may impact the findings from meta-analysis. Using a well-characterized genotypic (brain-derived neurotrophic factor: BDNF) and phenotypic constructs (working memory and reading comprehension), we employ an approach called Rosetta, which allows for the simultaneous examination of primary studies that employ related but incompletely overlapping data. We examined four studies of BDNF, working memory, and reading comprehension with a combined sample size of 1711 participants. Although the correlation between working memory and reading comprehension over all participants was high, as expected (ρ = 0.45), the correlation between working memory and reading comprehension was attenuated in the BDNF Met/Met genotype group (ρ = 0.18, n.s.) but not in the Val/Val (ρ = 0.44) or Val/Met (ρ = 0.41) groups. These findings indicate that Met/Met carriers may be a unique and robustly defined subgroup in terms of memory and reading comprehension. This study demonstrates the utility of the Rosetta method when examining complex phenotypes across multiple studies, including psychiatric genetic studies, as shown here, and also for the mega-analysis of cohorts generally.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

Trans-ancestry, Bayesian meta-analysis discovers 20 novel risk loci for inflammatory bowel disease in an African American, East Asian and European cohort

Inflammatory bowel disease (IBD) is an immune-mediated chronic intestinal disorder with major phenotypes: ulcerative colitis (UC) and Crohn's disease (CD). Multiple studies have identified over 240 IBD susceptibility loci. However, most studies have centered on European (EUR) and East Asian (EAS) populations. The prevalence of IBD in non-EUR, including African Americans (AAs), has risen in recent years. Here we present the first attempt to identify loci in AAs using a trans-ancestry Bayesian approach (MANTRA) accounting for heterogeneity between diverse ancestries while allowing for the similarity between closely related populations. We meta-analyzed genome-wide association studies (GWAS) and Immunochip data from a 2015 EUR meta-analysis of 38 155 IBD cases and 48 485 controls and EAS Immunochip study of 2824 IBD cases and 3719 controls, and our recent AA IBD GWAS of 2345 cases and 5002 controls. Across the major IBD phenotypes, we found significant evidence for 92% of 205 loci lead SNPs from the 2015 meta-analysis, but also for three IBD loci only established in latter studies. We detected 20 novel loci, all containing immunity-related genes or genes with other evidence for IBD or immune-mediated disease relevance: PLEKHG5;TNFSFR25 (encoding death receptor 3, receptor for TNFSF15 gene product TL1A), XKR6, ELMO1, BC021024;PI4KB;PSMD4 and APLP1 for IBD; AUTS2, XKR6, OSER1, TET2;AK094561, BCAP29 and APLP1 for CD; and GABBR1;MOG, DQ570892, SPDEF;ILRUN, SMARCE1;CCR7;KRT222;KRT24;KRT25, ANKS1A;TCP11, IL7, LRRC18;WDFY4, XKR6 and TNFSF4 for UC. Our study highlights the value of combining low-powered genomic studies from understudied populations of diverse ancestral backgrounds together with a high-powered study to enable novel locus discovery, including potentially important therapeutic IBD gene targets.

Common genetic risk factors in ASD and ADHD co-occurring families

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

Large-scale genome-wide association study of coronary artery disease in genetically diverse populations

We report a genome-wide association study (GWAS) of coronary artery disease (CAD) incorporating nearly a quarter of a million cases, in which existing studies are integrated with data from cohorts of white, Black and Hispanic individuals from the Million Veteran Program. We document near equivalent heritability of CAD across multiple ancestral groups, identify 95 novel loci, including nine on the X chromosome, detect eight loci of genome-wide significance in Black and Hispanic individuals, and demonstrate that two common haplotypes at the 9p21 locus are responsible for risk stratification in all populations except those of African origin, in which these haplotypes are virtually absent. Moreover, in the largest GWAS for angiographically derived coronary atherosclerosis performed to date, we find 15 loci of genome-wide significance that robustly overlap with established loci for clinical CAD. Phenome-wide association analyses of novel loci and polygenic risk scores (PRSs) augment signals related to insulin resistance, extend pleiotropic associations of these loci to include smoking and family history, and precisely document the markedly reduced transferability of existing PRSs to Black individuals. Downstream integrative analyses reinforce the critical roles of vascular endothelial, fibroblast, and smooth muscle cells in CAD susceptibility, but also point to a shared biology between atherosclerosis and oncogenesis. This study highlights the value of diverse populations in further characterizing the genetic architecture of CAD.

MicroRNA and MicroRNA-Target Variants Associated with Autism Spectrum Disorder and Related Disorders

Autism spectrum disorder (ASD) is a childhood neurodevelopmental disorder with a complex and heterogeneous genetic etiology. MicroRNA (miRNA), a class of small non-coding RNAs, could regulate ASD risk genes post-transcriptionally and affect broad molecular pathways related to ASD and associated disorders. Using whole-genome sequencing, we analyzed 272 samples in 73 families in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Families with at least one ASD patient were recruited and were further assessed for language impairment, reading impairment, and other associated phenotypes. A total of 5104 miRNA variants and 1,181,148 3′ untranslated region (3′ UTR) variants were identified in the dataset. After applying several filtering criteria, including population allele frequency, brain expression, miRNA functional regions, and inheritance patterns, we identified high-confidence variants in five brain-expressed miRNAs (targeting 326 genes) and 3′ UTR miRNA target regions of 152 genes. Some genes, such as SCP2 and UCGC, were identified in multiple families. Using Gene Ontology overrepresentation analysis and protein–protein interaction network analysis, we identified clusters of genes and pathways that are important for neurodevelopment. The miRNAs and miRNA target genes identified in this study are potentially involved in neurodevelopmental disorders and should be considered for further functional studies.

Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program

BACKGROUND

The Centers for Disease Control and Prevention contracted with laboratories to sequence the SARS-CoV-2 genome from positive samples across the United States to enable public health officials to investigate the impact of variants on disease severity as well as the effectiveness of vaccines and treatment. Herein we present the initial results correlating RT-PCR quality control metrics with sample collection and sequencing methods from full SARS-CoV-2 viral genomic sequencing of 24,441 positive patient samples between April and June 2021.

METHODS

RT-PCR confirmed (N Gene Ct value < 30) positive patient samples, with nucleic acid extracted from saliva, nasopharyngeal and oropharyngeal swabs were selected for viral whole genome SARS-CoV-2 sequencing. Sequencing was performed using Illumina COVIDSeq™ protocol on either the NextSeq550 or NovaSeq6000 systems. Informatic variant calling, and lineage analysis were performed using DRAGEN COVID Lineage applications on Illumina's Basespace cloud analytical system. All sequence data and variant calls were uploaded to NCBI and GISAID.

RESULTS

An association was observed between higher sequencing coverage, quality, and samples with a lower Ct value, with < 27 being optimal, across both sequencing platforms and sample collection methods. Both nasopharyngeal swabs and saliva samples were found to be optimal samples of choice for SARS-CoV-2 surveillance sequencing studies, both in terms of strain identification and sequencing depth of coverage, with NovaSeq 6000 providing higher coverage than the NextSeq 550. The most frequent variants identified were the B.1.617.2 Delta (India) and P.1 Gamma (Brazil) variants in the samples sequenced between April 2021 and June 2021. At the time of submission, the most common variant > 99% of positives sequenced was Omicron.

CONCLUSION

These initial analyses highlight the importance of sequencing platform, sample collection methods, and RT-PCR Ct values in guiding surveillance efforts. These surveillance studies evaluating genetic changes of SARS-CoV-2 have been identified as critical by the CDC that can affect many aspects of public health including transmission, disease severity, diagnostics, therapeutics, and vaccines.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Centers for Mendelian Genomics: A decade of facilitating gene discovery

PURPOSE

Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration.

METHODS

Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution.

RESULTS

We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher.

CONCLUSION

The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study

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 ).

Genetic discovery and risk characterization in type 2 diabetes across diverse populations

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⁻⁸), which replicated in independent studies of African ancestry (p = 6.26 × 10⁻²³). We identified a multiethnic risk variant in the BACE2 gene (rs13052926, RAF = 14.1%, OR = 1.08, p = 5.75 × 10⁻⁹), which also replicated in independent studies (p = 3.45 × 10⁻⁴). We also observed a significant difference in the performance of a multiethnic genetic risk score (GRS) across population groups (p_{heterogeneity} = 3.85 × 10⁻²⁰). 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.

Mapping and characterization of structural variation in 17,795 human genomes

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Translational derepression of Elavl4 isoforms at their alternative 5' UTRs determines neuronal development

Neurodevelopment requires precise regulation of gene expression, including post-transcriptional regulatory events such as alternative splicing and mRNA translation. However, translational regulation of specific isoforms during neurodevelopment and the mechanisms behind it remain unknown. Using RNA-seq analysis of mouse neocortical polysomes, here we report translationally repressed and derepressed mRNA isoforms during neocortical neurogenesis whose orthologs include risk genes for neurodevelopmental disorders. We demonstrate that the translation of distinct mRNA isoforms of the RNA binding protein (RBP), Elavl4, in radial glia progenitors and early neurons depends on its alternative 5' UTRs. Furthermore, 5' UTR-driven Elavl4 isoform-specific translation depends on upstream control by another RBP, Celf1. Celf1 regulation of Elavl4 translation dictates development of glutamatergic neurons. Our findings reveal a dynamic interplay between distinct RBPs and alternative 5' UTRs in neuronal development and underscore the risk of post-transcriptional dysregulation in co-occurring neurodevelopmental disorders.

Minority-centric meta-analyses of blood lipid levels identify novel loci in the Population Architecture using Genomics and Epidemiology (PAGE) study

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.

Blood lipid levels are closely linked to cardio-metabolic diseases, and genetic factors play an important role in their metabolism and regulation. Although over 400 loci have been identified through genetic association studies, the genetic architecture of lipid levels is not fully characterized. The lack of representation of diverse populations in previous studies resulted in a large gap in understanding the genetic background of lipid traits between European and minority populations, including African Americans, Hispanics, Hawaiians, and Native Americans. In our current analyses which included ancestrally diverse populations, we identified nine novel loci, one novel gene, and minority-specific independent signals at eight known loci, and pinpointed potential functional variants at two known loci. We further observed smaller effect sizes of reported lipids-associated loci in African Americans and Hispanics than those in Europeans, and better performance of polygenic risk scores using minority-specific instead of European-derived effect sizes when estimating genetic predisposition in minority populations. Our findings showed the benefits of including multi-ethnic studies in identification and refinement of lipids-associated loci, which will help to reduce the existing disparities and to pave the road to precision medicine.

A phenome-wide association study (PheWAS) in the Population Architecture using Genomics and Epidemiology (PAGE) study reveals potential pleiotropy in African Americans

We performed a hypothesis-generating phenome-wide association study (PheWAS) to identify and characterize cross-phenotype associations, where one SNP is associated with two or more phenotypes, between thousands of genetic variants assayed on the Metabochip and hundreds of phenotypes in 5,897 African Americans as part of the Population Architecture using Genomics and Epidemiology (PAGE) I study. The PAGE I study was a National Human Genome Research Institute-funded collaboration of four study sites accessing diverse epidemiologic studies genotyped on the Metabochip, a custom genotyping chip that has dense coverage of regions in the genome previously associated with cardio-metabolic traits and outcomes in mostly European-descent populations. Here we focus on identifying novel phenome-genome relationships, where SNPs are associated with more than one phenotype. To do this, we performed a PheWAS, testing each SNP on the Metabochip for an association with up to 273 phenotypes in the participating PAGE I study sites. We identified 133 putative pleiotropic variants, defined as SNPs associated at an empirically derived p-value threshold of p<0.01 in two or more PAGE study sites for two or more phenotype classes. We further annotated these PheWAS-identified variants using publicly available functional data and local genetic ancestry. Amongst our novel findings is SPARC rs4958487, associated with increased glucose levels and hypertension. SPARC has been implicated in the pathogenesis of diabetes and is also known to have a potential role in fibrosis, a common consequence of multiple conditions including hypertension. The SPARC example and others highlight the potential that PheWAS approaches have in improving our understanding of complex disease architecture by identifying novel relationships between genetic variants and an array of common human phenotypes.

Forming Big Datasets through Latent Class Concatenation of Imperfectly Matched Databases Features

Informatics researchers often need to combine data from many different sources to increase statistical power and study subtle or complicated effects. Perfect overlap of measurements across academic studies is rare since virtually every dataset is collected for a unique purpose and without coordination across parties not-at-hand (i.e., informatics researchers in the future). Thus, incomplete concordance of measurements across datasets poses a major challenge for researchers seeking to combine public databases. In any given field, some measurements are fairly standard, but every organization collecting data makes unique decisions on instruments, protocols, and methods of processing the data. This typically denies literal concatenation of the raw data since constituent cohorts do not have the same measurements (i.e., columns of data). When measurements across datasets are similar prima facie, there is a desire to combine the data to increase power, but mixing non-identical measurements could greatly reduce the sensitivity of the downstream analysis. Here, we discuss a statistical method that is applicable when certain patterns of missing data are found; namely, it is possible to combine datasets that measure the same underlying constructs (or latent traits) when there is only partial overlap of measurements across the constituent datasets. Our method, ROSETTA empirically derives a set of common latent trait metrics for each related measurement domain using a novel variation of factor analysis to ensure equivalence across the constituent datasets. The advantage of combining datasets this way is the simplicity, statistical power, and modeling flexibility of a single joint analysis of all the data. Three simulation studies show the performance of ROSETTA on datasets with only partially overlapping measurements (i.e., systematically missing information), benchmarked to a condition of perfectly overlapped data (i.e., full information). The first study examined a range of correlations, while the second study was modeled after the observed correlations in a well-characterized clinical, behavioral cohort. Both studies consistently show significant correlations >0.94, often >0.96, indicating the robustness of the method and validating the general approach. The third study varied within and between domain correlations and compared ROSETTA to multiple imputation and meta-analysis as two commonly used methods that ostensibly solve the same data integration problem. We provide one alternative to meta-analysis and multiple imputation by developing a method that statistically equates similar but distinct manifest metrics into a set of empirically derived metrics that can be used for analysis across all datasets.

Genetic analyses of diverse populations improves discovery for complex traits

Genome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development and clinical guidelines. However, the majority of discovery efforts are based on data from populations of European ancestry^1-3. In light of the differential genetic architecture that is known to exist between populations, bias in representation can exacerbate existing disease and healthcare disparities. Critical variants may be missed if they have a low frequency or are completely absent in European populations, especially as the field shifts its attention towards rare variants, which are more likely to be population-specific^4-10. Additionally, effect sizes and their derived risk prediction scores derived in one population may not accurately extrapolate to other populations^11,12. Here we demonstrate the value of diverse, multi-ethnic participants in large-scale genomic studies. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioural phenotypes in 49,839 non-European individuals. Using strategies tailored for analysis of multi-ethnic and admixed populations, we describe a framework for analysing diverse populations, identify 27 novel loci and 38 secondary signals at known loci, as well as replicate 1,444 GWAS catalogue associations across these traits. Our data show evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts and insights into clinical implications. In the United States-where minority populations have a disproportionately higher burden of chronic conditions¹³-the lack of representation of diverse populations in genetic research will result in inequitable access to precision medicine for those with the highest burden of disease. We strongly advocate for continued, large genome-wide efforts in diverse populations to maximize genetic discovery and reduce health disparities.

Genetics of Chronic Kidney Disease Stages Across Ancestries: The PAGE Study

Background

Chronic kidney disease (CKD) is common and disproportionally burdens United States ethnic minorities. Its genetic determinants may differ by disease severity and clinical stages. To uncover genetic factors associated CKD severity among high-risk ethnic groups, we performed genome-wide association studies (GWAS) in diverse populations within the Population Architecture using Genomics and Epidemiology (PAGE) study.

Methods

We assembled multi-ethnic genome-wide imputed data on CKD non-overlapping cases [4,150 mild to moderate CKD, 1,105 end-stage kidney disease (ESKD)] and non-CKD controls for up to 41,041 PAGE participants (African Americans, Hispanics/Latinos, East Asian, Native Hawaiian, and American Indians). We implemented a generalized estimating equation approach for GWAS using ancestry combined data while adjusting for age, sex, principal components, study, and ethnicity.

Results

The GWAS identified a novel genome-wide associated locus for mild to moderate CKD nearby NMT2 (rs10906850, p = 3.7 × 10^-8) that replicated in the United Kingdom Biobank white British (p = 0.008). Several variants at the APOL1 locus were associated with ESKD including the APOL1 G1 rs73885319 (p = 1.2 × 10^-9). There was no overlap among associated loci for CKD and ESKD traits, even at the previously reported APOL1 locus (p = 0.76 for CKD). Several additional loci were associated with CKD or ESKD at p-values below the genome-wide threshold. These loci were often driven by variants more common in non-European ancestry.

Conclusion

Our genetic study identified a novel association at NMT2 for CKD and showed for the first time strong associations of the APOL1 variants with ESKD across multi-ethnic populations. Our findings suggest differences in genetic effects across CKD severity and provide information for study design of genetic studies of CKD in diverse populations.

Insights into genetics, human biology and disease gleaned from family based genomic studies

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

The genetic underpinnings of variation in ages at menarche and natural menopause among women from the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) Study: A trans-ethnic meta-analysis

Current knowledge of the genetic architecture of key reproductive events across the female life course is largely based on association studies of European descent women. The relevance of known loci for age at menarche (AAM) and age at natural menopause (ANM) in diverse populations remains unclear. We investigated 32 AAM and 14 ANM previously-identified loci and sought to identify novel loci in a trans-ethnic array-wide study of 196,483 SNPs on the MetaboChip (Illumina, Inc.). A total of 45,364 women of diverse ancestries (African, Hispanic/Latina, Asian American and American Indian/Alaskan Native) in the Population Architecture using Genomics and Epidemiology (PAGE) Study were included in cross-sectional analyses of AAM and ANM. Within each study we conducted a linear regression of SNP associations with self-reported or medical record-derived AAM or ANM (in years), adjusting for birth year, population stratification, and center/region, as appropriate, and meta-analyzed results across studies using multiple meta-analytic techniques. For both AAM and ANM, we observed more directionally consistent associations with the previously reported risk alleles than expected by chance (p-values_binomial≤0.01). Eight densely genotyped reproductive loci generalized significantly to at least one non-European population. We identified one trans-ethnic array-wide SNP association with AAM and two significant associations with ANM, which have not been described previously. Additionally, we observed evidence of independent secondary signals at three of six AAM trans-ethnic loci. Our findings support the transferability of reproductive trait loci discovered in European women to women of other race/ethnicities and indicate the presence of additional trans-ethnic associations both at both novel and established loci. These findings suggest the benefit of including diverse populations in future studies of the genetic architecture of female growth and development.

Validation of a microRNA target site polymorphism in H3F3B that is potentially associated with a broad schizophrenia phenotype

Despite much progress, few genetic findings for schizophrenia have been assessed by functional validation experiments at the molecular level. We previously reported evidence for genetic linkage of broadly defined schizophrenia to chromosome 17q25 in a sample of 24 multiplex families. 2,002 SNPs under this linkage peak were analyzed for evidence of linkage disequilibrium using the posterior probability of linkage (PPL) framework. SNP rs1060120 produced the strongest evidence for association, with a PPLD|L score of 0.21. This SNP is located within the 3'UTR of the histone gene H3F3B and colocalizes with potential gene target miR-616. A custom miRNA target prediction program predicted that the binding of miR-616 to H3F3B transcripts would be altered by the allelic variants of rs1060120. We used dual luciferase assays to experimentally validate this interaction. The rs1060120 A allele significantly reduced luciferase expression, indicating a stronger interaction with miR-616 than the G allele (p = 0.000412). These results provide functional validation that this SNP could alter schizophrenia epigenetic mechanisms thereby contributing to schizophrenia-related disease risk.

Trans-ethnic analysis of metabochip data identifies two new loci associated with BMI

OBJECTIVE

Body mass index (BMI) is commonly used to assess obesity, which is associated with numerous diseases and negative health outcomes. BMI has been shown to be a heritable, polygenic trait, with close to 100 loci previously identified and replicated in multiple populations. We aim to replicate known BMI loci and identify novel associations in a trans-ethnic study population.

SUBJECTS

Using eligible participants from the Population Architecture using Genomics and Epidemiology consortium, we conducted a trans-ethnic meta-analysis of 102 514 African Americans, Hispanics, Asian/Native Hawaiian, Native Americans and European Americans. Participants were genotyped on over 200 000 SNPs on the Illumina Metabochip custom array, or imputed into the 1000 Genomes Project (Phase I). Linear regression of the natural log of BMI, adjusting for age, sex, study site (if applicable), and ancestry principal components, was conducted for each race/ethnicity within each study cohort. Race/ethnicity-specific, and combined meta-analyses used fixed-effects models.

RESULTS

We replicated 15 of 21 BMI loci included on the Metabochip, and identified two novel BMI loci at 1q41 (rs2820436) and 2q31.1 (rs10930502) at the Metabochip-wide significance threshold (P<2.5 × 10-7). Bioinformatic functional investigation of SNPs at these loci suggests a possible impact on pathways that regulate metabolism and adipose tissue.

CONCLUSION

Conducting studies in genetically diverse populations continues to be a valuable strategy for replicating known loci and uncovering novel BMI associations.

Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium

AIMS/HYPOTHESIS

Elevated levels of fasting glucose and fasting insulin in non-diabetic individuals are markers of dysregulation of glucose metabolism and are strong risk factors for type 2 diabetes. Genome-wide association studies have discovered over 50 SNPs associated with these traits. Most of these loci were discovered in European populations and have not been tested in a well-powered multi-ethnic study. We hypothesised that a large, ancestrally diverse, fine-mapping genetic study of glycaemic traits would identify novel and population-specific associations that were previously undetectable by European-centric studies.

METHODS

A multiethnic study of up to 26,760 unrelated individuals without diabetes, of predominantly Hispanic/Latino and African ancestries, were genotyped using the Metabochip. Transethnic meta-analysis of racial/ethnic-specific linear regression analyses were performed for fasting glucose and fasting insulin. We attempted to replicate 39 fasting glucose and 17 fasting insulin loci. Genetic fine-mapping was performed through sequential conditional analyses in 15 regions that included both the initially reported SNP association(s) and denser coverage of SNP markers. In addition, Metabochip-wide analyses were performed to discover novel fasting glucose and fasting insulin loci. The most significant SNP associations were further examined using bioinformatic functional annotation.

RESULTS

Previously reported SNP associations were significantly replicated (p ≤ 0.05) in 31/39 fasting glucose loci and 14/17 fasting insulin loci. Eleven glycaemic trait loci were refined to a smaller list of potentially causal variants through transethnic meta-analysis. Stepwise conditional analysis identified two loci with independent secondary signals (G6PC2-rs477224 and GCK-rs2908290), which had not previously been reported. Population-specific conditional analyses identified an independent signal in G6PC2 tagged by the rare variant rs77719485 in African ancestry. Further Metabochip-wide analysis uncovered one novel fasting insulin locus at SLC17A2-rs75862513.

CONCLUSIONS/INTERPRETATION

These findings suggest that while glycaemic trait loci often have generalisable effects across the studied populations, transethnic genetic studies help to prioritise likely functional SNPs, identify novel associations that may be population-specific and in turn have the potential to influence screening efforts or therapeutic discoveries.

DATA AVAILABILITY

The summary statistics from each of the ancestry-specific and transethnic (combined ancestry) results can be found under the PAGE study on dbGaP here: https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000356.v1.p1.

Correction to: Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium

The authors regret that Yinchang Lu's name incorrectly included a middle initial in the author list. The details given in this erratum are correct.

Fine-mapping of lipid regions in global populations discovers ethnic-specific signals and refines previously identified lipid loci

Genome-wide association studies have identified over 150 loci associated with lipid traits, however, no large-scale studies exist for Hispanics and other minority populations. Additionally, the genetic architecture of lipid-influencing loci remains largely unknown. We performed one of the most racially/ethnically diverse fine-mapping genetic studies of HDL-C, LDL-C, and triglycerides to-date using SNPs on the MetaboChip array on 54,119 individuals: 21,304 African Americans, 19,829 Hispanic Americans, 12,456 Asians, and 530 American Indians. The majority of signals found in these groups generalize to European Americans. While we uncovered signals unique to racial/ethnic populations, we also observed systematically consistent lipid associations across these groups. In African Americans, we identified three novel signals associated with HDL-C (LPL, APOA5, LCAT) and two associated with LDL-C (ABCG8, DHODH). In addition, using this population, we refined the location for 16 out of the 58 known MetaboChip lipid loci. These results can guide tailored screening efforts, reveal population-specific responses to lipid-lowering medications, and aid in the development of new targeted drug therapies.

Trans-ethnic fine-mapping of genetic loci for body mass index in the diverse ancestral populations of the Population Architecture using Genomics and Epidemiology (PAGE) Study reveals evidence for multiple signals at established loci

Most body mass index (BMI) genetic loci have been identified in studies of primarily European ancestries. The effect of these loci in other racial/ethnic groups is less clear. Thus, we aimed to characterize the generalizability of 170 established BMI variants, or their proxies, to diverse US populations and trans-ethnically fine-map 36 BMI loci using a sample of >102,000 adults of African, Hispanic/Latino, Asian, European and American Indian/Alaskan Native descent from the Population Architecture using Genomics and Epidemiology Study. We performed linear regression of the natural log of BMI (18.5-70 kg/m²) on the additive single nucleotide polymorphisms (SNPs) at BMI loci on the MetaboChip (Illumina, Inc.), adjusting for age, sex, population stratification, study site, or relatedness. We then performed fixed-effect meta-analyses and a Bayesian trans-ethnic meta-analysis to empirically cluster by allele frequency differences. Finally, we approximated conditional and joint associations to test for the presence of secondary signals. We noted directional consistency with the previously reported risk alleles beyond what would have been expected by chance (binomial p < 0.05). Nearly, a quarter of the previously described BMI index SNPs and 29 of 36 densely-genotyped BMI loci on the MetaboChip replicated/generalized in trans-ethnic analyses. We observed multiple signals at nine loci, including the description of seven loci with novel multiple signals. This study supports the generalization of most common genetic loci to diverse ancestral populations and emphasizes the importance of dense multiethnic genomic data in refining the functional variation at genetic loci of interest and describing several loci with multiple underlying genetic variants.

Fine mapping of QT interval regions in global populations refines previously identified QT interval loci and identifies signals unique to African and Hispanic descent populations

BACKGROUND

The electrocardiographically measured QT interval (QT) is heritable and its prolongation is an established risk factor for several cardiovascular diseases. Yet, most QT genetic studies have been performed in European ancestral populations, possibly reducing their global relevance.

OBJECTIVE

To leverage diversity and improve biological insight, we fine mapped 16 of the 35 previously identified QT loci (46%) in populations of African American (n = 12,410) and Hispanic/Latino (n = 14,837) ancestry.

METHODS

Racial/ethnic-specific multiple linear regression analyses adjusted for heart rate and clinical covariates were examined separately and in combination after inverse-variance weighted trans-ethnic meta-analysis.

RESULTS

The 16 fine-mapped QT loci included on the Illumina Metabochip represented 21 independent signals, of which 16 (76%) were significantly (P-value≤9.1×10-5) associated with QT. Through sequential conditional analysis we also identified three trans-ethnic novel SNPs at ATP1B1, SCN5A-SCN10A, and KCNQ1 and three Hispanic/Latino-specific novel SNPs at NOS1AP and SCN5A-SCN10A (two novel SNPs) with evidence of associations with QT independent of previous identified GWAS lead SNPs. Linkage disequilibrium patterns helped to narrow the region likely to contain the functional variants at several loci, including NOS1AP, USP50-TRPM7, and PRKCA, although intervals surrounding SLC35F1-PLN and CNOT1 remained broad in size (>100 kb). Finally, bioinformatics-based functional characterization suggested a regulatory function in cardiac tissues for the majority of independent signals that generalized and the novel SNPs.

CONCLUSION

Our findings suggest that a majority of identified SNPs implicate gene regulatory dysfunction in QT prolongation, that the same loci influence variation in QT across global populations, and that additional, novel, population-specific QT signals exist.

Generalization and fine mapping of European ancestry-based central adiposity variants in African ancestry populations

BACKGROUND/OBJECTIVES

Central adiposity measures such as waist circumference (WC) and waist-to-hip ratio (WHR) are associated with cardiometabolic disorders independently of body mass index (BMI) and are gaining clinically utility. Several studies report genetic variants associated with central adiposity, but most utilize only European ancestry populations. Understanding whether the genetic associations discovered among mainly European descendants are shared with African ancestry populations will help elucidate the biological underpinnings of abdominal fat deposition.

SUBJECTS/METHODS

To identify the underlying functional genetic determinants of body fat distribution, we conducted an array-wide association meta-analysis among persons of African ancestry across seven studies/consortia participating in the Population Architecture using Genomics and Epidemiology (PAGE) consortium. We used the Metabochip array, designed for fine-mapping cardiovascular-associated loci, to explore novel array-wide associations with WC and WHR among 15 945 African descendants using all and sex-stratified groups. We further interrogated 17 known WHR regions for African ancestry-specific variants.

RESULTS

Of the 17 WHR loci, eight single-nucleotide polymorphisms (SNPs) located in four loci were replicated in the sex-combined or sex-stratified meta-analyses. Two of these eight independently associated with WHR after conditioning on the known variant in European descendants (rs12096179 in TBX15-WARS2 and rs2059092 in ADAMTS9). In the fine-mapping assessment, the putative functional region was reduced across all four loci but to varying degrees (average 40% drop in number of putative SNPs and 20% drop in genomic region). Similar to previous studies, the significant SNPs in the female-stratified analysis were stronger than the significant SNPs from the sex-combined analysis. No novel associations were detected in the array-wide analyses.

CONCLUSIONS

Of 17 previously identified loci, four loci replicated in the African ancestry populations of this study. Utilizing different linkage disequilibrium patterns observed between European and African ancestries, we narrowed the suggestive region containing causative variants for all four loci.

Strategies for Enriching Variant Coverage in Candidate Disease Loci on a Multiethnic Genotyping Array

Investigating genetic architecture of complex traits in ancestrally diverse populations is imperative to understand the etiology of disease. However, the current paucity of genetic research in people of African and Latin American ancestry, Hispanic and indigenous peoples in the United States is likely to exacerbate existing health disparities for many common diseases. The Population Architecture using Genomics and Epidemiology, Phase II (PAGE II), Study was initiated in 2013 by the National Human Genome Research Institute to expand our understanding of complex trait loci in ethnically diverse and well characterized study populations. To meet this goal, the Multi-Ethnic Genotyping Array (MEGA) was designed to substantially improve fine-mapping and functional discovery by increasing variant coverage across multiple ethnicities at known loci for metabolic, cardiovascular, renal, inflammatory, anthropometric, and a variety of lifestyle traits. Studying the frequency distribution of clinically relevant mutations, putative risk alleles, and known functional variants across multiple populations will provide important insight into the genetic architecture of complex diseases and facilitate the discovery of novel, sometimes population-specific, disease associations. DNA samples from 51,650 self-identified African ancestry (17,328), Hispanic/Latino (22,379), Asian/Pacific Islander (8,640), and American Indian (653) and an additional 2,650 participants of either South Asian or European ancestry, and other reference panels have been genotyped on MEGA by PAGE II. MEGA was designed as a new resource for studying ancestrally diverse populations. Here, we describe the methodology for selecting trait-specific content for use in multi-ethnic populations and how enriching MEGA for this content may contribute to deeper biological understanding of the genetic etiology of complex disease.

Behavioral and Molecular Genetics of Reading-Related AM and FM Detection Thresholds

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h ²  = 0.20) and FM (h ²  = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.

Variant Discovery and Fine Mapping of Genetic Loci Associated with Blood Pressure Traits in Hispanics and African Americans

Despite the substantial burden of hypertension in US minority populations, few genetic studies of blood pressure have been conducted in Hispanics and African Americans, and it is unclear whether many of the established loci identified in European-descent populations contribute to blood pressure variation in non-European descent populations. Using the Metabochip array, we sought to characterize the genetic architecture of previously identified blood pressure loci, and identify novel cardiometabolic variants related to systolic and diastolic blood pressure in a multi-ethnic US population including Hispanics (n = 19,706) and African Americans (n = 18,744). Several known blood pressure loci replicated in African Americans and Hispanics. Fourteen variants in three loci (KCNK3, FGF5, ATXN2-SH2B3) were significantly associated with blood pressure in Hispanics. The most significant diastolic blood pressure variant identified in our analysis, rs2586886/KCNK3 (P = 5.2 x 10⁻⁹), also replicated in independent Hispanic and European-descent samples. African American and trans-ethnic meta-analysis data identified novel variants in the FGF5, ULK4 and HOXA-EVX1 loci, which have not been previously associated with blood pressure traits. Our identification and independent replication of variants in KCNK3, a gene implicated in primary hyperaldosteronism, as well as a variant in HOTTIP (HOXA-EVX1) suggest that further work to clarify the roles of these genes may be warranted. Overall, our findings suggest that loci identified in European descent populations also contribute to blood pressure variation in diverse populations including Hispanics and African Americans—populations that are understudied for hypertension genetic risk factors.

Fine-mapping, novel loci identification, and SNP association transferability in a genome-wide association study of QRS duration in African Americans

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

Predictive cytokine biomarkers of clinical response to glatiramer acetate therapy in multiple sclerosis

A prospective study of 62 patients with relapsing-remitting multiple sclerosis (RRMS) treated with Glatiramer acetate (GA) was conducted to evaluate the value of baseline and treatment-modulated cytokines in predicting the clinical response to the drug after 2years of therapy. There were 32 responders and 30 non-responders. GA upregulated Th2/regulatory cytokines and inhibited Th1 cytokines in sera or PBMC supernatants 3 and 6months into treatment. We found two prognostic models with clinical utility. A model based on IL-18 at baseline, the change in TNFa from baseline to 3months, the change in IL-4 from baseline to 6months, and the change in the log of the ratio of TNFa/IL-4 from baseline to 6months had an area under the curve (AUC) of 0.80. A high IL-18 level at baseline and a reduction of TNF-alpha over time are associated with a response to GA. Although the study identified predictive biomarkers of clinical response to GA, the results will need to be validated in other data sets.

Cefoxitin Plasma and Subcutaneous Adipose Tissue Concentration in Patients Undergoing Sleeve Gastrectomy

PURPOSE

Antibiotic dosing in obese surgical patients has not been adequately evaluated. The objective of this study was to identify whether currently prescribed doses of cefoxitin achieve adequate and sustained plasma and tissue concentrations in obese patients undergoing sleeve gastrectomy.

METHODS

A prospective evaluation of plasma and tissue cefoxitin concentrations in patients undergoing sleeve gastrectomy was performed. On the day of the surgical procedure, venous blood samples (5 mL) were collected just before cefoxitin administration and then at 5, 30, 60, 120, and 240 minutes after dose administration. In addition, subcutaneous adipose tissue was collected from the surgical site at the time of surgical incision and at closure. Cefoxitin concentrations in the collected samples were quantified by using an HPLC-ultraviolet method. A standard noncompartmental analysis was performed for each individual cefoxitin plasma concentration-time profile. In addition, the ratio of tissue to plasma concentration was calculated for all patients.

FINDINGS

Plasma and tissue pharmacokinetics of cefoxitin were evaluated in 6 patients undergoing sleeve gastrectomy. The mean age and BMI were 48.7 (6.2) years and 42.8 (7.1) kg/m(2), respectively. At the time of surgical closure, subcutaneous adipose tissue concentrations of cefoxitin were subtherapeutic (<8 µg/mL) in all evaluated patients.

IMPLICATIONS

Current dosing strategies for cefoxitin in obese surgical patients may be inadequate, and there is an urgent need to define the appropriate dosage.

New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10⁻⁸). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Prospective associations of coronary heart disease loci in African Americans using the MetaboChip: the PAGE study

BACKGROUND

Coronary heart disease (CHD) is a leading cause of morbidity and mortality in African Americans. However, there is a paucity of studies assessing genetic determinants of CHD in African Americans. We examined the association of published variants in CHD loci with incident CHD, attempted to fine map these loci, and characterize novel variants influencing CHD risk in African Americans.

METHODS AND RESULTS

Up to 8,201 African Americans (including 546 first CHD events) were genotyped using the MetaboChip array in the Atherosclerosis Risk in Communities (ARIC) study and Women's Health Initiative (WHI). We tested associations using Cox proportional hazard models in sex- and study-stratified analyses and combined results using meta-analysis. Among 44 validated CHD loci available in the array, we replicated and fine-mapped the SORT1 locus, and showed same direction of effects as reported in studies of individuals of European ancestry for SNPs in 22 additional published loci. We also identified a SNP achieving array wide significance (MYC: rs2070583, allele frequency 0.02, P = 8.1 × 10(-8)), but the association did not replicate in an additional 8,059 African Americans (577 events) from the WHI, HealthABC and GeneSTAR studies, and in a meta-analysis of 5 cohort studies of European ancestry (24,024 individuals including 1,570 cases of MI and 2,406 cases of CHD) from the CHARGE Consortium.

CONCLUSIONS

Our findings suggest that some CHD loci previously identified in individuals of European ancestry may be relevant to incident CHD in African Americans.

Differential behavior patterns in cynomolgus monkey Macaca fascicularis in home cage in response to human gaze

BACKGROUND

Non-human primates, when encountering human beings, show wariness and alertness. These behaviors differ when there is direct human gaze vs. when human averts his gaze.

METHODS

We observed cynomolgus monkey in their home cage and studied their behaviors in response to human gaze. Four behaviors were analyzed: opening mouth, staring at observer, agitated activity, and approaching observer.

RESULTS

Three behaviors appeared to be sensitive to human gaze between when the human observer gazed at the monkey and when the human observer looked away. Individual animals also displayed subpatterns of responses to human gaze.

CONCLUSIONS

These results indicate that, even in their home cage, monkeys display a heightened level of awareness when gazed upon by a human observer, suggesting that human gaze may elicit emotional reactions. Further, under the human gaze, distinct behavioral subpatterns were apparent within the monkey cohort in our study, indicative of subgroups within the cohort.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Pleiotropic associations of risk variants identified for other cancers with lung cancer risk: the PAGE and TRICL consortia

BACKGROUND

Genome-wide association studies have identified hundreds of genetic variants associated with specific cancers. A few of these risk regions have been associated with more than one cancer site; however, a systematic evaluation of the associations between risk variants for other cancers and lung cancer risk has yet to be performed.

METHODS

We included 18023 patients with lung cancer and 60543 control subjects from two consortia, Population Architecture using Genomics and Epidemiology (PAGE) and Transdisciplinary Research in Cancer of the Lung (TRICL). We examined 165 single-nucleotide polymorphisms (SNPs) that were previously associated with at least one of 16 non-lung cancer sites. Study-specific logistic regression results underwent meta-analysis, and associations were also examined by race/ethnicity, histological cell type, sex, and smoking status. A Bonferroni-corrected P value of 2.5×10(-5) was used to assign statistical significance.

RESULTS

The breast cancer SNP LSP1 rs3817198 was associated with an increased risk of lung cancer (odds ratio [OR] = 1.10; 95% confidence interval [CI] = 1.05 to 1.14; P = 2.8×10(-6)). This association was strongest for women with adenocarcinoma (P = 1.2×10(-4)) and not statistically significant in men (P = .14) with this cell type (P het by sex = .10). Two glioma risk variants, TERT rs2853676 and CDKN2BAS1 rs4977756, which are located in regions previously associated with lung cancer, were associated with increased risk of adenocarcinoma (OR = 1.16; 95% CI = 1.10 to 1.22; P = 1.1×10(-8)) and squamous cell carcinoma (OR = 1.13; CI = 1.07 to 1.19; P = 2.5×10(-5)), respectively.

CONCLUSIONS

Our findings demonstrate a novel pleiotropic association between the breast cancer LSP1 risk region marked by variant rs3817198 and lung cancer risk.

A genome scan for loci shared by autism spectrum disorder and language impairment

OBJECTIVE

The authors conducted a genetic linkage study of families that have both autism spectrum disorder (ASD) and language-impaired probands to find common communication impairment loci. The hypothesis was that these families have a high genetic loading for impairments in language ability, thus influencing the language and communication deficits of the family members with ASD. Comprehensive behavioral phenotyping of the families also enabled linkage analysis of quantitative measures, including normal, subclinical, and disordered variation in all family members for the three general autism symptom domains: social, communication, and compulsive behaviors.

METHOD

The primary linkage analysis coded persons with either ASD or specific language impairment as "affected." The secondary linkage analysis consisted of quantitative metrics of autism-associated behaviors capturing normal to clinically severe variation, measured in all family members.

RESULTS

Linkage to language phenotypes was established at two novel chromosomal loci, 15q23-26 and 16p12. The secondary analysis of normal and disordered quantitative variation in social and compulsive behaviors established linkage to two loci for social behaviors (at 14q and 15q) and one locus for repetitive behaviors (at 13q).

CONCLUSION

These data indicate shared etiology of ASD and specific language impairment at two novel loci. Additionally, nonlanguage phenotypes based on social aloofness and rigid personality traits showed compelling evidence for linkage in this study group. Further genetic mapping is warranted at these loci.

Fine Mapping and Identification of BMI Loci in African Americans

Genome-wide association studies (GWASs) primarily performed in European-ancestry (EA) populations have identified numerous loci associated with body mass index (BMI). However, it is still unclear whether these GWAS loci can be generalized to other ethnic groups, such as African Americans (AAs). Furthermore, the putative functional variant or variants in these loci mostly remain under investigation. The overall lower linkage disequilibrium in AA compared to EA populations provides the opportunity to narrow in or fine-map these BMI-related loci. Therefore, we used the Metabochip to densely genotype and evaluate 21 BMI GWAS loci identified in EA studies in 29,151 AAs from the Population Architecture using Genomics and Epidemiology (PAGE) study. Eight of the 21 loci (SEC16B, TMEM18, ETV5, GNPDA2, TFAP2B, BDNF, FTO, and MC4R) were found to be associated with BMI in AAs at 5.8 × 10(-5). Within seven out of these eight loci, we found that, on average, a substantially smaller number of variants was correlated (r(2) > 0.5) with the most significant SNP in AA than in EA populations (16 versus 55). Conditional analyses revealed GNPDA2 harboring a potential additional independent signal. Moreover, Metabochip-wide discovery analyses revealed two BMI-related loci, BRE (rs116612809, p = 3.6 × 10(-8)) and DHX34 (rs4802349, p = 1.2 × 10(-7)), which were significant when adjustment was made for the total number of SNPs tested across the chip. These results demonstrate that fine mapping in AAs is a powerful approach for both narrowing in on the underlying causal variants in known loci and discovering BMI-related loci.

Genetic variants associated with fasting glucose and insulin concentrations in an ethnically diverse population: results from the Population Architecture using Genomics and Epidemiology (PAGE) study

BACKGROUND

Multiple genome-wide association studies (GWAS) within European populations have implicated common genetic variants associated with insulin and glucose concentrations. In contrast, few studies have been conducted within minority groups, which carry the highest burden of impaired glucose homeostasis and type 2 diabetes in the U.S.

METHODS

As part of the 'Population Architecture using Genomics and Epidemiology (PAGE) Consortium, we investigated the association of up to 10 GWAS-identified single nucleotide polymorphisms (SNPs) in 8 genetic regions with glucose or insulin concentrations in up to 36,579 non-diabetic subjects including 23,323 European Americans (EA) and 7,526 African Americans (AA), 3,140 Hispanics, 1,779 American Indians (AI), and 811 Asians. We estimated the association between each SNP and fasting glucose or log-transformed fasting insulin, followed by meta-analysis to combine results across PAGE sites.

RESULTS

Overall, our results show that 9/9 GWAS SNPs are associated with glucose in EA (p = 0.04 to 9 × 10-15), versus 3/9 in AA (p= 0.03 to 6 × 10-5), 3/4 SNPs in Hispanics, 2/4 SNPs in AI, and 1/2 SNPs in Asians. For insulin we observed a significant association with rs780094/GCKR in EA, Hispanics and AI only.

CONCLUSIONS

Generalization of results across multiple racial/ethnic groups helps confirm the relevance of some of these loci for glucose and insulin metabolism. Lack of association in non-EA groups may be due to insufficient power, or to unique patterns of linkage disequilibrium.

Generalization and dilution of association results from European GWAS in populations of non-European ancestry: the PAGE study

The vast majority of genome-wide association study (GWAS) findings reported to date are from populations with European Ancestry (EA), and it is not yet clear how broadly the genetic associations described will generalize to populations of diverse ancestry. The Population Architecture Using Genomics and Epidemiology (PAGE) study is a consortium of multi-ancestry, population-based studies formed with the objective of refining our understanding of the genetic architecture of common traits emerging from GWAS. In the present analysis of five common diseases and traits, including body mass index, type 2 diabetes, and lipid levels, we compare direction and magnitude of effects for GWAS-identified variants in multiple non-EA populations against EA findings. We demonstrate that, in all populations analyzed, a significant majority of GWAS-identified variants have allelic associations in the same direction as in EA, with none showing a statistically significant effect in the opposite direction, after adjustment for multiple testing. However, 25% of tagSNPs identified in EA GWAS have significantly different effect sizes in at least one non-EA population, and these differential effects were most frequent in African Americans where all differential effects were diluted toward the null. We demonstrate that differential LD between tagSNPs and functional variants within populations contributes significantly to dilute effect sizes in this population. Although most variants identified from GWAS in EA populations generalize to all non-EA populations assessed, genetic models derived from GWAS findings in EA may generate spurious results in non-EA populations due to differential effect sizes. Regardless of the origin of the differential effects, caution should be exercised in applying any genetic risk prediction model based on tagSNPs outside of the ancestry group in which it was derived. Models based directly on functional variation may generalize more robustly, but the identification of functional variants remains challenging.

Imputation of coding variants in African Americans: better performance using data from the exome sequencing project

SUMMARY

Although the 1000 Genomes haplotypes are the most commonly used reference panel for imputation, medical sequencing projects are generating large alternate sets of sequenced samples. Imputation in African Americans using 3384 haplotypes from the Exome Sequencing Project, compared with 2184 haplotypes from 1000 Genomes Project, increased effective sample size by 8.3-11.4% for coding variants with minor allele frequency <1%. No loss of imputation quality was observed using a panel built from phenotypic extremes. We recommend using haplotypes from Exome Sequencing Project alone or concatenation of the two panels over quality score-based post-imputation selection or IMPUTE2's two-panel combination.

CONTACT

yunli@med.unc.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.