The role of accelerometer-derived sleep traits on glycated haemoglobin and glucose levels: a Mendelian randomization study.

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.

Associations between diabetes-related genetic risk scores and residual beta cell function in type 1 diabetes: the GUTDM1 study.

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], HbA1c 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.

Chromosome 20p11.2 deletions cause congenital hyperinsulinism via the loss of FOXA2 or its regulatory elements.

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.

Genetic modifiers of rare variants in monogenic developmental disorder loci.

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.

Genome-wide association analysis of composite sleep health scores in 413,904 individuals.

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Applying a genetic risk score model to enhance prediction of future multiple sclerosis diagnosis at first presentation with optic neuritis.

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.

Hyperglycaemia is a causal risk factor for upper limb pathologies.

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.

Penetrance and expressivity of mitochondrial variants in a large clinically unselected population.

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.

Influence of family history on penetrance of hereditary cancers in a population setting.

Background: We sought to investigate how penetrance of familial cancer syndromes varies with family history using a population-based cohort. Methods: We analysed 454,712 UK Biobank participants with exome sequence and clinical data (data collected between March 2006 and June 2021). We identified participants with a self-reported family history of breast or colorectal cancer and a pathogenic/likely pathogenic variant in the major genes responsible for hereditary breast cancer or Lynch syndrome. We calculated survival to cancer diagnosis (controlled for sex, death, recruitment centre, screening and prophylactic surgery). Findings: Women with a pathogenic BRCA1 or BRCA2 variant had an increased risk of breast cancer that was higher in those with a first-degree family history (relative hazard 10.3 and 7.8, respectively) than those without (7.2 and 4.7). Penetrance to age 60 was also higher in those with a family history (44.7%, CI 32.2-59.3 and 24.1%, CI 17.5-32.6) versus those without (22.8%, CI 15.9-32.0 and 17.9%, CI 13.8-23.0). A similar pattern was seen in Lynch syndrome: individuals with a pathogenic MLH1, MSH2 or MSH6 variant had an increased risk of colorectal cancer that was significantly higher in those with a family history (relative hazard 35.6, 48.0 and 9.9) than those without (13.0, 15.4 and 7.2). Penetrance to age 60 was also higher for carriers of a pathogenic MLH1 or MSH2 variant in those with a family history (30.9%, CI 18.1-49.3 and 38.3%, CI 21.5-61.8) versus those without (20.5% CI 9.6-40.5 and 8.3% CI 2.1-30.4), but not for MSH6 (6.5% CI 2.7-15.1 with family history versus 8.3%, CI 5.1-13.2). Relative risk increases were also observed both within and across conditions. Interpretation: Individuals with pathogenic cancer syndrome variants may be at a less elevated risk of cancer in the absence of a first-degree family history, so in the context of results return, family history should be considered when counselling patients on the risks and benefits of potential follow-up care. Funding: The current work is supported by the MRC (grant no MR/T00200X/1). The MRC had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Correlations in sleeping patterns and circadian preference between spouses.

Spouses may affect each other's sleeping behaviour. In 47,420 spouse-pairs from the UK Biobank, we found a weak positive phenotypic correlation between spouses for self-reported sleep duration (r = 0.11; 95% CI = 0.10, 0.12) and a weak inverse correlation for chronotype (diurnal preference) (r = -0.11; -0.12, -0.10), which replicated in up to 127,035 23andMe spouse-pairs. Using accelerometer data on 3454 UK Biobank spouse-pairs, the correlation for derived sleep duration was similar to self-report (r = 0.12; 0.09, 0.15). Timing of diurnal activity was positively correlated (r = 0.24; 0.21, 0.27) in contrast to the inverse correlation for chronotype. In Mendelian randomization analysis, positive effects of sleep duration (mean difference=0.13; 0.04, 0.23 SD per SD) and diurnal activity (0.49; 0.03, 0.94) were observed, as were inverse effects of chronotype (-0.15; -0.26, -0.04) and snoring (-0.15; -0.27, -0.04). Findings support the notion that an individual's sleep may impact that of their partner, promoting opportunities for sleep interventions at the family-level.

Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records.

BACKGROUND: Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records. METHODS: We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time. RESULTS: Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after. CONCLUSIONS: We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

Identification and analysis of individuals who deviate from their genetically-predicted phenotype.

Findings from genome-wide association studies have facilitated the generation of genetic predictors for many common human phenotypes. Stratifying individuals misaligned to a genetic predictor based on common variants may be important for follow-up studies that aim to identify alternative causal factors. Using genome-wide imputed genetic data, we aimed to classify 158,951 unrelated individuals from the UK Biobank as either concordant or deviating from two well-measured phenotypes. We first applied our methods to standing height: our primary analysis classified 244 individuals (0.15%) as misaligned to their genetically predicted height. We show that these individuals are enriched for self-reporting being shorter or taller than average at age 10, diagnosed congenital malformations, and rare loss-of-function variants in genes previously catalogued as causal for growth disorders. Secondly, we apply our methods to LDL cholesterol (LDL-C). We classified 156 (0.12%) individuals as misaligned to their genetically predicted LDL-C and show that these individuals were enriched for both clinically actionable cardiovascular risk factors and rare genetic variants in genes previously shown to be involved in metabolic processes. Individuals whose LDL-C was higher than expected based on the genetic predictor were also at higher risk of developing coronary artery disease and type-two diabetes, even after adjustment for measured LDL-C, BMI and age, suggesting upward deviation from genetically predicted LDL-C is indicative of generally poor health. Our results remained broadly consistent when performing sensitivity analysis based on a variety of parametric and non-parametric methods to define individuals deviating from polygenic expectation. Our analyses demonstrate the potential importance of quantitatively identifying individuals for further follow-up based on deviation from genetic predictions.

The Role of ONECUT1 Variants in Monogenic and Type 2 Diabetes Mellitus.

ONECUT1 (also known as HNF6) is a transcription factor involved in pancreatic development and ß-cell function. Recently, biallelic variants in ONECUT1 were reported as a cause of neonatal diabetes mellitus (NDM) in two subjects, and missense monoallelic variants were associated with type 2 diabetes and possibly maturity-onset diabetes of the young (MODY). Here we examine the role of ONECUT1 variants in NDM, MODY, and type 2 diabetes in large international cohorts of subjects with monogenic diabetes and >400,000 subjects from UK Biobank. We identified a biallelic frameshift ONECUT1 variant as the cause of NDM in one individual. However, we found no enrichment of missense or null ONECUT1 variants among 484 individuals clinically suspected of MODY, in whom all known genes had been excluded. Finally, using a rare variant burden test in the UK Biobank European cohort, we identified a significant association between heterozygous ONECUT1 null variants and type 2 diabetes (P = 0.006) but did not find an association between missense variants and type 2 diabetes. Our results confirm biallelic ONECUT1 variants as a cause of NDM and highlight monoallelic null variants as a risk factor for type 2 diabetes. These findings confirm the critical role of ONECUT1 in human ß-cell function.

Utility of genetic risk scores in type 1 diabetes.

Iterative advances in understanding of the genetics of type 1 diabetes have identified >70 genetic regions associated with risk of the disease, including strong associations across the HLA class II region that account for >50% of heritability. The increased availability of genetic data combined with the decreased costs of generating these data, have facilitated the development of polygenic scores that aggregate risk variants from associated loci into a single number: either a genetic risk score (GRS) or a polygenic risk score (PRS). PRSs incorporate the risk of many possibly correlated variants from across the genome, even if they do not reach genome-wide significance, whereas GRSs estimate the cumulative contribution of a smaller subset of genetic variants that reach genome-wide significance. Type 1 diabetes GRSs have utility in diabetes classification, aiding discrimination between type 1 diabetes, type 2 diabetes and MODY. Type 1 diabetes GRSs are also being used in newborn screening studies to identify infants at risk of future presentation of the disease. Most early studies of type 1 diabetes genetics have been conducted in European ancestry populations, but, to develop accurate GRSs across diverse ancestries, large case-control cohorts from non-European populations are still needed. The current barriers to GRS implementation within healthcare are mainly related to a lack of guidance and knowledge on integration with other biomarkers and clinical variables. Once these limitations are addressed, there is huge potential for 'test and treat' approaches to be used to tailor care for individuals with type 1 diabetes.

Associations between sleep health and grey matter volume in the UK Biobank cohort (n = 33 356).

As suggested by previous research, sleep health is assumed to be a key determinant of future morbidity and mortality. In line with this, recent studies have found that poor sleep is associated with impaired cognitive function. However, to date, little is known about brain structural abnormalities underlying this association. Although recent findings link sleep health deficits to specific alterations in grey matter volume, evidence remains inconsistent and reliant on small sample sizes. Addressing this problem, the current preregistered study investigated associations between sleep health and grey matter volume (139 imaging-derived phenotypes) in the UK Biobank cohort (33 356 participants). Drawing on a large sample size and consistent data acquisition, sleep duration, insomnia symptoms, daytime sleepiness, chronotype, sleep medication and sleep apnoea were examined. Our main analyses revealed that long sleep duration was systematically associated with larger grey matter volume of basal ganglia substructures. Insomnia symptoms, sleep medication and sleep apnoea were not associated with any of the 139 imaging-derived phenotypes. Short sleep duration, daytime sleepiness as well as late and early chronotype were associated with solitary imaging-derived phenotypes (no recognizable pattern, small effect sizes). To our knowledge, this is the largest study to test associations between sleep health and grey matter volume. Clinical implications of the association between long sleep duration and larger grey matter volume of basal ganglia are discussed. Insomnia symptoms as operationalized in the UK Biobank do not translate into grey matter volume findings.

Type 1 Diabetes Genetic Risk Score Differentiates Subgroups of Ketosis-Prone Diabetes.

OBJECTIVE: To determine whether genetic risk for type 1 diabetes (T1D) differentiates the four Aß subgroups of ketosis-prone diabetes (KPD), where A+ and A- define the presence or absence of islet autoantibodies and ß+ and ß- define the presence or absence of ß-cell function. RESEARCH DESIGN AND METHODS: We compared T1D genetic risk scores (GRS) of patients with KPD across subgroups, race/ethnicity, ß-cell function, and glycemia. RESULTS: Among 426 patients with KPD (54% Hispanic, 31% African American, 11% White), rank order of GRS was A+ß- > A+ß+ = A-ß- > A-ß+. GRS of A+ß- KPD was lower than that of a T1D cohort, and GRS of A-ß+ KPD was higher than that of a type 2 diabetes cohort. GRS was lowest among African American patients, with a similar distribution across KPD subgroups. CONCLUSIONS: T1D genetic risk delineates etiologic differences among KPD subgroups. Patients with A+ß- KPD have the highest and those with A-ß+ KPD the lowest GRS.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.

Identification and analysis of individuals who deviate from their genetically-predicted phenotype.

Findings from genome-wide association studies have facilitated the generation of genetic predictors for many common human phenotypes. Stratifying individuals misaligned to a genetic predictor based on common variants may be important for follow-up studies that aim to identify alternative causal factors. Using genome-wide imputed genetic data, we aimed to classify 158,951 unrelated individuals from the UK Biobank as either concordant or deviating from two well-measured phenotypes. We first applied our methods to standing height: our primary analysis classified 244 individuals (0.15%) as misaligned to their genetically predicted height. We show that these individuals are enriched for self-reporting being shorter or taller than average at age 10, diagnosed congenital malformations, and rare loss-of-function variants in genes previously catalogued as causal for growth disorders. Secondly, we apply our methods to LDL cholesterol. We classified 156 (0.12%) individuals as misaligned to their genetically predicted LDL cholesterol and show that these individuals were enriched for both clinically actionable cardiovascular risk factors and rare genetic variants in genes previously shown to be involved in metabolic processes. Individuals whose LDL-C was higher than expected based on the genetic predictor were also at higher risk of developing coronary artery disease and type-two diabetes, even after adjustment for measured LDL-C, BMI and age, suggesting upward deviation from genetically predicted LDL-C is indicative of generally poor health. Our results remained broadly consistent when performing sensitivity analysis based on a variety of parametric and non-parametric methods to define individuals deviating from polygenic expectation. Our analyses demonstrate the potential importance of quantitatively identifying individuals for further follow-up based on deviation from genetic predictions. Author Summary: Human genetics is becoming increasingly useful to help predict human traits across a population owing to findings from large-scale genetic association studies and advances in the power of genetic predictors. This provides an opportunity to potentially identify individuals that deviate from genetic predictions for a common phenotype under investigation. For example, an individual may be genetically predicted to be tall, but be shorter than expected. It is potentially important to identify individuals who deviate from genetic predictions as this can facilitate further follow-up to assess likely causes. Using 158,951 unrelated individuals from the UK Biobank, with height and LDL cholesterol, as exemplar traits, we demonstrate that approximately 0.15% & 0.12% of individuals deviate from their genetically predicted phenotypes respectively. We observed these individuals to be enriched for a range of rare clinical diagnoses, as well as rare genetic factors that may be causal. Our analyses also demonstrate several methods for detecting individuals who deviate from genetic predictions that can be applied to a range of continuous human phenotypes.