Bi-allelic variants in MYH3 cause recessively-inherited arthrogryposis.

Arthrogryposis is a clinical feature defined by congenital joint contractures in two or more different body areas which occurs in between 1/3000 and 1/5000 live births. Variants in multiple genes have been associated with distal arthrogryposis syndromes. Heterozygous variants in MYH3 have been identified to cause the dominantly-inherited distal arthrogryposis conditions, Freeman-Sheldon syndrome, Sheldon-Hall syndrome, and multiple pterygium syndrome. In contrast, MYH3 variants underlie both dominantly and recessively inherited Contractures, Pterygia, and Spondylocarpotarsal Fusion syndromes (CPSFS) which are characterized by extensive bony abnormalities in addition to congenital contractures. Here we report two affected sibs with distal arthrogryposis born to unaffected, distantly related parents. Sequencing revealed that both sibs were homozygous for two ultra-rare MYH3 variants, c.3445G>A (p.Glu1149Lys) and c.4760T>C (p.Leu1587Pro). Sequencing and deletion/duplication analysis of 169 other arthrogryposis genes yielded no other compelling candidate variants. This is the first report of biallelic variants in MYH3 being implicated in a distal arthrogryposis phenotype without the additional features of CPSFS. Thus, akin to CPSFS, both dominant and recessively inherited distal arthrogryposis can be caused by variants in MYH3.

A corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery.

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

Biallelic variants in CSMD1 are implicated in a neurodevelopmental disorder with intellectual disability and variable cortical malformations.

CSMD1 (Cub and Sushi Multiple Domains 1) is a well-recognized regulator of the complement cascade, an important component of the innate immune response. CSMD1 is highly expressed in the central nervous system (CNS) where emergent functions of the complement pathway modulate neural development and synaptic activity. While a genetic risk factor for neuropsychiatric disorders, the role of CSMD1 in neurodevelopmental disorders is unclear. Through international variant sharing, we identified inherited biallelic CSMD1 variants in eight individuals from six families of diverse ancestry who present with global developmental delay, intellectual disability, microcephaly, and polymicrogyria. We modeled CSMD1 loss-of-function (LOF) pathogenesis in early-stage forebrain organoids differentiated from CSMD1 knockout human embryonic stem cells (hESCs). We show that CSMD1 is necessary for neuroepithelial cytoarchitecture and synchronous differentiation. In summary, we identified a critical role for CSMD1 in brain development and biallelic CSMD1 variants as the molecular basis of a previously undefined neurodevelopmental disorder.

STR mutations on chromosome 15q cause thyrotropin resistance by activating a primate-specific enhancer of MIR7-2/MIR1179.

Thyrotropin (TSH) is the master regulator of thyroid gland growth and function. Resistance to TSH (RTSH) describes conditions with reduced sensitivity to TSH. Dominantly inherited RTSH has been linked to a locus on chromosome 15q, but its genetic basis has remained elusive. Here we show that non-coding mutations in a (TTTG)4 short tandem repeat (STR) underlie dominantly inherited RTSH in all 82 affected participants from 12 unrelated families. The STR is contained in a primate-specific Alu retrotransposon with thyroid-specific cis-regulatory chromatin features. Fiber-seq and RNA-seq studies revealed that the mutant STR activates a thyroid-specific enhancer cluster, leading to haplotype-specific upregulation of the bicistronic MIR7-2/MIR1179 locus 35 kb downstream and overexpression of its microRNA products in the participants' thyrocytes. An imbalance in signaling pathways targeted by these micro-RNAs provides a working model for this cause of RTSH. This finding broadens our current knowledge of genetic defects altering pituitary-thyroid feedback regulation.

The clinical utility of sequencing the entirety of CFTR.

BACKGROUND: Cystic fibrosis (CF) is caused by deleterious variants in each CFTR gene. We investigated the utility of whole-gene CFTR sequencing when fewer than two pathogenic or likely pathogenic (P/LP) variants were detected by conventional testing (sequencing of exons and flanking introns) of CFTR. METHODS: Individuals with features of CF and a CF-diagnostic sweat chloride concentration with zero or one P/LP variants identified by conventional testing enrolled in the CF Mutation Analysis Program (MAP) underwent whole-gene CFTR sequencing. Replication was performed on individuals enrolled in the CF Genome Project (CFGP), followed by phenotype review and interrogation of other genes. RESULTS: Whole-gene sequencing identified a second P/LP variant in 20/43 MAP enrollees (47 %) and 10/22 CFGP enrollees (45 %) who had one P/LP variant after conventional testing. No P/LP variants were detected when conventional testing was negative (MAP: n = 43; CFGP: n = 13). Genome-wide analysis was unable to find an alternative etiology in CFGP participants with fewer than two P/LP CFTR variants and CF could not be confirmed in 91 % following phenotype re-review. CONCLUSIONS: Whole-gene CFTR analysis is beneficial in individuals with one previously-identified P/LP variant and a CF-diagnostic sweat chloride. Negative conventional CFTR testing indicates that the phenotype should be re-evaluated.

CFAP47 is a novel causative gene implicated in X-linked polycystic kidney disease.

Autosomal dominant polycystic kidney disease (ADPKD) is a well-described condition in which ~80% of cases have a genetic explanation, while the genetic basis of sporadic cystic kidney disease in adults remains unclear in ~30% of cases. This study aimed to identify novel genes associated with polycystic kidney disease (PKD) in patients with sporadic cystic kidney disease in which a clear genetic change was not identified in established genes. A next-generation sequencing panel analyzed known genes related to renal cysts in 118 sporadic cases, followed by whole-genome sequencing on 47 unrelated individuals without identified candidate variants. Three male patients were found to have rare missense variants in the X-linked gene Cilia And Flagella Associated Protein 47 (CFAP47). CFAP47 was expressed in primary cilia of human renal tubules, and knockout mice exhibited vacuolation of tubular cells and tubular dilation, providing evidence that CFAP47 is a causative gene involved in cyst formation. This discovery of CFAP47 as a newly identified gene associated with PKD, displaying X-linked inheritance, emphasizes the need for further cases to understand the role of CFAP47 in PKD.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Whole Exome Sequencing Uncovers the Genetic Complexity of Bicuspid Aortic Valve in Families with Early Onset Complications.

Bicuspid Aortic Valve (BAV) is the most common adult congenital heart lesion with an estimated population prevalence of 1%. We hypothesize that early onset complications of BAV (EBAV) are driven by specific impactful genetic variants. We analyzed whole exome sequences (WES) to identify rare coding variants that contribute to BAV disease in 215 EBAV families. Predicted pathogenic variants of causal genes were present in 111 EBAV families (51% of total), including genes that cause BAV (8%) or heritable thoracic aortic disease (HTAD, 17%). After appropriate filtration, we also identified 93 variants in 26 novel genes that are associated with autosomal dominant congenital heart phenotypes, including recurrent deleterious variation of FBN2, MYH6, channelopathy genes, and type 1 and 5 collagen genes. These findings confirm our hypothesis that unique rare genetic variants contribute to early onset complications of BAV disease.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing.

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing (ES), the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare disease. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery which should in turn increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints, and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks like Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Lisch Epithelial Corneal Dystrophy Is Caused by Heterozygous Loss-of-Function Variants in MCOLN1.

PURPOSE: To report the genetic etiology of Lisch epithelial corneal dystrophy (LECD). DESIGN: Multicenter cohort study. METHODS: A discovery cohort of 27 individuals with LECD from 17 families, including 7 affected members from the original LECD family, 6 patients from 2 new families and 14 simplex cases, was recruited. A cohort of 6 individuals carrying a pathogenic MCOLN1 (mucolipin 1) variant was reviewed for signs of LECD. Next-generation sequencing or targeted Sanger sequencing were used in all patients to identify pathogenic or likely pathogenic variants and penetrance of variants. RESULTS: Nine rare heterozygous MCOLN1 variants were identified in 23 of 27 affected individuals from 13 families. The truncating nature of 7 variants and functional testing of 1 missense variant indicated that they result in MCOLN1 haploinsufficiency. Importantly, in the homozygous and compound-heterozygous state, 4 of 9 LECD-associated variants cause the rare lysosomal storage disorder mucolipidosis IV (MLIV). Autosomal recessive MLIV is a systemic disease and comprises neurodegeneration as well as corneal opacity of infantile-onset with epithelial autofluorescent lysosomal inclusions. However, the 6 parents of 3 patients with MLIV confirmed to carry pathogenic MCOLN1 variants did not have the LECD phenotype, suggesting MCOLN1 haploinsufficiency may be associated with reduced penetrance and variable expressivity. CONCLUSIONS: MCOLN1 haploinsufficiency is the major cause of LECD. Based on the overlapping clinical features of corneal epithelial cells with autofluorescent inclusions reported in both LECD and MLIV, it is concluded that some carriers of MCOLN1 haploinsufficiency-causing variants present with LECD.

SAMHD1 compound heterozygous rare variants associated with moyamoya and mitral valve disease in the absence of other features of Aicardi-Goutières syndrome.

Aicardi-Goutières syndrome (AGS) is an autosomal recessive inflammatory syndrome that manifests as an early-onset encephalopathy with both neurologic and extraneurologic clinical findings. AGS has been associated with pathogenic variants in nine genes: TREX1, RNASEH2B, RNASEH2C, RNASEH2A, SAMHD1, ADAR, IFIH1, LSM11, and RNU7-1. Diagnosis is established by clinical findings (encephalopathy and acquired microcephaly, intellectual and physical impairments, dystonia, hepatosplenomegaly, sterile pyrexia, and/or chilblains), characteristic abnormalities on cranial CT (calcification of the basal ganglia and white matter) and MRI (leukodystrophic changes), or the identification of pathogenic/likely pathogenic variants in the known genes. One of the genes associated with AGS, SAMHD1, has also been associated with a spectrum of cerebrovascular diseases, including moyamoya disease (MMD). In this report, we describe a 31-year-old male referred to genetics for MMD since childhood who lacked the hallmark features of AGS patients but was found to have compound heterozygous SAMHD1 variants. He later developed mitral valve insufficiency due to recurrent chordal rupture and ultimately underwent a heart transplant at 37 years of age. Thus, these data suggest that SAMHD1 pathogenic variants can cause MMD without typical AGS symptoms and support that SAMHD1 should be assessed in MMD patients even in the absence of AGS features.

A loss of function variant in AGPAT3 underlies intellectual disability and retinitis pigmentosa (IDRP) syndrome.

Intellectual disability (ID) and retinal dystrophy (RD) are the frequently found features of multiple syndromes involving additional systemic manifestations. Here, we studied a family with four members presenting severe ID and retinitis pigmentosa (RP). Using genome wide genotyping and exome sequencing, we identified a nonsense variant c.747 C > A (p.Tyr249Ter) in exon 7 of AGPAT3 which co-segregates with the disease phenotype. Western blot analysis of overexpressed WT and mutant AGPAT3 in HEK293T cells showed the absence of AGPAT3, suggesting instability of the truncated protein. Knockdown of Agpat3 in the embryonic mouse brain caused marked deficits in neuronal migration, strongly suggesting that reduced expression of AGPAT3 affects neuronal function. Altogether, our data indicates that AGPAT3 activity is essential for neuronal functioning and loss of its activity probably causes intellectual disability and retinitis pigmentosa (IDRP) syndrome.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition.

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.