Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition.

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

Ancient Clostridium DNA and variants of tetanus neurotoxins associated with human archaeological remains.

The analysis of microbial genomes from human archaeological samples offers a historic snapshot of ancient pathogens and provides insights into the origins of modern infectious diseases. Here, we analyze metagenomic datasets from 38 human archaeological samples and identify bacterial genomic sequences related to modern-day Clostridium tetani, which produces the tetanus neurotoxin (TeNT) and causes the disease tetanus. These genomic assemblies had varying levels of completeness, and a subset of them displayed hallmarks of ancient DNA damage. Phylogenetic analyses revealed known C. tetani clades as well as potentially new Clostridium lineages closely related to C. tetani. The genomic assemblies encode 13 TeNT variants with unique substitution profiles, including a subgroup of TeNT variants found exclusively in ancient samples from South America. We experimentally tested a TeNT variant selected from an ancient Chilean mummy sample and found that it induced tetanus muscle paralysis in mice, with potency comparable to modern TeNT. Thus, our ancient DNA analysis identifies DNA from neurotoxigenic C. tetani in archaeological human samples, and a novel variant of TeNT that can cause disease in mammals.

Mechanisms of epigenetic regulation of transcription by lncRNAs in plants.

Long noncoding RNAs (lncRNAs) are an ubiquitous feature of eukaryotic genomes and in recent decades have been shown to be highly abundant and varied. Many prominent examples have been described as having essential roles in regulating the expression of genes in different developmental and environmental contexts. As a result, much work has been done on elucidating the mechanisms by which they modulate the expression of protein coding genes. In this review, we focus on those which have been characterized in plants. We specifically examine common epigenetic mechanisms that regulate gene expression at the level of transcription. In this regard, we focus on the lncRNAs in plants that have primarily been associated with controlling the chromatin environment of genes at the level of modifications, RNA POLYMERASE II (RNAPII) processivity and efficiency of transcription, and mediating the formation of transcriptionally activating and repressive chromatin loops. We discuss open questions in plant lncRNA epigenetic regulation and opportunities for future study of functionally significant lncRNAs with yet-unknown epigenetic mechanisms.

Enhancing diversity analysis by repeatedly rarefying next generation sequencing data describing microbial communities.

Amplicon sequencing has revolutionized our ability to study DNA collected from environmental samples by providing a rapid and sensitive technique for microbial community analysis that eliminates the challenges associated with lab cultivation and taxonomic identification through microscopy. In water resources management, it can be especially useful to evaluate ecosystem shifts in response to natural and anthropogenic landscape disturbances to signal potential water quality concerns, such as the detection of toxic cyanobacteria or pathogenic bacteria. Amplicon sequencing data consist of discrete counts of sequence reads, the sum of which is the library size. Groups of samples typically have different library sizes that are not representative of biological variation; library size normalization is required to meaningfully compare diversity between them. Rarefaction is a widely used normalization technique that involves the random subsampling of sequences from the initial sample library to a selected normalized library size. This process is often dismissed as statistically invalid because subsampling effectively discards a portion of the observed sequences, yet it remains prevalent in practice and the suitability of rarefying, relative to many other normalization approaches, for diversity analysis has been argued. Here, repeated rarefying is proposed as a tool to normalize library sizes for diversity analyses. This enables (i) proportionate representation of all observed sequences and (ii) characterization of the random variation introduced to diversity analyses by rarefying to a smaller library size shared by all samples. While many deterministic data transformations are not tailored to produce equal library sizes, repeatedly rarefying reflects the probabilistic process by which amplicon sequencing data are obtained as a representation of the amplified source microbial community. Specifically, it evaluates which data might have been obtained if a particular sample's library size had been smaller and allows graphical representation of the effects of this library size normalization process upon diversity analysis results.

Experimental and natural evidence of SARS-CoV-2-infection-induced activation of type I interferon responses.

Type I interferons (IFNs) are our first line of defense against virus infection. Recent studies have suggested the ability of SARS-CoV-2 proteins to inhibit IFN responses. Emerging data also suggest that timing and extent of IFN production is associated with manifestation of COVID-19 severity. In spite of progress in understanding how SARS-CoV-2 activates antiviral responses, mechanistic studies into wild-type SARS-CoV-2-mediated induction and inhibition of human type I IFN responses are scarce. Here we demonstrate that SARS-CoV-2 infection induces a type I IFN response in vitro and in moderate cases of COVID-19. In vitro stimulation of type I IFN expression and signaling in human airway epithelial cells is associated with activation of canonical transcriptions factors, and SARS-CoV-2 is unable to inhibit exogenous induction of these responses. Furthermore, we show that physiological levels of IFNα detected in patients with moderate COVID-19 is sufficient to suppress SARS-CoV-2 replication in human airway cells.

PhyloCorrelate: inferring bacterial gene-gene functional associations through large-scale phylogenetic profiling.

MOTIVATION: Statistical detection of co-occurring genes across genomes, known as 'phylogenetic profiling', is a powerful bioinformatic technique for inferring gene-gene functional associations. However, this can be a challenging task given the size and complexity of phylogenomic databases, difficulty in accounting for phylogenetic structure, inconsistencies in genome annotation and substantial computational requirements. RESULTS: We introduce PhyloCorrelate-a computational framework for gene co-occurrence analysis across large phylogenomic datasets. PhyloCorrelate implements a variety of co-occurrence metrics including standard correlation metrics and model-based metrics that account for phylogenetic history. By combining multiple metrics, we developed an optimized score that exhibits a superior ability to link genes with overlapping GO terms and KEGG pathways, enabling gene function prediction. Using genomic and functional annotation data from the Genome Taxonomy Database and AnnoTree, we performed all-by-all comparisons of gene occurrence profiles across the bacterial tree of life, totaling 154 217 052 comparisons for 28 315 genes across 27 372 bacterial genomes. All predictions are available in an online database, which instantaneously returns the top correlated genes for any PFAM, TIGRFAM or KEGG query. In total, PhyloCorrelate detected 29 762 high confidence associations between bacterial gene/protein pairs, and generated functional predictions for 834 DUFs and proteins of unknown function. AVAILABILITYAND IMPLEMENTATION: PhyloCorrelate is available as a web-server at phylocorrelate.uwaterloo.ca as well as an R package for analysis of custom datasets. We anticipate that PhyloCorrelate will be broadly useful as a tool for predicting function and interactions for gene families. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Endogenous CRISPR-Cas Systems in Group I Clostridium botulinum and Clostridium sporogenes Do Not Directly Target the Botulinum Neurotoxin Gene Cluster.

Most strains of proteolytic group I Clostridium botulinum (G1 C. botulinum) and some strains of Clostridium sporogenes possess genes encoding botulinum neurotoxin (BoNT), a potent neuroparalytic agent. Within G1 C. botulinum, conserved bont gene clusters of three major toxin serotypes (bont/A/B/F) can be found on conjugative plasmids and/or within chromosomal pathogenicity islands. CRISPR-Cas systems enable site-specific targeting of previously encountered mobile genetic elements (MGE) such as plasmids and bacteriophage through the creation of a spacer library complementary to protospacers within the MGEs. To examine whether endogenous CRISPR-Cas systems restrict the transfer of bont gene clusters across strains we conducted a bioinformatic analysis profiling endogenous CRISPR-Cas systems from 241 G1 C. botulinum and C. sporogenes strains. Approximately 6,200 CRISPR spacers were identified across the strains and Type I-B, III-A/B/D cas genes and CRISPR array features were identified in 83% of the strains. Mapping the predicted spacers against the masked strain and RefSeq plasmid dataset identified 56,000 spacer-protospacer matches. While spacers mapped heavily to targets within bont(+) plasmids, no protospacers were identified within the bont gene clusters. These results indicate the toxin is not a direct target of CRISPR-Cas but the plasmids predominantly responsible for its mobilization are. Finally, while the presence of a CRISPR-Cas system did not reliably indicate the presence or absence of a bont gene cluster, comparative genomics across strains indicates they often occupy the same hypervariable loci common to both species, potentially suggesting similar mechanisms are involved in the acquisition and curation of both genomic features.

Expression of endocannabinoid system components in human airway epithelial cells: impact of sex and chronic respiratory disease status.

Cannabis smoking is the dominant route of delivery, with the airway epithelium functioning as the site of first contact. The endocannabinoid system is responsible for mediating the physiological effects of inhaled phytocannabinoids. The expression of the endocannabinoid system in the airway epithelium and contribution to normal physiological responses remains to be defined. To begin to address this knowledge gap, a curated dataset of 1090 unique human bronchial brushing gene expression profiles was created. The dataset included 616 healthy subjects, 136 subjects with asthma, and 338 subjects with COPD. A 32-gene endocannabinoid signature was analysed across all samples with sex and disease-specific analyses performed. Immunohistochemistry and immunoblots were performed to probe in situ and in vitro protein expression. CB1, CB2, and TRPV1 protein signal is detectable in human airway epithelial cells in situ and in vitro, justifying examining the downstream endocannabinoid pathway. Sex status was associated with differential expression of 7 of 32 genes. In contrast, disease status was associated with differential expression of 21 of 32 genes in people with asthma and 26 of 32 genes in people with COPD. We confirm at the protein level that TRPV1, the most differentially expressed candidate in our analyses, was upregulated in airway epithelial cells from people with asthma relative to healthy subjects. Our data demonstrate that the endocannabinoid system is expressed in human airway epithelial cells with expression impacted by disease status and minimally by sex. The data suggest that cannabis consumers may have differential physiological responses in the respiratory mucosa.

Phylogenomics of 8,839 Clostridioides difficile genomes reveals recombination-driven evolution and diversification of toxin A and B.

Clostridioides difficile is the major worldwide cause of antibiotic-associated gastrointestinal infection. A pathogenicity locus (PaLoc) encoding one or two homologous toxins, toxin A (TcdA) and toxin B (TcdB), is essential for C. difficile pathogenicity. However, toxin sequence variation poses major challenges for the development of diagnostic assays, therapeutics, and vaccines. Here, we present a comprehensive phylogenomic analysis of 8,839 C. difficile strains and their toxins including 6,492 genomes that we assembled from the NCBI short read archive. A total of 5,175 tcdA and 8,022 tcdB genes clustered into 7 (A1-A7) and 12 (B1-B12) distinct subtypes, which form the basis of a new method for toxin-based subtyping of C. difficile. We developed a haplotype coloring algorithm to visualize amino acid variation across all toxin sequences, which revealed that TcdB has diversified through extensive homologous recombination throughout its entire sequence, and formed new subtypes through distinct recombination events. In contrast, TcdA varies mainly in the number of repeats in its C-terminal repetitive region, suggesting that recombination-mediated diversification of TcdB provides a selective advantage in C. difficile evolution. The application of toxin subtyping is then validated by classifying 351 C. difficile clinical isolates from Brigham and Women's Hospital in Boston, demonstrating its clinical utility. Subtyping partitions TcdB into binary functional and antigenic groups generated by intragenic recombinations, including two distinct cell-rounding phenotypes, whether recognizing frizzled proteins as receptors, and whether it can be efficiently neutralized by monoclonal antibody bezlotoxumab, the only FDA-approved therapeutic antibody. Our analysis also identifies eight universally conserved surface patches across the TcdB structure, representing ideal targets for developing broad-spectrum therapeutics. Finally, we established an open online database (DiffBase) as a central hub for collection and classification of C. difficile toxins, which will help clinicians decide on therapeutic strategies targeting specific toxin variants, and allow researchers to monitor the ongoing evolution and diversification of C. difficile.

ABCF1 Regulates dsDNA-induced Immune Responses in Human Airway Epithelial Cells.

Background: The airway epithelium represents a critical component of the human lung that helps orchestrate defenses against respiratory tract viral infections, which are responsible for more than 2.5 million deaths/year globally. Innate immune activities of the airway epithelium rely on Toll-like receptors (TLRs), nucleotide binding and leucine-rich-repeat pyrin domain containing (NLRP) receptors, and cytosolic nucleic acid sensors. ATP Binding Cassette (ABC) transporters are ubiquitous across all three domains of life-Archaea, Bacteria, and Eukarya-and expressed in the human airway epithelium. ABCF1, a unique ABC family member that lacks a transmembrane domain, has been defined as a cytosolic nucleic acid sensor that regulates CXCL10, interferon-ß expression, and downstream type I interferon responses. We tested the hypothesis that ABCF1 functions as a dsDNA nucleic acid sensor in human airway epithelial cells important in regulating antiviral responses. Methods: Expression and localization experiments were performed using in situ hybridization and immunohistochemistry in human lung tissue, while confirmatory transcript and protein expression was performed in human airway epithelial cells. Functional experiments were performed with siRNA methods in a human airway epithelial cell line. Complementary transcriptomic analyses were performed to explore the contributions of ABCF1 to gene expression patterns. Results: Using archived human lung and human airway epithelial cells, we confirm expression of ABCF1 gene and protein expression in these tissue samples, with a role for mediating CXCL10 production in response to dsDNA viral mimic challenge. Although, ABCF1 knockdown was associated with an attenuation of select genes involved in the antiviral responses, Gene Ontology analyses revealed a greater interaction of ABCF1 with TLR signaling suggesting a multifactorial role for ABCF1 in innate immunity in human airway epithelial cells. Conclusion: ABCF1 is a candidate cytosolic nucleic acid sensor and modulator of TLR signaling that is expressed at gene and protein levels in human airway epithelial cells. The precise level where ABCF1 protein functions to modulate immune responses to pathogens remains to be determined but is anticipated to involve IRF-3 and CXCL10 production.

Predicting the recombination potential of severe acute respiratory syndrome coronavirus 2 and Middle East respiratory syndrome coronavirus.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) recently emerged to cause widespread infections in humans. SARS-CoV-2 infections have been reported in the Kingdom of Saudi Arabia, where Middle East respiratory syndrome coronavirus (MERS-CoV) causes seasonal outbreaks with a case fatality rate of ~37 %. Here we show that there exists a theoretical possibility of future recombination events between SARS-CoV-2 and MERS-CoV RNA. Through computational analyses, we have identified homologous genomic regions within the ORF1ab and S genes that could facilitate recombination, and have analysed co-expression patterns of the cellular receptors for SARS-CoV-2 and MERS-CoV, ACE2 and DPP4, respectively, to identify human anatomical sites that could facilitate co-infection. Furthermore, we have investigated the likely susceptibility of various animal species to MERS-CoV and SARS-CoV-2 infection by comparing known virus spike protein-receptor interacting residues. In conclusion, we suggest that a recombination between SARS-CoV-2 and MERS-CoV RNA is possible and urge public health laboratories in high-risk areas to develop diagnostic capability for the detection of recombined coronaviruses in patient samples.

Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue.

In December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged, causing the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV, the agent responsible for the 2003 SARS outbreak, utilises angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) host molecules for viral entry. ACE2 and TMPRSS2 have recently been implicated in SARS-CoV-2 viral infection. Additional host molecules including ADAM17, cathepsin L, CD147 and GRP78 may also function as receptors for SARS-CoV-2.To determine the expression and in situ localisation of candidate SARS-CoV-2 receptors in the respiratory mucosa, we analysed gene expression datasets from airway epithelial cells of 515 healthy subjects, gene promoter activity analysis using the FANTOM5 dataset containing 120 distinct sample types, single cell RNA sequencing (scRNAseq) of 10 healthy subjects, proteomic datasets, immunoblots on multiple airway epithelial cell types, and immunohistochemistry on 98 human lung samples.We demonstrate absent to low ACE2 promoter activity in a variety of lung epithelial cell samples and low ACE2 gene expression in both microarray and scRNAseq datasets of epithelial cell populations. Consistent with gene expression, rare ACE2 protein expression was observed in the airway epithelium and alveoli of human lung, confirmed with proteomics. We present confirmatory evidence for the presence of TMPRSS2, CD147 and GRP78 protein in vitro in airway epithelial cells and confirm broad in situ protein expression of CD147 and GRP78 in the respiratory mucosa.Collectively, our data suggest the presence of a mechanism dynamically regulating ACE2 expression in human lung, perhaps in periods of SARS-CoV-2 infection, and also suggest that alternative receptors for SARS-CoV-2 exist to facilitate initial host cell infection.