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.

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.

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.

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.

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.

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.

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.

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci.

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.