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.

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.

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.

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.

New observations on maternal age effect on germline de novo mutations.

Germline mutations are the source of evolution and contribute substantially to many health-related processes. Here we use whole-genome deep sequencing data from 693 parents-offspring trios to examine the de novo point mutations (DNMs) in the offspring. Our estimate for the mutation rate per base pair per generation is 1.05 × 10(-8), well within the range of previous studies. We show that maternal age has a small but significant correlation with the total number of DNMs in the offspring after controlling for paternal age (0.51 additional mutations per year, 95% CI: 0.29, 0.73), which was not detectable in the smaller and younger parental cohorts of earlier studies. Furthermore, while the total number of DNMs increases at a constant rate for paternal age, the contribution from the mother increases at an accelerated rate with age.These observations have implications related to the incidence of de novo mutations relating to maternal age.

Utility of whole-genome sequencing for detection of newborn screening disorders in a population cohort of 1,696 neonates.

PURPOSE: To assess the potential of whole-genome sequencing (WGS) to replicate and augment results from conventional blood-based newborn screening (NBS). METHODS: Research-generated WGS data from an ancestrally diverse cohort of 1,696 infants and both parents of each infant were analyzed for variants in 163 genes involved in disorders included or under discussion for inclusion in US NBS programs. WGS results were compared with results from state NBS and related follow-up testing. RESULTS: NBS genes are generally well covered by WGS. There is a median of one (range: 0-6) database-annotated pathogenic variant in the NBS genes per infant. Results of WGS and NBS in detecting 28 state-screened disorders and four hemoglobin traits were concordant for 88.6% of true positives (n = 35) and 98.9% of true negatives (n = 45,757). Of the five infants affected with a state-screened disorder, WGS identified two whereas NBS detected four. WGS yielded fewer false positives than NBS (0.037 vs. 0.17%) but more results of uncertain significance (0.90 vs. 0.013%). CONCLUSION: WGS may help rule in and rule out NBS disorders, pinpoint molecular diagnoses, and detect conditions not amenable to current NBS assays.

Systematic evaluation of underlying defects in DNA repair as an approach to case-only assessment of familial prostate cancer.

Risk assessment for prostate cancer is challenging due to its genetic heterogeneity. In this study, our goal was to develop an operational framework to select and evaluate gene variants that may contribute to familial prostate cancer risk. Drawing on orthogonal sources, we developed a candidate list of genes relevant to prostate cancer, then analyzed germline exomes from 12 case-only prostate cancer patients from high-risk families to identify patterns of protein-damaging gene variants. We described an average of 5 potentially disruptive variants in each individual and annotated them in the context of public databases representing human variation. Novel damaging variants were found in several genes of relevance to prostate cancer. Almost all patients had variants associated with defects in DNA damage response. Many also had variants linked to androgen signaling. Treatment of primary T-lymphocytes from these prostate cancer patients versus controls with DNA damaging agents showed elevated levels of the DNA double strand break (DSB) marker γH2AX (p < 0.05), supporting the idea of an underlying defect in DNA repair. This work suggests the value of focusing on underlying defects in DNA damage in familial prostate cancer risk assessment and demonstrates an operational framework for exome sequencing in case-only prostate cancer genetic evaluation.

Genetic Variants That Predispose to DNA Double-Strand Breaks in Lymphocytes From a Subset of Patients With Familial Colorectal Carcinomas.

BACKGROUND & AIMS: DNA structural lesions are prevalent in sporadic colorectal cancer. Therefore, we proposed that gene variants that predispose to DNA double-strand breaks (DSBs) would be found in patients with familial colorectal carcinomas of an undefined genetic basis (UFCRC). METHODS: We collected primary T cells from 25 patients with UFCRC and matched patients without colorectal cancer (controls) and assayed for DSBs. We performed exome sequence analyses of germline DNA from 20 patients with UFCRC and 5 undiagnosed patients with polyposis. The prevalence of identified variants in genes linked to DNA integrity was compared with that of individuals without a family history of cancer. The effects of representative variants found to be associated with UFCRC was confirmed in functional assays with HCT116 cells. RESULTS: Primary T cells from most patients with UFCRC had increased levels of the DSB marker γ(phosphorylated)histone2AX (γH2AX) after treatment with DNA damaging agents, compared with T cells from controls (P < .001). Exome sequence analysis identified a mean 1.4 rare variants per patient that were predicted to disrupt functions of genes relevant to DSBs. Controls (from public databases) had a much lower frequency of variants in the same genes (P < .001). Knockdown of representative variant genes in HCT116 CRC cells increased γH2AX. A detailed analysis of immortalized patient-derived B cells that contained variants in the Werner syndrome, RecQ helicase-like gene (WRN, encoding T705I), and excision repair cross-complementation group 6 (ERCC6, encoding N180Y) showed reduced levels of these proteins and increased DSBs, compared with B cells from controls. This phenotype was rescued by exogenous expression of WRN or ERCC6. Direct analysis of the recombinant variant proteins confirmed defective enzymatic activities. CONCLUSIONS: These results provide evidence that defects in suppression of DSBs underlie some cases of UFCRC; these can be identified by assays of circulating lymphocytes. We specifically associated UFCRC with variants in WRN and ERCC6 that reduce the capacity for repair of DNA DSBs. These observations could lead to a simple screening strategy for UFCRC, and provide insight into the pathogenic mechanisms of colorectal carcinogenesis.