Elucidating the Role of Wildtype and Variant FGFR2 Structural Dynamics in (Dys)Function and Disorder.

The fibroblast growth factor receptor 2 (FGFR2) gene is one of the most extensively studied genes with many known mutations implicated in several human disorders, including oncogenic ones. Most FGFR2 disease-associated gene mutations are missense mutations that result in constitutive activation of the FGFR2 protein and downstream molecular pathways. Many tertiary structures of the FGFR2 kinase domain are publicly available in the wildtype and mutated forms and in the inactive and activated state of the receptor. The current literature suggests a molecular brake inhibiting the ATP-binding A loop from adopting the activated state. Mutations relieve this brake, triggering allosteric changes between active and inactive states. However, the existing analysis relies on static structures and fails to account for the intrinsic structural dynamics. In this study, we utilize experimentally resolved structures of the FGFR2 tyrosine kinase domain and machine learning to capture the intrinsic structural dynamics, correlate it with functional regions and disease types, and enrich it with predicted structures of variants with currently no experimentally resolved structures. Our findings demonstrate the value of machine learning-enabled characterizations of structure dynamics in revealing the impact of mutations on (dys)function and disorder in FGFR2.

Estimated Prevalence and Clinical Manifestations of UBA1 Variants Associated With VEXAS Syndrome in a Clinical Population.

Importance: VEXAS (vacuoles, E1-ubiquitin-activating enzyme, X-linked, autoinflammatory, somatic) syndrome is a disease with rheumatologic and hematologic features caused by somatic variants in UBA1. Pathogenic variants are associated with a broad spectrum of clinical manifestations. Knowledge of prevalence, penetrance, and clinical characteristics of this disease have been limited by ascertainment biases based on known phenotypes. Objective: To determine the prevalence of pathogenic variants in UBA1 and associated clinical manifestations in an unselected population using a genomic ascertainment approach. Design, Setting, and Participants: This retrospective observational study evaluated UBA1 variants in exome data from 163 096 participants within the Geisinger MyCode Community Health Initiative. Clinical phenotypes were determined from Geisinger electronic health record data from January 1, 1996, to January 1, 2022. Exposures: Exome sequencing was performed. Main Outcomes and Measures: Outcome measures included prevalence of somatic UBA1 variation; presence of rheumatologic, hematologic, pulmonary, dermatologic, and other findings in individuals with somatic UBA1 variation on review of the electronic health record; review of laboratory data; bone marrow biopsy pathology analysis; and in vitro enzymatic assays. Results: In 163 096 participants (mean age, 52.8 years; 94% White; 61% women), 11 individuals harbored likely somatic variants at known pathogenic UBA1 positions, with 11 of 11 (100%) having clinical manifestations consistent with VEXAS syndrome (9 male, 2 female). A total of 5 of 11 individuals (45%) did not meet criteria for rheumatologic and/or hematologic diagnoses previously associated with VEXAS syndrome; however, all individuals had anemia (hemoglobin: mean, 7.8 g/dL; median, 7.5 g/dL), which was mostly macrocytic (10/11 [91%]) with concomitant thrombocytopenia (10/11 [91%]). Among the 11 patients identified, there was a pathogenic variant in 1 male participant prior to onset of VEXAS-related signs or symptoms and 2 female participants had disease with heterozygous variants. A previously unreported UBA1 variant (c.1861A>T; p.Ser621Cys) was found in a symptomatic patient, with in vitro data supporting a catalytic defect and pathogenicity. Together, disease-causing UBA1 variants were found in 1 in 13 591 unrelated individuals (95% CI, 1:7775-1:23 758), 1 in 4269 men older than 50 years (95% CI, 1:2319-1:7859), and 1 in 26 238 women older than 50 years (95% CI, 1:7196-1:147 669). Conclusions and Relevance: This study provides an estimate of the prevalence and a description of the clinical manifestations of UBA1 variants associated with VEXAS syndrome within a single regional health system in the US. Additional studies are needed in unselected and genetically diverse populations to better define general population prevalence and phenotypic spectrum.

Leveraging population-based exome screening to impact clinical care: The evolution of variant assessment in the Geisinger MyCode research project.

Exome and genome sequencing are increasingly utilized in research studies and clinical care and can provide clinically relevant information beyond the initial intent for sequencing, including medically actionable secondary findings. Despite ongoing debate about sharing this information with patients and participants, a growing number of clinical laboratories and research programs routinely report secondary findings that increase the risk for selected diseases. Recently, there has been a push to maximize the potential benefit of this practice by implementing proactive genomic screening at the population level irrespective of medical history, but the feasibility of deploying population-scale proactive genomic screening requires scaling key elements of the genomic data evaluation process. Herein, we describe the motivation, development, and implementation of a population-scale variant-first screening pipeline combining bioinformatics-based filtering with a manual review process to screen for clinically relevant findings in research exomes generated through the DiscovEHR collaboration within Geisinger's MyCode® research project. Consistent with other studies, this pipeline yields a screen-positive detection rate between 2.1 and 2.6% (depending on inclusion of those with prior indication-based testing) in 130,048 adult MyCode patient-participants screened for clinically relevant findings in 60 genes. Our variant-first pipeline affords cost and time savings by filtering out negative cases, thereby avoiding analysis of each exome one-by-one, as typically employed in the diagnostic setting. While research is still needed to fully appreciate the benefits of population genomic screening, MyCode provides the first demonstration of a program at scale to help shape how population genomic screening is integrated into routine clinical care.

Genotype-first analysis of a generally healthy population cohort supports genetic testing for diagnosis of hereditary angioedema of unknown cause.

BACKGROUND: Hereditary angioedema (HAE) is a potentially life-threatening group of conditions that is often underdiagnosed or misdiagnosed. As HAE is typically diagnosed by detecting C1 inhibitor deficiency, there is a critical need for methods that can identify affected individuals with normal C1 inhibitor. The recent discovery of associations between PLG K330E and ANGPT1 A119S and HAE of unknown genetic cause (HAE-U), has raised the possibility that genetic evaluation could be used to diagnose HAE-U in patients with unexplained angioedema or non-confirmatory laboratory testing. CASE PRESENTATION: We analyzed genome sequences from a generally healthy population cohort of 2820 adults and identified PLG K330E in one individual. Subsequent review of this participant's medical history revealed symptoms clinically attributed to allergy of unknown etiology but that are consistent with published descriptions of HAE patients carrying the PLG K330E variant. The participant, a 31 year old female, reported lip and tongue angioedema, without wheals, which did not respond to treatment with steroids or antihistamines. CONCLUSIONS: The genotype-first approach demonstrated that detection of PLG K330E in undiagnosed or misdiagnosed individuals can identify patients actually affected with HAE-U. The genetic diagnosis will facilitate selection of appropriate treatment, discontinuation of therapies ineffective for this condition, and timely diagnosis of affected family members. The results support a role of PLG K330E in the pathogenesis of HAE and suggest that genetic testing be considered as an approach to diagnose patients with unexplained angioedema.

Mendelian Inconsistent Signatures from 1314 Ancestrally Diverse Family Trios Distinguish Biological Variation from Sequencing Error.

Next-generation sequencing enables advances in the clinical application of genomics by providing high-throughput detection of genomic variation. However, next-generation sequencing technologies, especially whole-genome sequencing (WGS), are often associated with a high false-positive rate. Trio-based WGS can contribute significantly towards improved quality control methods. Mendelian-inconsistent calls (MIC) in parent-child trios are commonly attributed to erroneous sequencing calls, as the true de novo mutation rate is extremely low compared with MIC incidence. Here, we analyzed WGS data from 1314 mother, father, and child trios across ethnically diverse populations with the goal of characterizing MIC. Genotype calls in a trio can be used to assign different signatures to MIC. MIC occur more frequently within repeats but show varying distribution and error mechanisms across repeat types. MIC are enriched within poly-A/T runs in short interspersed nuclear elements. Alignability scores, allele balance, and relative parental read depth vary among MIC signatures and these differences should be considered when designing filters for MIC reduction. MIC cluster in germline deletions and these MIC also segregate with population. Our results provide a basis for making decisions on how each MIC type should be evaluated before discarding them as errors or including them in alternative applications. With the reduction of sequencing cost, family trio whole genome and exome analysis are being performed more routinely in clinical practice. We provide a reference that can be used for annotating MIC with their frequencies in a larger population to aid in the filtering of candidate de novo mutations.

Genomic and molecular characterization of preterm birth.

Preterm birth (PTB) complications are the leading cause of long-term morbidity and mortality in children. By using whole blood samples, we integrated whole-genome sequencing (WGS), RNA sequencing (RNA-seq), and DNA methylation data for 270 PTB and 521 control families. We analyzed this combined dataset to identify genomic variants associated with PTB and secondary analyses to identify variants associated with very early PTB (VEPTB) as well as other subcategories of disease that may contribute to PTB. We identified differentially expressed genes (DEGs) and methylated genomic loci and performed expression and methylation quantitative trait loci analyses to link genomic variants to these expression and methylation changes. We performed enrichment tests to identify overlaps between new and known PTB candidate gene systems. We identified 160 significant genomic variants associated with PTB-related phenotypes. The most significant variants, DEGs, and differentially methylated loci were associated with VEPTB. Integration of all data types identified a set of 72 candidate biomarker genes for VEPTB, encompassing genes and those previously associated with PTB. Notably, PTB-associated genes RAB31 and RBPJ were identified by all three data types (WGS, RNA-seq, and methylation). Pathways associated with VEPTB include EGFR and prolactin signaling pathways, inflammation- and immunity-related pathways, chemokine signaling, IFN-γ signaling, and Notch1 signaling. Progress in identifying molecular components of a complex disease is aided by integrated analyses of multiple molecular data types and clinical data. With these data, and by stratifying PTB by subphenotype, we have identified associations between VEPTB and the underlying biology.

Pitfalls of clinical exome and gene panel testing: alternative transcripts.

PURPOSE: Clinical exome and gene panel testing can provide molecular diagnoses for patients with rare Mendelian disorders, but for many patients these tests are nonexplanatory. We investigated whether interrogation of alternative transcripts in known disease genes could provide answers for additional patients. METHODS: We integrated alternative transcripts for known neonatal epilepsy genes with RNA-Seq data to identify brain-expressed coding regions that are not evaluated by popular neonatal epilepsy clinical gene panel and exome tests. RESULTS: We found brain-expressed alternative coding regions in 89 (30%) of 292 neonatal epilepsy genes. The 147 regions encompass 15,713 bases that are noncoding in the primary transcripts analyzed by the clinical tests. Alternative coding regions from at least 5 genes carry reported pathogenic variants. Three candidate variants in these regions were identified in public exome data from 337 epilepsy patients. Incorporating alternative transcripts into the analysis of neonatal epilepsy genes in 44 patient genomes identified the pathogenic variant for the epilepsy case and 2 variants of uncertain significance (VUS) among the 43 control cases. CONCLUSION: Assessment of alternative transcripts in exon-based clinical genetic tests, including gene panel, exome, and genome sequencing, may provide diagnoses for patients for whom standard testing is unrevealing, without introducing many VUS.

Author Correction: Parent-of-origin-specific signatures of de novo mutations.

In the version of this article published, the P values for the enrichment of single mutation categories were inadvertently not corrected for multiple testing. After multiple-testing correction, only two of the six mutation categories mentioned are still statistically significant. To reflect this, the text "More specifically, paternally derived DNMs are enriched in transitions in A[.]G contexts, especially ACG>ATG and ATG>ACG (Bonferroni-corrected P = 1.3 × 10-2 and P = 1 × 10-3, respectively). Additionally, we observed overrepresentation of ATA>ACA mutations (Bonferroni-corrected P = 4.28 × 10-2) for DNMs of paternal origin. Among maternally derived DNMs, CCA>CTA, GCA>GTA and TCT>TGT mutations were significantly overrepresented (Bonferroni-corrected P = 4 × 10-4, P = 5 × 10-4, P = 1 × 10-3, respectively)" should read "More specifically, CCA>CTA and GCA>GTA mutations were significantly overenriched on the maternal allele (Bonferroni-corrected P = 0.0192 and P = 0.048, respectively)." Additionally, the last sentence to the legend for Fig. 3b should read "Green boxes highlight the mutation categories that differ significantly" instead of "Green boxes highlight the mutation categories that differ more than 1% of mutation load with a bootstrapping P value <0.05." Corrected versions of Fig. 3b and Supplementary Table 25 appear with the Author Correction.

Gene-Centric Analysis of Preeclampsia Identifies Maternal Association at PLEKHG1.

The genetic susceptibility to preeclampsia, a pregnancy-specific complication with significant maternal and fetal morbidity, has been poorly characterized. To identify maternal genes associated with preeclampsia risk, we assembled 498 cases and 1864 controls of European ancestry from preeclampsia case-control collections in 5 different US sites (with additional matched population controls), genotyped samples on a cardiovascular gene-centric array composed of variants from ≈2000 genes selected based on prior genetic studies of cardiovascular and metabolic diseases and performed case-control genetic association analysis on 27 429 variants passing quality control. In silico replication testing of 9 lead signals with P<10-4 was performed in independent European samples from the SOPHIA (Study of Pregnancy Hypertension in Iowa) and Inova cohorts (212 cases, 456 controls). Multiethnic assessment of lead signals was then performed in samples of black (26 cases, 136 controls), Hispanic (132 cases, 468 controls), and East Asian (9 cases, 80 controls) ancestry. Multiethnic meta-analysis (877 cases, 3004 controls) revealed a study-wide statistically significant association of the rs9478812 variant in the pleiotropic PLEKHG1 gene (odds ratio, 1.40 [1.23-1.60]; Pmeta=5.90×10-7). The rs9478812 effect was even stronger in the subset of European cases with known early-onset preeclampsia (236 cases diagnosed <37 weeks, 1864 controls; odds ratio, 1.59 [1.27-1.98]; P=4.01×10-5). PLEKHG1 variants have previously been implicated in genome-wide association studies of blood pressure, body weight, and neurological disorders. Although larger studies are required to further define maternal preeclampsia heritability, this study identifies a novel maternal risk locus for further investigation.

Germline de novo mutation clusters arise during oocyte aging in genomic regions with high double-strand-break incidence.

Clustering of mutations has been observed in cancer genomes as well as for germline de novo mutations (DNMs). We identified 1,796 clustered DNMs (cDNMs) within whole-genome-sequencing data from 1,291 parent-offspring trios to investigate their patterns and infer a mutational mechanism. We found that the number of clusters on the maternal allele was positively correlated with maternal age and that these clusters consisted of more individual mutations with larger intermutational distances than those of paternal clusters. More than 50% of maternal clusters were located on chromosomes 8, 9 and 16, in previously identified regions with accelerated maternal mutation rates. Maternal clusters in these regions showed a distinct mutation signature characterized by C>G transversions. Finally, we found that maternal clusters were associated with processes involving double-strand-breaks (DSBs), such as meiotic gene conversions and de novo deletion events. This result suggested accumulation of DSB-induced mutations throughout oocyte aging as the mechanism underlying the formation of maternal mutation clusters.

Mutation in an alternative transcript of CDKL5 in a boy with early-onset seizures.

Infantile-onset epilepsies are a set of severe, heterogeneous disorders for which clinical genetic testing yields causative mutations in ∼20%-50% of affected individuals. We report the case of a boy presenting with intractable seizures at 2 wk of age, for whom gene panel testing was unrevealing. Research-based whole-genome sequencing of the proband and four unaffected family members identified a de novo mutation, NM_001323289.1:c.2828_2829delGA in CDKL5, a gene associated with X-linked early infantile epileptic encephalopathy 2. CDKL5 has multiple alternative transcripts, and the mutation lies in an exon in the brain-expressed forms. The mutation was undetected by gene panel sequencing because of its intronic location in the CDKL5 transcript typically used to define the exons of this gene for clinical exon-based tests (NM_003159). This is the first report of a patient with a mutation in an alternative transcript of CDKL5 This finding suggests that incorporating alternative transcripts into the design and variant interpretation of exon-based tests, including gene panel and exome sequencing, could improve the diagnostic yield.

A genome-wide association study identifies only two ancestry specific variants associated with spontaneous preterm birth.

Preterm birth (PTB), or the delivery prior to 37 weeks of gestation, is a significant cause of infant morbidity and mortality. Although twin studies estimate that maternal genetic contributions account for approximately 30% of the incidence of PTB, and other studies reported fetal gene polymorphism association, to date no consistent associations have been identified. In this study, we performed the largest reported genome-wide association study analysis on 1,349 cases of PTB and 12,595 ancestry-matched controls from the focusing on genomic fetal signals. We tested over 2 million single nucleotide polymorphisms (SNPs) for associations with PTB across five subpopulations: African (AFR), the Americas (AMR), European, South Asian, and East Asian. We identified only two intergenic loci associated with PTB at a genome-wide level of significance: rs17591250 (P = 4.55E-09) on chromosome 1 in the AFR population and rs1979081 (P = 3.72E-08) on chromosome 8 in the AMR group. We have queried several existing replication cohorts and found no support of these associations. We conclude that the fetal genetic contribution to PTB is unlikely due to single common genetic variant, but could be explained by interactions of multiple common variants, or of rare variants affected by environmental influences, all not detectable using a GWAS alone.

Experience with genomic sequencing in pediatric patients with congenital cardiac defects in a large community hospital.

BACKGROUND: Congenital cardiac defects, whether isolated or as part of a larger syndrome, are the most common type of human birth defect occurring on average in about 1% of live births depending on the malformation. As there is an expanding understanding of the underlying molecular mechanisms by which a cardiac defect may occur, there is a need to assess the current rates of diagnosis of cardiac defects by molecular sequencing in a clinical setting. METHODS AND RESULTS: In this report, we evaluated 34 neonatal and pediatric patients born with a cardiac defect and their parents using exomized preexisting whole genome sequencing (WGS) data to model clinically available exon-based tests. Overall, we identified candidate variants in previously reported cardiac-related genes in 35% (12/34) of the probands. These include clearly pathogenic variants in two of 34 patients (6%) and variants of uncertain significance in relevant genes in 10 patients (26%), of these latter 10, 2 segregated with clinically apparent findings in the family trios. CONCLUSIONS: These findings suggest that with current knowledge of the proteins underlying CHD, genomic sequencing can identify the underlying genetic etiology in certain patients; however, this technology currently does not have a high enough yield to be of routine clinical use in the screening of pediatric congenital cardiac defects.

Tools for annotation and comparison of structural variation.

The impact of structural variants (SVs) on a variety of organisms and diseases like cancer has become increasingly evident. Methods for SV detection when studying genomic differences across cells, individuals or populations are being actively developed. Currently, just a few methods are available to compare different SVs callsets, and no specialized methods are available to annotate SVs that account for the unique characteristics of these variant types. Here, we introduce SURVIVOR_ant, a tool that compares types and breakpoints for candidate SVs from different callsets and enables fast comparison of SVs to genomic features such as genes and repetitive regions, as well as to previously established SV datasets such as from the 1000 Genomes Project. As proof of concept we compared 16 SV callsets generated by different SV calling methods on a single genome, the Genome in a Bottle sample HG002 (Ashkenazi son), and annotated the SVs with gene annotations, 1000 Genomes Project SV calls, and four different types of repetitive regions. Computation time to annotate 134,528 SVs with 33,954 of annotations was 22 seconds on a laptop.

Genomic analysis of an infant with intractable diarrhea and dilated cardiomyopathy.

We describe a case of an infant presenting with intractable diarrhea who subsequently developed dilated cardiomyopathy, for whom a diagnosis was not initially achieved despite extensive clinical testing, including panel-based genetic testing. Research-based whole-genome sequences of the proband and both parents were analyzed by the SAVANNA pipeline, a variant prioritization strategy integrating features of variants, genes, and phenotypes, which was implemented using publicly available tools. Although the intestinal morphological abnormalities characteristic of congenital tufting enteropathy (CTE) were not observed in the initial clinical gastrointestinal tract biopsies of the proband, an intronic variant, EPCAM c.556-14A>G, previously identified as pathogenic for CTE, was found in the homozygous state. A newborn cousin of the proband also presenting with intractable diarrhea was found to carry the same homozygous EPCAM variant, and clinical testing revealed intestinal tufting and loss of EPCAM staining. This variant, however, was considered nonexplanatory for the proband's dilated cardiomyopathy, which could be a sequela of the child's condition and/or related to other genetic variants, which include de novo mutations in the genes NEDD4L and GSK3A and a maternally inherited SCN5A variant. This study illustrates three ways in which genomic sequencing can aid in the diagnosis of clinically challenging patients: differential diagnosis despite atypical clinical presentation, distinguishing the possibilities of a syndromic condition versus multiple conditions, and generating hypotheses for novel contributory genes.

Utilization of genomic sequencing for population screening of immunodeficiencies in the newborn.

PurposeImmunodeficiency screening has been added to many state-directed newborn screening programs. The current methodology is limited to screening for severe T-cell lymphopenia disorders. We evaluated the potential of genomic sequencing to augment current newborn screening for immunodeficiency, including identification of non-T cell disorders.MethodsWe analyzed whole-genome sequencing (WGS) and clinical data from a cohort of 1,349 newborn-parent trios by genotype-first and phenotype-first approaches. For the genotype-first approach, we analyzed predicted protein-impacting variants in 329 immunodeficiency-related genes in the WGS data. As a phenotype-first approach, electronic health records were used to identify children with clinical features suggestive of immunodeficiency. Genomes of these children and their parents were analyzed using a separate pipeline for identification of candidate pathogenic variants for rare Mendelian disorders.ResultsWGS provides adequate coverage for most known immunodeficiency-related genes. 13,476 distinct variants and 8,502 distinct predicted protein-impacting variants were identified in this cohort; five individuals carried potentially pathogenic variants requiring expert clinical correlation. One clinically asymptomatic individual was found genomically to have complement component 9 deficiency. Of the symptomatic children, one was molecularly identified as having an immunodeficiency condition and two were found to have other molecular diagnoses.ConclusionNeonatal genomic sequencing can potentially augment newborn screening for immunodeficiency.

Parent-of-origin-specific signatures of de novo mutations.

De novo mutations (DNMs) originating in gametogenesis are an important source of genetic variation. We use a data set of 7,216 autosomal DNMs with resolved parent of origin from whole-genome sequencing of 816 parent-offspring trios to investigate differences between maternally and paternally derived DNMs and study the underlying mutational mechanisms. Our results show that the number of DNMs in offspring increases not only with paternal age, but also with maternal age, and that some genome regions show enrichment for maternally derived DNMs. We identify parent-of-origin-specific mutation signatures that become more pronounced with increased parental age, pointing to different mutational mechanisms in spermatogenesis and oogenesis. Moreover, we find DNMs that are spatially clustered to have a unique mutational signature with no significant differences between parental alleles, suggesting a different mutational mechanism. Our findings provide insights into the molecular mechanisms that underlie mutagenesis and are relevant to disease and evolution in humans.

Whole-Genome Sequencing of a Healthy Aging Cohort.

Studies of long-lived individuals have revealed few genetic mechanisms for protection against age-associated disease. Therefore, we pursued genome sequencing of a related phenotype-healthy aging-to understand the genetics of disease-free aging without medical intervention. In contrast with studies of exceptional longevity, usually focused on centenarians, healthy aging is not associated with known longevity variants, but is associated with reduced genetic susceptibility to Alzheimer and coronary artery disease. Additionally, healthy aging is not associated with a decreased rate of rare pathogenic variants, potentially indicating the presence of disease-resistance factors. In keeping with this possibility, we identify suggestive common and rare variant genetic associations implying that protection against cognitive decline is a genetic component of healthy aging. These findings, based on a relatively small cohort, require independent replication. Overall, our results suggest healthy aging is an overlapping but distinct phenotype from exceptional longevity that may be enriched with disease-protective genetic factors. VIDEO ABSTRACT.