Comparison of DNA sequences with protein sequences

Community differences and potential function along the particle size spectrum of microbes in the twilight zone

BACKGROUND

The twilight zone, which extends from the base of the euphotic zone to a depth of 1000 m, is the major area of particulate organic carbon (POC) remineralization in the ocean. However, little is known about the microbial community and metabolic activity that are directly associated with POC remineralization in this consistently underexplored realm. Here, we utilized a large-volume in situ water transfer system to collect the microbes on different-sized particles from the twilight zone in three regions and analyzed their composition and metabolic function by metagenomic analysis.

RESULTS

Distinct prokaryotic communities with significantly lower diversity and less endemic species were detected on particles in the South East Asian Time-series Study (SEATS) compared with the other two regions, perhaps due to the in situ physicochemical conditions and low labile nutrient availability in this region. Observable transitions in community composition and function at the upper and lower boundaries of the twilight zone suggest that microbes respond differently to (and potentially drive the transformation of) POC through this zone. Substantial variations among different particle sizes were observed, with smaller particles typically exhibiting lower diversity but harboring a greater abundance of carbon degradation-associated genes than the larger particles. Such a pattern might arise due to the relatively larger surface area of the smaller particles relative to their volume, which likely provides more sites for microbial colonization, increasing their chance of being remineralized. This makes them less likely to be transferred to the deep ocean, and thus, they contribute more to carbon recycling than to long-term sequestration. Both contig-based and metagenome-assembled genome-(MAG-) based analyses revealed a high diversity of the Carbohydrate-Active enZymes (CAZy) family. This indicates the versatile carbohydrate metabolisms of the microbial communities associated with sinking particles that modulate the remineralization and export of POC in the twilight zone.

CONCLUSION

Our study reveals significant shifts in microbial community composition and function in the twilight zone, with clear differences among the three particle sizes. Microbes with diverse metabolic potential exhibited different responses to the POC entering the twilight zone and also collectively drove the transformation of POC through this zone. These findings provided insights into the diversity of prokaryotes in sinking particles and their roles in POC remineralization and export in marine ecosystems. Video Abstract.

Functionally characterizing obesity-susceptibility genes using CRISPR/Cas9, in vivo imaging and deep learning

Hundreds of loci have been robustly associated with obesity-related traits, but functional characterization of candidate genes remains a bottleneck. Aiming to systematically characterize candidate genes for a role in accumulation of lipids in adipocytes and other cardiometabolic traits, we developed a pipeline using CRISPR/Cas9, non-invasive, semi-automated fluorescence imaging and deep learning-based image analysis in live zebrafish larvae. Results from a dietary intervention show that 5 days of overfeeding is sufficient to increase the odds of lipid accumulation in adipocytes by 10 days post-fertilization (dpf, n = 275). However, subsequent experiments show that across 12 to 16 established obesity genes, 10 dpf is too early to detect an effect of CRISPR/Cas9-induced mutations on lipid accumulation in adipocytes (n = 1014), and effects on food intake at 8 dpf (n = 1127) are inconsistent with earlier results from mammals. Despite this, we observe effects of CRISPR/Cas9-induced mutations on ectopic accumulation of lipids in the vasculature (sh2b1 and sim1b) and liver (bdnf); as well as on body size (pcsk1, pomca, irs1); whole-body LDLc and/or total cholesterol content (irs2b and sh2b1); and pancreatic beta cell traits and/or glucose content (pcsk1, pomca, and sim1a). Taken together, our results illustrate that CRISPR/Cas9- and image-based experiments in zebrafish larvae can highlight direct effects of obesity genes on cardiometabolic traits, unconfounded by their - not yet apparent - effect on excess adiposity.

Co-translational protein aggregation and ribosome stalling as a broad-spectrum antibacterial mechanism

Drug-resistant bacteria pose an urgent global health threat, necessitating the development of antibacterial compounds with novel modes of action. Protein biosynthesis accounts for up to half of the energy expenditure of bacterial cells, and consequently inhibiting the efficiency or fidelity of the bacterial ribosome is a major target of existing antibiotics. Here, we describe an alternative mode of action that affects the same process: allowing translation to proceed but causing co-translational aggregation of the nascent peptidic chain. We show that treatment with an aggregation-prone peptide induces formation of polar inclusion bodies and activates the SsrA ribosome rescue pathway in bacteria. The inclusion bodies contain ribosomal proteins and ribosome hibernation factors, as well as mRNAs and cognate nascent chains of many proteins in amyloid-like structures, with a bias for membrane proteins with a fold rich in long-range beta-sheet interactions. The peptide is bactericidal against a wide range of pathogenic bacteria in planktonic growth and in biofilms, and reduces bacterial loads in mouse models of Escherichia coli and Acinetobacter baumannii infections. Our results indicate that disrupting protein homeostasis via co-translational aggregation constitutes a promising strategy for development of broad-spectrum antibacterials.

Identification and Analysis of Circular RNAs in Mammary Gland from Yaks Between Lactation and Dry Period

Lactation is a complex physiological process regulated by numerous genes and factors. Circular RNA (circRNA), a non-coding RNA, acts as a molecular sponge that sequesters microRNAs (miRNAs) to regulate target gene expression. Although circRNA has been linked to mammary gland lactation, its specific role in yaks remains underexplored. This study employed circular RNA sequencing (circRNA-seq) to examine the differential expression of circRNAs in yak mammary tissues during lactation and the dry period. Additionally, an enrichment analysis of the differentially expressed circRNAs (DECs) was performed. A competing endogenous RNA (ceRNA) network was then constructed to explore the potential of their roles in lactation and mammary gland development. We detected 18,905 circRNAs in yak mammary tissue, among which 302 showed differential expression. The host genes of these DECs were enriched in functions and pathways associated with yak milk synthesis and composition. Through the construction of a ceRNA network and the enrichment analysis of associated mRNAs, this study identified ceRNAs potentially involved in regulating lactation and mammary gland development. In conclusion, circRNAs in yak mammary tissues were identified and analyzed across lactation and dry periods, establishing a ceRNA network related to lactation regulation. These findings provide novel insights into the regulatory mechanisms governing lactation in yaks (Bos grunniens).

Gra-CRC-miRTar: The pre-trained nucleotide-to-graph neural networks to identify potential miRNA targets in colorectal cancer

Colorectal cancer (CRC) is the third most diagnosed cancer and the second deadliest cancer worldwide representing a major public health problem. In recent years, increasing evidence has shown that microRNA (miRNA) can control the expression of targeted human messenger RNA (mRNA) by reducing their abundance or translation, acting as oncogenes or tumor suppressors in various cancers, including CRC. Due to the significant up-regulation of oncogenic miRNAs in CRC, elucidating the underlying mechanism and identifying dysregulated miRNA targets may provide a basis for improving current therapeutic interventions. In this paper, we proposed Gra-CRC-miRTar, a pre-trained nucleotide-to-graph neural network framework, for identifying potential miRNA targets in CRC. Different from previous studies, we constructed two pre-trained models to encode RNA sequences and transformed them into de Bruijn graphs. We employed different graph neural networks to learn the latent representations. The embeddings generated from de Bruijn graphs were then fed into a Multilayer Perceptron (MLP) to perform the prediction tasks. Our extensive experiments show that Gra-CRC-miRTar achieves better performance than other deep learning algorithms and existing predictors. In addition, our analyses also successfully revealed 172 out of 201 functional interactions through experimentally validated miRNA-mRNA pairs in CRC. Collectively, our effort provides an accurate and efficient framework to identify potential miRNA targets in CRC, which can also be used to reveal miRNA target interactions in other malignancies, facilitating the development of novel therapeutics. The Gra-CRC-miRTar web server can be found at: http://gra-crc-mirtar.com/.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

Summary

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

Availability and implementation

The software is available at https://github.com/TravisWheelerLab/BATH.

Gra-CRC-miRTar: The pre-trained nucleotide-to-graph neural networks to identify potential miRNA targets in colorectal cancer

Colorectal cancer (CRC) is the third most diagnosed cancer and the second deadliest cancer worldwide representing a major public health problem. In recent years, increasing evidence has shown that microRNA (miRNA) can control the expression of targeted human messenger RNA (mRNA) by reducing their abundance or translation, acting as oncogenes or tumor suppressors in various cancers, including CRC. Due to the significant up-regulation of oncogenic miRNAs in CRC, elucidating the underlying mechanism and identifying dysregulated miRNA targets may provide a basis for improving current therapeutic interventions. In this paper, we proposed Gra-CRC-miRTar, a pre-trained nucleotide-to-graph neural network framework, for identifying potential miRNA targets in CRC. Different from previous studies, we constructed two pre-trained models to encode RNA sequences and transformed them into de Bruijn graphs. We employed different graph neural networks to learn the latent representations. The embeddings generated from de Bruijn graphs were then fed into a Multilayer Perceptron (MLP) to perform the prediction tasks. Our extensive experiments show that Gra-CRC-miRTar achieves better performance than other deep learning algorithms and existing predictors. In addition, our analyses also successfully revealed 172 out of 201 functional interactions through experimentally validated miRNA-mRNA pairs in CRC. Collectively, our effort provides an accurate and efficient framework to identify potential miRNA targets in CRC, which can also be used to reveal miRNA target interactions in other malignancies, facilitating the development of novel therapeutics.

Identification of potential pathogenic targets and survival strategies of Vibrio vulnificus through population genomics

Vibrio vulnificus, a foodborne pathogen, has a high mortality rate. Despite its relevance to public health, the identification of virulence genes associated with the pathogenicity of currently known clinical isolates of V. vulnificus is incomplete and its synergistic pathogenesis remains unclear. Here, we integrate whole genome sequencing (WGS), genome-wide association studies (GWAS), and genome-wide epistasis studies (GWES), along with phenotype characterization to investigate the pathogenesis and survival strategies of V. vulnificus. GWAS and GWES identified a total of six genes (purH, gmr, yiaV, dsbD, ramA, and wbpA) associated with the pathogenicity of clinical isolates related to nucleotide/amino acid transport and metabolism, cell membrane biogenesis, signal transduction mechanisms, and protein turnover. Of these, five were newly discovered potential specific virulence genes of V. vulnificus in this study. Furthermore, GWES combined with phenotype experiments indicated that V. vulnificus isolates were clustered into two ecological groups (EGs) that shared distinct biotic and abiotic factors, and ecological strategies. Our study reveals pathogenic mechanisms and their evolution in V. vulnificus to provide a solid foundation for designing new vaccines and therapeutic targets.

The hidden diversity of ancient bornaviral sequences from X and P genes in vertebrate genomes

Endogenous bornavirus-like elements (EBLs) are heritable sequences derived from bornaviruses in vertebrate genomes that originate from transcripts of ancient bornaviruses. EBLs have been detected using sequence similarity searches such as tBLASTn, whose technical limitations may hinder the detection of EBLs derived from small and/or rapidly evolving viral X and P genes. Indeed, no EBLs derived from the X and P genes of orthobornaviruses have been detected to date in vertebrate genomes. Here, we aimed to develop a novel strategy to detect such 'hidden' EBLs. To this aim, we focused on the 1.9-kb read-through transcript of orthobornaviruses, which encodes a well-conserved N gene and small and rapidly evolving X and P genes. We show a series of evidence supporting the existence of EBLs derived from orthobornaviral X and P genes (EBLX/Ps) in mammalian genomes. Furthermore, we found that an EBLX/P is expressed as a fusion transcript with the cellular gene, ZNF451, which potentially encodes the ZNF451/EBLP fusion protein in miniopterid bat cells. This study contributes to a deeper understanding of ancient bornaviruses and co-evolution between bornaviruses and their hosts. Furthermore, our data suggest that endogenous viral elements are more abundant than those previously appreciated using BLAST searches alone, and further studies are required to understand ancient viruses more accurately.

Iva xanthiifolia leaf extract reduced the diversity of indigenous plant rhizosphere bacteria

BACKGROUND

Iva xanthiifolia, native to North America, is now widely distributed in northeastern China and has become a vicious invasive plant. This article aims to probe the role of leaf extract in the invasion of I. xanthiifolia.

METHODS

We collected the rhizosphere soil of Amaranthus tricolor and Setaria viridis in the invasive zone, the noninvasive zone and the noninvasive zone treated with extract from I. xanthiifolia leaf, and obtained I. xanthiifolia rhizosphere soil in the invasive zone. All wild plants were identified by Xu Yongqing. I. xanthiifolia (collection number: RQSB04100), A. tricolor (collection number: 831,030) and S. viridis (collection number: CF-0002-034) are all included in Chinese Virtual Herbarium ( https://www.cvh.ac.cn/index.php ). The soil bacterial diversity was analyzed based on the Illumina HiSeq sequencing platform. Subsequently, taxonomic analysis and Faprotax functional prediction were performed.

RESULTS

The results showed that the leaf extract significantly reduced the diversity of indigenous plant rhizosphere bacteria. A. tricolor and S. viridis rhizobacterial phylum and genus abundances were significantly reduced under the influence of I. xanthiifolia or its leaf extract. The results of functional prediction showed that bacterial abundance changes induced by leaf extracts could potentially hinder nutrient cycling in native plants and increased bacterial abundance in the A. tricolor rhizosphere related to aromatic compound degradation. In addition, the greatest number of sensitive Operational Taxonomic Units (OTUs) appeared in the rhizosphere when S. viridis was in response to the invasion of I. xanthiifolia. It can be seen that A. tricolor and S. viridis have different mechanisms in response to the invasion of I. xanthiifolia.

CONCLUSION

I. xanthiifolia leaves material has potential role in invasion by altering indigenous plant rhizosphere bacteria.

Comparative oxidation proteomics analyses suggest redox regulation of cytosolic translation in rice leaves upon Magnaporthe oryzae infection

Pathogen attack can increase plant levels of reactive oxygen species (ROS), which act as signaling molecules to activate plant defense mechanisms. Elucidating these processes is crucial for understanding redox signaling pathways in plant defense responses. Using an iodo-tandem mass tag (TMT)-based quantitative proteomics approach, we mapped 3362 oxidized cysteine sites in 2275 proteins in rice leaves. Oxidized proteins were involved in gene expression, peptide biosynthetic processes, stress responses, ROS metabolic processes, and translation pathways. Magnaporthe oryzae infection led to increased oxidative modification levels of 512 cysteine sites in 438 proteins, including many transcriptional regulators and ribosomal proteins. Ribosome profiling (Ribo-seq) analysis revealed that the oxidative modification of ribosomal proteins promoted the translational efficiency of many mRNAs involved in defense response pathways, thereby affecting rice immunity. Our results suggest that increased oxidative modification of ribosomal proteins in rice leaves promotes cytosolic translation, thus revealing a novel function of post-translational modifications. Furthermore, the oxidation-sensitive proteins identified here provide a valuable resource for research on protein redox regulation and can guide future mechanistic studies.

Improved DNA-Versus-Protein Homology Search for Protein Fossils

Protein fossils, i.e., noncoding DNA descended from coding DNA, arise frequently from transposable elements (TEs), decayed genes, and viral integrations. They can reveal, and mislead about, evolutionary history and relationships. They have been detected by comparing DNA to protein sequences, but current methods are not optimized for this task. We describe a powerful DNA-protein homology search method. We use a 64×21 substitution matrix, which is fitted to sequence data, automatically learning the genetic code. We detect subtly homologous regions by considering alternative possible alignments between them, and calculate significance (probability of occurring by chance between random sequences). Our method detects TE protein fossils much more sensitively than blastx, and faster. Of the  ∼ 7 major categories of eukaryotic TE, three were long thought absent in mammals: we find two of them in the human genome, polinton and DIRS/Ngaro. This method increases our power to find ancient fossils, and perhaps to detect non-standard genetic codes. The alternative-alignments and significance paradigm is not specific to DNA-protein comparison, and could benefit homology search generally. This is an extended version of a conference paper (Yao & Frith, 2021).

Detection of Burkholderia pseudomallei with CRISPR-Cas12a based on specific sequence tags

Melioidosis is a bacterial infection caused by Burkholderia pseudomallei (B. pseudomallei), posing a significant threat to public health. Rapid and accurate detection of B. pseudomallei is crucial for preventing and controlling melioidosis. However, identifying B. pseudomallei is challenging due to its high similarity to other species in the same genus. To address this issue, this study proposed a dual-target method that can specifically identify B. pseudomallei in less than 40 min. We analyzed 1722 B. pseudomallei genomes to construct large-scale pan-genomes and selected specific sequence tags in their core genomes that effectively distinguish B. pseudomallei from its closely related species. Specifically, we selected two specific tags, LC1 and LC2, which we combined with the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR associated proteins (Cas12a) system and recombinase polymerase amplification (RPA) pre-amplification. Our analysis showed that the dual-target RPA-CRISPR/Cas12a assay has a sensitivity of approximately 0.2 copies/reaction and 10 fg genomic DNA for LC1, and 2 copies/reaction and 20 fg genomic DNA for LC2. Additionally, our method can accurately and rapidly detect B. pseudomallei in human blood and moist soil samples using the specific sequence tags mentioned above. In conclusion, the dual-target RPA-CRISPR/Cas12a method is a valuable tool for the rapid and accurate identification of B. pseudomallei in clinical and environmental samples, aiding in the prevention and control of melioidosis.

Staphylococcus aureus Infection Induces the Production of the Neutrophil Chemoattractants CXCL1, CXCL2, CXCL3, CXCL5, CCL3, and CCL7 by Murine Osteoblasts

Staphylococcus aureus is the principal causative agent of osteomyelitis, a serious bacterial infection of bone that is associated with progressive inflammatory damage. Bone-forming osteoblasts have increasingly been recognized to play an important role in the initiation and progression of detrimental inflammation at sites of infection and have been demonstrated to release an array of inflammatory mediators and factors that promote osteoclastogenesis and leukocyte recruitment following bacterial challenge. In the present study, we describe elevated bone tissue levels of the potent neutrophil-attracting chemokines CXCL1, CXCL2, CXCL3, CXCL5, CCL3, and CCL7 in a murine model of posttraumatic staphylococcal osteomyelitis. RNA sequencing (RNA-Seq) gene ontology analysis of isolated primary murine osteoblasts showed enrichment in differentially expressed genes involved in cell migration and chemokine receptor binding and chemokine activity following S. aureus infection, and a rapid increase in the expression of mRNA encoding CXCL1, CXCL2, CXCL3, CXCL5, CCL3, and CCL7, in these cells. Importantly, we have confirmed that such upregulated gene expression results in protein production with the demonstration that S. aureus challenge elicits the rapid and robust release of these chemokines by osteoblasts and does so in a bacterial dose-dependent manner. Furthermore, we have confirmed the ability of soluble osteoblast-derived chemokines to elicit the migration of a neutrophil-like cell line. As such, these studies demonstrate the robust production of CXCL1, CXCL2, CXCL3, CXCL5, CCL3, and CCL7 by osteoblasts in response to S. aureus infection, and the release of such neutrophil-attracting chemokines provides an additional mechanism by which osteoblasts could drive the inflammatory bone loss associated with staphylococcal osteomyelitis.

Rhizosphere Priestia species altered cowpea root transcriptome and enhanced growth under drought and nutrient deficiency

Priestia species isolated from the cowpea rhizosphere altered the transcriptome of cowpea roots by colonization and enhanced nutrient uptake, antioxidant mechanisms, and photosynthesis, protecting cowpea from drought and nutrient deficiency. Cowpea is a significant grain legume crop primarily grown in sub-Saharan Africa, Asia, and South America. Drought and nutrient deficiency affect the growth and yield of cowpea. To address this challenge, we studied the phyto-beneficial effects of stress-tolerant rhizobacteria on the biomass yield of cowpea under water- and nutrient-deficit conditions. Among the bacteria isolated, two rhizobacillus genotypes, C8 (Priestia filamentosa; basonym: Bacillus filamentosus) and C29 (Priestia aryabhattai; basonym: Bacillus aryabhattai) were evaluated for the improvement of seed germination and growth of cowpea under stress. Our study revealed that C8 protected cowpea from stress by facilitating phosphorus and potassium uptake, protecting it from oxidative damage, reducing transpiration, and enhancing CO₂ assimilation. A 17% increase in root biomass upon C8 inoculation was concomitant with the induction of stress tolerance genes in cowpea roots predominantly involved in growth and metabolic processes, cell wall organization, ion homeostasis, and cellular responses to phosphate starvation. Our results indicate a metabolic alteration in cowpea root triggered by P. filamentosa, leading to efficient nutrient reallocation in the host plant. We propose inoculation with P. filamentosa as an effective strategy for improving the yield of cowpea in low-input agriculture, where chemical fertilization and irrigation are less accessible to resource-poor farmers.

Induction of protective interferon-β responses in murine osteoblasts following Staphylococcus aureus infection

Introduction

The refractory and recurrent nature of chronic staphylococcal osteomyelitis may be due, at least in part, to the ability of Staphylococcus aureus to invade and persist within bone-forming osteoblasts. However, osteoblasts are now recognized to respond to S. aureus infection and produce numerous immune mediators and bone regulatory factors that can shape the host response. Type I interferons (IFNs) are best known for their antiviral effects, but it is becoming apparent that they impact host susceptibility to a wide range of pathogens including S. aureus.

Methods

Here, we have assessed the local expression of IFN-β by specific capture ELISA in an established in vivo mouse model of staphylococcal osteomyelitis. RNA Tag-Seq analysis, specific capture ELISAs, and/or immunoblot analyses, were then used to assess the expression of type I IFNs and select IFN stimulated genes (ISGs) in S. aureus infected primary murine osteoblasts. The effect of IFN-β on intracellular S. aureus burden was assessed in vitro following recombinant cytokine treatment by serial colony counts of liberated bacteria.

Results

We report the presence of markedly elevated IFN-β levels in infected bone tissue in a mouse model of staphylococcal osteomyelitis. RNA Tag-Seq analysis of S. aureus infected osteoblasts showed enrichment of genes associated with type I IFN signaling and ISGs, and elevated expression of mRNA encoding IFN-β and ISG products. IFN-β production was confirmed with the demonstration that S. aureus induces its rapid and robust release by osteoblasts in a dose-dependent manner. Furthermore, we showed increased protein expression of the ISG products IFIT1 and IFIT3 by infected osteoblasts and demonstrate that this occurs secondary to the release of IFN-β by these cells. Finally, we have determined that exposure of S. aureus-infected osteoblasts to IFN-β markedly reduces the number of viable bacteria harbored by these cells.

Discussion

Together, these findings indicate an ability of osteoblasts to respond to bacteria by producing IFN-β that can act in an autocrine and/or paracrine manner to elicit ISG expression and mitigate S. aureus infection.

Building the Chordata Olfactory Receptor Database using more than 400,000 receptors annotated by Genome2OR

Olfactory receptors are poorly annotated for most genome-sequenced chordates. To address this deficiency, we developed a nhmmer-based olfactory receptor annotation tool Genome2OR ( https://github.com/ToHanwei/Genome2OR.git ), and used it to process 1,695 sequenced chordate genomes in the NCBI Assembly database as of January, 2021. In total, 765,248 olfactory receptor genes were annotated, with 404,426 functional genes and 360,822 pseudogenes, which represents a four-fold increase in the number of annotated olfactory receptors. Based on the annotation data, we built a database called Chordata Olfactory Receptor Database (CORD, https://cord.ihuman.shanghaitech.edu.cn ) for archiving, analysing and disseminating the data. Beyond the primary data, we offer derivative information, including pictures of species, cross references to public databases, structural models, sequence similarity networks and sequence profiles in the CORD. Furthermore, we did brief analyses on these receptors, including building a huge protein sequence similarity network covering all receptors in the database, and clustering them into 20 communities, classifying the 20 communities into three categories based on their presences/absences in ray-finned fish and/or lobe-finned fish. We infer that olfactory receptors should have unique activation and desensitization mechanisms by analysing their sequences and structural models. We believe the CORD can benefit the researchers and the general public who are interested in olfaction.

Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes

The integration of mitochondrial genome fragments into the nuclear genome is well documented, and the transfer of these mitochondrial nuclear pseudogenes (numts) is thought to be an ongoing evolutionary process. With the increasing number of eukaryotic genomes available, genome-wide distributions of numts are often surveyed. However, inconsistencies in genome quality can reduce the accuracy of numt estimates, and methods used for identification can be complicated by the diverse sizes and ages of numts. Numts have been previously characterized in rodent genomes and it was postulated that they might be more prevalent in a group of voles with rapidly evolving karyotypes. Here, we examine 37 rodent genomes, and an additional 26 vertebrate genomes, while also considering numt detection methods. We identify numts using DNA:DNA and protein:translated-DNA similarity searches and compare numt distributions among rodent and vertebrate taxa to assess whether some groups are more susceptible to transfer. A combination of protein sequence comparisons (protein:translated-DNA) and BLASTN genomic DNA searches detect 50% more numts than genomic DNA:DNA searches alone. In addition, higher-quality RefSeq genomes produce lower estimates of numts than GenBank genomes, suggesting that lower quality genome assemblies can overestimate numts abundance. Phylogenetic analysis shows that mitochondrial transfers are not associated with karyotypic diversity among rodents. Surprisingly, we did not find a strong correlation between numt counts and genome size. Estimates using DNA: DNA analyses can underestimate the amount of mitochondrial DNA that is transferred to the nucleus.

Uncovering the characteristics of the gut microbiota in patients with acute ischemic stroke and phlegm-heat syndrome

Growing evidence has indicated that the characteristics of gut microbiota are associated with acute ischemic stroke (AIS). Phlegm-heat syndrome (PHS), a specific pathological state of the AIS, is one of the common traditional Chinese syndromes of stroke. The long duration of PHS in patients with AIS could lead to poor clinical outcomes. Gut microbiota characteristics in patients with both AIS and PHS, and their relationship remains unknown. This study was designed to investigate the alterations in gut microbiota in patients with AIS and PHS through a cross-sectional study. Fecal samples were collected from 10 patients with AIS and non-PHS (ntAIS), 7 patients with AIS and PHS (tAIS), and 10 healthy controls (HC). Samples were profiled via Illumina sequencing of the 16S rRNA V3-V4. Stroke severity was assessed at admission by the National Institutes of Health Stroke Scale (NIHSS) and modified Rankin scale (mRS); their correlation with gut microbiota was investigated. The alpha-diversity of the bacterial communities was significantly higher in the fecal samples of patients with tAIS than in patients with ntAIS (Shannon index, P = 0.037). In addition, the combined tAIS and ntAIS group (tntAIS) exhibited higher microbiotic diversity when compared with HC (chao1, P = 0.019). The structure of intestinal microbiota was effectively distinguished between the tAIS and ntAIS group (ANOSIM, r = 0.337, P = 0.007). Additionally, the gut microbiota structure was significantly different between the tntAIS and HC groups (ANOSIM, r = 0.217, P = 0.005). The genera, Ruminococcaceae_ UCG_002 and Christensenellaceae_R-7_group, were implicated in the discrimination of PHS from non-PHS. The order Lactobacillales and family Lachnospiraceae were significantly negatively correlated with NIHSS and mRS at admission (P < 0.05). By contrast, the order Desulfovibrionales, families Christensenellaceae and Desulfovibrionaceae, and genera Ruminococcaceae UCG-014 and Ruminococcaceae UCG-002 were significantly positively correlated with NIHSS and mRS at admission (P < 0.05). This study is the first to profile the characteristics of gut microbiota in patients with AIS and PHS, compared with those with non-PHS. The genera, Ruminococcaceae_ UCG_002 and Christensenellaceae_R-7_group, may be objective indicators of this traditional Chinese medicine (TCM) syndrome in AIS. Furthermore, it provides a microbe-inspired biological basis for TCM syndrome differentiation.

Characteristics of Wild Bird Resistomes and Dissemination of Antibiotic Resistance Genes in Interconnected Bird-Habitat Systems Revealed by Similarity of blaTEM Polymorphic Sequences

Wild birds are known to harbor and discharge antibiotic-resistant bacteria (ARB) and their associated antibiotic resistance genes (ARGs). However, assessments of their contribution to the dissemination of antibiotic resistance in the environment are limited to culture-dependent bacterial snapshots. Here, we present a high-throughput sequencing study that corroborates extensive ARG exchange between wild bird feces and their habitats and implies the need to scrutinize high-mobility birds as potential vectors for global propagation of ARGs. We characterized the resistome (281 ARGs) and microbiome of seven wild bird species and their terrestrial and aquatic habitats. The resistomes of bird feces were influenced by the microbial community structure, mobile genetic elements (MGEs), and residual antibiotics. We designated 33 ARGs found in more than 90% of the bird fecal samples as core ARGs of wild bird feces, among which 16 ARGs were shared as core ARGs in both wild bird feces and their habitats; these genes represent a large proportion of both the bird feces (35.0 ± 15.9%) and the environmental resistome (29.9 ± 21.4%). One of the most detected β-lactam resistance genes (bla_TEM, commonly harbored by multidrug resistant "superbugs") was used as molecular marker to demonstrate the high interconnectivity of ARGs between the microbiomes of wild birds and their habitats. Overall, this work provides a comprehensive analysis of the wild bird resistome and underscores the importance to consider genetic exchange between animals and the environment in the One Health approach.

COLONOMICS - integrative omics data of one hundred paired normal-tumoral samples from colon cancer patients

Colonomics is a multi-omics dataset that includes 250 samples: 50 samples from healthy colon mucosa donors and 100 paired samples from colon cancer patients (tumor/adjacent). From these samples, Colonomics project includes data from genotyping, DNA methylation, gene expression, whole exome sequencing and micro-RNAs (miRNAs) expression. It also includes data from copy number variation (CNV) from tumoral samples. In addition, clinical data from all these samples is available. The aims of the project were to explore and integrate these datasets to describe colon cancer at molecular level and to compare normal and tumoral tissues. Also, to improve screening by finding biomarkers for the diagnosis and prognosis of colon cancer. This project has its own website including four browsers allowing users to explore Colonomics datasets. Since generated data could be reuse for the scientific community for exploratory or validation purposes, here we describe omics datasets included in the Colonomics project as well as results from multi-omics layers integration.

LI-Detector: a Method for Curating Ordered Gene-Replacement Libraries

In recent years the availability of genome sequence information has grown logarithmically resulting in the identification of a plethora of uncharacterized genes. To address this gap in functional annotation, many high-throughput screens have been devised to uncover novel gene functions. Gene-replacement libraries are one such tool that can be screened in a high-throughput way to link genotype and phenotype and are key community resources. However, for a phenotype to be attributed to a specific gene, there needs to be confidence in the genotype. Construction of large libraries can be laborious and occasionally errors will arise. Here, we present a rapid and accurate method for the validation of any ordered library where a gene has been replaced or disrupted by a uniform linear insertion (LI). We applied our method (LI-detector) to the well-known Keio library of Escherichia coli gene-deletion mutants. Our method identified 3,718 constructed mutants out of a total of 3,728 confirmed isolates, with a success rate of 99.7% for identifying the correct kanamycin cassette position. This data set provides a benchmark for the purity of the Keio mutants and a screening method for mapping the position of any linear insertion, such as an antibiotic resistance cassette in any ordered library. IMPORTANCE The construction of ordered gene replacement libraries requires significant investment of time and resources to create a valuable community resource. During construction, technical errors may result in a limited number of incorrect mutants being made. Such mutants may confound the output of subsequent experiments. Here, using the remarkable E. coli Keio knockout library, we describe a method to rapidly validate the construction of every mutant.

Light attenuation in enriched purple phototrophic bacteria cultures: Implications for modelling and reactor design

Light attenuation in enriched purple phototrophic bacteria (PPB) cultures has not been studied, and its understanding is critical for proper process modelling and reactor design, especially for scaled systems. This work evaluated the effect of different biomass concentrations, reactor configurations, wastewater matrices, and growth conditions, on the attenuation extent of near infra-red (NIR) and ultraviolet-visible (UV-VIS) light spectra. The results show that increased biomass concentrations lead to higher light attenuation, and that PPB absorb both VIS and NIR wavelengths, with both fractions of the spectrum being equally absorbed at biomass concentrations above 1,000 g COD·m^-3. A flat plate configuration showed less attenuation compared with cylindrical reactors illuminated from the top, representative for open ponds. Neither a complex wastewater matrix nor the presence of polyhydroxyalkanoates (under nutrient limited conditions) affected light attenuation significantly. The pigment concentration (both bacteriochlorophyll and carotenoids) however, had a strong effect, with significant attenuation in the presence of pigments. Attenuation predictions using the Lambert-Beer law (excluding scattering) and the Schuster model (including scattering) indicated that light scattering had a minimal effect. A proposed mathematical model, based on the Lambert-Beer law and a Monod function for light requirements, allowed effective prediction of the kinetics of photoheterotrophic growth. This resulted in a half saturation coefficient of 4.6 W·m^-2. Finally, the results showed that in dense outdoor PPB cultures (≥1,000 g COD·m^-3), effective light penetration is only 5 cm, which biases design away from horizontal lagoons, and towards non-incident multi-panel systems such as flat plate reactors.

Composition and Ecological Roles of the Core Microbiome along the Abyssal-Hadal Transition Zone Sediments of the Mariana Trench

The unique geological features of hadal trenches are known to influence both the structure and ecological function of microbial communities. It is also well known that heterotrophs and chemoautotrophs dominate the hadal and abyssal pelagic zones, respectively. Here, a metagenomic investigation was conducted on sediment samples obtained from the abyssal-hadal transition zone in the Mariana Trench to gain a better understanding of the general diversity and potential function of the core microbiome in this zone. A high level of cosmopolitanism existed in the core microbiome referred from a high community similarity among different stations. Niche differentiation along the fine-scale of different sediment layers was observed, especially for major archaeal groups, largely due to sediment depth and the source of organic matter. A prevalence of nitrogen biogeochemical cycles driven by various nitrifying groups with the capability of dark carbon fixation in the abyssal-hadal biosphere was also demonstrated. The predominance of heterotrophic over chemolithoautotrophic pathways in this transition zone was found, and a high abundance of genes related to respiration and carbon fixation (i.e., the intact Calvin and rTCA cycles) were detected as well, which might reflect the intensive microbial activities known to occur in this deep biosphere. The presence of those metabolic processes and associated microbes were reflected by functional and genetic markers generated from the metagenomic data in the current study. However, their roles and contributions to the nitrogen/carbon biogeochemical cycles and flux in the abyssal-hadal transition zone still need further analysis. IMPORTANCE The Mariana Trench is the deepest oceanic region on earth, its microbial ecological exploration has become feasible with the rapid progress of submersible and metagenomic sequencing. We investigated the community compositions and metabolic functions of the core microbiome along the abyssal-hadal transition zone of the Mariana Trench, although most studies by far were focused on the pelagic zone. We found a predominance of heterotrophic groups and related metabolic pathways, which were closely associated with nitrogen biogeochemical cycles driven by various nitrifying groups with the capability of dark carbon fixation.

Outdoor demonstration-scale flat plate photobioreactor for resource recovery with purple phototrophic bacteria

To make purple phototrophic bacteria (PPB)-based technologies a reality for resource recovery, research must be demonstrated outdoors, using scaled reactors. In this study, a 10 m long PPB-enriched flat plate photobioreactor (FPPBR) with a volume of 0.95 m³ was operated for 253 days, fed with poultry processing wastewater. Different operational strategies were tested, including varying influent types, retention times, feeding strategies, and anaerobic/aerobic conditions in a novel mixed metabolic mode concept. The overall results show that regardless of the fermented wastewater fed (raw or after solid removal via dissolved air flotation) and the varying environmental conditions (e.g., light exposure and temperatures), the FPPBR provided effective volatile fatty acids (VFAs), N, and P removals (average efficiencies of >90%, 34-77%, and 28-45%, respectively). The removal of N and P was limited by the availability of biodegradable COD. Biomass (C, N and P) could be harvested at ∼90% VS/TS ratio, 58% crude protein content and a suitable amino acid profile for potential feed applications. During fully anaerobic operation with semicontinuous/day-only feeding, the FPPBR showed biomass productivities between 25 and 84 g VS m^-2 d^-1 (high due to solid influx; the productivities estimated from COD removal rates were 6.0-24 g VS•m^-2•d^-1 (conservative values)), and soluble COD removal rates of up to 1.0 g•L^-1•d^-1 (overall average of 0.34 ± 0.16 g•L^-1•d^-1). Under these conditions, the relative abundance of PPB in the harvested biomass was up to 56%. A minimum overall HRT of 2-2.4 d (1.0-1.2 d when only fed during the day) is recommended to avoid PPB washout, assuming no biomass retention. A combined daily-illuminated-anaerobic/night-aerobic operation (supplying air during night-time) exploiting photoheterotrophy during the day and aerobic chemoheterotrophy of the same bacteria at night improved the overall removal performance, avoiding VFA accumulation during the night. However, while enabling enhanced treatment, this resulted in a lower relative abundance of PPB and reduced biomass productivities, highlighting the need to balance resource recovery and treatment goals.

Broad-spectrum fungal resistance in sorghum is conferred through the complex regulation of an immune receptor gene embedded in a natural antisense transcript

Sorghum (Sorghum bicolor), the fifth most widely grown cereal crop globally, provides food security for millions of people. Anthracnose caused by the fungus Colletotrichum sublineola is a major disease of sorghum worldwide. We discovered a major fungal resistance locus in sorghum composed of the nucleotide-binding leucine-rich repeat receptor gene ANTHRACNOSE RESISTANCE GENE1 (ARG1) that is completely nested in an intron of a cis-natural antisense transcript (NAT) gene designated CARRIER OF ARG1 (CARG). Susceptible genotypes express CARG and two alternatively spliced ARG1 transcripts encoding truncated proteins lacking the leucine-rich repeat domains. In resistant genotypes, elevated expression of an intact allele of ARG1, attributed to the loss of CARG transcription and the presence of miniature inverted-repeat transposable element sequences, resulted in broad-spectrum resistance to fungal pathogens with distinct virulence strategies. Increased ARG1 expression in resistant genotypes is also associated with higher histone H3K4 and H3K36 methylation. In susceptible genotypes, lower ARG1 expression is associated with reduced H3K4 and H3K36 methylation and increased expression of NATs of CARG. The repressive chromatin state associated with H3K9me2 is low in CARG-expressing genotypes within the CARG exon and higher in genotypes with low CARG expression. Thus, ARG1 is regulated by multiple mechanisms and confers broad-spectrum, strong resistance to fungal pathogens.

Insights into carbon-fixation pathways through metagonomics in the sediments of deep-sea cold seeps

Carbon fixation by chemoautotrophic microorganisms in the dark ocean has a major impact on global carbon cycling and ecological relationships in the ocean's interior. At present, six pathways of autotrophic carbon fixation have been found: the Calvin cycle, the reductive Acetyl-CoA or Wood-Ljungdahl pathway (rAcCoA), the reductive tricarboxylic acid cycle (rTCA), the 3-hydroxypropionate bicycle (3HP), the 3-hydroxypropionate/4-hydroxybutyrate cycle (3HP/4HB), and the dicarboxylate/4-hydroxybutyrate cycle (DC/4HB). Although our knowledge about carbon fixation pathways in the ocean has increased significantly, carbon fixation pathways in the cold seeps are still unknown. In this study, we collected sediment samples from two cold seeps and one trough in the south China sea (SCS), and investigated with metagenomic and metagenome assembled genomes (MAGs). We found that six autotrophic carbon fixation pathways present in the cold seeps and trough with rTCA cycle was the most common pathway, whose genes were particularly high in the cold seeps and increased with sediment depths; the rAcCoA cycle mainly occurred in the cold seep regions, and the abundance of module genes increased with sediment depths. We also elucidated members of chemoautotrophic microorganisms involved in these six carbon-fixation pathways. The rAcCoA, rTCA and DC/4-HB cycles required significantly less energy probably play an important role in the deep-sea environments, especially in the cold seeps. This study provided metabolic insights into the carbon fixation pathways in the cold seeps, and laid the foundation for future detailed study on processes and rates of carbon fixation in the deep-sea ecosystems.

Insulin-like peptide 3 stimulates hemocytes to proliferate in anautogenous and facultatively autogenous mosquitoes

Most mosquito species are anautogenous, which means they must blood feed on a vertebrate host to produce eggs, while a few are autogenous and can produce eggs without blood feeding. Egg formation is best understood in the anautogenous mosquito Aedes aegypti, where insulin-like peptides (ILPs), ovary ecdysteroidogenic hormone (OEH) and 20-hydroxyecdysone (20E) interact to regulate gonadotrophic cycles. Circulating hemocytes also approximately double in abundance in conjunction with a gonadotrophic cycle, but the factors responsible for stimulating this increase remain unclear. Focusing on Ae. aegypti, we determined that hemocyte abundance similarly increased in intact blood-fed females and decapitated blood-fed females that were injected with ILP3, whereas OEH, 20E or heat-killed bacteria had no stimulatory activity. ILP3 upregulated insulin-insulin growth factor signaling in hemocytes, but few genes - including almost no transcripts for immune factors - were differentially expressed. ILP3 also stimulated circulating hemocytes to increase in two other anautogenous (Anopheles gambiae and Culex quinquefasciatus) and two facultatively autogenous mosquitoes (Aedes atropalpus and Culex pipiens molestus), but had no stimulatory activity in the obligately autogenous mosquito Toxorhynchites amboinensis. Altogether, our results identify ILPs as the primary regulators of hemocyte proliferation in association with egg formation, but also suggest this response has been lost in the evolution of obligate autogeny.

Loss of YhcB results in dysregulation of coordinated peptidoglycan, LPS and phospholipid synthesis during Escherichia coli cell growth

The cell envelope is essential for viability in all domains of life. It retains enzymes and substrates within a confined space while providing a protective barrier to the external environment. Destabilising the envelope of bacterial pathogens is a common strategy employed by antimicrobial treatment. However, even in one of the best studied organisms, Escherichia coli, there remain gaps in our understanding of how the synthesis of the successive layers of the cell envelope are coordinated during growth and cell division. Here, we used a whole-genome phenotypic screen to identify mutants with a defective cell envelope. We report that loss of yhcB, a conserved gene of unknown function, results in loss of envelope stability, increased cell permeability and dysregulated control of cell size. Using whole genome transposon mutagenesis strategies, we report the comprehensive genetic interaction network of yhcB, revealing all genes with a synthetic negative and a synthetic positive relationship. These genes include those previously reported to have a role in cell envelope biogenesis. Surprisingly, we identified genes previously annotated as essential that became non-essential in a ΔyhcB background. Subsequent analyses suggest that YhcB functions at the junction of several envelope biosynthetic pathways coordinating the spatiotemporal growth of the cell, highlighting YhcB as an as yet unexplored antimicrobial target.

All life depends on a cell envelope to enclose the chemical reactions that make life possible. But how do cell envelopes grow? How each component of the cell envelope is incorporated into the envelope at the correct amount, in the correct place, and at the correct time, to prevent cell death, has been a long-standing question in bacteriology. Using a unique combination of high throughput chemical genetic screens we identified yhcB, a conserved gene of unknown function, required for the maintenance of cell envelope integrity in Escherichia coli. Loss of YhcB results in aberrant cell size driven by the production of excess membrane phospholipids. Subsequent molecular and biochemical analyses suggest YhcB influences the spatiotemporal biogenesis of LPS, peptidoglycan and membrane phospholipids. Our data indicate YhcB is a key regulator of cell envelope growth in Gram-negative bacteria playing a crucial role in coordinating cell width, elongation, and division to maintain cell envelope integrity.

Tissue-specific gene expression shows a cynipid wasp repurposes oak host gene networks to create a complex and novel parasite-specific organ

Every organism on Earth depends on interactions with other organisms to survive. In each of these interactions, an organism must utilize the limited toolbox of genes and proteins it possesses to successfully manipulate or cooperate with another species, but it can also co-opt the genome machinery of its partner to expand its available tools. Insect-induced plant galls are an extreme example of this, wherein an insect hijacks the plant's genome to direct the initiation and development of galls consisting of plant tissue. However, previous transcriptomic studies have not evaluated individual tissues within a gall to determine the full extent to which a galling insect manipulates its host plant. Here we demonstrate that the cynipid wasp Dryocosmus quercuspalustris creates a complex parasite-specific organ from red oak tissue via massive changes in host gene expression. Our results show that the gall wasp is not merely modifying oak leaf tissue but creating extensive changes in gene expression between galled and ungalled tissue (differential expression in 28% of genes) and distinct gall tissue types (20% of genes). The outer gall tissue shows increases in various plant defence systems, which is consistent with its predicted functional role of protecting the wasp larva. The inner larval capsule shows suppression of large parts of the plant innate immune system and evidence for the wasp utilizing the plant's RNA interference mechanisms, which may be a potential mechanism for the wasp's control on gall growth.

DpaA Detaches Braun's Lipoprotein from Peptidoglycan

Gram-negative bacteria have a unique cell envelope with a lipopolysaccharide-containing outer membrane that is tightly connected to a thin layer of peptidoglycan. The tight connection between the outer membrane and peptidoglycan is needed to maintain the outer membrane as an impermeable barrier for many toxic molecules and antibiotics. Enterobacteriaceae such as Escherichia coli covalently attach the abundant outer membrane-anchored lipoprotein Lpp (Braun’s lipoprotein) to tripeptides in peptidoglycan, mediated by the transpeptidases LdtA, LdtB, and LdtC. LdtD and LdtE are members of the same family of ld-transpeptidases but they catalyze a different reaction, the formation of 3-3 cross-links in the peptidoglycan. The function of the sixth homologue in E. coli, LdtF, remains unclear, although it has been shown to become essential in cells with inhibited lipopolysaccharide export to the outer membrane. We now show that LdtF hydrolyzes the Lpp-peptidoglycan linkage, detaching Lpp from peptidoglycan, and have renamed LdtF to peptidoglycan meso-diaminopimelic acid protein amidase A (DpaA). We show that the detachment of Lpp from peptidoglycan is beneficial for the cell under certain stress conditions and that the deletion of dpaA allows frequent transposon inactivation in the lapB (yciM) gene, whose product downregulates lipopolysaccharide biosynthesis. DpaA-like proteins have characteristic sequence motifs and are present in many Gram-negative bacteria, of which some have no Lpp, raising the possibility that DpaA has other substrates in these species. Overall, our data show that the Lpp-peptidoglycan linkage in E. coli is more dynamic than previously appreciated.

Purple phototrophic bacteria are outcompeted by aerobic heterotrophs in the presence of oxygen

There is an ongoing debate around the effect of microaerobic/aerobic conditions on the wastewater treatment performance and stability of enriched purple phototrophic bacteria (PPB) cultures. It is well known that oxygen-induced oxidative conditions inhibit the synthesis of light harvesting complexes, required for photoheterotrophy. However, in applied research, several publications have reported efficient wastewater treatment at high dissolved oxygen (DO) levels. This study evaluated the impact of different DO concentrations (0-0.25 mg·L^-1, 0-0.5 mg·L^-1 and 0-4.5 mg·L^-1) on the COD, nitrogen and phosphorus removal performances, the biomass yields, and the final microbial communities of PPB-enriched cultures, treating real wastewaters (domestic and poultry processing wastewater). The results show that the presence of oxygen suppressed photoheterotrophic growth, which led to a complete pigment and colour loss in a matter of 20-30 h after starting the batch. Under aerobic conditions, chemoheterotrophy was the dominant catabolic pathway, with wastewater treatment performances similar to those achieved in common aerobic reactors, rather than those corresponding to phototrophic systems (i.e. considerable total COD decrease (45-57% aerobically vs. ± 10% anaerobically). This includes faster consumption of COD and nutrients, lower nutrient removal efficiencies (50-58% vs. 72-99% for NH₄⁺-N), lower COD:N:P substrate ratios (100:4.5-5.0:0.4-0.8 vs. 100:6.7-12:0.9-1.2), and lower apparent biomass yields (0.15-0.31 vs. 0.8-1.2 g COD_biomass·g COD_removed^-1)). The suppression of photoheterotrophy inevitably resulted in a reduction of the relative PPB abundances in all the aerated tests (below 20% at the end of the tests), as PPB lost their main competitive advantage against competing aerobic heterotrophic microbes. This was explained by the lower aerobic PPB growth rates (2.4 d^-1 at 35 °C) when compared to common growth rates for aerobic heterotrophs (6.0 d^-1 at 20 °C). Therefore, PPB effectively outcompete other microbes under illuminated-anaerobic conditions, but not under aerobic or even micro-aerobic conditions, as shown by continuously aerated tests controlled at undetectable DO levels. While their aerobic heterotrophic capabilities provide some resilience, at non-sterile conditions PPB cannot dominate when growing chemoheterotrophically, and will be outcompeted.

Mussismilia braziliensis White Plague Disease Is Characterized by an Affected Coral Immune System and Dysbiosis

Infectious diseases are one of the major drivers of coral reef decline worldwide. White plague-like disease (WPL) is a widespread disease with a complex etiology that infects several coral species, including the Brazilian endemic species Mussismilia braziliensis. Gene expression profiles of healthy and WPL-affected M. braziliensis were analyzed in winter and summer seasons. The de novo assembly of the M. braziliensis transcriptome from healthy and white plague samples produced a reference transcriptome containing 119,088 transcripts. WPL-diseased samples were characterized by repression of immune system and cellular defense processes. Autophagy and cellular adhesion transcripts were also repressed in WPL samples, suggesting exhaustion of the coral host defenses. Seasonal variation leads to plasticity in transcription with upregulation of intracellular signal transduction, apoptosis regulation, and oocyte development in the summer. Analysis of the active bacterial rRNA indicated that Pantoea bacteria were more abundant in WPL corals, while Tistlia, Fulvivirga, and Gammaproteobacteria Ga0077536 were more abundant in healthy samples. Cyanobacteria proliferation was also observed in WPL, mostly in the winter. These results indicate a scenario of dysbiosis in WPL-affected M. braziliensis, with the loss of potentially symbiotic bacteria and proliferation of opportunistic microbes after the start of the infection process.

A first insight into the genome of Prototheca wickerhamii, a major causative agent of human protothecosis

BACKGROUND

Colourless microalgae of the Prototheca genus are the only known plants that have consistently been implicated in a range of clinically relevant opportunistic infections in both animals and humans. The Prototheca algae are emerging pathogens, whose incidence has increased importantly over the past two decades. Prototheca wickerhamii is a major human pathogen, responsible for at least 115 cases worldwide. Although the algae are receiving more attention nowadays, there is still a substantial knowledge gap regarding their biology, and pathogenicity in particular. Here we report, for the first time, the complete nuclear genome, organelle genomes, and transcriptome of the P. wickerhamii type strain ATCC 16529.

RESULTS

The assembled genome size was of 16.7 Mbp, making it the smallest and most compact genome sequenced so far among the protothecans. Key features of the genome included a high overall GC content (64.5%), a high number (6081) and proportion (45.9%) of protein-coding genes, and a low repetitive sequence content (2.2%). The vast majority (90.6%) of the predicted genes were confirmed with the corresponding transcripts upon RNA-sequencing analysis. Most (93.2%) of the genes had their putative function assigned when searched against the InterProScan database. A fourth (23.3%) of the genes were annotated with an enzymatic activity possibly associated with the adaptation to the human host environment. The P. wickerhamii genome encoded a wide array of possible virulence factors, including those already identified in two model opportunistic fungal pathogens, i.e. Candida albicans and Trichophyton rubrum, and thought to be involved in invasion of the host or elicitation of the adaptive stress response. Approximately 6% of the P. wickerhamii genes matched a Pathogen-Host Interaction Database entry and had a previously experimentally proven role in the disease development. Furthermore, genes coding for proteins (e.g. ATPase, malate dehydrogenase) hitherto considered as potential virulence factors of Prototheca spp. were demonstrated in the P. wickerhamii genome.

CONCLUSIONS

Overall, this study is the first to describe the genetic make-up of P. wickerhamii and discovers proteins possibly involved in the development of protothecosis.

Characterization of three bacterial glycoside hydrolase family 9 endoglucanases with different modular architectures isolated from a compost metagenome

BACKGROUND

Environmental bacteria express a wide diversity of glycoside hydrolases (GH). Screening and characterization of GH from metagenomic sources provides an insight into biomass degradation strategies of non-cultivated prokaryotes.

METHODS

In the present report, we screened a compost metagenome for lignocellulolytic activities and identified six genes encoding enzymes belonging to family GH9 (GH9a-f). Three of these enzymes (GH9b, GH9d and GH9e) were successfully expressed and characterized.

RESULTS

A phylogenetic analysis of the catalytic domain of pro- and eukaryotic GH9 enzymes suggested the existence of two major subgroups. Bacterial GH9s displayed a wide variety of modular architectures and those harboring an N-terminal Ig-like domain, such as GH9b and GH9d, segregated from the remainder. We purified and characterized GH9 endoglucanases from both subgroups and examined their stabilities, substrate specificities and product profiles. GH9e exhibited an original hydrolysis pattern, liberating an elevated proportion of oligosaccharides longer than cellobiose. All of the enzymes exhibited processive behavior and a synergistic action on crystalline cellulose. Synergy was also evidenced between GH9d and a GH48 enzyme identified from the same metagenome.

CONCLUSIONS

The characterized GH9 enzymes displayed different modular architectures and distinct substrate and product profiles. The presence of a cellulose binding domain was shown to be necessary for binding and digestion of insoluble cellulosic substrates, but not for processivity.

GENERAL SIGNIFICANCE

The identification of six GH9 enzymes from a compost metagenome and the functional variety of three characterized members highlight the importance of this enzyme family in bacterial biomass deconstruction.

Nutrient-imbalanced conditions shift the interplay between zooplankton and gut microbiota

Background

Nutrient stoichiometry of phytoplankton frequently changes with aquatic ambient nutrient concentrations, which is mainly influenced by anthropogenic water treatment and the ecosystem dynamics. Consequently, the stoichiometry of phytoplankton can markedly alter the metabolism and growth of zooplankton. However, the effects of nutrient-imbalanced prey on the interplay between zooplankton and their gut microbiota remain unknown. Using metatranscriptome, a 16 s rRNA amplicon-based neutral community model (NCM) and experimental validation, we investigated the interactions between Daphnia magna and its gut microbiota in a nutrient-imbalanced algal diet.

Results

Our results showed that in nutrient-depleted water, the nutrient-enriched zooplankton gut stimulated the accumulation of microbial polyphosphate in fecal pellets under phosphorus limitation and the microbial assimilation of ammonia under nitrogen limitation. Compared with the nutrient replete group, both N and P limitation markedly promoted the gene expression of the gut microbiome for organic matter degradation but repressed that for anaerobic metabolisms. In the nutrient limited diet, the gut microbial community exhibited a higher fit to NCM (R² = 0.624 and 0.781, for N- and P-limitation, respectively) when compared with the Control group (R² = 0.542), suggesting increased ambient-gut exchange process favored by compensatory feeding. Further, an additional axenic grazing experiment revealed that the growth of D. magna can still benefit from gut microbiota under a nutrient-imbalanced diet.

Conclusions

Together, these results demonstrated that under a nutrient-imbalanced diet, the microbes not only benefit themselves by absorbing excess nutrients inside the zooplankton gut but also help zooplankton to survive during nutrient limitation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07333-z.

Anaerobic methane oxidation coupled to denitrification is an important potential methane sink in deep-sea cold seeps

Microbes play a crucial role in mediating the methane flux in deep-sea cold seep ecosystems, where only methane-related microbes have been well studied, while the whole microbial community and their ecological functions were still largely unknown. Here, we utilized metagenomic data to investigate how the structure and metabolism of microbial community shift in the reduced sediment habitats along the spatial scales. Microbial communities in cold seeps and troughs formed two distinct clades likely driven by environmental factors, such as total sulfur, total phosphate and NO₃^-, rather than geographical proximity. The predominance of Methanosarcinales reflected a high potential for methane production. In addition to the already well-reported ANME-1/SRB consortia, prevalence of bacterial Methylomirabilis and archaeal Methanoperedens as important performers in the n-damo process with respective of nitrite and nitrate as respective electron acceptor was observed in deep-sea hydrate-bearing regions as well. Aerobic methane oxidization was conducted mainly by type I methanotrophs at Site F (Formosa Ridge), but also via the n-damo process by Methanoperedens and Methylomirabilis in the Haima seep and Xisha Trough, respectively. Based on the high abundance of those denitrifying-dependent methane oxidizers and their related functional genes, we concluded that the previously overlooked n-damo process might be a major methane sink in cold seeps or in gas hydrate-bearing sediments if nitrate is available in the anoxic zones. The signature of isotopic labeling would be essential to confirm the contribution of different anaerobic methane oxidizing pathways in deep-sea cold seep ecosystems.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Dickeya poaceiphila sp. nov., a plant-pathogenic bacterium isolated from sugar cane (Saccharum officinarum)

The genus Dickeya is an important group of plant pathogens that currently comprises 10 recognized species. Although most Dickeya isolates originated from infected cultivated plants, they are also isolated from water. The genomic sequence of the Australian strain NCPPB 569^T clearly established its separation from the previously characterized Dickeya species. The average nucleotide identity and digital DNA-DNA hybridization values obtained by comparing strain NCPPB 569^T with strains of characterized Dickeya species were lower than 87 and 32 %, respectively, supporting the delineation of a new species. The name Dickeya poaceiphila sp. nov. is proposed for this taxon with the type strain NCPPB 569^T (=CFBP 8731^T). Two other strains isolated in Australia, CFBP 1537 and CFBP 2040, also belong to this species. Phenotypic and genomic comparisons enabled the identification of traits distinguishing D. poaceiphila isolates from strains of other Dickeya species.

Sequencing and Analysis of the Complete Organellar Genomes of Prototheca wickerhamii

Of the Prototheca genus, Prototheca wickerhamii has the highest clinical significance in humans. However, neither nuclear nor organellar genomes of this species were sequenced until now. The hitherto determined and analyzed mitochondrial and plastid genomes of the alleged P. wickerhamii species belong in fact to another species, recently named Prototheca xanthoriae. This study provides a first insight into the organellar genomes of a true P. wickerhamii (type strain ATCC 16529). The P. wickerhamii mitochondrion had a 53.8-kb genome, which was considerably larger than that of Prototheca ciferrii (formerly Prototheca zopfii gen. 1) and Prototheca bovis (formerly Prototheca zopfii gen. 2), yet similarly functional, with the differences in size attributable to a higher number of introns and the presence of extra unique putative genes. The 48-kb plastid genome of P. wickerhamii, compared to autotrophic Trebouxiophyceae, was highly reduced due to the elimination of the photosynthesis-related genes. The gene content of the plastid genome of P. wickerhamii was, however, very similar to other colorless Prototheca species. Plastid genome-based phylogeny reinforced the polyphyly of the genus Prototheca, with Helicosporidium and Auxenochlorella branching within clades of Prototheca species. Phylogenetic reconstruction also confirmed the close relationship of P. wickerhamii and P. xanthoriae, which is reflected in the synteny of their organellar genomes. Interestingly, the entire set of atp genes was lost in P. wickerhamii plastid genome while being preserved in P. xanthoriae.

Translating GWAS-identified loci for cardiac rhythm and rate using an in vivo image- and CRISPR/Cas9-based approach

A meta-analysis of genome-wide association studies (GWAS) identified eight loci that are associated with heart rate variability (HRV), but candidate genes in these loci remain uncharacterized. We developed an image- and CRISPR/Cas9-based pipeline to systematically characterize candidate genes for HRV in live zebrafish embryos. Nine zebrafish orthologues of six human candidate genes were targeted simultaneously in eggs from fish that transgenically express GFP on smooth muscle cells (Tg[acta2:GFP]), to visualize the beating heart. An automated analysis of repeated 30 s recordings of beating atria in 381 live, intact zebrafish embryos at 2 and 5 days post-fertilization highlighted genes that influence HRV (hcn4 and si:dkey-65j6.2 [KIAA1755]); heart rate (rgs6 and hcn4); and the risk of sinoatrial pauses and arrests (hcn4). Exposure to 10 or 25 µM ivabradine—an open channel blocker of HCNs—for 24 h resulted in a dose-dependent higher HRV and lower heart rate at 5 days post-fertilization. Hence, our screen confirmed the role of established genes for heart rate and rhythm (RGS6 and HCN4); showed that ivabradine reduces heart rate and increases HRV in zebrafish embryos, as it does in humans; and highlighted a novel gene that plays a role in HRV (KIAA1755).

Ross River Virus Provokes Differentially Expressed MicroRNA and RNA Interference Responses in Aedes aegypti Mosquitoes

Alphaviruses are globally distributed and predominately transmitted by mosquitoes. Aedes species are common vectors for the clinically important alphaviruses-Chikungunya, Sindbis, and Ross River (RRV) viruses-with Aedes aegypti also being a vector for the flaviviruses dengue, Yellow Fever, and Zika viruses. Ae. aegypti was putatively implicated in the large 1979-1980 South Pacific Islands outbreak of RRV-the leading cause of arboviral disease in Australia today. The RNA interference (RNAi) defense response in mosquitoes involves a number of small RNAs, with their kinetics induced by alphaviruses being poorly understood, particularly at the tissue level. We compared the small RNA profiles between RRV-infected and noninfected Ae. aegypti midgut and fat body tissues at 2, 6, and 12 days post-inoculation (dpi). RRV induced an incremental RNAi response, yielding short interfering and P-element-induced-wimpy-testis (PIWI)-interacting RNAs. Fourteen host microRNAs were differentially expressed due to RRV with the majority in the fat body at 2 dpi. The largely congruent pattern of microRNA regulation with previous reports for alphaviruses and divergence from those for flaviviruses suggests a degree of conservation, whereas patterns of microRNA expression unique to this study provide novel insights into the tissuespecific hostvirus attributes of Ae. aegypti responses to this previously unexplored oldworld alphavirus.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

Whole-transcriptome analysis of the toxic effects of zearalenone exposure on ceRNA networks in porcine granulosa cells

Zearalenone (ZEA), an estrogen-like mycotoxin, is commonly detected in animal feeds including improperly stored grains. It has been well demonstrated that ovarian granulosa cells (GCs) perform vital roles during follicular development, however, the competing endogenous RNA (ceRNA) network in GCs after ZEA exposure remains to be well described. Here, for the first time, we adopted whole-transcriptome sequence technology to explore the molecular mechanism of ZEA toxicology on porcine GCs. The results provide evidence that the cell cycle of porcine GCs is arrested in the G2/M phase after exposure to ZEA. Furthermore, bioinformation analysis found that cell cycle arrest related genes were perturbed, including CDK1, CCNB1, CDC25A, and CDC25C, which was consistent with the results of RT-qPCR, immunofluorescence, and Western Blotting. Based on the whole-transcriptome sequence data, by constructing ceRNA networks related to cell cycle arrest, we observed that ZEA exposure arrested cell cycle progression at the G2/M phase in porcine GCs, and non-coding RNAs (ncRNAs) played an important role in this process via regulating the expressions of cell cycle arrest related genes. Taken together, our data here provides strong data to support that the toxicological mechanism regarding the widely distributed toxicant ZEA acts through ceRNA networks in porcine granulosa cells.

Decoupling behavioral and transcriptional responses to color in an eyeless cnidarian

Background

Animals have specific molecular, physiological, and behavioral responses to light that are influenced by wavelength and intensity. Predictable environmental changes – predominantly solar and lunar cycles – drive endogenous daily oscillations by setting internal pacemakers, otherwise known as the circadian clock. Cnidarians have been a focal group to discern the evolution of light responsiveness due to their phylogenetic position as a sister phylum to bilaterians and broad range of light-responsive behaviors and physiology. Marine species that occupy a range of depths will experience different ranges of wavelengths and light intensities, which may result in variable phenotypic responses. Here, we utilize the eyeless sea anemone Nematostella vectensis, an estuarine anemone that typically resides in shallow water habitats, to compare behavioral and molecular responses when exposed to different light conditions.

Results

Quantitative measures of locomotion clearly showed that this species responds to light in the blue and green spectral range with a circadian activity profile, in contrast to a circatidal activity profile in the red spectral range and in constant darkness. Differences in average day/night locomotion was significant in each condition, with overall peak activity during the dark period. Comparative analyses of 96 transcriptomes from individuals sampled every 4 h in each lighting treatment revealed complex differences in gene expression between colors, including in many of the genes likely involved in the cnidarian circadian clock. Transcriptional profiling showed the majority of genes are differentially expressed when comparing mid-day with mid-night, and mostly in red light. Gene expression profiles were largely unique in each color, although animals in blue and green were overall more similar to each other than to red light.

Conclusions

Together, these analyses support the hypothesis that cnidarians are sensitive to red light, and this perception results in a rich transcriptional and divergent behavioral response. Future work determining the specific molecular mechanisms driving the circadian and potential circatidal rhythms measured here would be impactful to connect gene expression variation with behavioral variation in this eyeless species.

New insights into the pelagic microorganisms involved in the methane cycle in the meromictic Lake Pavin through metagenomics

Advances in metagenomics have given rise to the possibility of obtaining genome sequences from uncultured microorganisms, even for those poorly represented in the microbial community, thereby providing an important means to study their ecology and evolution. In this study, metagenomic sequencing was carried out at four sampling depths having different oxygen concentrations or environmental conditions in the water column of Lake Pavin. By analyzing the sequenced reads and matching the contigs to the proxy genomes of the closest cultivated relatives, we evaluated the metabolic potential of the dominant planktonic species involved in the methane cycle. We demonstrated that methane-producing communities were dominated by the genus Methanoregula while methane-consuming communities were dominated by the genus Methylobacter, thus confirming prior observations. Our work allowed the reconstruction of a draft of their core metabolic pathways. Hydrogenotrophs, the genes required for acetate activation in the methanogen genome, were also detected. Regarding methanotrophy, Methylobacter was present in the same areas as the non-methanotrophic, methylotrophic Methylotenera, which could suggest a relationship between these two groups. Furthermore, the presence of a large gene inventory for nitrogen metabolism (nitrate transport, denitrification, nitrite assimilation and nitrogen fixation, for instance) was detected in the Methylobacter genome.

Metabolic response of prokaryotic microbes to sporadic hypoxia in a eutrophic subtropical estuary

Coastal eutrophication and consequent oxygen depletion (hypoxia) occurs worldwide due to increased human activity. The paucity of genomic information of microbes in hypoxia prone coastal waters have hindered our understanding of microorganism related causation and adaption to the environment. Here, using metagenomic approach, we investigated microbial metabolic capability in heavily polluted Pearl River estuary. Our results highlighted the possible roles of microbial metabolic activity in the formation of bottom water hypoxia by revealing enriched organic degradation related microbial genes in the bottom layer beneath surface phytoplankton bloom. Microbial nitrate reduction in hypoxia layer was low, possibly due to the low pH and fluctuating oxygen level. On contrary, high abundance of sulfate-reducing, and antibiotic and metal resistance related genes were detected in bottom and surface layers, respectively, indicating microbial adaptation to oxygen depletion and pollution. Our study provides gene level information on the interactive relations between microbial functions and environmental stress.