Inborn errors of immunity: an expanding universe of disease and genetic architecture.

Inborn errors of immunity (IEIs) are generally considered to be rare monogenic disorders of the immune system that cause immunodeficiency, autoinflammation, autoimmunity, allergy and/or cancer. Here, we discuss evidence that IEIs need not be rare disorders or exclusively affect the immune system. Namely, an increasing number of patients with IEIs present with severe dysregulations of the central nervous, digestive, renal or pulmonary systems. Current challenges in the diagnosis of IEIs that result from the segregated practice of specialized medicine could thus be mitigated, in part, by immunogenetic approaches. Starting with a brief historical overview of IEIs, we then discuss the technological advances that are facilitating the immunogenetic study of IEIs, progress in understanding disease penetrance in IEIs, the expanding universe of IEIs affecting distal organ systems and the future of genetic, biochemical and medical discoveries in this field.

Biochemical Basis for Dominant Inheritance, Variable Penetrance, and Maternal Effects in RBP4 Congenital Eye Disease.

Gestational vitamin A (retinol) deficiency poses a risk for ocular birth defects and blindness. We identified missense mutations in RBP4, encoding serum retinol binding protein, in three families with eye malformations of differing severity, including bilateral anophthalmia. The mutant phenotypes exhibit dominant inheritance, but incomplete penetrance. Maternal transmission significantly increases the probability of phenotypic expression. RBP normally delivers retinol from hepatic stores to peripheral tissues, including the placenta and fetal eye. The disease mutations greatly reduce retinol binding to RBP, yet paradoxically increase the affinity of RBP for its cell surface receptor, STRA6. By occupying STRA6 nonproductively, the dominant-negative proteins disrupt vitamin A delivery from wild-type proteins within the fetus, but also, in the case of maternal transmission, at the placenta. These findings establish a previously uncharacterized mode of maternal inheritance, distinct from imprinting and oocyte-derived mRNA, and define a group of hereditary disorders plausibly modulated by dietary vitamin A.

Jagged 1 Rescues the Duchenne Muscular Dystrophy Phenotype.

Duchenne muscular dystrophy (DMD), caused by mutations at the dystrophin gene, is the most common form of muscular dystrophy. There is no cure for DMD and current therapeutic approaches to restore dystrophin expression are only partially effective. The absence of dystrophin in muscle results in dysregulation of signaling pathways, which could be targets for disease therapy and drug discovery. Previously, we identified two exceptional Golden Retriever muscular dystrophy (GRMD) dogs that are mildly affected, have functional muscle, and normal lifespan despite the complete absence of dystrophin. Now, our data on linkage, whole-genome sequencing, and transcriptome analyses of these dogs compared to severely affected GRMD and control animals reveals that increased expression of Jagged1 gene, a known regulator of the Notch signaling pathway, is a hallmark of the mild phenotype. Functional analyses demonstrate that Jagged1 overexpression ameliorates the dystrophic phenotype, suggesting that Jagged1 may represent a target for DMD therapy in a dystrophin-independent manner. PAPERCLIP.

Affinity-optimizing enhancer variants disrupt development.

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease1-3. The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb4. Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly4-6. However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies.

A nondegenerate code of deleterious variants in Mendelian loci contributes to complex disease risk.

Although countless highly penetrant variants have been associated with Mendelian disorders, the genetic etiologies underlying complex diseases remain largely unresolved. By mining the medical records of over 110 million patients, we examine the extent to which Mendelian variation contributes to complex disease risk. We detect thousands of associations between Mendelian and complex diseases, revealing a nondegenerate, phenotypic code that links each complex disorder to a unique collection of Mendelian loci. Using genome-wide association results, we demonstrate that common variants associated with complex diseases are enriched in the genes indicated by this "Mendelian code." Finally, we detect hundreds of comorbidity associations among Mendelian disorders, and we use probabilistic genetic modeling to demonstrate that Mendelian variants likely contribute nonadditively to the risk for a subset of complex diseases. Overall, this study illustrates a complementary approach for mapping complex disease loci and provides unique predictions concerning the etiologies of specific diseases.

A mechanism for gene-environment interaction in the etiology of congenital scoliosis.

Congenital scoliosis, a lateral curvature of the spine caused by vertebral defects, occurs in approximately 1 in 1,000 live births. Here we demonstrate that haploinsufficiency of Notch signaling pathway genes in humans can cause this congenital abnormality. We also show that in a mouse model, the combination of this genetic risk factor with an environmental condition (short-term gestational hypoxia) significantly increases the penetrance and severity of vertebral defects. We demonstrate that hypoxia disrupts FGF signaling, leading to a temporary failure of embryonic somitogenesis. Our results potentially provide a mechanism for the genesis of a host of common sporadic congenital abnormalities through gene-environment interaction.

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.

Low disease risk and penetrance in Leber hereditary optic neuropathy.

The risk of Leber hereditary optic neuropathy (LHON) has largely been extrapolated from disease cohorts, which underestimate the population prevalence of pathogenic primary LHON variants as a result of incomplete disease penetrance. Understanding the true population prevalence of primary LHON variants, alongside the rate of clinical disease, provides a better understanding of disease risk and variant penetrance. We identified pathogenic primary LHON variants in whole-genome sequencing data of a well-characterized population-based control cohort and found that the prevalence is far greater than previously estimated, as it occurs in approximately 1 in 800 individuals. Accordingly, we were able to more accurately estimate population risk and disease penetrance in LHON variant carriers, validating our findings by using other large control datasets. These findings will inform accurate counseling in relation to the risk of vision loss in LHON variant carriers and disease manifestation in their family. This Matters Arising paper is in response to Lopez Sanchez et al. (2021), published in The American Journal of Human Genetics. See also the response by Mackey et al. (2022), published in this issue.

Is the disease risk and penetrance in Leber hereditary optic neuropathy actually low?

Pedigree analysis showed that a large proportion of Leber hereditary optic neuropathy (LHON) family members who carry a mitochondrial risk variant never lose vision. Mitochondrial haplotype appears to be a major factor influencing the risk of vision loss from LHON. Mitochondrial variants, including m.14484T>C and m.11778G>A, have been added to gene arrays, and thus many patients and research participants are tested for LHON mutations. Analysis of the UK Biobank and Australian cohort studies found more than 1 in 1,000 people in the general population carry either the m.14484T>C or the m.11778G>A LHON variant. None of the subset of carriers examined had visual acuity at 20/200 or worse, suggesting a very low penetrance of LHON. Haplogroup analysis of m.14484T>C carriers showed a high rate of haplogroup U subclades, previously shown to have low penetrance in pedigrees. Penetrance calculations of the general population are lower than pedigree calculations, most likely because of modifier genetic factors. This Matters Arising Response paper addresses the Watson et al. (2022) Matters Arising paper, published concurrently in The American Journal of Human Genetics.

The penetrance of rare variants in cardiomyopathy-associated genes: A cross-sectional approach to estimating penetrance for secondary findings.

Understanding the penetrance of pathogenic variants identified as secondary findings (SFs) is of paramount importance with the growing availability of genetic testing. We estimated penetrance through large-scale analyses of individuals referred for diagnostic sequencing for hypertrophic cardiomyopathy (HCM; 10,400 affected individuals, 1,332 variants) and dilated cardiomyopathy (DCM; 2,564 affected individuals, 663 variants), using a cross-sectional approach comparing allele frequencies against reference populations (293,226 participants from UK Biobank and gnomAD). We generated updated prevalence estimates for HCM (1:543) and DCM (1:220). In aggregate, the penetrance by late adulthood of rare, pathogenic variants (23% for HCM, 35% for DCM) and likely pathogenic variants (7% for HCM, 10% for DCM) was substantial for dominant cardiomyopathy (CM). Penetrance was significantly higher for variant subgroups annotated as loss of function or ultra-rare and for males compared to females for variants in HCM-associated genes. We estimated variant-specific penetrance for 316 recurrent variants most likely to be identified as SFs (found in 51% of HCM- and 17% of DCM-affected individuals). 49 variants were observed at least ten times (14% of affected individuals) in HCM-associated genes. Median penetrance was 14.6% (±14.4% SD). We explore estimates of penetrance by age, sex, and ancestry and simulate the impact of including future cohorts. This dataset reports penetrance of individual variants at scale and will inform the management of individuals undergoing genetic screening for SFs. While most variants had low penetrance and the costs and harms of screening are unclear, some individuals with highly penetrant variants may benefit from SFs.

Inherited mutations affecting the SRCAP complex are central in moderate-penetrance predisposition to uterine leiomyomas.

Uterine leiomyomas (ULs) are benign smooth muscle tumors that are common in premenopausal women. Somatic alterations in MED12, HMGA2, FH, genes encoding subunits of the SRCAP complex, and genes involved in Cullin 3-RING E3 ligase neddylation are mutually exclusive UL drivers. Established predisposition genes explain only partially the estimated heritability of leiomyomas. Here, we examined loss-of-function variants across 18,899 genes in a cohort of 233,614 White European women, revealing variants in four genes encoding SRCAP complex subunits (YEATS4, ZNHIT1, DMAP1, and ACTL6A) with a significant association to ULs, and YEATS4 and ZNHIT1 strikingly rank first and second, respectively. Positive mutation status was also associated with younger age at diagnosis and hysterectomy. Moderate-penetrance UL risk was largely attributed to rare non-synonymous mutations affecting the SRCAP complex. To examine this disease phenotype more closely, we set out to identify inherited mutations affecting the SRCAP complex in our in-house sample collection of Finnish individuals with ULs (n = 860). We detected one individual with an ACTL6A splice-site mutation, two individuals with a YEATS4 missense mutation, and four individuals with DMAP1 mutations: one splice-site, one nonsense, and two missense variants. These individuals had large and/or multiple ULs, were often diagnosed at an early age, and many had family history of ULs. When a somatic second hit was found, ACTL6A and DMAP1 were silenced in tumors by somatic mutation and YEATS4 by promoter hypermethylation. Decreased H2A.Z staining was observed in the tumors, providing further evidence for the pathogenic nature of the germline mutations. Our results establish inactivation of genes encoding SRCAP complex subunits as a central contributor to moderate-penetrance UL predisposition.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Evidence for 28 genetic disorders discovered by combining healthcare and research data.

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

A theoretical perspective on Waddington's genetic assimilation experiments.

Genetic assimilation is the process by which a phenotype that is initially induced by an environmental stimulus becomes stably inherited in the absence of the stimulus after a few generations of selection. While the concept has attracted much debate after being introduced by C. H. Waddington 70 y ago, there have been few experiments to quantitatively characterize the phenomenon. Here, we revisit and organize the results of Waddington's original experiments and follow-up studies that attempted to replicate his results. We then present a theoretical model to illustrate the process of genetic assimilation and highlight several aspects that we think require further quantitative studies, including the gradual increase of penetrance, the statistics of delay in assimilation, and the frequency of unviability during selection. Our model captures Waddington's picture of developmental paths in a canalized landscape using a stochastic dynamical system with alternative trajectories that can be controlled by either external signals or internal variables. It also reconciles two descriptions of the phenomenon-Waddington's, expressed in terms of an individual organism's developmental paths, and that of Bateman in terms of the population distribution crossing a hypothetical threshold. Our results provide theoretical insight into the concepts of canalization, phenotypic plasticity, and genetic assimilation.

Meta-Analysis of Penetrance and Systematic Review on Transition to Disease in Genetic Hypertrophic Cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is characterized by unexplained left ventricular hypertrophy and is classically caused by pathogenic or likely pathogenic variants (P/LP) in genes encoding sarcomere proteins. Not all subclinical variant carriers will manifest clinically overt disease because penetrance (proportion of sarcomere or sarcomere-related P/LP variant carriers who develop disease) is variable, age dependent, and not reliably predicted. METHODS: A systematic search of the literature was performed. We used random-effects generalized linear mixed model meta-analyses to contrast the cross-sectional prevalence and penetrance of sarcomere or sarcomere-related genes in 2 different contexts: clinically-based studies on patients and families with HCM versus population or community-based studies. Longitudinal family/clinical studies were additionally analyzed to investigate the rate of phenotypic conversion from subclinical to overt HCM during follow-up. RESULTS: In total, 455 full-text manuscripts and articles were assessed. In family/clinical studies, the prevalence of sarcomere variants in patients diagnosed with HCM was 34%. The penetrance across all genes in nonproband relatives carrying P/LP variants identified during cascade screening was 57% (95% CI, 52%-63%), and the mean age at HCM diagnosis was 38 years (95% CI, 36%-40%). Penetrance varied from ≈32% for MYL3 (myosin light chain 3) to ≈55% for MYBPC3 (myosin-binding protein C3), ≈60% for TNNT2 (troponin T2) and TNNI3 (troponin I3), and ≈65% for MYH7 (myosin heavy chain 7). Population-based genetic studies demonstrate that P/LP sarcomere variants are present in the background population but at a low prevalence of <1%. The penetrance of HCM in incidentally identified P/LP variant carriers was also substantially lower at ≈11%, ranging from 0% in Atherosclerosis Risk in Communities to 18% in UK Biobank. In longitudinal family studies, the pooled phenotypic conversion across all genes was 15% over an average of ≈8 years of follow-up, starting from a mean of ≈16 years of age. However, short-term gene-specific phenotypic conversion varied between ≈12% for MYBPC3 and ≈23% for MYH7. CONCLUSIONS: The penetrance of P/LP variants is highly variable and influenced by currently undefined and context-dependent genetic and environmental factors. Additional longitudinal studies are needed to improve our understanding of true lifetime penetrance in families and in the community and to identify drivers of the transition from subclinical to overt HCM.

Reclassification of a likely pathogenic Dutch founder variant in KCNH2; implications of reduced penetrance.

BACKGROUND: Variants in KCNH2, encoding the human ether a-go-go (hERG) channel that is responsible for the rapid component of the cardiac delayed rectifier K+ current (IKr), are causal to long QT syndrome type 2 (LQTS2). We identified eight index patients with a new variant of unknown significance (VUS), KCNH2:c.2717C > T:p.(Ser906Leu). We aimed to elucidate the biophysiological effect of this variant, to enable reclassification and consequent clinical decision-making. METHODS: A genotype-phenotype overview of the patients and relatives was created. The biophysiological effects were assessed independently by manual-, and automated calibrated patch clamp. HEK293a cells expressing (i) wild-type (WT) KCNH2, (ii) KCNH2-p.S906L alone (homozygous, Hm) or (iii) KCNH2-p.S906L in combination with WT (1:1) (heterozygous, Hz) were used for manual patching. Automated patch clamp measured the variants function against known benign and pathogenic variants, using Flp-In T-rex HEK293 KCNH2-variant cell lines. RESULTS: Incomplete penetrance of LQTS2 in KCNH2:p.(Ser906Leu) carriers was observed. In addition, some patients were heterozygous for other VUSs in CACNA1C, PKP2, RYR2 or AKAP9. The phenotype of carriers of KCNH2:p.(Ser906Leu) ranged from asymptomatic to life-threatening arrhythmic events. Manual patch clamp showed a reduced current density by 69.8 and 60.4% in KCNH2-p.S906L-Hm and KCNH2-p.S906L-Hz, respectively. The time constant of activation was significantly increased with 80.1% in KCNH2-p.S906L-Hm compared with KCNH2-WT. Assessment of KCNH2-p.S906L-Hz by calibrated automatic patch clamp assay showed a reduction in current density by 35.6%. CONCLUSION: The reduced current density in the KCNH2-p.S906L-Hz indicates a moderate loss-of-function. Combined with the reduced penetrance and variable phenotype, we conclude that KCNH2:p.(Ser906Leu) is a low penetrant likely pathogenic variant for LQTS2.

Rare genetic variants in genes and loci linked to dominant monogenic developmental disorders cause milder related phenotypes in the general population.

Many rare monogenic diseases are known to be caused by deleterious variants in thousands of genes, however the same variants can also be found in people without the associated clinical phenotypes. The penetrance of these monogenic variants is generally unknown in the wider population, as they are typically identified in small clinical cohorts of affected individuals and families with highly penetrant variants. Here, we investigated the phenotypic effect of rare, potentially deleterious variants in genes and loci where similar variants are known to cause monogenic developmental disorders (DDs) in a large population cohort. We used UK Biobank to investigate phenotypes associated with rare protein-truncating and missense variants in 599 monoallelic DDG2P genes by using whole-exome-sequencing data from ∼200,000 individuals and rare copy-number variants overlapping known DD loci by using SNP-array data from ∼500,000 individuals. We found that individuals with these likely deleterious variants had a mild DD-related phenotype, including lower fluid intelligence, slower reaction times, lower numeric memory scores, and longer pairs matching times compared to the rest of the UK Biobank cohort. They were also shorter, had a higher BMI, and had significant socioeconomic disadvantages: they were less likely to be employed or be able to work and had a lower income and higher deprivation index. Our findings suggest that many genes routinely tested within pediatric genetics have deleterious variants with intermediate penetrance that may cause lifelong sub-clinical phenotypes in the general adult population.

Reduced penetrance of MODY-associated HNF1A/HNF4A variants but not GCK variants in clinically unselected cohorts.

The true prevalence and penetrance of monogenic disease variants are often not known because of clinical-referral ascertainment bias. We comprehensively assess the penetrance and prevalence of pathogenic variants in HNF1A, HNF4A, and GCK that account for >80% of monogenic diabetes. We analyzed clinical and genetic data from 1,742 clinically referred probands, 2,194 family members, clinically unselected individuals from a US health system-based cohort (n = 132,194), and a UK population-based cohort (n = 198,748). We show that one in 1,500 individuals harbor a pathogenic variant in one of these genes. The penetrance of diabetes for HNF1A and HNF4A pathogenic variants was substantially lower in the clinically unselected individuals compared to clinically referred probands and was dependent on the setting (32% in the population, 49% in the health system cohort, 86% in a family member, and 98% in probands for HNF1A). The relative risk of diabetes was similar across the clinically unselected cohorts highlighting the role of environment/other genetic factors. Surprisingly, the penetrance of pathogenic GCK variants was similar across all cohorts (89%-97%). We highlight that pathogenic variants in HNF1A, HNF4A, and GCK are not ultra-rare in the population. For HNF1A and HNF4A, we need to tailor genetic interpretation and counseling based on the setting in which a pathogenic monogenic variant was identified. GCK is an exception with near-complete penetrance in all settings. This along with the clinical implication of diagnosis makes it an excellent candidate for the American College of Medical Genetics secondary gene list.

Segregation analysis of 17,425 population-based breast cancer families: Evidence for genetic susceptibility and risk prediction.

Rare pathogenic variants in known breast cancer-susceptibility genes and known common susceptibility variants do not fully explain the familial aggregation of breast cancer. To investigate plausible genetic models for the residual familial aggregation, we studied 17,425 families ascertained through population-based probands, 86% of whom were screened for pathogenic variants in BRCA1, BRCA2, PALB2, CHEK2, ATM, and TP53 via gene-panel sequencing. We conducted complex segregation analyses and fitted genetic models in which breast cancer incidence depended on the effects of known susceptibility genes and other unidentified major genes and a normally distributed polygenic component. The proportion of familial variance explained by the six genes was 46% at age 20-29 years and decreased steadily with age thereafter. After allowing for these genes, the best fitting model for the residual familial variance included a recessive risk component with a combined genotype frequency of 1.7% (95% CI: 0.3%-5.4%) and a penetrance to age 80 years of 69% (95% CI: 38%-95%) for homozygotes, which may reflect the combined effects of multiple variants acting in a recessive manner, and a polygenic variance of 1.27 (95% CI: 0.94%-1.65), which did not vary with age. The proportion of the residual familial variance explained by the recessive risk component was 40% at age 20-29 years and decreased with age thereafter. The model predicted age-specific familial relative risks consistent with those observed by large epidemiological studies. The findings have implications for strategies to identify new breast cancer-susceptibility genes and improve disease-risk prediction, especially at a young age.

A spectrum of recessiveness among Mendelian disease variants in UK Biobank.

Recent work has found increasing evidence of mitigated, incompletely penetrant phenotypes in heterozygous carriers of recessive Mendelian disease variants. We leveraged whole-exome imputation within the full UK Biobank cohort (n ∼ 500K) to extend such analyses to 3,475 rare variants curated from ClinVar and OMIM. Testing these variants for association with 58 quantitative traits yielded 102 significant associations involving variants previously implicated in 34 different diseases. Notable examples included a POR missense variant implicated in Antley-Bixler syndrome that associated with a 1.76 (SE 0.27) cm increase in height and an ABCA3 missense variant implicated in interstitial lung disease that associated with reduced FEV1/FVC ratio. Association analyses with 1,134 disease traits yielded five additional variant-disease associations. We also observed contrasting levels of recessiveness between two more-common, classical Mendelian diseases. Carriers of cystic fibrosis variants exhibited increased risk of several mitigated disease phenotypes, whereas carriers of spinal muscular atrophy alleles showed no evidence of altered phenotypes. Incomplete penetrance of cystic fibrosis carrier phenotypes did not appear to be mediated by common allelic variation on the functional haplotype. Our results show that many disease-associated recessive variants can produce mitigated phenotypes in heterozygous carriers and motivate further work exploring penetrance mechanisms.