Leveraging health systems data to characterize a large effect variant conferring risk for liver disease in Puerto Ricans.

The integration of genomic data into health systems offers opportunities to identify genomic factors underlying the continuum of rare and common disease. We applied a population-scale haplotype association approach based on identity-by-descent (IBD) in a large multi-ethnic biobank to a spectrum of disease outcomes derived from electronic health records (EHRs) and uncovered a risk locus for liver disease. We used genome sequencing and in silico approaches to fine-map the signal to a non-coding variant (c.2784-12T>C) in the gene ABCB4. In vitro analysis confirmed the variant disrupted splicing of the ABCB4 pre-mRNA. Four of five homozygotes had evidence of advanced liver disease, and there was a significant association with liver disease among heterozygotes, suggesting the variant is linked to increased risk of liver disease in an allele dose-dependent manner. Population-level screening revealed the variant to be at a carrier rate of 1.95% in Puerto Rican individuals, likely as the result of a Puerto Rican founder effect. This work demonstrates that integrating EHR and genomic data at a population scale can facilitate strategies for understanding the continuum of genomic risk for common diseases, particularly in populations underrepresented in genomic medicine.

Bayesian multivariate reanalysis of large genetic studies identifies many new associations.

Genome-wide association studies (GWAS) have now been conducted for hundreds of phenotypes of relevance to human health. Many such GWAS involve multiple closely-related phenotypes collected on the same samples. However, the vast majority of these GWAS have been analyzed using simple univariate analyses, which consider one phenotype at a time. This is despite the fact that, at least in simulation experiments, multivariate analyses have been shown to be more powerful at detecting associations. Here, we conduct multivariate association analyses on 13 different publicly-available GWAS datasets that involve multiple closely-related phenotypes. These data include large studies of anthropometric traits (GIANT), plasma lipid traits (GlobalLipids), and red blood cell traits (HaemgenRBC). Our analyses identify many new associations (433 in total across the 13 studies), many of which replicate when follow-up samples are available. Overall, our results demonstrate that multivariate analyses can help make more effective use of data from both existing and future GWAS.

Genome-wide association of copy-number variation reveals an association between short stature and the presence of low-frequency genomic deletions.

Height is a model polygenic trait that is highly heritable. Genome-wide association studies have identified hundreds of single-nucleotide polymorphisms associated with stature, but the role of structural variation in determining height is largely unknown. We performed a genome-wide association study of copy-number variation and stature in a clinical cohort of children who had undergone comparative genomic hybridization (CGH) microarray analysis for clinical indications. We found that subjects with short stature had a greater global burden of copy-number variants (CNVs) and a greater average CNV length than did controls (p < 0.002). These associations were present for lower-frequency (<5%) and rare (<1%) deletions, but there were no significant associations seen for duplications. Known gene-deletion syndromes did not account for our findings, and we saw no significant associations with tall stature. We then extended our findings into a population-based cohort and found that, in agreement with the clinical cohort study, an increased burden of lower-frequency deletions was associated with shorter stature (p = 0.015). Our results suggest that in individuals undergoing copy-number analysis for clinical indications, short stature increases the odds that a low-frequency deletion will be found. Additionally, copy-number variation might contribute to genetic variation in stature in the general population.

Gencrypt: one-way cryptographic hashes to detect overlapping individuals across samples.

Meta-analysis across genome-wide association studies is a common approach for discovering genetic associations. However, in some meta-analysis efforts, individual-level data cannot be broadly shared by study investigators due to privacy and Institutional Review Board concerns. In such cases, researchers cannot confirm that each study represents a unique group of people, leading to potentially inflated test statistics and false positives. To resolve this problem, we created a software tool, Gencrypt, which utilizes a security protocol known as one-way cryptographic hashes to allow overlapping participants to be identified without sharing individual-level data.

Genome-wide polygenic score to predict chronic kidney disease across ancestries.

Chronic kidney disease (CKD) is a common complex condition associated with high morbidity and mortality. Polygenic prediction could enhance CKD screening and prevention; however, this approach has not been optimized for ancestrally diverse populations. By combining APOL1 risk genotypes with genome-wide association studies (GWAS) of kidney function, we designed, optimized and validated a genome-wide polygenic score (GPS) for CKD. The new GPS was tested in 15 independent cohorts, including 3 cohorts of European ancestry (n = 97,050), 6 cohorts of African ancestry (n = 14,544), 4 cohorts of Asian ancestry (n = 8,625) and 2 admixed Latinx cohorts (n = 3,625). We demonstrated score transferability with reproducible performance across all tested cohorts. The top 2% of the GPS was associated with nearly threefold increased risk of CKD across ancestries. In African ancestry cohorts, the APOL1 risk genotype and polygenic component of the GPS had additive effects on the risk of CKD.

Evidence for positive selection on Drosophila melanogaster seminal fluid protease homologs.

Proteins present in the seminal fluid of Drosophila melanogaster (accessory gland proteins Acps) contribute to female postmating behavioral changes, sperm storage, sperm competition, and immunity. Consequently, male-female coevolution and host-pathogen interactions are thought to underlie the rapid, adaptive evolution that characterizes several Acp-encoding genes. We propose that seminal fluid proteases are likely targets of selection due to their demonstrated or potential roles in between-sex interactions and immune processes. We use within- and between-species sequence data for 5 predicted protease-encoding Acp loci to test this hypothesis. Our polymorphism-based analyses find evidence for positive selection at 2 genes, both of which encode predicted serine protease homologs. One of these genes, CG6069, also shows evidence for consistent selection on a subset of codons over a deeper evolutionary time scale. The second gene, CG9997, was previously shown to be essential for normal sperm usage, suggesting that sexual selection may underlie its history of adaptation.

Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies.

Genetic predictions of height differ among human populations and these differences have been interpreted as evidence of polygenic adaptation. These differences were first detected using SNPs genome-wide significantly associated with height, and shown to grow stronger when large numbers of sub-significant SNPs were included, leading to excitement about the prospect of analyzing large fractions of the genome to detect polygenic adaptation for multiple traits. Previous studies of height have been based on SNP effect size measurements in the GIANT Consortium meta-analysis. Here we repeat the analyses in the UK Biobank, a much more homogeneously designed study. We show that polygenic adaptation signals based on large numbers of SNPs below genome-wide significance are extremely sensitive to biases due to uncorrected population stratification. More generally, our results imply that typical constructions of polygenic scores are sensitive to population stratification and that population-level differences should be interpreted with caution. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

The deleterious mutation load is insensitive to recent population history.

Human populations have undergone major changes in population size in the past 100,000 years, including recent rapid growth. How these demographic events have affected the burden of deleterious mutations in individuals and the frequencies of disease mutations in populations remains unclear. We use population genetic models to show that recent human demography has probably had little impact on the average burden of deleterious mutations. This prediction is supported by two exome sequence data sets showing that individuals of west African and European ancestry carry very similar burdens of damaging mutations. We further show that for many diseases, rare alleles are unlikely to contribute a large fraction of the heritable variation, and therefore the impact of recent growth is likely to be modest. However, for those diseases that have a direct impact on fitness, strongly deleterious rare mutations probably do have an important role, and recent growth will have increased their impact.

Evidence of widespread selection on standing variation in Europe at height-associated SNPs.

Strong signatures of positive selection at newly arising genetic variants are well documented in humans(1-8), but this form of selection may not be widespread in recent human evolution(9). Because many human traits are highly polygenic and partly determined by common, ancient genetic variation, an alternative model for rapid genetic adaptation has been proposed: weak selection acting on many pre-existing (standing) genetic variants, or polygenic adaptation(10-12). By studying height, a classic polygenic trait, we demonstrate the first human signature of widespread selection on standing variation. We show that frequencies of alleles associated with increased height, both at known loci and genome wide, are systematically elevated in Northern Europeans compared with Southern Europeans (P < 4.3 × 10(-4)). This pattern mirrors intra-European height differences and is not confounded by ancestry or other ascertainment biases. The systematic frequency differences are consistent with the presence of widespread weak selection (selection coefficients ∼10(-3)-10(-5) per allele) rather than genetic drift alone (P < 10(-15)).