Search

The role of rare non-coding variation in complex human phenotypes is still largely unknown. To elucidate the impact of rare variants in regulatory elements, we performed a whole-genome sequencing association analysis for height using 333,100 individuals from three datasets: UK Biobank (N = 200,003), TOPMed (N = 87,652) and All of Us (N = 45,445). We performed rare ( < 0.1% minor-allele-frequency) single-variant and aggregate testing of non-coding variants in regulatory regions based on proximal-regulatory, intergenic-regulatory and deep-intronic annotation. We observed 29 independent variants associated with height at P <  6 × 10 - 10 after conditioning on previously reported variants, with effect sizes ranging from -7cm to +4.7 cm. We also identified and replicated non-coding aggregate-based associations proximal to HMGA1 containing variants associated with a 5 cm taller height and of highly-conserved variants in MIR497HG on chromosome 17. We have developed an approach for identifying non-coding rare variants in regulatory regions with large effects from whole-genome sequencing data associated with complex traits., (© 2024. The Author(s).)

Human genetic studies of common variants have provided substantial insight into the biological mechanisms that govern ovarian ageing 1 . Here we report analyses of rare protein-coding variants in 106,973 women from the UK Biobank study, implicating genes with effects around five times larger than previously found for common variants (ETAA1, ZNF518A, PNPLA8, PALB2 and SAMHD1). The SAMHD1 association reinforces the link between ovarian ageing and cancer susceptibility 1 , with damaging germline variants being associated with extended reproductive lifespan and increased all-cause cancer risk in both men and women. Protein-truncating variants in ZNF518A are associated with shorter reproductive lifespan-that is, earlier age at menopause (by 5.61 years) and later age at menarche (by 0.56 years). Finally, using 8,089 sequenced trios from the 100,000 Genomes Project (100kGP), we observe that common genetic variants associated with earlier ovarian ageing associate with an increased rate of maternally derived de novo mutations. Although we were unable to replicate the finding in independent samples from the deCODE study, it is consistent with the expected role of DNA damage response genes in maintaining the genetic integrity of germ cells. This study provides evidence of genetic links between age of menopause and cancer risk., (© 2024. The Author(s).)

Aims/hypothesis: Use of genetic risk scores (GRS) may help to distinguish between type 1 diabetes and type 2 diabetes, but less is known about whether GRS are associated with disease severity or progression after diagnosis. Therefore, we tested whether GRS are associated with residual beta cell function and glycaemic control in individuals with type 1 diabetes., Methods: Immunochip arrays and TOPMed were used to genotype a cross-sectional cohort (n=479, age 41.7 ± 14.9 years, duration of diabetes 16.0 years [IQR 6.0-29.0], HbA 1c 55.6 ± 12.2 mmol/mol). Several GRS, which were originally developed to assess genetic risk of type 1 diabetes (GRS-1, GRS-2) and type 2 diabetes (GRS-T2D), were calculated. GRS-C1 and GRS-C2 were based on SNPs that have previously been shown to be associated with residual beta cell function. Regression models were used to investigate the association between GRS and residual beta cell function, assessed using the urinary C-peptide/creatinine ratio, and the association between GRS and continuous glucose monitor metrics., Results: Higher GRS-1 and higher GRS-2 both showed a significant association with undetectable UCPCR (OR 0.78; 95% CI 0.69, 0.89 and OR 0.84: 95% CI 0.75, 0.93, respectively), which were attenuated after correction for sex and age of onset (GRS-2) and disease duration (GRS-1). Higher GRS-C2 was associated with detectable urinary C-peptide/creatinine ratio (≥0.01 nmol/mmol) after correction for sex and age of onset (OR 6.95; 95% CI 1.19, 40.75). A higher GRS-T2D was associated with less time below range (TBR) (OR for TBR<4% 1.41; 95% CI 1.01 to 1.96) and lower glucose coefficient of variance (β -1.53; 95% CI -2.76, -0.29)., Conclusions/interpretation: Diabetes-related GRS are associated with residual beta cell function in individuals with type 1 diabetes. These findings suggest some genetic contribution to preservation of beta cell function., (© 2024. The Author(s).)

Penetrance is the probability that an individual with a pathogenic genetic variant develops a specific disease. Knowing the penetrance of variants for monogenic disorders is important for counseling of individuals. Until recently, estimates of penetrance have largely relied on affected individuals and their at-risk family members being clinically referred for genetic testing, a 'phenotype-first' approach. This approach substantially overestimates the penetrance of variants because of ascertainment bias. The recent availability of whole-genome sequencing data in individuals from very-large-scale population-based cohorts now allows 'genotype-first' estimates of penetrance for many conditions. Although this type of population-based study can underestimate penetrance owing to recruitment biases, it provides more accurate estimates of penetrance for secondary or incidental findings. Here, we provide guidance for the conduct of penetrance studies to ensure that robust genotypes and phenotypes are used to accurately estimate penetrance of variants and groups of similarly annotated variants from population-based studies., (© 2024. Springer Nature America, Inc.)

Colorectal cancer (CRC) is a leading cause of cancer mortality worldwide. Accurate cancer risk assessment approaches could increase rates of early CRC diagnosis, improve health outcomes for patients and reduce pressure on diagnostic services. The faecal immunochemical test (FIT) for blood in stool is widely used in primary care to identify symptomatic patients with likely CRC. However, there is a 6-16% noncompliance rate with FIT in clinic and ~90% of patients over the symptomatic 10 µg/g test threshold do not have CRC. A polygenic risk score (PRS) quantifies an individual's genetic risk of a condition based on many common variants. Existing PRS for CRC have so far been used to stratify asymptomatic populations. We conducted a retrospective cohort study of 50,387 UK Biobank participants with a CRC symptom in their primary care record at age 40+. A PRS based on 201 variants, 5 genetic principal components and 22 other risk factors and markers for CRC were assessed for association with CRC diagnosis within 2 years of first symptom presentation using logistic regression. Associated variables were included in an integrated risk model and trained in 80% of the cohort to predict CRC diagnosis within 2 years. An integrated risk model combining PRS, age, sex, and patient-reported symptoms was predictive of CRC development in a testing cohort (receiver operating characteristic area under the curve, ROCAUC: 0.76, 95% confidence interval: 0.71-0.81). This model has the potential to improve early diagnosis of CRC, particularly in cases of patient noncompliance with FIT., (© 2024. The Author(s).)

Persistent congenital hyperinsulinism (HI) is a rare genetically heterogeneous condition characterised by dysregulated insulin secretion leading to life-threatening hypoglycaemia. For up to 50% of affected individuals screening of the known HI genes does not identify a disease-causing variant. Large deletions have previously been used to identify novel regulatory regions causing HI. Here, we used genome sequencing to search for novel large (>1 Mb) deletions in 180 probands with HI of unknown cause and replicated our findings in a large cohort of 883 genetically unsolved individuals with HI using off-target copy number variant calling from targeted gene panels. We identified overlapping heterozygous deletions in five individuals (range 3-8 Mb) spanning chromosome 20p11.2. The pancreatic beta-cell transcription factor gene, FOXA2, a known cause of HI was deleted in two of the five individuals. In the remaining three, we found a minimal deleted region of 2.4 Mb adjacent to FOXA2 that encompasses multiple non-coding regulatory elements that are in conformational contact with FOXA2. Our data suggests that the deletions in these three children may cause disease through the dysregulation of FOXA2 expression. These findings provide new insights into the regulation of FOXA2 in the beta-cell and confirm an aetiological role for chromosome 20p11.2 deletions in syndromic HI., (© 2024. The Author(s).)

Self-reported shorter/longer sleep duration, insomnia, and evening preference are associated with hyperglycaemia in observational analyses, with similar observations in small studies using accelerometer-derived sleep traits. Mendelian randomization (MR) studies support an effect of self-reported insomnia, but not others, on glycated haemoglobin (HbA1c). To explore potential effects, we used MR methods to assess effects of accelerometer-derived sleep traits (duration, mid-point least active 5-h, mid-point most active 10-h, sleep fragmentation, and efficiency) on HbA1c/glucose in European adults from the UK Biobank (UKB) (n = 73,797) and the MAGIC consortium (n = 146,806). Cross-trait linkage disequilibrium score regression was applied to determine genetic correlations across accelerometer-derived, self-reported sleep traits, and HbA1c/glucose. We found no causal effect of any accelerometer-derived sleep trait on HbA1c or glucose. Similar MR results for self-reported sleep traits in the UKB sub-sample with accelerometer-derived measures suggested our results were not explained by selection bias. Phenotypic and genetic correlation analyses suggested complex relationships between self-reported and accelerometer-derived traits indicating that they may reflect different types of exposure. These findings suggested accelerometer-derived sleep traits do not affect HbA1c. Accelerometer-derived measures of sleep duration and quality might not simply be 'objective' measures of self-reported sleep duration and insomnia, but rather captured different sleep characteristics., (© 2024. The Author(s).)

Rare damaging variants in a large number of genes are known to cause monogenic developmental disorders (DDs) and have also been shown to cause milder subclinical phenotypes in population cohorts. Here, we show that carrying multiple (2-5) rare damaging variants across 599 dominant DD genes has an additive adverse effect on numerous cognitive and socioeconomic traits in UK Biobank, which can be partially counterbalanced by a higher educational attainment polygenic score (EA-PGS). Phenotypic deviators from expected EA-PGS could be partly explained by the enrichment or depletion of rare DD variants. Among carriers of rare DD variants, those with a DD-related clinical diagnosis had a substantially lower EA-PGS and more severe phenotype than those without a clinical diagnosis. Our results suggest that the overall burden of both rare and common variants can modify the expressivity of a phenotype, which may then influence whether an individual reaches the threshold for clinical disease., (© 2024. The Author(s).)

Optic neuritis (ON) is associated with numerous immune-mediated inflammatory diseases, but 50% patients are ultimately diagnosed with multiple sclerosis (MS). Differentiating MS-ON from non-MS-ON acutely is challenging but important; non-MS ON often requires urgent immunosuppression to preserve vision. Using data from the United Kingdom Biobank we showed that combining an MS-genetic risk score (GRS) with demographic risk factors (age, sex) significantly improved MS prediction in undifferentiated ON; one standard deviation of MS-GRS increased the Hazard of MS 1.3-fold (95% confidence interval 1.07-1.55, P < 0.01). Participants stratified into quartiles of predicted risk developed incident MS at rates varying from 4% (95%CI 0.5-7%, lowest risk quartile) to 41% (95%CI 33-49%, highest risk quartile). The model replicated across two cohorts (Geisinger, USA, and FinnGen, Finland). This study indicates that a combined model might enhance individual MS risk stratification, paving the way for precision-based ON treatment and earlier MS disease-modifying therapy., (© 2024. The Author(s).)

Whole genome sequencing (WGS) from large clinically unselected cohorts provides a unique opportunity to assess the penetrance and expressivity of rare and/or known pathogenic mitochondrial variants in population. Using WGS from 179 862 clinically unselected individuals from the UK Biobank, we performed extensive single and rare variant aggregation association analyses of 15 881 mtDNA variants and 73 known pathogenic variants with 15 mitochondrial disease-relevant phenotypes. We identified 12 homoplasmic and one heteroplasmic variant (m.3243A>G) with genome-wide significant associations in our clinically unselected cohort. Heteroplasmic m.3243A>G (MAF = 0.0002, a known pathogenic variant) was associated with diabetes, deafness and heart failure and 12 homoplasmic variants increased aspartate aminotransferase levels including three low-frequency variants (MAF ~0.002 and beta~0.3 SD). Most pathogenic mitochondrial disease variants (n = 66/74) were rare in the population (<1:9000). Aggregated or single variant analysis of pathogenic variants showed low penetrance in unselected settings for the relevant phenotypes, except m.3243A>G. Multi-system disease risk and penetrance of diabetes, deafness and heart failure greatly increased with m.3243A>G level ≥ 10%. The odds ratio of these traits increased from 5.61, 12.3 and 10.1 to 25.1, 55.0 and 39.5, respectively. Diabetes risk with m.3243A>G was further influenced by type 2 diabetes genetic risk. Our study of mitochondrial variation in a large-unselected population identified novel associations and demonstrated that pathogenic mitochondrial variants have lower penetrance in clinically unselected settings. m.3243A>G was an exception at higher heteroplasmy showing a significant impact on health making it a good candidate for incidental reporting., (© The Author(s) 2023. Published by Oxford University Press.)

Background: Diabetes (regardless of type) and obesity are associated with a range of musculoskeletal disorders. The causal mechanisms driving these associations are unknown for many upper limb pathologies. We used genetic techniques to test the causal link between glycemia, obesity and musculoskeletal conditions., Methods: In the UK Biobank's unrelated European cohort (N = 379 708) we performed mendelian randomisation (MR) analyses to test for a causal effect of long-term high glycaemia and adiposity on four musculoskeletal pathologies: frozen shoulder, Dupuytren's disease, carpal tunnel syndrome and trigger finger. We also performed single-gene MR using rare variants in the GCK gene., Results: Using MR, we found evidence that long-term high glycaemia has a causal role in the aetiology of upper limb conditions. A 10-mmol/mol increase in genetically predicted haemoglobin A1C (HbA1c) was associated with frozen shoulder: odds ratio (OR) = 1.50 [95% confidence interval (CI), 1.20-1.88], Dupuytren's disease: OR = 1.17 (95% CI, 1.01-1.35), trigger finger: OR = 1.30 (95% CI, 1.09-1.55) and carpal tunnel syndrome: OR = 1.20 (95% CI, 1.09-1.33). Carriers of GCK mutations have increased odds of frozen shoulder: OR = 7.16 (95% CI, 2.93-17.51) and carpal tunnel syndrome: OR = 2.86 (95% CI, 1.50-5.44) but not Dupuytren's disease or trigger finger. We found evidence that an increase in genetically predicted body mass index (BMI) of 5 kg/m2 was associated with carpal tunnel syndrome: OR = 1.13 (95% CI, 1.10-1.16) and associated negatively with Dupuytren's disease: OR = 0.94 (95% CI, 0.90-0.98), but no evidence of association with frozen shoulder or trigger finger. Trigger finger (OR 1.96 (95% CI, 1.42-2.69) P = 3.6e-05) and carpal tunnel syndrome [OR 1.63 (95% CI, 1.36-1.95) P = 8.5e-08] are associated with genetically predicted unfavourable adiposity increase of one standard deviation of body fat., Conclusions: Our study consistently demonstrates a causal role of long-term high glycaemia in the aetiology of upper limb musculoskeletal conditions. Clinicians treating diabetes patients should be aware of these complications in clinic, specifically those managing the care of GCK mutation carriers. Upper limb musculoskeletal conditions should be considered diabetes complications., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Epidemiological Association.)