Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

An FBN1 deep intronic variant is associated with pseudoexon formation and a variable Marfan phenotype in a five generation family

Exome sequencing of genes associated with heritable thoracic aortic disease (HTAD) failed to identify a pathogenic variant in a large family with Marfan syndrome (MFS). A genome-wide linkage analysis for thoracic aortic disease identified a peak at 15q21.1, and genome sequencing identified a novel deep intronic FBN1 variant that segregated with thoracic aortic disease in the family (LOD score 2.7) and was predicted to alter splicing. RT-PCR and bulk RNA sequencing of RNA harvested from fibroblasts explanted from the affected proband revealed an insertion of a pseudoexon between exons 13 and 14 of the FBN1 transcript, predicted to lead to nonsense mediated decay (NMD). Treating the fibroblasts with an NMD inhibitor, cycloheximide, greatly improved the detection of the pseudoexon-containing transcript. Family members with the FBN1 variant had later onset aortic events and fewer MFS systemic features than typical for individuals with haploinsufficiency of FBN1. Variable penetrance of the phenotype and negative genetic testing in MFS families should raise the possibility of deep intronic FBN1 variants and the need for additional molecular studies.

Dominant-negative variant in SLC1A4 causes an autosomal dominant epilepsy syndrome

SLC1A4 is a trimeric neutral amino acid transporter essential for shuttling L-serine from astrocytes into neurons. Individuals with biallelic variants in SLC1A4 are known to have spastic tetraplegia, thin corpus callosum, and progressive microcephaly (SPATCCM) syndrome, but individuals with heterozygous variants are not thought to have disease. We identify an 8-year-old patient with global developmental delay, spasticity, epilepsy, and microcephaly who has a de novo heterozygous three amino acid duplication in SLC1A4 (L86_M88dup). We demonstrate that L86_M88dup causes a dominant-negative N-glycosylation defect of SLC1A4, which in turn reduces the plasma membrane localization of SLC1A4 and the transport rate of SLC1A4 for L-serine.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

The Genetic Landscape of Familial Pulmonary Fibrosis

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

Guidelines for genetic ancestry inference created through roundtable discussions

The use of genetic and genomic technology to infer ancestry is commonplace in a variety of contexts, particularly in biomedical research and for direct-to-consumer genetic testing. In 2013 and 2015, two roundtables engaged a diverse group of stakeholders toward the development of guidelines for inferring genetic ancestry in academia and industry. This report shares the stakeholder groups' work and provides an analysis of, commentary on, and views from the groundbreaking and sustained dialogue. We describe the engagement processes and the stakeholder groups' resulting statements and proposed guidelines. The guidelines focus on five key areas: application of genetic ancestry inference, assumptions and confidence/laboratory and statistical methods, terminology and population identifiers, impact on individuals and groups, and communication or translation of genetic ancestry inferences. We delineate the terms and limitations of the guidelines and discuss their critical role in advancing the development and implementation of best practices for inferring genetic ancestry and reporting the results. These efforts should inform both governmental regulation and self-regulation.

PRDM1 DNA-binding zinc finger domain is required for normal limb development and is disrupted in split hand/foot malformation

Split hand/foot malformation (SHFM) is a rare limb abnormality with clefting of the fingers and/or toes. For many individuals, the genetic etiology is unknown. Through whole-exome and targeted sequencing, we detected three novel variants in a gene encoding a transcription factor, PRDM1, that arose de novo in families with SHFM or segregated with the phenotype. PRDM1 is required for limb development; however, its role is not well understood and it is unclear how the PRDM1 variants affect protein function. Using transient and stable overexpression rescue experiments in zebrafish, we show that the variants disrupt the proline/serine-rich and DNA-binding zinc finger domains, resulting in a dominant-negative effect. Through gene expression assays, RNA sequencing, and CUT&RUN in isolated pectoral fin cells, we demonstrate that Prdm1a directly binds to and regulates genes required for fin induction, outgrowth and anterior/posterior patterning, such as fgfr1a, dlx5a, dlx6a and smo. Taken together, these results improve our understanding of the role of PRDM1 in the limb gene regulatory network and identified novel PRDM1 variants that link to SHFM in humans.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7 , TPM1 , and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes MYH2 , TPM2 , and TNNI2 , that encode parts of the skeletal muscle sarcomere, cause muscle diseases affecting skeletal muscle, such as the distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7 ) encoding sarcomeric proteins in which the same pathogenic variant affects both skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain mutations that also cause cardiac abnormalities. We report five families with DA due to heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 ( ACTC1 ). ACTC1 encodes a highly conserved actin that binds to myosin in both cardiac and skeletal muscle. Mutations in ACTC1 have previously been found to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition due to mutations in ACTC1 and suggests that some functions of actin, alpha, cardiac muscle 1 are shared in cardiac and skeletal muscle.

Full-length isoform sequencing for resolving the molecular basis of Charcot-Marie-Tooth 2A

Objectives

Transcript sequencing of patient derived samples has been shown to improve the diagnostic yield for solving cases of likely Mendelian disorders, yet the added benefit of full-length long-read transcript sequencing is largely unexplored.

Methods

We applied short-read and full-length isoform cDNA sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis.

Results

We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts a branch point critical for intron 6 spicing. Full-length long-read isoform cDNA sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates five distinct altered splicing transcripts. All five altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 protein levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A).

Discussion

This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Precision medicine for developmental and epileptic encephalopathies in Africa-strategies for a resource-limited setting

PURPOSE

Sub-Saharan Africa bears the highest burden of epilepsy worldwide. A presumed proportion is genetic, but this etiology is buried under the burden of infections and perinatal insults in a setting of limited awareness and few options for testing. Children with developmental and epileptic encephalopathies (DEEs) are most severely affected by this diagnostic gap in Africa, because the rate of actionable findings is highest in DEE-associated genes.

METHODS

We tested 234 genetically naive South African children diagnosed with/possible DEE using gene panels, exome sequencing, and chromosomal microarray. Statistical comparison of electroclinical features in children with and children without candidate variants was performed to identify characteristics most likely predictive of a positive genetic finding.

RESULTS

Of the 41 (of 234) children with likely/pathogenic variants, 26 had variants supporting precision therapy. Multivariate regression modeling highlighted neonatal or infantile-onset seizures and movement abnormalities as predictive of a positive genetic finding. We used this, coupled with an emphasis on precision medicine outcomes, to propose the pragmatic "Think-Genetics" strategy for early recognition of a possible genetic etiology.

CONCLUSION

Our findings emphasize the importance of an early genetic diagnosis in DEE. We designed the Think-Genetics strategy for early recognition, appropriate interim management, and genetic testing for DEE in resource-constrained settings.

Beyond the exome: what's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order and emerging technologies, such as optical genome mapping and long-read DNA or RNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to a consortium such as GREGoR, which is focused on elucidating the underlying cause of rare unsolved genetic disorders.

Erratum: Variants in PHF8 cause a spectrum of X-linked neurodevelopmental disorders and facial dysmorphology

[This corrects the article DOI: 10.1016/j.xhgg.2022.100102.].

MagicalRsq: Machine-learning-based genotype imputation quality calibration

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

Club cell secretory protein and lung function in children with cystic fibrosis

BACKGROUND

Club cell secretory protein (CC16) exerts anti-inflammatory functions in lung disease. We sought to determine the relation of serum CC16 deficits and genetic variants that control serum CC16 to lung function among children with cystic fibrosis (CF).

METHODS

We used longitudinal data from CF children (EPIC Study) with no positive cultures for Pseudomonas aeruginosa prior to enrollment. Circulating levels of CC16 and an inflammatory score (generated from CRP, SAA, calprotectin, G-CSF) were compared between participants with the lowest and highest FEV1 levels in adolescence (LLF and HLF groups, respectively; N = 130-per-group). Single nucleotide variants (SNVs) in the SCGB1A1, EHF-APIP loci were tested for association with circulating CC16 and with decline of FEV1 and FEV1/FVC% predicted levels between ages 7-16 using mixed models.

RESULTS

Compared with the HLF group, the LLF group had lower levels of CC16 (geometric means: 8.2 vs 6.5 ng/ml, respectively; p = 0.0002) and higher levels of the normalized inflammatory score (-0.21 vs 0.21, p = 0.0007). Participants in the lowest CC16 and highest inflammation tertile had the highest odds for having LLF (p<0.0001 for comparison with participants in the highest CC16 and lowest inflammation tertile). Among seven SNVs associated with circulating CC16, the top SNV rs3741240 was associated with decline of FEV1/FVC and, marginally, FEV1 (p = 0.003 and 0.025, respectively; N = 611 participants, 20,801 lung function observations).

CONCLUSIONS

Serum CC16 deficits are strongly associated with severity of CF lung disease and their effects are additive with systemic inflammation. The rs3741240 A allele is associated with low circulating CC16 and, possibly, accelerated lung function decline in CF.

Clinical, biochemical and genetic characteristics of MOGS-CDG: a rare congenital disorder of glycosylation

PURPOSE

To summarise the clinical, molecular and biochemical phenotype of mannosyl-oligosaccharide glucosidase-related congenital disorders of glycosylation (MOGS-CDG), which presents with variable clinical manifestations, and to analyse which clinical biochemical assay consistently supports diagnosis in individuals with bi-allelic variants in MOGS.

METHODS

Phenotypic characterisation was performed through an international and multicentre collaboration. Genetic testing was done by exome sequencing and targeted arrays. Biochemical assays on serum and urine were performed to delineate the biochemical signature of MOGS-CDG.

RESULTS

Clinical phenotyping revealed heterogeneity in MOGS-CDG, including neurological, immunological and skeletal phenotypes. Bi-allelic variants in MOGS were identified in 12 individuals from 11 families. The severity in each organ system was variable, without definite genotype correlation. Urine oligosaccharide analysis was consistently abnormal for all affected probands, whereas other biochemical analyses such as serum transferrin analysis was not consistently abnormal.

CONCLUSION

The clinical phenotype of MOGS-CDG includes multisystemic involvement with variable severity. Molecular analysis, combined with biochemical testing, is important for diagnosis. In MOGS-CDG, urine oligosaccharide analysis via matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry can be used as a reliable biochemical test for screening and confirmation of disease.

Variant-level matching for diagnosis and discovery: Challenges and opportunities

Here we describe MyGene2, Geno2MP, VariantMatcher, and Franklin; databases that provide variant-level information and phenotypic features to researchers, clinicians, healthcare providers and patients. Following the footsteps of the Matchmaker Exchange project that connects exome, genome, and phenotype databases at the gene level, these databases have as one goal to facilitate connection to one another using Data Connect, a standard for discovery and search of biomedical data from the Global Alliance for Genomics and Health (GA4GH).

Expanding ACTA2 genotypes with corresponding phenotypes overlapping with smooth muscle dysfunction syndrome

Pathogenic variants in ACTA2, encoding smooth muscle α-actin, predispose to thoracic aortic aneurysms and dissections. ACTA2 variants altering arginine 179 predispose to a more severe, multisystemic disease termed smooth muscle dysfunction syndrome (SMDS; OMIM 613834). Vascular complications of SMDS include patent ductus arteriosus (PDA) or aortopulmonary window, early-onset thoracic aortic disease (TAD), moyamoya-like cerebrovascular disease, and primary pulmonary hypertension. Patients also have dysfunction of other smooth muscle-dependent systems, including congenital mydriasis, hypotonic bladder, and gut hypoperistalsis. Here, we describe five patients with novel heterozygous ACTA2 missense variants, p.Arg179Gly, p.Met46Arg, p.Thr204Ile, p.Arg39Cys, and p.Ile66Asn, who have clinical complications that align or overlap with SMDS. Patients with the ACTA2 p.Arg179Gly and p.Thr204Ile variants display classic features of SMDS. The patient with the ACTA2 p.Met46Arg variant exhibits exclusively vascular complications of SMDS, including early-onset TAD, PDA, and moyamoya-like cerebrovascular disease. The patient with the ACTA2 p.Ile66Asn variant has an unusual vascular complication, a large fusiform internal carotid artery aneurysm. The patient with the ACTA2 p.Arg39Cys variant has pulmonary, gastrointestinal, and genitourinary complications of SMDS but no vascular manifestations. Identifying pathogenic ACTA2 variants associated with features of SMDS is critical for aggressive surveillance and management of vascular and nonvascular complications and delineating the molecular pathogenesis of SMDS.

SCN1A Variants as the Underlying Cause of Genetic Epilepsy with Febrile Seizures Plus in Two Multi-Generational Colombian Families

Genetic epilepsy with febrile seizures plus (GEFS+) is an autosomal dominant disorder with febrile or afebrile seizures that exhibits phenotypic variability. Only a few variants in SCN1A have been previously characterized for GEFS+, in Latin American populations where studies on the genetic and phenotypic spectrum of GEFS+ are scarce. We evaluated members in two multi-generational Colombian Paisa families whose affected members present with classic GEFS+. Exome and Sanger sequencing were used to detect the causal variants in these families. In each of these families, we identified variants in SCN1A causing GEFS+ with incomplete penetrance. In Family 047, we identified a heterozygous variant (c.3530C > G; p.(Pro1177Arg)) that segregates with GEFS+ in 15 affected individuals. In Family 167, we identified a previously unreported variant (c.725A > G; p.(Gln242Arg)) that segregates with the disease in a family with four affected members. Both variants are located in a cytoplasmic loop region in SCN1A and based on our findings the variants are classified as pathogenic and likely pathogenic, respectively. Our results expand the genotypic and phenotypic spectrum associated with SCN1A variants and will aid in improving molecular diagnostics and counseling in Latin American and other populations.

Leveraging TOPMed imputation server and constructing a cohort-specific imputation reference panel to enhance genotype imputation among cystic fibrosis patients

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Centers for Mendelian Genomics: A decade of facilitating gene discovery

PURPOSE

Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration.

METHODS

Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution.

RESULTS

We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher.

CONCLUSION

The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Exome sequencing identifies variants in infants with sacral agenesis

BACKGROUND

Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA.

METHODS

Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced.

RESULTS

Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children.

CONCLUSIONS

To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies.

Variants in PHF8 cause a spectrum of X-linked neurodevelopmental disorders and facial dysmorphology

Loss-of-function variants in PHD Finger Protein 8 (PHF8) cause Siderius X-linked intellectual disability (ID) syndrome, hereafter called PHF8-XLID. PHF8 is a histone demethylase that is important for epigenetic regulation of gene expression. PHF8-XLID is an under-characterized disorder with only five previous reports describing different PHF8 predicted loss-of-function variants in eight individuals. Features of PHF8-XLID include ID and craniofacial dysmorphology. In this report we present 16 additional individuals with PHF8-XLID from 11 different families of diverse ancestry. We also present five individuals from four different families who have ID and a variant of unknown significance in PHF8 with no other explanatory variant in another gene. All affected individuals exhibited developmental delay and all but two had borderline to severe ID. Of the two who did not have ID, one had dyscalculia and the other had mild learning difficulties. Craniofacial findings such as hypertelorism, microcephaly, elongated face, ptosis, and mild facial asymmetry were found in some affected individuals. Orofacial clefting was seen in three individuals from our cohort, suggesting that this feature is less common than previously reported. Autism spectrum disorder and attention deficit hyperactivity disorder, which were not previously emphasized in PHF8-XLID, were frequently observed in affected individuals. This series expands the clinical phenotype of this rare ID syndrome caused by loss of PHF8 function.