Genetic insights into the age-specific biological mechanisms governing human ovarian aging.

There is currently little evidence that the genetic basis of human phenotype varies significantly across the lifespan. However, time-to-event phenotypes are understudied and can be thought of as reflecting an underlying hazard, which is unlikely to be constant through life when values take a broad range. Here, we find that 74% of 245 genome-wide significant genetic associations with age at natural menopause (ANM) in the UK Biobank show a form of age-specific effect. Nineteen of these replicated discoveries are identified only by our modeling framework, which determines the time dependency of DNA-variant age-at-onset associations without a significant multiple-testing burden. Across the range of early to late menopause, we find evidence for significantly different underlying biological pathways, changes in the signs of genetic correlations of ANM to health indicators and outcomes, and differences in inferred causal relationships. We find that DNA damage response processes only act to shape ovarian reserve and depletion for women of early ANM. Genetically mediated delays in ANM were associated with increased relative risk of breast cancer and leiomyoma at all ages and with high cholesterol and heart failure for late-ANM women. These findings suggest that a better understanding of the age dependency of genetic risk factor relationships among health indicators and outcomes is achievable through appropriate statistical modeling of large-scale biobank data.

Liability-scale heritability estimation for biobank studies of low-prevalence disease.

Theory for liability-scale models of the underlying genetic basis of complex disease provides an important way to interpret, compare, and understand results generated from biological studies. In particular, through estimation of the liability-scale heritability (LSH), liability models facilitate an understanding and comparison of the relative importance of genetic and environmental risk factors that shape different clinically important disease outcomes. Increasingly, large-scale biobank studies that link genetic information to electronic health records, containing hundreds of disease diagnosis indicators that mostly occur infrequently within the sample, are becoming available. Here, we propose an extension of the existing liability-scale model theory suitable for estimating LSH in biobank studies of low-prevalence disease. In a simulation study, we find that our derived expression yields lower mean square error (MSE) and is less sensitive to prevalence misspecification as compared to previous transformations for diseases with ≤2% population prevalence and LSH of ≤0.45, especially if the biobank sample prevalence is less than that of the wider population. Applying our expression to 13 diagnostic outcomes of ≤3% prevalence in the UK Biobank study revealed important differences in LSH obtained from the different theoretical expressions that impact the conclusions made when comparing LSH across disease outcomes. This demonstrates the importance of careful consideration for estimation and prediction of low-prevalence disease outcomes and facilitates improved inference of the underlying genetic basis of ≤2% population prevalence diseases, especially where biobank sample ascertainment results in a healthier sample population.

Improving GWAS discovery and genomic prediction accuracy in biobank data.

Genetically informed, deep-phenotyped biobanks are an important research resource and it is imperative that the most powerful, versatile, and efficient analysis approaches are used. Here, we apply our recently developed Bayesian grouped mixture of regressions model (GMRM) in the UK and Estonian Biobanks and obtain the highest genomic prediction accuracy reported to date across 21 heritable traits. When compared to other approaches, GMRM accuracy was greater than annotation prediction models run in the LDAK or LDPred-funct software by 15% (SE 7%) and 14% (SE 2%), respectively, and was 18% (SE 3%) greater than a baseline BayesR model without single-nucleotide polymorphism (SNP) markers grouped into minor allele frequency-linkage disequilibrium (MAF-LD) annotation categories. For height, the prediction accuracy R2 was 47% in a UK Biobank holdout sample, which was 76% of the estimated [Formula: see text]. We then extend our GMRM prediction model to provide mixed-linear model association (MLMA) SNP marker estimates for genome-wide association (GWAS) discovery, which increased the independent loci detected to 16,162 in unrelated UK Biobank individuals, compared to 10,550 from BoltLMM and 10,095 from Regenie, a 62 and 65% increase, respectively. The average [Formula: see text] value of the leading markers increased by 15.24 (SE 0.41) for every 1% increase in prediction accuracy gained over a baseline BayesR model across the traits. Thus, we show that modeling genetic associations accounting for MAF and LD differences among SNP markers, and incorporating prior knowledge of genomic function, is important for both genomic prediction and discovery in large-scale individual-level studies.

Incarceration history is associated with HIV infection among community-recruited people who inject drugs in Europe: A propensity-score matched analysis of cross-sectional studies.

AIMS: We measured the association between a history of incarceration and HIV positivity among people who inject drugs (PWID) across Europe. DESIGN, SETTING AND PARTICIPANTS: This was a cross-sectional, multi-site, multi-year propensity-score matched analysis conducted in Europe. Participants comprised community-recruited PWID who reported a recent injection (within the last 12 months). MEASUREMENTS: Data on incarceration history, demographics, substance use, sexual behavior and harm reduction service use originated from cross-sectional studies among PWID in Europe. Our primary outcome was HIV status. Generalized linear mixed models and propensity-score matching were used to compare HIV status between ever- and never-incarcerated PWID. FINDINGS: Among 43 807 PWID from 82 studies surveyed (in 22 sites and 13 countries), 58.7% reported having ever been in prison and 7.16% (n = 3099) tested HIV-positive. Incarceration was associated with 30% higher odds of HIV infection [adjusted odds ratio (aOR) = 1.32, 95% confidence interval (CI) = 1.09-1.59]; the association between a history of incarceration and HIV infection was strongest among PWID, with the lowest estimated propensity-score for having a history of incarceration (aOR = 1.78, 95% CI = 1.47-2.16). Additionally, mainly injecting cocaine and/or opioids (aOR = 2.16, 95% CI = 1.33-3.53), increased duration of injecting drugs (per 8 years aOR = 1.31, 95% CI = 1.16-1.48), ever sharing needles/syringes (aOR = 1.91, 95% CI = 1.59-2.28) and increased income inequality among the general population (measured by the Gini index, aOR = 1.34, 95% CI = 1.18-1.51) were associated with a higher odds of HIV infection. Older age (per 8 years aOR = 0.84, 95% CI = 0.76-0.94), male sex (aOR = 0.77, 95% CI = 0.65-0.91) and reporting pharmacies as the main source of clean syringes (aOR = 0.72, 95% CI = 0.59-0.88) were associated with lower odds of HIV positivity. CONCLUSIONS: A history of incarceration appears to be independently associated with HIV infection among people who inject drugs (PWID) in Europe, with a stronger effect among PWID with lower probability of incarceration.

Probabilistic inference of the genetic architecture underlying functional enrichment of complex traits.

We develop a Bayesian model (BayesRR-RC) that provides robust SNP-heritability estimation, an alternative to marker discovery, and accurate genomic prediction, taking 22 seconds per iteration to estimate 8.4 million SNP-effects and 78 SNP-heritability parameters in the UK Biobank. We find that only ≤10% of the genetic variation captured for height, body mass index, cardiovascular disease, and type 2 diabetes is attributable to proximal regulatory regions within 10kb upstream of genes, while 12-25% is attributed to coding regions, 32-44% to introns, and 22-28% to distal 10-500kb upstream regions. Up to 24% of all cis and coding regions of each chromosome are associated with each trait, with over 3,100 independent exonic and intronic regions and over 5,400 independent regulatory regions having ≥95% probability of contributing ≥0.001% to the genetic variance of these four traits. Our open-source software (GMRM) provides a scalable alternative to current approaches for biobank data.

Genomic architecture and prediction of censored time-to-event phenotypes with a Bayesian genome-wide analysis.

While recent advancements in computation and modelling have improved the analysis of complex traits, our understanding of the genetic basis of the time at symptom onset remains limited. Here, we develop a Bayesian approach (BayesW) that provides probabilistic inference of the genetic architecture of age-at-onset phenotypes in a sampling scheme that facilitates biobank-scale time-to-event analyses. We show in extensive simulation work the benefits BayesW provides in terms of number of discoveries, model performance and genomic prediction. In the UK Biobank, we find many thousands of common genomic regions underlying the age-at-onset of high blood pressure (HBP), cardiac disease (CAD), and type-2 diabetes (T2D), and for the genetic basis of onset reflecting the underlying genetic liability to disease. Age-at-menopause and age-at-menarche are also highly polygenic, but with higher variance contributed by low frequency variants. Genomic prediction into the Estonian Biobank data shows that BayesW gives higher prediction accuracy than other approaches.