Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

High Clinical Exome Sequencing Diagnostic Rates and Novel Phenotypic Expansions for Nonisolated Microphthalmia, Anophthalmia, and Coloboma

Purpose

A molecular diagnosis is only made in a subset of individuals with nonisolated microphthalmia, anophthalmia, and coloboma (MAC). This may be due to underutilization of clinical (whole) exome sequencing (cES) and an incomplete understanding of the genes that cause MAC. The purpose of this study is to determine the efficacy of cES in cases of nonisolated MAC and to identify new MAC phenotypic expansions.

Methods

We determined the efficacy of cES in 189 individuals with nonisolated MAC. We then used cES data, a validated machine learning algorithm, and previously published expression data, case reports, and animal models to determine which candidate genes were most likely to contribute to the development of MAC.

Results

We found the efficacy of cES in nonisolated MAC to be between 32.3% (61/189) and 48.1% (91/189). Most genes affected in our cohort were not among genes currently screened in clinically available ophthalmologic gene panels. A subset of the genes implicated in our cohort had not been clearly associated with MAC. Our analyses revealed sufficient evidence to support low-penetrance MAC phenotypic expansions involving nine of these human disease genes.

Conclusions

We conclude that cES is an effective means of identifying a molecular diagnosis in individuals with nonisolated MAC and may identify putatively damaging variants that would be missed if only a clinically available ophthalmologic gene panel was obtained. Our data also suggest that deleterious variants in BRCA2, BRIP1, KAT6A, KAT6B, NSF, RAC1, SMARCA4, SMC1A, and TUBA1A can contribute to the development of MAC.

Clinical exome sequencing efficacy and phenotypic expansions involving anomalous pulmonary venous return

Anomalous pulmonary venous return (APVR) frequently occurs with other congenital heart defects (CHDs) or extra-cardiac anomalies. While some genetic causes have been identified, the optimal approach to genetic testing in individuals with APVR remains uncertain, and the etiology of most cases of APVR is unclear. Here, we analyzed molecular data from 49 individuals to determine the diagnostic yield of clinical exome sequencing (ES) for non-isolated APVR. A definitive or probable diagnosis was made for 8 of those individuals yielding a diagnostic efficacy rate of 16.3%. We then analyzed molecular data from 62 individuals with APVR accrued from three databases to identify novel APVR genes. Based on data from this analysis, published case reports, mouse models, and/or similarity to known APVR genes as revealed by a machine learning algorithm, we identified 3 genes-EFTUD2, NAA15, and NKX2-1-for which there is sufficient evidence to support phenotypic expansion to include APVR. We also provide evidence that 3 recurrent copy number variants contribute to the development of APVR: proximal 1q21.1 microdeletions involving RBM8A and PDZK1, recurrent BP1-BP2 15q11.2 deletions, and central 22q11.2 deletions involving CRKL. Our results suggest that ES and chromosomal microarray analysis (or genome sequencing) should be considered for individuals with non-isolated APVR for whom a genetic etiology has not been identified, and that genetic testing to identify an independent genetic etiology of APVR is not warranted in individuals with EFTUD2-, NAA15-, and NKX2-1-related disorders.

Modulating a prebiotic food source influences inflammation and immune-regulating gut microbes and metabolites: insights from the BE GONE trial

BACKGROUND

Accessible prebiotic foods hold strong potential to jointly target gut health and metabolic health in high-risk patients. The BE GONE trial targeted the gut microbiota of obese surveillance patients with a history of colorectal neoplasia through a straightforward bean intervention.

METHODS

This low-risk, non-invasive dietary intervention trial was conducted at MD Anderson Cancer Center (Houston, TX, USA). Following a 4-week equilibration, patients were randomized to continue their usual diet without beans (control) or to add a daily cup of study beans to their usual diet (intervention) with immediate crossover at 8-weeks. Stool and fasting blood were collected every 4 weeks to assess the primary outcome of intra and inter-individual changes in the gut microbiome and in circulating markers and metabolites within 8 weeks. This study was registered on ClinicalTrials.gov as NCT02843425, recruitment is complete and long-term follow-up continues.

FINDINGS

Of the 55 patients randomized by intervention sequence, 87% completed the 16-week trial, demonstrating an increase on-intervention in diversity [n = 48; linear mixed effect and 95% CI for inverse Simpson index: 0.16 (0.02, 0.30); p = 0.02] and shifts in multiple bacteria indicative of prebiotic efficacy, including increased Faecalibacterium, Eubacterium and Bifidobacterium (all p < 0.05). The circulating metabolome showed parallel shifts in nutrient and microbiome-derived metabolites, including increased pipecolic acid and decreased indole (all p < 0.002) that regressed upon returning to the usual diet. No significant changes were observed in circulating lipoproteins within 8 weeks; however, proteomic biomarkers of intestinal and systemic inflammatory response, fibroblast-growth factor-19 increased, and interleukin-10 receptor-α decreased (p = 0.01).

INTERPRETATION

These findings underscore the prebiotic and potential therapeutic role of beans to enhance the gut microbiome and to regulate host markers associated with metabolic obesity and colorectal cancer, while further emphasizing the need for consistent and sustainable dietary adjustments in high-risk patients.

FUNDING

This study was funded by the American Cancer Society.

Hematopoietic stem and progenitor cells confer cross-protective trained immunity in mouse models

Recent studies suggest that infection reprograms hematopoietic stem and progenitor cells (HSPCs) to enhance innate immune responses upon secondary infectious challenge, a process called "trained immunity." However, the specificity and cell types responsible for this response remain poorly defined. We established a model of trained immunity in mice in response to Mycobacterium avium infection. scRNA-seq analysis revealed that HSPCs activate interferon gamma-response genes heterogeneously upon primary challenge, while rare cell populations expand. Macrophages derived from trained HSPCs demonstrated enhanced bacterial killing and metabolism, and a single dose of recombinant interferon gamma exposure was sufficient to induce similar training. Mice transplanted with influenza-trained HSPCs displayed enhanced immunity against M. avium challenge and vice versa, demonstrating cross protection against antigenically distinct pathogens. Together, these results indicate that heterogeneous responses to infection by HSPCs can lead to long-term production of bone marrow derived macrophages with enhanced function and confer cross-protection against alternative pathogens.

High molecular diagnostic yields and novel phenotypic expansions involving syndromic anorectal malformations

Evidence suggests that genetic factors contribute to the development of anorectal malformations (ARMs). However, the etiology of the majority of ARMs cases remains unclear. Exome sequencing (ES) may be underutilized in the diagnostic workup of ARMs due to uncertainty regarding its diagnostic yield. In a clinical database of ~17,000 individuals referred for ES, we identified 130 individuals with syndromic ARMs. A definitive or probable diagnosis was made in 45 of these individuals for a diagnostic yield of 34.6% (45/130). The molecular diagnostic yield of individuals who initially met criteria for VACTERL association was lower than those who did not (26.8% vs 44.1%; p = 0.0437), suggesting that non-genetic factors may play an important role in this subset of syndromic ARM cases. Within this cohort, we identified two individuals who carried de novo pathogenic frameshift variants in ADNP, two individuals who were homozygous for pathogenic variants in BBS1, and single individuals who carried pathogenic or likely pathogenic variants in CREBBP, EP300, FANCC, KDM6A, SETD2, and SMARCA4. The association of these genes with ARMs was supported by previously published cases, and their similarity to known ARM genes as demonstrated using a machine learning algorithm. These data suggest that ES should be considered for all individuals with syndromic ARMs in whom a molecular diagnosis has not been made, and that ARMs represent a low penetrance phenotype associated with Helsmoortel-van der Aa syndrome, Bardet-Biedl syndrome 1, Rubinstein-Taybi syndromes 1 and 2, Fanconi anemia group C, Kabuki syndrome 2, SETD2-related disorders, and Coffin-Siris syndrome 4.

Exome sequencing efficacy and phenotypic expansions involving esophageal atresia/tracheoesophageal fistula plus

Esophageal atresia/tracheoesophageal fistula (EA/TEF) is a life-threatening birth defect that often occurs with other major birth defects (EA/TEF+). Despite advances in genetic testing, a molecular diagnosis can only be made in a minority of EA/TEF+ cases. Here, we analyzed clinical exome sequencing data and data from the DECIPHER database to determine the efficacy of exome sequencing in cases of EA/TEF+ and to identify phenotypic expansions involving EA/TEF. Among 67 individuals with EA/TEF+ referred for clinical exome sequencing, a definitive or probable diagnosis was made in 11 cases for an efficacy rate of 16% (11/67). This efficacy rate is significantly lower than that reported for other major birth defects, suggesting that polygenic, multifactorial, epigenetic, and/or environmental factors may play a particularly important role in EA/TEF pathogenesis. Our cohort included individuals with pathogenic or likely pathogenic variants that affect TCF4 and its downstream target NRXN1, and FANCA, FANCB, and FANCC, which are associated with Fanconi anemia. These cases, previously published case reports, and comparisons to other EA/TEF genes made using a machine learning algorithm, provide evidence in support of a potential pathogenic role for these genes in the development of EA/TEF.

Selection for or against escape from nonsense mediated decay is a novel signature for the detection of cancer genes

Escape from nonsense mediated mRNA decay (NMD-) can produce activated or inactivated gene products, and bias in rates of escape can identify functionally important genes in germline disease. We hypothesized that the same would be true of cancer genes, and tested for NMD- bias within The Cancer Genome Atlas pan-cancer somatic mutation dataset. We identify 29 genes that show significantly elevated or suppressed rates of NMD-. This novel approach to cancer gene discovery reveals genes not previously cataloged as potentially tumorigenic, and identifies many potential driver mutations in known cancer genes for functional characterization.

A joint modeling approach for longitudinal microbiome data improves ability to detect microbiome associations with disease

Changes in the composition of the microbiome over time are associated with myriad human illnesses. Unfortunately, the lack of analytic techniques has hindered researchers' ability to quantify the association between longitudinal microbial composition and time-to-event outcomes. Prior methodological work developed the joint model for longitudinal and time-to-event data to incorporate time-dependent biomarker covariates into the hazard regression approach to disease outcomes. The original implementation of this joint modeling approach employed a linear mixed effects model to represent the time-dependent covariates. However, when the distribution of the time-dependent covariate is non-Gaussian, as is the case with microbial abundances, researchers require different statistical methodology. We present a joint modeling framework that uses a negative binomial mixed effects model to determine longitudinal taxon abundances. We incorporate these modeled microbial abundances into a hazard function with a parameterization that not only accounts for the proportional nature of microbiome data, but also generates biologically interpretable results. Herein we demonstrate the performance improvements of our approach over existing alternatives via simulation as well as a previously published longitudinal dataset studying the microbiome during pregnancy. The results demonstrate that our joint modeling framework for longitudinal microbiome count data provides a powerful methodology to uncover associations between changes in microbial abundances over time and the onset of disease. This method offers the potential to equip researchers with a deeper understanding of the associations between longitudinal microbial composition changes and disease outcomes. This new approach could potentially lead to new diagnostic biomarkers or inform clinical interventions to help prevent or treat disease.

Increased Moraxella and Streptococcus species abundance after severe bronchiolitis is associated with recurrent wheezing

BACKGROUND

The role of the airway microbiome in the development of recurrent wheezing and asthma remains uncertain, particularly in the high-risk group of infants hospitalized for bronchiolitis.

OBJECTIVE

We sought to examine the relation of the nasal microbiota at bronchiolitis-related hospitalization and 3 later points to the risk of recurrent wheezing by age 3 years.

METHODS

In 17 US centers researchers collected clinical data and nasal swabs from infants hospitalized for bronchiolitis. Trained parents collected nasal swabs 3 weeks after hospitalization and, when healthy, during the summer and 1 year after hospitalization. We applied 16S rRNA gene sequencing to all nasal swabs. We used joint modeling to examine the relation of longitudinal nasal microbiota abundances to the risk of recurrent wheezing.

RESULTS

Among 842 infants hospitalized for bronchiolitis, there was 88% follow-up at 3 years, and 31% had recurrent wheezing. The median age at enrollment was 3.2 months (interquartile range, 1.7-5.8 months). In joint modeling analyses adjusting for 16 covariates, including viral cause, a 10% increase in relative abundance of Moraxella or Streptococcus species 3 weeks after day 1 of hospitalization was associated with an increased risk of recurrent wheezing (hazard ratio [HR] of 1.38 and 95% high-density interval [HDI] of 1.11-1.85 and HR of 1.76 and 95% HDI of 1.13-3.19, respectively). Increased Streptococcus species abundance the summer after hospitalization was also associated with a greater risk of recurrent wheezing (HR, 1.76; 95% HDI, 1.15-3.27).

CONCLUSIONS

Enrichment of Moraxella or Streptococcus species after bronchiolitis hospitalization was associated with recurrent wheezing by age 3 years, possibly providing new avenues to ameliorate the long-term respiratory outcomes of infants with severe bronchiolitis.

The association between anterior nares and nasopharyngeal microbiota in infants hospitalized for bronchiolitis

BACKGROUND

The airway microbiome is a subject of great interest for the study of respiratory disease. Anterior nare samples are more accessible than samples from deeper within the nasopharynx. However, the correlation between the microbiota found in the anterior nares and the microbiota found within the nasopharynx is unknown. We assessed the anterior nares and nasopharyngeal microbiota to determine (1) the relation of the microbiota from these two upper airway sites and (2) if associations were maintained between the microbiota from these two sites and two bronchiolitis severity outcomes.

RESULTS

Among 815 infants hospitalized at 17 US centers for bronchiolitis with optimal 16S rRNA gene sequence reads from both nasal swab and nasopharyngeal aspirate samples, there were strong intra-individual correlations in the microbial communities between the two sample types, especially relating to Haemophilus and Moraxella genera. By contrast, we found a high abundance of Staphylococcus genus in the nasal swabs-a pattern not found in the nasopharyngeal samples and not informative when predicting the dominant nasopharyngeal genera. While these disparities may have been due to sample processing differences (i.e., nasal swabs were mailed at ambient temperature to emulate processing of future parent collected swabs while nasopharyngeal aspirates were mailed on dry ice), a previously reported association between Haemophilus-dominant nasopharyngeal microbiota and the increased severity of bronchiolitis was replicated utilizing the nasal swab microbiota and the same outcome measures: intensive care use (adjusted OR 6.43; 95% CI 2.25-20.51; P < 0.001) and hospital length-of-stay (adjusted OR 4.31; 95% CI, 1.73-11.11; P = 0.002). Additionally, Moraxella-dominant nasopharyngeal microbiota was previously identified as protective against intensive care use, a result that was replicated when analyzing the nasal swab microbiota (adjusted OR 0.30; 95% CI, 0.11-0.64; P = 0.01).

CONCLUSIONS

While the microbiota of the anterior nares and the nasopharynx are distinct, there is considerable overlap between the bacterial community compositions from these two anatomic sites. Despite processing differences between the samples, these results indicate that microbiota severity associations from the nasopharynx are recapitulated in the anterior nares, suggesting that nasal swab samples not only are effective sample types, but also can be used to detect microbial risk markers.

Cesarean or Vaginal Birth Does Not Impact the Longitudinal Development of the Gut Microbiome in a Cohort of Exclusively Preterm Infants

The short and long-term impact of birth mode on the developing gut microbiome in neonates has potential implications for the health of infants. In term infants, the microbiome immediately following birth across multiple body sites corresponds to birth mode, with increased Bacteroides in vaginally delivered infants. We aimed to determine the impact of birth mode of the preterm gut microbiome over the first 100 days of life and following neonatal intensive care unit (NICU) discharge. In total, 867 stool samples from 46 preterm infants (21 cesarean and 25 vaginal), median gestational age 27 weeks, were sequenced (V4 region 16S rRNA gene, Illumina MiSeq). Of these, 776 samples passed quality filtering and were included in the analysis. The overall longitudinal alpha-diversity and within infant beta-diversity was comparable between cesarean and vaginally delivered infants. Vaginally delivered infants kept significantly more OTUs from 2 months of life and following NICU discharge, but OTUs lost, gained, and regained were not different based on birth mode. Furthermore, the temporal progression of dominant genera was comparable between birth modes and no significant difference was found for any genera following adjustment for covariates. Lastly, preterm gut community types (PGCTs) showed some moderate differences in very early life, but progressed toward a comparable pattern by week 5. No PGCT was significantly associated with cesarean or vaginal birth. Unlike term infants, birth mode was not significantly associated with changes in microbial diversity, composition, specific taxa, or overall microbial development in preterm infants. This may result from the dominating effects of NICU exposures including the universal use of antibiotics immediately following birth and/or the lack of Bacteroides colonizing preterm infants.

Cesarean or Vaginal Birth Does Not Impact the Longitudinal Development of the Gut Microbiome in a Cohort of Exclusively Preterm Infants

The short and long-term impact of birth mode on the developing gut microbiome in neonates has potential implications for the health of infants. In term infants, the microbiome immediately following birth across multiple body sites corresponds to birth mode, with increased Bacteroides in vaginally delivered infants. We aimed to determine the impact of birth mode of the preterm gut microbiome over the first 100 days of life and following neonatal intensive care unit (NICU) discharge. In total, 867 stool samples from 46 preterm infants (21 cesarean and 25 vaginal), median gestational age 27 weeks, were sequenced (V4 region 16S rRNA gene, Illumina MiSeq). Of these, 776 samples passed quality filtering and were included in the analysis. The overall longitudinal alpha-diversity and within infant beta-diversity was comparable between cesarean and vaginally delivered infants. Vaginally delivered infants kept significantly more OTUs from 2 months of life and following NICU discharge, but OTUs lost, gained, and regained were not different based on birth mode. Furthermore, the temporal progression of dominant genera was comparable between birth modes and no significant difference was found for any genera following adjustment for covariates. Lastly, preterm gut community types (PGCTs) showed some moderate differences in very early life, but progressed toward a comparable pattern by week 5. No PGCT was significantly associated with cesarean or vaginal birth. Unlike term infants, birth mode was not significantly associated with changes in microbial diversity, composition, specific taxa, or overall microbial development in preterm infants. This may result from the dominating effects of NICU exposures including the universal use of antibiotics immediately following birth and/or the lack of Bacteroides colonizing preterm infants.