FastqCleaner: an interactive Bioconductor application for quality-control, filtering and trimming of FASTQ files

Microbiota of pest insect Nezara viridula mediate detoxification and plant defense repression

The Southern green shield bug, Nezara viridula, is an invasive piercing and sucking pest insect that feeds on crop plants and poses a threat to global food production. Given that insects are known to live in a close relationship with microorganisms, our study provides insights into the community composition and function of the N. viridula-associated microbiota and its effect on host-plant interactions. We discovered that N. viridula hosts both vertically and horizontally transmitted microbiota throughout different developmental stages and their salivary glands harbor a thriving microbial community that is transmitted to the plant while feeding. The N. viridula microbiota was shown to aid its host with the detoxification of a plant metabolite, namely 3-nitropropionic acid, and repression of host plant defenses. Our results demonstrate that the N. viridula-associated microbiota plays an important role in interactions between insects and plants and could therefore be considered a valuable target for the development of sustainable pest control strategies.

Comparative Analysis and Identification of Terpene Synthase Genes in Convallaria keiskei Leaf, Flower and Root Using RNA-Sequencing Profiling

The 'Lilly of the Valley' species, Convallaria, is renowned for its fragrant white flowers and distinctive fresh and green floral scent, attributed to a rich composition of volatile organic compounds (VOCs). However, the molecular mechanisms underlying the biosynthesis of this floral scent remain poorly understood due to a lack of transcriptomic data. In this study, we conducted the first comparative transcriptome analysis of C. keiskei, encompassing the leaf, flower, and root tissues. Our aim was to investigate the terpene synthase (TPS) genes and differential gene expression (DEG) patterns associated with essential oil biosynthesis. Through de novo assembly, we generated a substantial number of unigenes, with the highest count in the root (146,550), followed by the flower (116,434) and the leaf (72,044). Among the identified unigenes, we focused on fifteen putative ckTPS genes, which are involved in the synthesis of mono- and sesquiterpenes, the key aromatic compounds responsible for the essential oil biosynthesis in C. keiskei. The expression of these genes was validated using quantitative PCR analysis. Both DEG and qPCR analyses revealed the presence of ckTPS genes in the flower transcriptome, responsible for the synthesis of various compounds such as geraniol, germacrene, kaurene, linalool, nerolidol, trans-ocimene and valencene. The leaf transcriptome exhibited genes related to the biosynthesis of kaurene and trans-ocimene. In the root, the identified unigenes were associated with synthesizing kaurene, trans-ocimene and valencene. Both analyses indicated that the genes involved in mono- and sesquiterpene biosynthesis are more highly expressed in the flower compared to the leaf and root. This comprehensive study provides valuable resources for future investigations aiming to unravel the essential oil-biosynthesis-related genes in the Convallaria genus.

Transcriptome profile analysis of Indian mustard (Brassica juncea L.) during seed germination reveals the drought stress-induced genes associated with energy, hormone, and phenylpropanoid pathways

Indian mustard (Brassica juncea L. Czern and Coss) is an important oil and vegetable crop frequently affected by seasonal drought stress during seed germination, which retards plant growth and causes yield loss considerably. However, the gene networks regulating responses to drought stress in leafy Indian mustard remain elusive. Here, we elucidated the underlying gene networks and pathways of drought response in leafy Indian mustard using next-generation transcriptomic techniques. Phenotypic analysis showed that the drought-tolerant leafy Indian mustard cv. 'WeiLiang' (WL) had a higher germination rate, antioxidant capacity, and better growth performance than the drought-sensitive cv. 'ShuiDong' (SD). Transcriptome analysis identified differentially expressed genes (DEGs) in both cultivars under drought stress during four germination time points (i.e., 0, 12, 24, and 36 h); most of which were classified as drought-responsive, seed germination, and dormancy-related genes. In the Kyoto Encyclopedia of Genes and Genome (KEGG) analyses, three main pathways (i.e., starch and sucrose metabolism, phenylpropanoid biosynthesis, and plant hormone signal transduction) were unveiled involved in response to drought stress during seed germination. Furthermore, Weighted Gene Co-expression Network Analysis (WGCNA) identified several hub genes (novel.12726, novel.1856, BjuB027900, BjuA003402, BjuA021578, BjuA005565, BjuB006596, novel.12977, and BjuA033308) associated with seed germination and drought stress in leafy Indian mustard. Taken together, these findings deepen our understanding of the gene networks for drought responses during seed germination in leafy Indian mustard and provide potential target genes for the genetic improvement of drought tolerance in this crop.

Endogenous Retrovirus Elements Are Co-Expressed with IFN Stimulation Genes in the JAK-STAT Pathway

Background: Endogenous retrovirus (ERV) elements can act as proximal regulatory elements in promoting interferon (IFN) responses. Previous relevant studies have mainly focused on IFN-stimulated genes (ISGs). However, the role of ERV elements as cis-regulatory motifs in regulating genes of the JAK-STAT pathway remains poorly understood. In our study, we analyzed the changes in ERV elements and genes under both IFN stimulation and blockade of the signaling pathway. Methods: The effects of interferon on cells under normal conditions and knockout of the receptor were compared based on the THP1_IFNAR1_KO and THP1_IFNAR2_mutant cell lines. The correlation between differentially expressed ERVs (DHERVs) and differentially expressed genes (DEGs) as DEHERV-G pairs was explored with construction of gene regulatory networks related to ERV and induced by proinflammatory cytokines. Results: A total of 430 DEHERV loci and 190 DEGs were identified in 842 DEHERV-G pairs that are common to the three groups. More than 87% of DEHERV-G pairs demonstrated a consistent expression pattern. ISGs such as AIM2, IFIT1, IFIT2, IFIT3, STAT1, and IRF were activated via the JAK-STAT pathway in response to interferon stimulation. Thus, STAT1, STAT2, and IRF1 appear to play core roles in regulatory networks and are closely associated with ERVs. Conclusions: The RNA expression of ISGs and ERV elements is correlated, indicating that ERV elements are closely linked to host innate immune responses.

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

LABRADOR-A Computational Workflow for Virus Detection in High-Throughput Sequencing Data

High-throughput sequencing (HTS) allows detection of known and unknown viruses in samples of broad origin. This makes HTS a perfect technology to determine whether or not the biological products, such as vaccines are free from the adventitious agents, which could support or replace extensive testing using various in vitro and in vivo assays. Due to bioinformatics complexities, there is a need for standardized and reliable methods to manage HTS generated data in this field. Thus, we developed LABRADOR—an analysis pipeline for adventitious virus detection. The pipeline consists of several third-party programs and is divided into two major parts: (i) direct reads classification based on the comparison of characteristic profiles between reads and sequences deposited in the database supported with alignment of to the best matching reference sequence and (ii) de novo assembly of contigs and their classification on nucleotide and amino acid levels. To meet the requirements published in guidelines for biologicals’ safety we generated a custom nucleotide database with viral sequences. We tested our pipeline on publicly available HTS datasets and showed that LABRADOR can reliably detect viruses in mixtures of model viruses, vaccines and clinical samples.

Transcriptomic analysis to reveal the differentially expressed miRNA targets and their miRNAs in response to Ralstonia solanacearum in ginger species

BACKGROUND

Bacterial wilt is the most devastating disease in ginger caused by Ralstonia solanacearum. Even though ginger (Zingiber officinale) and mango ginger (Curcuma amada) are from the same family Zingiberaceae, the latter is resistant to R. solanacearum infection. MicroRNAs have been identified in many crops which regulates plant-pathogen interaction, either through silencing genes or by blocking mRNA translation. However, miRNA's vital role and its targets in mango ginger in protecting bacterial wilt is not yet studied extensively. In the present study, using the "psRNATarget" server, we analyzed available ginger (susceptible) and mango ginger (resistant) transcriptome to delineate and compare the microRNAs (miRNA) and their target genes (miRTGs).

RESULTS

A total of 4736 and 4485 differential expressed miRTGs (DEmiRTGs) were identified in ginger and mango ginger, respectively, in response to R. solanacearum. Functional annotation results showed that mango ginger had higher enrichment than ginger in top enriched GO terms. Among the DEmiRTGs, 2105 were common in ginger and mango ginger. However, 2337 miRTGs were expressed only in mango ginger which includes 62 defence related and upregulated miRTGs. We also identified 213 miRTGs upregulated in mango ginger but downregulated in ginger, out of which 23 DEmiRTGS were defence response related. We selected nine miRNA/miRTGs pairs from the data set of common miRTGs of ginger and mango ginger and validated using qPCR.

CONCLUSIONS

Our data covered the expression information of 9221 miRTGs. We identified nine miRNA/miRTGs key candidate pairs in response to R. solanacearum infection in ginger. This is the first report of the integrated analysis of miRTGs and miRNAs in response to R. solanacearum infection among ginger species. This study is expected to deliver several insights in understanding the miRNA regulatory network in ginger and mango ginger response to bacterial wilt.

Transcriptome-Wide Identification and Quantification of Caffeoylquinic Acid Biosynthesis Pathway and Prediction of Its Putative BAHDs Gene Complex in A. spathulifolius

The phenylpropanoid pathway is a major secondary metabolite pathway that helps plants overcome biotic and abiotic stress and produces various byproducts that promote human health. Its byproduct caffeoylquinic acid is a soluble phenolic compound present in many angiosperms. Hydroxycinnamate-CoA shikimate/quinate transferase is a significant enzyme that plays a role in accumulating CQA biosynthesis. This study analyzed transcriptome-wide identification of the phenylpropanoid to caffeoylquinic acid biosynthesis candidate genes in A. spathulifolius flowers and leaves. Transcriptomic analyses of the flowers and leaves showed a differential expression of the PPP and CQA biosynthesis regulated unigenes. An analysis of PPP-captive unigenes revealed a major duplication in the following genes: PAL, 120 unigenes in leaves and 76 in flowers; C3′H, 169 unigenes in leaves and 140 in flowers; 4CL, 41 unigenes in leaves and 27 in flowers; and C4H, 12 unigenes in leaves and 4 in flowers. The phylogenetic analysis revealed 82 BAHDs superfamily members in leaves and 72 in flowers, among which five unigenes encode for HQT and three for HCT. The three HQT are common to both leaves and flowers, whereas the two HQT were specialized for leaves. The pattern of HQT synthesis was upregulated in flowers, whereas HCT was expressed strongly in the leaves of A. spathulifolius. Overall, 4CL, C4H, and HQT are expressed strongly in flowers and CAA and HCT show more expression in leaves. As a result, the quantification of HQT and HCT indicates that CQA biosynthesis is more abundant in the flowers and synthesis of caffeic acid in the leaves of A. spathulifolius.

Metagenome sequencing to unveil the resistome in a deep subtropical lake on the Yunnan-Guizhou Plateau, China

Erhai Lake in the Yunnan-Guizhou Plateau, a typical subtropical deep lake, experienced long-time input of nutrients and pollutants, which provides beneficial conditions for evolution and spread of resistance genes. In this study, the profile of bacterial community, antibiotic and metal resistance genes in sediments of Erhai Lake were analyzed via high-thought sequencing of 16S rRNA and metagenomic DNA. Proteobacteria, Firmicutes, Nitrospirae, and Bacteroidetes were found to be the main composition of the bacterial community. Multidrug-, bacitracin-, macrolide-lincosamide-streptogramin (MLS)- and tetracycline resistance genes were the primary antibiotic resistance gene (ARG) types with high relative abundance, whereas Multi-metal-, and arsenic- resistance genes were the primary metal resistance gene (MRG) types. The effects of nitrogen and phosphorus on the abundance of ARGs and MRGs depended on the different types. Some specific ARG (fosmidomycin resistance genes) and MRG types (multi-metal resistance genes) also showed significant geographic distribution. Bacterial community was the main driver shaping the resistome. Nutrients also played an essential role in structuring the bacterial community and resistome in the sediments of Erhai Lake. This study sheds light on the distribution and fate of resistome under a high load of nitrogen and phosphorus in a deep lake.