NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy.

BACKGROUND: NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-ß/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects. METHODS: We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. RESULTS: Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model. CONCLUSION: Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

Closing the gap: Solving complex medically relevant genes at scale.

Comprehending the mechanism behind human diseases with an established heritable component represents the forefront of personalized medicine. Nevertheless, numerous medically important genes are inaccurately represented in short-read sequencing data analysis due to their complexity and repetitiveness or the so-called 'dark regions' of the human genome. The advent of PacBio as a long-read platform has provided new insights, yet HiFi whole-genome sequencing (WGS) cost remains frequently prohibitive. We introduce a targeted sequencing and analysis framework, Twist Alliance Dark Genes Panel (TADGP), designed to offer phased variants across 389 medically important yet complex autosomal genes. We highlight TADGP accuracy across eleven control samples and compare it to WGS. This demonstrates that TADGP achieves variant calling accuracy comparable to HiFi-WGS data, but at a fraction of the cost. Thus, enabling scalability and broad applicability for studying rare diseases or complementing previously sequenced samples to gain insights into these complex genes. TADGP revealed several candidate variants across all cases and provided insight into LPA diversity when tested on samples from rare disease and cardiovascular disease cohorts. In both cohorts, we identified novel variants affecting individual disease-associated genes (e.g., IKZF1, KCNE1). Nevertheless, the annotation of the variants across these 389 medically important genes remains challenging due to their underrepresentation in ClinVar and gnomAD. Consequently, we also offer an annotation resource to enhance the evaluation and prioritization of these variants. Overall, we can demonstrate that TADGP offers a cost-efficient and scalable approach to routinely assess the dark regions of the human genome with clinical relevance.

Somatic mutations of esophageal adenocarcinoma: a comparison between Black and White patients.

Esophageal adenocarcinoma is the most common histological subtype of esophageal cancer in Western countries and shows poor prognosis with rapid growth. EAC is characterized by a strong male predominance and racial disparity. EAC is up to fivefold more common among Whites than Blacks, yet Black patients with EAC have poorer survival rates. The racial disparity remains largely unknown, and there is limited knowledge of mutations in EAC regarding racial disparities. We used whole-exome sequencing to show somatic mutation profiles derived from tumor samples from 18 EAC male patients. We identified three molecular subgroups based on the pre-defined esophageal cancer-specific mutational signatures. Group 1 is associated with age and NTHL1 deficiency-related signatures. Group 2 occurs primarily in Black patients and is associated with signatures related to DNA damage from oxidative stress and NTHL1 deficiency-related signatures. Group 3 is associated with defective homologous recombination-based DNA often caused by BRCA mutation in White patients. We observed significantly mutated race related genes (LCE2B in Black, SDR39U1 in White) were (q-value < 0.1). Our findings underscore the possibility of distinct molecular mutation patterns in EAC among different races. Further studies are needed to validate our findings, which could contribute to precision medicine in EAC.

Germline Genetic Testing and Survival Outcomes Among Children With Rhabdomyosarcoma: A Report From the Children's Oncology Group.

Importance: Determining the impact of germline cancer-predisposition variants (CPVs) on outcomes could inform novel approaches to testing and treating children with rhabdomyosarcoma. Objective: To assess whether CPVs are associated with outcome among children with rhabdomyosarcoma. Design, Setting, and Participants: In this cohort study, data were obtained for individuals, aged 0.01-23.23 years, newly diagnosed with rhabdomyosarcoma who were treated across 171 Children's Oncology Group sites from March 15, 1999, to December 8, 2017. Data analysis was performed from June 16, 2021, to May 15, 2023. Exposure: The presence of a CPV in 24 rhabdomyosarcoma-associated cancer-predisposition genes (CPGs) or an expanded set of 63 autosomal-dominant CPGs. Main Outcomes and Measures: Overall survival (OS) and event-free survival (EFS) were the main outcomes, using the Kaplan-Meier estimator to assess survival probabilities and the Cox proportional hazards regression model to adjust for clinical covariates. Analyses were stratified by tumor histology and the fusion status of PAX3 or PAX7 to the FOXO1 gene. Results: In this study of 580 individuals with rhabdomyosarcoma, the median patient age was 5.9 years (range, 0.01-23.23 years), and the male-to-female ratio was 1.5 to 1 (351 [60.5%] male). For patients with CPVs in rhabdomyosarcoma-associated CPGs, EFS was 48.4% compared with 57.8% for patients without a CPV (P = .10), and OS was 53.7% compared with 65.3% for patients without a CPV (P = .06). After adjustment, patients with CPVs had significantly worse OS (adjusted hazard ratio [AHR], 2.49 [95% CI, 1.39-4.45]; P = .002), and the outcomes were not better among patients with embryonal histology (EFS: AHR, 2.25 [95% CI, 1.25-4.06]; P = .007]; OS: AHR, 2.83 [95% CI, 1.47-5.43]; P = .002]). These associations were not due to the development of a second malignant neoplasm, and importantly, patients with fusion-negative rhabdomyosarcoma who harbored a CPV had similarly inferior outcomes as patients with fusion-positive rhabdomyosarcoma without CPVs (EFS: AHR, 1.35 [95% CI, 0.71-2.59]; P = .37; OS: AHR, 1.71 [95% CI, 0.84-3.47]; P = .14). There were no significant differences in outcome by CPV status of the 63 CPG set. Conclusions and Relevance: This cohort study identified a group of patients with embryonal rhabdomyosarcoma who had a particularly poor outcome. Other important clinical findings included that individuals with TP53 had poor outcomes independent of second malignant neoplasms and that patients with fusion-negative rhabdomyosarcoma who harbored a CPV had outcomes comparable to patients with fusion-positive rhabdomyosarcoma. These findings suggest that germline CPV testing may aid in clinical prognosis and should be considered in prospective risk-based clinical trials.

Genetic sex validation for sample tracking in next-generation sequencing clinical testing.

OBJECTIVE: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. RESULTS: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors (49.09%), samples from transgender participants (3.64%) and stem cell or bone marrow transplant patients (7.27%) along with undetermined sample mix-ups (40%) for which sample swaps occurred prior to arrival at genome centers, however the exact cause of the events at the sampling sites resulting in the mix-ups were not able to be determined.

Inter and intra-host diversity of RSV in hematopoietic stem cell transplant adults with normal and delayed viral clearance.

Respiratory syncytial virus (RSV) infection in immunocompromised individuals often leads to prolonged illness, progression to severe lower respiratory tract infection, and even death. How the host immune environment of the hematopoietic stem cell transplant (HCT) adults can affect viral genetic variation during an acute infection is not understood well. In the present study, we performed whole genome sequencing of RSV/A or RSV/B from samples collected longitudinally from HCT adults with normal (<14 days) and delayed (≥14 days) RSV clearance who were enrolled in a ribavirin trial. We determined the inter-host and intra-host genetic variation of RSV and the effect of mutations on putative glycosylation sites. The inter-host variation of RSV is centered in the attachment (G) and fusion (F) glycoprotein genes followed by polymerase (L) and matrix (M) genes. Interestingly, the overall genetic variation was constant between normal and delayed clearance groups for both RSV/A and RSV/B. Intra-host variation primarily occurred in the G gene followed by non-structural protein (NS1) and L genes; however, gain or loss of stop codons and frameshift mutations appeared only in the G gene and only in the delayed viral clearance group. Potential gain or loss of O-linked glycosylation sites in the G gene occurred both in RSV/A and RSV/B isolates. For RSV F gene, loss of N-linked glycosylation site occurred in three RSV/B isolates within an antigenic epitope. Both oral and aerosolized ribavirin did not cause any mutations in the L gene. In summary, prolonged viral shedding and immune deficiency resulted in RSV variation, especially in structural mutations in the G gene, possibly associated with immune evasion. Therefore, sequencing and monitoring of RSV isolates from immunocompromised patients are crucial as they can create escape mutants that can impact the effectiveness of upcoming vaccines and treatments.

Tryptophan Metabolites And Their Predicted Microbial Sources In Fecal Samples From Healthy Individuals.

Gut microbiota produce tryptophan metabolites (TMs) important to homeostasis. However, measuring TM levels in stool and determining their microbial sources can be difficult. Here, we measured TMs from the indole pathway in fecal samples from 21 healthy adults with the goal to: 1) determine fecal TM concentrations in healthy individuals; 2) link TM levels to bacterial abundance using 16S and whole genome shotgun (WGS) sequencing data; and 3) predict likely bacterial sources of TM production. Within our samples, we identified 151 genera (16S) and 592 bacterial species (WGS). Eight TMs were found in ≥17 fecal samples, including four in all persons. To our knowledge, we are the first to report fecal levels for indole-3-lactate, indole-3-propionate, and 3-indoleacrylate levels in healthy persons. Overall, indole, indole-3-acetate (IAA), and skatole accounted for 86% of the eight TMs measured. Significant correlations were found between seven TMs and 29 bacterial species.  Predicted multiple TM sources support the notion of a complex network of TM production and regulation. Further, the data suggest key roles for Collinsella aerofaciens and IAA, a metabolite reported to maintain intestinal homeostasis through enhanced barrier integrity and anti-inflammatory/antioxidant activities. These findings extend our understanding of TMs and their relationship to the microbial species that act as effectors and/or regulators in the healthy intestine and may lead to novel strategies designed to manipulate tryptophan metabolism to prevent disease and/or restore health to the dysbiotic gut.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data.

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

Dominant negative variants in KIF5B cause osteogenesis imperfecta via down regulation of mTOR signaling.

BACKGROUND: Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown. METHODS: To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy. RESULTS: C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells. CONCLUSION: We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development.

Genetic Sex Validation for Sample Tracking in Clinical Testing.

Objective: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. Results: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors, samples from transgender participants and stem cell or bone marrow transplant patients along with undetermined sample mix-ups.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures.

Background: The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes. Results: Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair. Conclusions: These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Functional Genomics of Gastrointestinal Escherichia coli Isolated from Patients with Cancer and Diarrhea.

We describe the epidemiology and clinical characteristics of 29 patients with cancer and diarrhea in whom Enteroaggregative Escherichia coli (EAEC) was initially identified by GI BioFire panel multiplex. E. coli strains were successfully isolated from fecal cultures in 14 of 29 patients. Six of the 14 strains were identified as EAEC and 8 belonged to other diverse E. coli groups of unknown pathogenesis. We investigated these strains by their adherence to human intestinal organoids, cytotoxic responses, antibiotic resistance profile, full sequencing of their genomes, and annotation of their functional virulome. Interestingly, we discovered novel and enhanced adherence and aggregative patterns for several diarrheagenic pathotypes that were not previously seen when co-cultured with immortalized cell lines. EAEC isolates displayed exceptional adherence and aggregation to human colonoids compared not only to diverse GI E. coli , but also compared to prototype strains of other diarrheagenic E. coli . Some of the diverse E. coli strains that could not be classified as a conventional pathotype also showed an enhanced aggregative and cytotoxic response. Notably, we found a high carriage rate of antibiotic resistance genes in both EAEC strains and diverse GI E. coli isolates and observed a positive correlation between adherence to colonoids and the number of metal acquisition genes carried in both EAEC and the diverse E. coli strains. This work indicates that E. coli from cancer patients constitute strains of remarkable pathotypic and genomic divergence, including strains of unknown disease etiology with unique virulomes. Future studies will allow for the opportunity to re-define E. coli pathotypes with greater diagnostic accuracy and into more clinically relevant groupings.

Longitudinal host transcriptional responses to SARS-CoV-2 infection in adults with extremely high viral load.

Current understanding of viral dynamics of SARS-CoV-2 and host responses driving the pathogenic mechanisms in COVID-19 is rapidly evolving. Here, we conducted a longitudinal study to investigate gene expression patterns during acute SARS-CoV-2 illness. Cases included SARS-CoV-2 infected individuals with extremely high viral loads early in their illness, individuals having low SARS-CoV-2 viral loads early in their infection, and individuals testing negative for SARS-CoV-2. We could identify widespread transcriptional host responses to SARS-CoV-2 infection that were initially most strongly manifested in patients with extremely high initial viral loads, then attenuating within the patient over time as viral loads decreased. Genes correlated with SARS-CoV-2 viral load over time were similarly differentially expressed across independent datasets of SARS-CoV-2 infected lung and upper airway cells, from both in vitro systems and patient samples. We also generated expression data on the human nose organoid model during SARS-CoV-2 infection. The human nose organoid-generated host transcriptional response captured many aspects of responses observed in the above patient samples, while suggesting the existence of distinct host responses to SARS-CoV-2 depending on the cellular context, involving both epithelial and cellular immune responses. Our findings provide a catalog of SARS-CoV-2 host response genes changing over time.

Identification of USP9X as a leukemia susceptibility gene.

We recently reported that children with multiple birth defects have a significantly higher risk of childhood cancer. We performed whole-genome sequencing on a cohort of probands from this study with birth defects and cancer and their parents. Structural variant analysis identified a novel 5 kb de novo heterozygous inframe deletion overlapping the catalytic domain of USP9X in a female proband with multiple birth defects, developmental delay, and B-cell acute lymphoblastic leukemia (B-ALL). Her phenotype was consistent with female-restricted X-linked syndromic intellectual developmental disorder-99 (MRXS99F). Genotype-phenotype analysis including previously reported female probands (n = 42) demonstrated that MRXS99F probands with B-ALL (n = 3) clustered with subjects with loss-of-function (LoF) USP9X variants and multiple anomalies. The cumulative incidence of B-ALL among these female probands (7.1%) was significantly higher than an age- and sex-matched cohort (0.003%) from the Surveillance, Epidemiology, and End Results database (P < .0001, log-rank test). There are no reports of LoF variants in males. Males with hypomorphic missense variants have neurodevelopmental disorders without birth defects or leukemia risk. In contrast, in sporadic B-ALL, somatic LoF USP9X mutations occur in both males and females, and expression levels are comparable in leukemia samples from both sexes (P = .54), with the highest expressors being female patients with extra copies of the X-chromosome. Overall, we describe USP9X as a novel female-specific leukemia predisposition gene associated with multiple congenital, neurodevelopmental anomalies, and B-ALL risk. In contrast, USP9X serves as a tumor suppressor in sporadic pediatric B-ALL in both sexes, with low expression associated with poorer survival in patients with high-risk B-ALL.

Longitudinal host transcriptional responses to SARS-CoV-2 infection in adults with extremely high viral load.

Current understanding of viral dynamics of SARS-CoV-2 and host responses driving the pathogenic mechanisms in COVID-19 is rapidly evolving. Here, we conducted a longitudinal study to investigate gene expression patterns during acute SARS-CoV-2 illness. Cases included SARS-CoV-2 infected individuals with extremely high viral loads early in their illness, individuals having low SARS-CoV-2 viral loads early in their infection, and individuals testing negative for SARS-CoV-2. We could identify widespread transcriptional host responses to SARS-CoV-2 infection that were initially most strongly manifested in patients with extremely high initial viral loads, then attenuating within the patient over time as viral loads decreased. Genes correlated with SARS-CoV-2 viral load over time were similarly differentially expressed across independent datasets of SARS-CoV-2 infected lung and upper airway cells, from both in vitro systems and patient samples. We also generated expression data on the human nose organoid model during SARS-CoV-2 infection. The human nose organoid-generated host transcriptional response captured many aspects of responses observed in the above patient samples, while suggesting the existence of distinct host responses to SARS-CoV-2 depending on the cellular context, involving both epithelial and cellular immune responses. Our findings provide a catalog of SARS-CoV-2 host response genes changing over time.

Circulating tumor DNA sequencing of pediatric solid and brain tumor patients: An institutional feasibility study.

The potential of circulating tumor DNA (ctDNA) analysis to serve as a real-time "liquid biopsy" for children with central nervous system (CNS) and non-CNS solid tumors remains to be fully elucidated. We conducted a study to investigate the feasibility and potential clinical utility of ctDNA sequencing in pediatric patients enrolled on an institutional clinical genomics trial. A total of 240 patients had tumor DNA profiling performed during the study period. Plasma samples were collected at study enrollment from 217 patients and then longitudinally from a subset of patients. Successful cell-free DNA extraction and quantification occurred in 216 of 217 (99.5%) of these initial samples. Twenty-four patients were identified whose tumors harbored 30 unique variants that were potentially detectable on a commercially-available ctDNA panel. Twenty of these 30 mutations (67%) were successfully detected by next-generation sequencing in the ctDNA from at least one plasma sample. The rate of ctDNA mutation detection was higher in patients with non-CNS solid tumors (7/9, 78%) compared to those with CNS tumors (9/15, 60%). A higher ctDNA mutation detection rate was also observed in patients with metastatic disease (9/10, 90%) compared to non-metastatic disease (7/14, 50%), although tumor-specific variants were detected in a few patients in the absence of radiographic evidence of disease. This study illustrates the feasibility of incorporating longitudinal ctDNA analysis into the management of relapsed or refractory patients with childhood CNS or non-CNS solid tumors.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

The genomic landscape of familial glioma.

Glioma is a rare brain tumor with a poor prognosis. Familial glioma is a subset of glioma with a strong genetic predisposition that accounts for approximately 5% of glioma cases. We performed whole-genome sequencing on an exploratory cohort of 203 individuals from 189 families with a history of familial glioma and an additional validation cohort of 122 individuals from 115 families. We found significant enrichment of rare deleterious variants of seven genes in both cohorts, and the most significantly enriched gene was HERC2 (P = 0.0006). Furthermore, we identified rare noncoding variants in both cohorts that were predicted to affect transcription factor binding sites or cause cryptic splicing. Last, we selected a subset of discovered genes for validation by CRISPR knockdown screening and found that DMBT1, HP1BP3, and ZCH7B3 have profound impacts on proliferation. This study performs comprehensive surveillance of the genomic landscape of familial glioma.

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.

A biallelic frameshift indel in PPP1R35 as a cause of primary microcephaly.

Protein phosphatase 1 regulatory subunit 35 (PPP1R35) encodes a centrosomal protein required for recruiting microtubule-binding elongation machinery. Several proteins in this centriole biogenesis pathway correspond to established primary microcephaly (MCPH) genes, and multiple model organism studies hypothesize PPP1R35 as a candidate MCPH gene. Here, using exome sequencing (ES) and family-based rare variant analyses, we report a homozygous, frameshifting indel deleting the canonical stop codon in the last exon of PPP1R35 [Chr7: c.753_*3delGGAAGCGTAGACCinsCG (p.Trp251Cysfs*22)]; the variant allele maps in a 3.7 Mb block of absence of heterozygosity (AOH) in a proband with severe MCPH (-4.3 SD at birth, -6.1 SD by 42 months), pachygyria, and global developmental delay from a consanguineous Turkish kindred. Droplet digital PCR (ddPCR) confirmed mutant mRNA expression in fibroblasts. In silico prediction of the translation of mutant PPP1R35 is expected to be elongated by 18 amino acids before encountering a downstream stop codon. This complex indel allele is absent in public databases (ClinVar, gnomAD, ARIC, 1000 genomes) and our in-house database of 14,000+ exomes including 1800+ Turkish exomes supporting predicted pathogenicity. Comprehensive literature searches for PPP1R35 variants yielded two probands affected with severe microcephaly (-15 SD and -12 SD) with the same homozygous indel from a single, consanguineous, Iranian family from a cohort of 404 predominantly Iranian families. The lack of heterozygous cases in two large cohorts representative of the genetic background of these two families decreased our suspicion of a founder allele and supports the contention of a recurrent mutation. We propose two potential secondary structure mutagenesis models for the origin of this variant allele mediated by hairpin formation between complementary GC rich segments flanking the stop codon via secondary structure mutagenesis.