Advancing RNA-Seq analysis

Selection and Validation of Reference Genes for Gene Expression Studies in Euonymus japonicus Based on RNA Sequencing

Euonymus japonicus is one of the most low-temperature-tolerant evergreen broad-leaved tree species in the world and is widely used in urban greening. However, there are very few molecular biology studies on its low-temperature tolerance mechanism. So far, no researcher has selected and reported on its reference genes. In this study, 21 candidate reference genes (12 traditional housekeeping genes and 9 other genes) were initially selected based on gene expression and coefficient of variation (CV) through RNA-Seq (unpublished data), and qRT-PCR was used to detect the expression levels of candidate reference genes in three different groups of samples (leaves under different temperature stresses, leaves of plants at different growth stages, and different organs). After further evaluating the expression stability of these genes using geNorm, NormFinder, Bestkeeper, and RefFind, the results show that the traditional housekeeping gene eIF5A and the new reference gene RTNLB1 have good stability in the three different groups of samples, so they are reference genes with universality. In addition, we used eIF5A and RTNLB1 as reference genes to calibrate the expression pattern of the target gene EjMAH1, which confirmed this view. This article is the first to select and report on the reference gene of E. japonicus, laying the foundation for its low-temperature tolerance mechanism and other molecular biology research.

Protein-coding potential of non-canonical open reading frames in human transcriptome

In recent years, proteogenomics and ribosome profiling studies have identified a large number of proteins encoded by noncoding regions in the human genome. They are encoded by small open reading frames (sORFs) in the untranslated regions (UTRs) of mRNAs and long non-coding RNAs (lncRNAs). These sORF encoded proteins (SEPs) are often <150AA and show poor evolutionary conservation. A subset of them have been functionally characterized and shown to play an important role in fundamental biological processes including cardiac and muscle function, DNA repair, embryonic development and various human diseases. How many novel protein-coding regions exist in the human genome and what fraction of them are functionally important remains a mystery. In this review, we discuss current progress in unraveling SEPs, approaches used for their identification, their limitations and reliability of these identifications. We also discuss functionally characterized SEPs and their involvement in various biological processes and diseases. Lastly, we provide insights into their distinctive features compared to canonical proteins and challenges associated with annotating these in protein reference databases.

De novo assembly of Iron-Heart Cunninghamia lanceolata transcriptome and EST-SSR marker development for genetic diversity analysis

Iron-Heart Cunninghamia lanceolata, a wild relative of Chinese fir with valuable genetic and breeding traits, has been limited in genetic studies due to a lack of genomic resources and markers. In this study, we conducted transcriptome sequencing of Iron-Heart C. lanceolata leaves using Illumina NovaSeq 6000 and performed assembly and analysis. We obtained 45,326,576 clean reads and 115,501 unigenes. Comparative analysis in five functional databases resulted in successful annotation of 26,278 unigenes, with 6,693 unigenes annotated in all databases (5.79% of the total). UniProt and Pfam databases provided annotations for 22,673 and 18,315 unigenes, respectively. Gene Ontology analysis categorized 23,962 unigenes into three categories. KEGG database alignment annotated 10,195 unigenes, classifying them into five categories: metabolism, genetic information, biological systems, cellular processes, and environmental information processing. From the unigenes, we identified 5,645 SSRs, with dinucleotides repeats being the most common (41.47%). We observed variations in repeat numbers and base compositions, with the majority of markers ranging from 12 to 29 bp in length. We randomly selected 200 primer pairs and successfully amplified 15 pairs of polymorphic SSR primers, which effectively distinguished Chinese fir plants of different origins. This study provides insights into the genetic characteristics of Iron-Heart C. lanceolata and offers a foundation for future molecular marker development, breeding programs, genetic diversity analysis, and conservation strategies.

Transcriptome analysis of the gene expression of M. iliotibialis lateralis affected by dietary methionine restriction

Introduction: Methionine (Met) is an important amino acid related to the development of skeletal muscle. This study investigated the effects of dietary Met restriction on the gene expression of M. iliotibialis lateralis. Methods: A total of 84 day-old broiler chicks (Zhuanghe Dagu) with a similar initial body weight (207.62 ± 8.54 g) were used in this study. All birds were divided into two groups (CON; L-Met) based on the initial body weight. Each group consisted of six replicates with seven birds per replicate. The experimental period was 63 days (phase 1, days 1-21; phase 2, days 22-63). According to the nutritional requirements of Zhuanghe Dagu chickens, we provided a basal diet (0.39% Met levels during phase 1 and 0.35% Met levels during phase 2, as-fed basis) to the birds in the CON group, while we provided a Met-restricted diet (0.31% Met levels during phase 1 and 0.28% Met levels during phase 2, as-fed basis) to the birds in the L-Met group. The growth performance of broiler chicks and their M. iliotibialis lateralis development parameters were measured on days 21 and 63. Results and Discussion: In this study, dietary Met restriction did not affect the growth performance of broiler chicks but hindered the development of M. iliotibialis lateralis at both sampling timepoints. On the final day, three birds selected from each group (three from CON and three from L-Met) were used to obtain M. iliotibialis lateralis samples from leg muscle for further transcriptome analysis. Transcriptome analysis revealed that dietary Met restriction significantly upregulated 247 differentially expressed genes (DEGs) and downregulated 173 DEGs. Additionally, DEGs were mainly enriched in 10 pathways. Among DEGs, we observed that dietary Met restriction downregulated the expression of CSRP3, KY, FHL1, LMCD1, and MYOZ2 in M. iliotibialis lateralis. Therefore, we considered that dietary Met restriction had negative effects on the development of M. iliotibialis lateralis, and CSRP3, KY, FHL1, LMCD1, and MYOZ2 may serve as potential functional genes involved in this process.

Physiological and transcriptomic analysis of OsLHCB3 knockdown lines in rice

The photosystem II (PSII) outer antenna LHCB3 protein plays critical roles in distributing the excitation energy and modulating the rate of state transition for photosynthesis. Here, OsLHCB3 knockdown mutants were produced using the RNAi system. Phenotypic analyses showed that OsLHCB3 knockdown led to pale green leaves and lower chlorophyll contents at both tillering and heading stages. In addition, mutant lines exhibited decreased non-photochemical quenching (NPQ) capacity and net photosynthetic rate (Pn) by downregulating the expression of PSII-related genes. Moreover, RNA-seq experiments were performed at both tillering and heading stages. The differentially expressed genes (DEGs) mainly involved in chlorophyll binding response to abscisic acid, photosystem II, response to chitin, and DNA-binding transcription factor. Besides, our transcriptomic and physiological data indicated that OsLHCB3 was essential for binding chlorophyll, but not for the metabolism of chlorophyll in rice. OsLHCB3 RNAi knockdown plants affected the expression of PS II-related genes, but not PS I-related genes. Overall, these results suggest that OsLHCB3 also plays vital roles in regulating photosynthesis and antenna proteins in rice as well as responses to environment stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01387-z.

Integrative analysis of transcriptome and metabolome reveals probiotic effects on cecal metabolism in broilers

BACKGROUND

Probiotics play an important role in the host and have attracted widespread attention as an alternative to antibiotics. Arbor Acres broilers were used in the present experiment and fed different doses of compound probiotics at 1, 5, and 10 g kg^-1 . The effects of compound probiotics on broiler growth performance and cecal transcriptome and metabolome were investigated.

RESULTS

We discovered 425 differentially expressed genes (DEGs; upregulated: 256; downregulated: 169) in the cecal transcriptome study. These DEGs were assigned to fat metabolic pathways, such as the peroxisome proliferator-activated receptor (PPAR) signaling pathway, according to KEGG analysis. Probiotics downregulated LPL and upregulated PPARα expression in the cecum. In metabolome analysis of the cecum of cecum, we screened 86 differential metabolites and performed KEGG enrichment analysis of these metabolites. The KEGG analysis showed that these differentially expressed metabolites were annotated to nucleotide metabolism-related pathways, such as purine metabolism. In the cecum, probiotics upregulated the content of guanine, AMP, 3'-AMP, adenylosuccinate, deoxyguanosine, and ADP-ribose, whereas they downregulated the content of 5-hydroxyisourate. Comprehensive transcriptome and metabolome analysis revealed that glycolysis, gluconeogenesis, and glycerophospholipid metabolism pathways were jointly enriched in cecum of broilers fed a probiotic-containing diet.

CONCLUSION

This study provides valuable information for studying the regulation and gene metabolism network of probiotics on cecal metabolism in broilers. © 2022 Society of Chemical Industry.

Łabaj P, Mangul S. RNA-seq data science: From raw data to effective interpretation

RNA sequencing (RNA-seq) has become an exemplary technology in modern biology and clinical science. Its immense popularity is due in large part to the continuous efforts of the bioinformatics community to develop accurate and scalable computational tools to analyze the enormous amounts of transcriptomic data that it produces. RNA-seq analysis enables genes and their corresponding transcripts to be probed for a variety of purposes, such as detecting novel exons or whole transcripts, assessing expression of genes and alternative transcripts, and studying alternative splicing structure. It can be a challenge, however, to obtain meaningful biological signals from raw RNA-seq data because of the enormous scale of the data as well as the inherent limitations of different sequencing technologies, such as amplification bias or biases of library preparation. The need to overcome these technical challenges has pushed the rapid development of novel computational tools, which have evolved and diversified in accordance with technological advancements, leading to the current myriad of RNA-seq tools. These tools, combined with the diverse computational skill sets of biomedical researchers, help to unlock the full potential of RNA-seq. The purpose of this review is to explain basic concepts in the computational analysis of RNA-seq data and define discipline-specific jargon.

Production performance, egg quality, and uterine gene expression for layers as affected by N-Carbamylglutamate supplementation

Introduction

Supplementation of exogenous additives is a strategy to improve laying performance of layers by regulating uterine function. N-Carbamylglutamate (NCG) as an activator for endogenous arginine synthesis has the potential to regulate the laying performance of layers, but its effects have not been fully understood.

Methods

This study investigated the effects of dietary supplementation of NCG on production performance, egg quality, and uterine gene expression in layers. A total of 360 45-week-old layers with a genetic line of Jinghong No. 1 were used in this study. The experimental period was 14 weeks. All birds were divided into 4 treatments with 6 replicates per treatment and 15 birds per replicate. Dietary treatments were based on a basal diet and supplemented with 0, 0.08, 0.12, or 0.16% NCG to form C, N1, N2, and N3 groups.

Results and discussion

We found that layers in group N1 had higher egg production rate than those in group C. Egg weight was significantly reduced, while eggshell thickness was significantly improved, by treatment. However, the albumen height and Haugh unit were the lowest in group N3. Based on the above results, groups C and N1 were selected for further transcriptomics analysis of uterine tissue by RNA-seq. More than 7.4 Gb clean reads and 19,882 tentative genes were obtained using the Gallus gallus genome as a reference. Transcriptomics analysis in uterus tissue revealed that 95 differentially expressed genes (DEGs) were upregulated and 127 DEGs were downregulated. Functional annotation and pathway enrichment analysis showed that DEGs in uterine tissue were mainly enriched in glutathione metabolism, cholesterol metabolism, and glycerolipid metabolism, etc. Vitamin A metabolism-related gene, RBP1, nutrient transport-related gene, ALB, protein synthesis-related gene, METTL21C, and calcium transport-related gene, RYR2, CACNB2, RAMP3, and STAC, were significantly regulated by 0.08% NCG supplementation. Therefore, we concluded that NCG supplementation at a dose of 0.08% improved production performance and egg quality of layers by regulating uterus function.

Expression and Functional Studies of Leaf, Floral, and Fruit Developmental Genes in Non-model Species

Researchers working on evolutionary developmental plant biology are inclined to choose non-model taxa to address how specific features have been acquired during ontogeny and fixed during phylogeny. In this chapter we describe methods to extract RNA, to assemble de-novo transcriptomes, to isolate orthologous genes within gene families, and to evaluate expression and function of target genes. We have successfully optimized these protocols for non-model plant species including ferns, gymnosperms, and a large assortment of angiosperms. In the latter, we have ranged a large number of families including Aristolochiaceae, Apodanthaceae, Chloranthaceae, Orchidaceae, Papaveraceae, Rubiaceae, Solanaceae, and Tropaeolaceae.

Regulation of mRNA and miRNA in the response to Salmonella enterica serovar Enteritidis infection in chicken cecum

BACKGROUND

Salmonella enterica, serovar Enteritidis (SE) is a food-borne pathogen, which can cause great threat to human health through consumption of the contaminated poultry products. Chicken is the main host of SE. The mRNA and microRNA (miRNA) expression profiles were analyzed on cecum of Shouguang chicken via next-generation sequencing and bioinformatics approaches. The treated group was inoculated SE, and the control group was inoculated with phosphate buffer saline (PBS).

RESULTS

There were 1760 differentially expressed mRNAs in the SE-infected group, of which 1046 were up-regulated mRNA, and 714 were down-regulated mRNA. In addition, a total of 821 miRNAs were identified, and 174 miRNAs were differentially expressed, of which 100 were up-regulated and 74 were down-regulated. Functional enrichment of differentially expressed mRNAs was similar to miRNA target genes. The functional analysis results of differentially expressed mRNAs and miRNAs were performed. Immune-related processes and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were enriched by up-regulated mRNA. The down-regulated mRNAs were enriched in tissue development and metabolic-related KEGG pathways. The functional analysis of up-regulated miRNA target genes was similar to the down-regulated mRNAs. The down-regulated miRNA target genes were enriched in metabolic-related GO (Gene Ontology) -BP (Biological process) terms and KEGG pathways. The overlap of the up-regulated mRNA and the up-regulated miRNA target genes (class I) was 325, and the overlap of the down-regulated miRNA target genes (class II) was 169. The class I enriched in the immune-related GO-BP terms and KEGG pathways. The class II mainly enriched in metabolic-related GO-BP terms and KEGG pathways. Then we detected the expression of mRNA and miRNA through qRT-PCR. The results shown that the expression of HHIP, PGM1, HTR2B, ITGB5, RELN, SFRP1, TCF7L2, SCNN1A, NEK7, miR-20b-5p, miR-1662, miR-15a, miR-16-1-3p was significantly different between two groups. Dual-luciferase reporter assay was used to detect the relationship between miR-20b-5p and SCNN1A. The result indicated that miR-20b-5p regulate immune or metabolic responses after SE infection in Shouguang chickens by directly targeting SCNN1A.

CONCLUSIONS

The findings here contribute to the further analysis of the mechanism of mRNA and miRNA defense against SE infection, and provide a theoretical foundation for the molecular disease-resistant breeding of chickens.

Current Scenario of Exogenously Induced RNAi for Lepidopteran Agricultural Pest Control: From dsRNA Design to Topical Application

Invasive insects cost the global economy around USD 70 billion per year. Moreover, increasing agricultural insect pests raise concerns about global food security constraining and infestation rising after climate changes. Current agricultural pest management largely relies on plant breeding-with or without transgenes-and chemical pesticides. Both approaches face serious technological obsolescence in the field due to plant resistance breakdown or development of insecticide resistance. The need for new modes of action (MoA) for managing crop health is growing each year, driven by market demands to reduce economic losses and by consumer demand for phytosanitary measures. The disabling of pest genes through sequence-specific expression silencing is a promising tool in the development of environmentally-friendly and safe biopesticides. The specificity conferred by long dsRNA-base solutions helps minimize effects on off-target genes in the insect pest genome and the target gene in non-target organisms (NTOs). In this review, we summarize the status of gene silencing by RNA interference (RNAi) for agricultural control. More specifically, we focus on the engineering, development and application of gene silencing to control Lepidoptera through non-transforming dsRNA technologies. Despite some delivery and stability drawbacks of topical applications, we reviewed works showing convincing proof-of-concept results that point to innovative solutions. Considerations about the regulation of the ongoing research on dsRNA-based pesticides to produce commercialized products for exogenous application are discussed. Academic and industry initiatives have revealed a worthy effort to control Lepidoptera pests with this new mode of action, which provides more sustainable and reliable technologies for field management. New data on the genomics of this taxon may contribute to a future customized target gene portfolio. As a case study, we illustrate how dsRNA and associated methodologies could be applied to control an important lepidopteran coffee pest.

Transcriptome and metabolome analyses reveal the regulatory effects of compound probiotics on cecal metabolism in heat-stressed broilers

The effect of compound probiotics on the caecum of broilers under heat stress was assessed in this study. A total of 400 twenty-eight-day-old AA male broilers were randomly divided into 4 treatment groups, where each group had 5 replicates of 20 broilers. The 4 treatment groups were a heat stress control group (broilers receiving a normal diet) and groups HP I, HP II, and HP Ⅲ, consisting of broilers receiving 1, 5, and 10 g of compound probiotics added to each kilogram of feed, respectively. Compound probiotics (L. casei, L. acidophilus, and B. lactis at a ratio of 1:1:2) were used to formulate a compound probiotic powder, with 1 × 10¹⁰ CFU/g of effective viable bacteria. Heat stress treatment was performed at 32 ± 1°C from 9:00 to 17:00 every day from 28 d to 42 d. In d 28 to 42, compared with the HC group, the ADG of broilers in the HP II and III groups was significantly increased (P < 0.05); the ADFI difference between groups was not significant (P > 0.05); the FCR of HP II and III broilers was significantly decreased (P < 0.05); and the FCR of the HP I group increased, but the difference was not significant (P > 0.05). Transcriptome results demonstrate that 665 differential genes were screened (DEGs; upregulated: 366, downregulated: 299). The DEGs were enriched in the B cell receptor signaling pathway, the intestinal immune network for IgA synthesis, the Fc epsilon RI signaling pathway, and other signaling pathways, according to KEGG enrichment analysis. Metabolome analysis identified 92 differential metabolites (DAMs; upregulated: 48, downregulated: 44). KEGG enrichment analysis indicated significant enrichment of Pantothenate and CoA biosynthesis and beta-Alanine metabolism. The combined transcriptome and metabolome analysis revealed that the DAMs and DEGs were mostly involved in beta-alanine metabolism, arginine biosynthesis, amino sugar and nucleotide sugar, and alanine, aspartate, and glutamate metabolism. The results of this study suggest that the addition of compound probiotics has a positive effect on intestinal metabolites, improving the growth performance and contributing to the overall health of broilers under heat stress.

The study of cancer cell in stromal environment through induced pluripotent stem cell-derived mesenchymal stem cells

BACKGROUND

The development of mesenchymal stem cells (MSCs) has gained reputation from its therapeutic potential in stem cell regeneration, anti-inflammation, tumor suppression, and drug delivery treatment. Previous studies have shown MSCs have both promoting and suppressing effects against cancer cells. While the limitation of obtaining a large quantity of homologous MSCs for studies and treatment remains a challenge, an alternative approach involving the production of MSCs derived from induced pluripotent stem cells (iPSCs; induced MSCs [iMSCs]) may be a promising prospect given its ability to undergo prolonged passage and with similar therapeutic profiles as that of their MSC counterparts. However, the influence of iMSC in the interaction of cancer cells remains to be explored as such studies are not well established. In this study, we aim to differentiate iPSCs into MSC-like cells as a potential substitute for adult MSCs and evaluate its effect on non-small-cell lung cancer (NSCLC).

METHODS

iMSCs were derived from iPSCs and validated with reference to the International Society of Cellular Therapy guidelines on MSC criteria. To create a stromal environment, the conditioned medium (CM) of iMSCs was harvested and applied for coculturing of NSCLC of H1975 at different concentrations. The H1975 was then harvested for RNA extraction and subjected to next-generation sequencing (NGS) for analysis.

RESULTS

The morphology of iMSCs-CM-treated H1975 was different from an untreated H1975. Our NGS data suggest the occurrence of apoptotic events and the presence of cytokines from H1975's RNA that are treated with iMSCs-CM.

CONCLUSION

Our results have shown that iMSCs may suppress the growth of H1975 by releasing proapoptotic cytokines into coculture media. Using iPSC-derived MSC models allows a deeper study of tumor cross talk between MSC and cancer cells that can be applied for potential future cancer therapy.

Integrated Transcriptome and Metabolomic Analysis Reveal Anti-Angiogenic Properties of Disarib, a Novel Bcl2-Specific Inhibitor

Transcriptomic profiling of several drugs in cancer cell lines has been utilised to obtain drug-specific signatures and guided combination therapy to combat drug resistance and toxicity. Global metabolomics reflects changes due to altered activity of enzymes, environmental factors, etc. Integrating transcriptomics and metabolomics can provide genotype-phenotype correlation, providing meaningful insights into alterations in gene expression and its outcome to understand differential metabolism and guide therapy. This study uses a multi-omics approach to understand the global gene expression and metabolite changes induced by Disarib, a novel Bcl2-specific inhibitor in the Ehrlich adenocarcinoma (EAC) breast cancer mouse model. RNAseq analysis was performed on EAC mouse tumours treated with Disarib and compared to the controls. The expression of 6 oncogenes and 101 tumour suppressor genes interacting with Bcl2 and Bak were modulated upon Disarib treatment. Cancer hallmark pathways like DNA repair, Cell cycle, angiogenesis, and mitochondrial metabolism were downregulated, and programmed cell death platelet-related pathways were upregulated. Global metabolomic profiling using LC-MS revealed that Oncometabolites like carnitine, oleic acid, glycine, and arginine were elevated in tumour mice compared to normal and were downregulated upon Disarib treatment. Integrated transcriptomic and metabolomic profiles identified arginine metabolism, histidine, and purine metabolism to be altered upon Disarib treatment. Pro-angiogenic metabolites, arginine, palmitic acid, oleic acid, and myristoleic acid were downregulated in Disarib-treated mice. We further validated the effect of Disarib on angiogenesis by qRT-PCR analysis of genes in the VEGF pathway. Disarib treatment led to the downregulation of pro-angiogenic markers. Furthermore, the chorioallantoic membrane assay displayed a reduction in the formation of the number of secondary blood vessels upon Disarib treatment. Disarib reduces tumours by reducing oncometabolite and activating apoptosis and downregulating angiogenesis.

MicroSugar: A database of comprehensive miRNA target prediction framework for sugarcane (Saccharum officinarum L.)

microRNA (miRNA) is a group of small non-coding RNA that plays important role in post-transcription of gene expression. With the studies about miRNA increase in sugarcane, the researchers lack an exhaustive resource to achieve the data. To fill this gap, we developed MicroSugar, a database that supported mRNA and miRNA annotation for sugarcane (http://suc.gene-db.com). MicroSugar is an integrated resource developed for 194,528 genes including 80,746 unigenes from long reads of Pacbio platform and 468 miRNAs from 72 samples. Internode elongation (jointing) is the key biological characteristic for the growth of sugarcane tillers into sugarcane stems. The present study combined the sequencing data from the different stages in internode elongation of stem and tiller. In total, the 14,300 3' untranslated region (UTR) sequences were extracted from the gene sequences and 3019 mRNAs as target of 327 miRNA were identified by miRanda algorithm and Spearman's Rho of expression levels. To determine the gene functions regulated by these miRNAs, the gene ontology enrichment analysis was performed and it confirmed that the over-represented Gene Ontology (GO) terms were associated with organism formation indicating the growth controlling function by miRNAs in sugarcane. Moreover, MicroSugar is a comprehensive and integrated database with a user-friendly responsive template. By browsing, searching and downloading of the nucleotide sequences, expression and miRNA targets, the user can retrieve information promptly. The database provides a valuable resource to facilitate the understanding of miRNA in sugarcane development and growth which will contribute to the study of sugarcane and other plants.

Transcriptional response of the azooxanthellate octocoral Scleronephthya gracillimum to seawater acidification and thermal stress

The stress responses to increased seawater temperature and marine acidification were investigated using a microarray to reveal transcriptional changes in S. gracillimum. For the study, corals were exposed to different stress experiments; high temperature only (26 °C, 28 °C and 30 °C), low-pH only (pH 7.5, pH 7.0 and pH 6.5) and dual stress experiments (28 °C + pH 7.8, 28 °C + pH 7.5 and 28 °C + pH 7.0), mortality and morphological changes in 24 h exposure experiments were investigated. The survival rates of each experimental group were observed. The gene expression changes in single and dual stress exposed coals were measured and the differentially expressed genes were classified with gene ontology analysis. The top three enriched gene ontology terms of DEGs in response to dual stress were metal ion binding (23.4%), extracellular region (17.2%), and calcium ion binding (12.8%). The gene showing the greatest increase in expression as a response to the dual stress was hemagglutinin/amebocyte aggregation factor, followed by interferon-inducible GTPase 5 and the gene showing the greatest decrease as a response to the dual stress was Fas-associating death domain-containing protein, followed by oxidase 2. These results represented the transcriptomic study focused on the stress responses of the temperate asymbiotic soft coral exposed to single and dual stresses. The combined effect of thermal and acidification stress on corals triggered the negative regulation of ion binding and extracellular matrix coding genes and these genes might serve as a basis for research into coral-specific adaptations to stress responses and global climate change.

The cwp66 Gene Affects Cell Adhesion, Stress Tolerance, and Antibiotic Resistance in Clostridioides difficile

Clostridioides difficile is a Gram-positive, spore-forming anaerobic bacteria that is one of the leading causes of antibiotic-associated diarrhea. The cell wall protein 66 gene (cwp66) encodes a cell wall protein, which is the second major cell surface antigen of C. difficile. Although immunological approaches, such as antibodies and purified recombinant proteins, have been implemented to study the role of Cwp66 in cell adhesion, no deletion mutant of the cwp66 gene has yet been characterized. We constructed a cwp66 gene deletion mutant using Clustered Regularly Interspaced Short Palindromic Repeats Cpf1 (CRISPR-Cpf1) system. The phenotypic and transcriptomic changes of the Δcwp66 mutant compared with the wild-type (WT) strain were studied. The deletion of the cwp66 gene led to the decrease of cell adhesive capacity, cell motility, and stresses tolerance (to Triton X-100, acidic environment, and oxidative stress). Interestingly, the Δcwp66 mutant is more sensitive than the WT strain to clindamycin, ampicillin, and erythromycin but more resistant than the latter to vancomycin and metronidazole. Moreover, mannitol utilization capability in the Δcwp66 mutant was lost. Comparative transcriptomic analyses indicated that (i) 22.90-fold upregulation of cwpV gene and unable to express gpr gene were prominent in the Δcwp66 mutant; (ii) the cwp66 gene was involved in vancomycin resistance of C. difficile by influencing the expression of d-Alanine-d-Alanine ligase; and (iii) the mannose/fructose/sorbose IIC and IID components were upregulated in Δcwp66 mutant. The present work deepens our understanding of the contribution of the cwp66 gene to cell adhesion, stress tolerance, antibiotic resistance, and mannitol transportation of C. difficile. IMPORTANCE The cell wall protein 66 gene (cwp66) encodes a cell wall protein, which is the second major cell surface antigen of C. difficile. Although immunological approaches, such as antibodies and purified recombinant proteins, have been implemented to study the role of Cwp66 in cell adhesion, no deletion mutant of the cwp66 gene has yet been characterized. The current study provides direct evidence that the cwp66 gene serves as a major adhesion in C. difficile, and also suggested that deletion of the cwp66 gene led to the decrease of cell adhesive capacity, cell motility, and stresses tolerance (to Triton X-100, acidic environment, and oxidative stress). Interestingly, the antibiotic resistance and carbon source utilization profiles of the Δcwp66 mutant were significantly changed. These phenotypes were detrimental to the survival and pathogenesis of C. difficile in the human gut and may shed light on preventing C. difficile infection.

Phosphoryl guanidine oligonucleotides as primers for RNA-dependent DNA synthesis using murine leukemia virus reverse transcriptase

Modern approaches to the detection and analysis of low-copy-number RNAs are often based on the use of RNA-dependent DNA polymerases, for example, in reverse-transcription PCR. The accuracy and eff iciency of cDNA synthesis in the reverse-transcription reaction catalyzed by reverse transcriptase (RNA-dependent DNA polymerase) signif icantly affect the correctness of the results of PCR diagnostic assays and/or RNA sequencing. In this regard, many studies are focused on the optimization of the reverse-transcription reaction, including the search for more perfect primers necessary to obtain a full-length DNA copy of RNA under study. The best-known completely uncharged analogs of oligonucleotides – morpholine oligonucleotides and peptide nucleic acids – cannot be substrates for enzymes that process nucleic acids. The aim of this work was to conduct a pilot study of uncharged phosphoryl guanidine oligodeoxyribonucleotides (PGOs) as primers for mouse leukemia virus reverse transcriptase (MMLV H-). Specif ic features of elongation of partially and completely uncharged PGO primers were investigated. It was demonstrated that PGOs can be elongated eff iciently, e. g., in the presence of a fragment of human ribosomal RNA having complex spatial structure. It was shown that the proportion (%) of abortive elongation products of a PGO primer depends on buffer ionic strength, nucleotide sequence of the primer, and the presence and location of phosphoryl guanidine groups in the primer. The results indicate the suitability of PGOs, including completely electroneutral ones, as primers for reverse-transcription PCR, thereby opening up new prospects for the creation of experimental models for the analysis of highly structured RNA.

Transcriptome Sequencing and Screening of Anthocyanin-Related Genes in the Leaves of Acer truncatum Bunge

Acer truncatum Bunge is generally used as an ornamental tree because of its autumn leaves, although the viewing period is short-approximately 7-15 days. Color improvement of ornamental trees has consistently been an important research topic because color partially determines the value of the commodity; however, a lack of genomic data have limited the progress of molecular breeding research in this area. The purposes of this study were to obtain a transcriptome database for A. truncatum, screen anthocyanin biosynthesis-related genes, and reveal the mechanisms underlying leaf color transformation to provide a basis for increasing the viewing period or breeding cultivars that display red leaves throughout the growing season via gene regulation. In this study, although the use of an Illumina HiSeq 2000 platform and systematic bioinformatics analysis using both young and mature leaves as experimental materials, 233,912,882 clean reads were generated and 121,287 unique transcripts were retrieved. We selected 16 color-related genes (from the transcriptome results) for qRT-PCR to validate the results, and the expression trends of the selected genes were largely consistent with the transcriptome analysis results, with a consistency of 0.875. According to the results of the transcriptome analysis, the validation, and previous studies, we obtained sequences of genes related to anthocyanins, including CHS, CHI, ANS, UFGT, UGT75c1, DFR, BZ1, F3H, F3'H, LAR, ANR, FLS, and those of several transcription factors, including MYB1, BHLH, and WD40. Verifying specific regulation by one or several of these genes in the control of leaf color requires further research. The acquisition of transcriptomic information, especially information concerning anthocyanin biosynthesis-related genes and their base sequences, can provide a theoretical basis for the study of the molecular mechanisms determining changes in leaf color in Acer and is of great importance to the breeding of new cultivars.

Identification and Quantification of Microexons Using Bulk and Single-Cell RNA-Seq Data

The analysis of RNA-seq has greatly improved the characterization and understanding of the transcriptome. In particular, RNA-seq experiments have extended catalogs of alternative splicing events. However, the analysis of RNAs-seq data for detection and quantification of microexons, extremely short exons of length up to 30 nt, require specialized computational workflows. Here, we describe MicroExonator, a reproducible computational workflow for microexon splicing analysis using bulk or single-cell RNA-seq data.

Targeted RNA-sequencing analysis for fusion transcripts detection in tumor diagnostics: assessment of bioinformatic tools reliability in FFPE samples

Aim:

Diagnostic laboratories are progressively introducing next-generation sequencing (NGS) technologies in the routine workflow to meet the increasing clinical need for comprehensive molecular characterization in cancer patients for diagnosis and precision medicine, including fusion-transcripts detection. Nevertheless, the low quality of messenger RNA (mRNA) extracted from formalin-fixed paraffin-embedded (FFPE) samples may affect the transition from traditional single-gene testing approaches [like fluorescence in situ hybridization (FISH), immunohistochemistry (IHC), or polymerase chain reaction (PCR)] to NGS. The present study is aimed at assessing the overall accuracy of RNA fusion transcripts detection by NGS analysis in FFPE samples in real-world diagnostics.

Methods:

Herein, NGS data from 190 soft tissue tumors (STTs) and carcinoma cases, discussed in the context of the institutional Molecular Tumor Board, are reported and analyzed by FusionPlex^© Solid tumor kit through the manufacturer’s pipeline and by two well-known fast and accurate open-source tools [Arriba (ARR) and spliced transcripts alignment to reference (STAR)-fusion (SFU)].

Results:

The combination of FusionPlex^© Solid tumor with ArcherDX^® Analysis suite (ADx) analysis package has been proven to be sensitive and specific in STT samples, while partial loss of sensitivity has been found in carcinoma specimens.

Conclusions:

Albeit ARR and SFU showed lower sensitivity, the use of additional fusion-detection tools can contribute to reinforcing or extending the output obtained by ADx, particularly in the case of low-quality input data. Overall, our results sustain the clinical use of NGS for the detection of fusion transcripts in FFPE material.

Genome annotation with long RNA reads reveals new patterns of gene expression and improves single-cell analyses in an ant brain

BACKGROUND

Functional genomic analyses rely on high-quality genome assemblies and annotations. Highly contiguous genome assemblies have become available for a variety of species, but accurate and complete annotation of gene models, inclusive of alternative splice isoforms and transcription start and termination sites, remains difficult with traditional approaches.

RESULTS

Here, we utilized full-length isoform sequencing (Iso-Seq), a long-read RNA sequencing technology, to obtain a comprehensive annotation of the transcriptome of the ant Harpegnathos saltator. The improved genome annotations include additional splice isoforms and extended 3' untranslated regions for more than 4000 genes. Reanalysis of RNA-seq experiments using these annotations revealed several genes with caste-specific differential expression and tissue- or caste-specific splicing patterns that were missed in previous analyses. The extended 3' untranslated regions afforded great improvements in the analysis of existing single-cell RNA-seq data, resulting in the recovery of the transcriptomes of 18% more cells. The deeper single-cell transcriptomes obtained with these new annotations allowed us to identify additional markers for several cell types in the ant brain, as well as genes differentially expressed across castes in specific cell types.

CONCLUSIONS

Our results demonstrate that Iso-Seq is an efficient and effective approach to improve genome annotations and maximize the amount of information that can be obtained from existing and future genomic datasets in Harpegnathos and other organisms.

Comparability of reference-based and reference-free transcriptome analysis approaches at the gene expression level

Background

Lately, high-throughput RNA sequencing has been extensively used to elucidate the transcriptome landscape and dynamics of cell types of different species. In particular, for most non-model organisms lacking complete reference genomes with high-quality annotation of genetic information, reference-free (RF) de novo transcriptome analyses, rather than reference-based (RB) approaches, are widely used, and RF analyses have substantially contributed toward understanding the mechanisms regulating key biological processes and functions. To date, numerous bioinformatics studies have been conducted for assessing the workflow, production rate, and completeness of transcriptome assemblies within and between RF and RB datasets. However, the degree of consistency and variability of results obtained by analyzing gene expression levels through these two different approaches have not been adequately documented.

Results

In the present study, we evaluated the differences in expression profiles obtained with RF and RB approaches and revealed that the former tends to be satisfactorily replaced by the latter with respect to transcriptome repertoires, as well as from a gene expression quantification perspective. In addition, we urge cautious interpretation of these findings. Several genes that are lowly expressed, have long coding sequences, or belong to large gene families must be validated carefully, whenever gene expression levels are calculated using the RF method.

Conclusions

Our empirical results indicate important contributions toward addressing transcriptome-related biological questions in non-model organisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04226-0.

Metabolite and transcriptome analyses revealed the modulation of fructo-oligosaccharide on ileum metabolism of Taiping chickens

AIM

The metabolic markers and differentially expressed genes (DEGs) related to fructo-oligosaccharide (FOS) were screened, and the response of FOS to the ileum metabolic pathway of Taiping chickens was analysed.

METHODS AND RESULTS

Prebiotic are widely used in agricultural breeding for care and maintenance of animal health, especially FOS. Metabonomics evaluation of ileum of Taiping chicken ultra-performance liquid chromatography-quadruple time of-flight high-sensitivity mass spectrometry showed that 93 differentially altered metabolites were identified and divided into eight categories, of which organic acids and derivatives was the most important one. Transcriptomic analysis showed that DEGs were mainly enriched in drug metabolism-cytochrome p450, metabolism of xenobiotics by cytochrome p450, retinol metabolism and fat digestion and absorption. Integrated analysis of metabolite profiles and gene expression revealed that the significantly up-regulated GSTT1 was significantly correlated with most of the different lipid metabolites, suggesting that GSTT1 may play an important role in FOS regulation of lipid metabolism.

CONCLUSIONS

The results of this study suggest that supplementation of FOS can have a positive effect on gut metabolites, which may contribute to the overall health with indigenous chickens.

SIGNIFICANCE AND IMPACT OF THE STUDY

Insight into the responses of intestinal prebiotics of Taiping chicken is helpful to understand the role of prebiotics in maintaining intestinal microflora balance and improving immune response and productivity of poultry from the molecular and metabolic levels.

Comparative assessment of long-read error correction software applied to Nanopore RNA-sequencing data

MOTIVATION

Nanopore long-read sequencing technology offers promising alternatives to high-throughput short read sequencing, especially in the context of RNA-sequencing. However this technology is currently hindered by high error rates in the output data that affect analyses such as the identification of isoforms, exon boundaries, open reading frames and creation of gene catalogues. Due to the novelty of such data, computational methods are still actively being developed and options for the error correction of Nanopore RNA-sequencing long reads remain limited.

RESULTS

In this article, we evaluate the extent to which existing long-read DNA error correction methods are capable of correcting cDNA Nanopore reads. We provide an automatic and extensive benchmark tool that not only reports classical error correction metrics but also the effect of correction on gene families, isoform diversity, bias toward the major isoform and splice site detection. We find that long read error correction tools that were originally developed for DNA are also suitable for the correction of Nanopore RNA-sequencing data, especially in terms of increasing base pair accuracy. Yet investigators should be warned that the correction process perturbs gene family sizes and isoform diversity. This work provides guidelines on which (or whether) error correction tools should be used, depending on the application type.

BENCHMARKING SOFTWARE

https://gitlab.com/leoisl/LR_EC_analyser.

RNA-Seq-Based Profiling of pl Mutant Reveals Transcriptional Regulation of Anthocyanin Biosynthesis in Rice (Oryza sativa L.)

Purple-colored leaves in plants attain much interest for their important biological functions and could be a potential source of phenotypic marker in selecting individuals in breeding. The transcriptional profiling helps to precisely identify mechanisms of leaf pigmentation in crop plants. In this study, two genetically unlike rice genotypes, the mutant purple leaf (pl) and wild (WT) were selected for RNA-sequencing and identifying the differentially expressed genes (DEGs) that are regulating purple leaf color. In total, 609 DEGs were identified, of which 513 and 96 genes were up- and down-regulated, respectively. The identified DEGs are categorized into metabolic process, carboxylic acid biosynthesis, phenylpropanoids, and phenylpropanoid biosynthesis process enrichment by GO analysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) confirmed their association with phenylpropanoid synthesis, flavonoid synthesis, and phenylalanine metabolism. To explore molecular mechanism of purple leaf color, a set of anthocyanin biosynthetic and regulatory gene expression patterns were checked by qPCR. We found that OsPAL (Os02g0626100, Os02g0626400, Os04g0518400, Os05g0427400 and Os02g0627100), OsF3H (Os03g0122300), OsC4HL (Os05g0320700), and Os4CL5 (Os08g0448000) are associated with anthocyanin biosynthesis, and they were up-regulated in pl leaves. Two members of regulatory MYB genes (OsMYB55; Os05g0553400 and Os08g0428200), two bHLH genes (Os01g0196300 and Os04g0300600), and two WD40 genes (Os11g0132700 and Os11g0610700) also showed up-regulation in pl mutant. These genes might have significant and vital roles in pl leaf coloration and could provide reference materials for further experimentation to confirm the molecular mechanisms of anthocyanin biosynthesis in rice.

Transcriptomics and Metabolomics Reveal Purine and Phenylpropanoid Metabolism Response to Drought Stress in Dendrobium sinense, an Endemic Orchid Species in Hainan Island

Drought stress is a bottleneck factor for plant growth and development, especially in epiphytic orchids that absorb moisture mainly from the air. Recent studies have suggested that there are complex transcriptional regulatory networks related to drought stress in Dendrobium sinense. In this study, the transcription and metabolite alterations involved in drought stress response in D. sinense were investigated through RNA-seq and metabolomics. A total of 856 metabolites were identified from stressed and control samples, with 391 metabolites showing significant differences. With PacBio and Illumina RNA sequencing, 72,969 genes were obtained with a mean length of 2,486 bp, and 622 differentially expressed genes (DEGs) were identified. Correlation analysis showed 7 differential genes, and 39 differential metabolites were involved in interaction networks. The network analysis of differential genes and metabolites suggested that the pathways of purine metabolism and phenylpropanoid biosynthesis may play an important role in drought response in D. sinense. These results provide new insights and reference data for culturally important medicinal plants and the protection of endangered orchids.

Antimicrobial Peptides Derived From Insects Offer a Novel Therapeutic Option to Combat Biofilm: A Review

Biofilms form a complex layer with defined structures, that attach on biotic or abiotic surfaces, are tough to eradicate and tend to cause some resistance against most antibiotics. Several studies confirmed that biofilm-producing bacteria exhibit higher resistance compared to the planktonic form of the same species. Antibiotic resistance factors are well understood in planktonic bacteria which is not so in case of biofilm producing forms. This may be due to the lack of available drugs with known resistance mechanisms for biofilms. Existing antibiotics cannot eradicate most biofilms, especially of ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species). Insects produce complex and diverse set of chemicals for survival and defense. Antimicrobial peptides (AMPs), produced by most insects, generally have a broad spectrum of activity and the potential to bypass the resistance mechanisms of classical antibiotics. Besides, AMPs may well act synergistically with classical antibiotics for a double-pronged attack on infections. Thus, AMPs could be promising alternatives to overcome medically important biofilms, decrease the possibility of acquired resistance and treatment of multidrug-resistant pathogens including ESKAPE. The present review focuses on insect-derived AMPs with special reference to anti-biofilm-based strategies. It covers the AMP composition, pathways and mechanisms of action, the formation of biofilms, impact of biofilms on human diseases, current strategies as well as therapeutic options to combat biofilm with antimicrobial peptides from insects. In addition, the review also illustrates the importance of bioinformatics tools and molecular docking studies to boost the importance of select bioactive peptides those can be developed as drugs, as well as suggestions for further basic and clinical research.

Study on the differences of gene expression between pear and apple wild cultivation materials based on RNA-seq technique

BACKGROUND

Pears and apples are both perennial deciduous trees of the Rosaceae family, and both are important economic fruit trees worldwide. The emergence of many varieties in the market has been mostly domesticated from wild to cultivated and regulated by the differential expression of genes. However, the molecular process and pathways underlying this phenomenon remain unclear. Four typical wild and cultivar pear and apple trees at three developmental stages were used in our study to investigate the molecular process at the transcriptome level.

RESULT

Physiological observations indicated the obvious differences of size, weight, sugar acid content and peel color in wild and cultivar fruit among each developmental stage. Using next-generation sequencing based RNA-seq expression profiling technology, we produced a transcriptome in procession of a large fraction of annotated pear and apple genes, and provided a molecular basis underlying the phenomenon of wild and cultivar fruit tree differences. 5921 and 5744 differential expression genes were identified in pear and apple at three developmental stages respectively. We performed temporal and spatial differential gene expression profiling in developing fruits. Several key pathways such as signal transduction, photosynthesis, translation and many metabolisms were identified as involved in the differentiation of wild and cultivar fruits.

CONCLUSION

In this study, we reported on the next-generation sequencing study of the temporal and spatial mRNA expression profiling of pear and apple fruit trees. Also, we demonstrated that the integrated analysis of pear and apple transcriptome, which strongly revealed the consistent process of domestication in Rosaceae fruit trees. The results will be great influence to the improvement of cultivar species and the utilization of wild resources.

Effect of flutianil on the morphology and gene expression of powdery mildew

Flutianil, a fungicide effective only on powdery mildew, was previously reported to affect the host cell's haustorial formation and nutrient absorption. Studies were conducted to investigate flutianil's primary site of action on Blumeria graminis morphology using transmission electron microscope (TEM) observation and RNA sequencing (RAN-seq) techniques. TEM observation revealed that flutianil caused the extra-haustorial matrix and fungal cell wall to be obscured, without remarkable changes of other fungal organelles. RNA-seq analysis indicated that, unlike other powdery-mildew fungicides, flutianil did not significantly affect the constantly expressed genes for the survival of B. graminis. Genes whose expression is up- or downregulated by flutianil were found; these are the three sugar transporter genes and various effector genes, mainly expressed in haustoria. These findings indicate that the primary site of action of flutianil might be in the haustoria.

Combination of RNA-Seq transcriptomics and iTRAQ proteomics reveal the mechanism involved in fresh-cut yam yellowing

The aim of this study was to examine the regulation of transcriptomics and proteomics related to the yellowing of fresh-cut yams after storage. The comparison of yellow fresh-cut yam (YFY) vs. white fresh-cut yam (control) revealed 6894 upregulated and 6800 downregulated differentially expressed genes along with 1277 upregulated and 677 downregulated differentially expressed proteins. The results showed that the total carotenoids, flavonoids, and bisdemethoxycurcumin in YFY were higher than in the control due to the significant up-regulation of critical genes in the carotenoid biosynthesis pathway, flavonoid biosynthesis pathway, and stilbenoid, diarylheptanoid, and gingerol biosynthesis pathway. In addition, the tricarboxylic acid cycle and phenylpropanoid biosynthesis were both enhanced in YFY compared to the control, providing energy and precursors for the formation of yellow pigments. The results suggest that the synthesis of yellow pigments is regulated by critical genes, which might explain the yellowing of fresh-cut yam after storage.

Response of the ileum transcriptome to probiotic and fructo-oligosaccharides in Taiping chicken

Taiping chicken is indigenous chickens (Gallus gallus domesticus), which was one of China's excellent poultry species, is an excellent chicken in Gansu Province. As the problems caused by the overuse of antibiotics become more and more severe, people begin to look for ways to replace them. Among them, probiotics and fructo-oligosaccharides are the research hotspot to replace antibiotics. Probiotics and fructo-oligosaccharides can promote the absorption of nutrients, improve the ability to resist and prevent diseases, and improve the intestinal tissue morphology. In this study, we used RNA-Seq analysis to study the gene expression in ileum tissue after Taiping chicken was given probiotics and fructo-oligosaccharides. In total, 67 genes were differentially expressed in the ileum. Ten of the differently expressed genes were further validated by RT-qPCR. In addition, these differentially expressed genes were mainly enriched to tyrosine metabolism, AGE-RAGE signaling pathway in diabetic complications, phenylalanine metabolism, and pyrimidine metabolism. The results which this study provides contribute to our understanding application of probiotics and fructo-oligosaccharides in indigenous chickens production and provide a theoretical basis for the genetic development of indigenous chickens.

Response of the ileum transcriptome to fructo-oligosaccharides in Taiping chickens

The aim of this study was to investigate the effects of fructo-oligosaccharide (FOS) supplementation intake of Taiping chickens (Gallus gallus domesticus) and its stimulating effects on ileum. 120 healthy chickens were randomly divided into two groups; control group (CT) and fructo-oligosaccharides group (FOS). At the 60th day of age, ileum mucosa of three chickens per group were collected and performed transcriptome profiling of Taiping chicken ileum mucosa using the Hiseq™ 2500 sequencing platform. Compared with CT group, 50 genes were differentially expressed in the FOS group. Ten of the differently expressed genes were further validated by RT-qPCR. In addition, gene ontology and Kyoto encyclopedia of genes and genomes analyses revealed that these differentially expressed genes were mainly enriched to drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, retinol metabolism, fat digestion and absorption, herpes simplex infection and valine, leucine and isoleucine biosynthesis. The results of this study provided the help to our understanding application of fructo-oligosaccharides in indigenous chicken production and provide a theoretical basis for the genetic development of indigenous chickens.

RNA-Seq transcriptome analysis of ileum in Taiping chicken supplemented with the dietary probiotic

Taiping chicken is indigenous chickens (Gallus gallus domesticus), which was one of China's excellent poultry species. As the problems caused by the overuse of antibiotics become more and more concern, people begin to look for ways to replace them. Among them, probiotics and their preparations are the research hotspot to replace antibiotics. Probiotics can promote the absorption of nutrients, improve the ratio of meat to meat, resist and prevent diseases, and improve the intestinal tissue morphology. Here, we performed transcriptome profiling of Taiping chicken ileum which was given probiotics by using the Hiseq™ 2500 sequencing platform. A total of 18 genes were differentially expressed in the ileum under control group and probiotics group. Thirteen genes were upregulated with a range of fold change from 1.02 to 8.61, and 5 were downregulated with a range of fold change from - 1.06 to - 2.29. Ten of the differently expressed genes were further validated by qRT-PCR. In addition, Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that these differentially expressed genes were mainly enriched to bile secretion, alpha-linolenic acid metabolism, linoleic acid metabolism, ether lipid metabolism, and arachidonic acid metabolism. This study will help us to understand the role of probiotics in indigenous chicken production and provide theoretical basis for the genetic development of indigenous chickens.

Guidelines for Setting Up a mRNA Sequencing Experiment and Best Practices for Bioinformatic Data Analysis

RNA-sequencing, commonly referred to as RNA-seq, is the most recently developed method for the analysis of transcriptomes. It uses high-throughput next-generation sequencing technologies and has revolutionized our understanding of the complexity and dynamics of whole transcriptomes.In this chapter, we recall the key developments in transcriptome analysis and dissect the different steps of the general workflow that can be run by users to design and perform a mRNA-seq experiment as well as to process mRNA-seq data obtained by the Illumina technology. The chapter proposes guidelines for completing a mRNA-seq study properly and makes available recommendations for best practices based on recent literature and on the latest developments in technology and algorithms. We also remark the large number of choices available (especially for bioinformatic data analysis) in front of which the scientist may be in trouble.In the last part of the chapter we discuss the new frontiers of single-cell RNA-seq and isoform sequencing by long read technology.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef

RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through the identification of differentially expressed (DE) genes and pathways. A total of 80 steers with extreme phenotypes were selected from the University of Florida multibreed Angus-Brahman herd. The average slaughter age was 12.91±8.69 months. Tenderness, juiciness and connective tissue assessed by sensory panel, along with marbling, Warner-Bratzler Shear Force (WBSF) and cooking loss, were measured in longissimus dorsi muscle. Total RNA was extracted from muscle and one RNA-seq library per sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq-3000 platform to generate 2×101 bp paired-end reads. The overall read mapping rate using the Btau_4.6.1 reference genome was 63%. A total of 8,799 genes were analyzed using two different methodologies, an expression association and a DE analysis. A gene and exon expression association analysis was carried out using a meat quality index on all 80 samples as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). A gene and isoform DE evaluation was performed analyzing two groups with extreme WBSF, tenderness and marbling. A total of 676 (adjusted p-value≤0.05), 70 (adjusted p-value≤0.1) and 198 (adjusted p-value≤0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR≤0.1) were DE. Cytoskeletal and transmembrane anchoring genes and pathways were identified in the expression association, DE and the gene enrichment analyses; these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal interaction with myocyte and organelle membranes, contributing to cytoskeletal structure and architecture maintenance postmortem.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

Accuracy assessment of fusion transcript detection via read-mapping and de novo fusion transcript assembly-based methods

BACKGROUND

Accurate fusion transcript detection is essential for comprehensive characterization of cancer transcriptomes. Over the last decade, multiple bioinformatic tools have been developed to predict fusions from RNA-seq, based on either read mapping or de novo fusion transcript assembly.

RESULTS

We benchmark 23 different methods including applications we develop, STAR-Fusion and TrinityFusion, leveraging both simulated and real RNA-seq. Overall, STAR-Fusion, Arriba, and STAR-SEQR are the most accurate and fastest for fusion detection on cancer transcriptomes.

CONCLUSION

The lower accuracy of de novo assembly-based methods notwithstanding, they are useful for reconstructing fusion isoforms and tumor viruses, both of which are important in cancer research.

Transcriptomic profiling and analysis of differentially expressed genes in asparagus bean (Vigna unguiculata ssp. sesquipedalis) under salt stress

Asparagus bean (Vigna unguiculata ssp. sesquipedalis) is a warm season legume which is widely distributed over subtropical regions and semiarid areas. It is mainly grown as a significant protein source in developing countries. Salinity, as one of the main abiotic stress factors, constrains the normal growth and yield of asparagus bean. This study used two cultivars (a salt-sensitive genotype and a salt-tolerant genotype) under salt stress vs. control to identify salt-stress-induced genes in asparagus bean using RNA sequencing. A total of 692,086,838 high-quality clean reads, assigned to 121,138 unigenes, were obtained from control and salt-treated libraries. Then, 216 root-derived DEGs (differentially expressed genes) and 127 leaf-derived DEGs were identified under salt stress between the two cultivars. Of these DEGs, thirteen were assigned to six transcription factors (TFs), including AP2/EREBP, CCHC(Zn), C2H2, WRKY, WD40-like and LIM. GO analysis indicated four DEGs might take effects on the "oxidation reduction", "transport" and "signal transduction" process. Moreover, expression of nine randomly-chosen DEGs was verified by quantitative real-time-PCR (qRT-PCR) analysis. Predicted function of the nine tested DEGs was mainly involved in the KEGG pathway of cation transport, response to osmotic stress, and phosphorelay signal transduction system. A salt-stress-related pathway of "SNARE interactions in vesicular transport" was concerned. As byproducts, 15, 321 microsatellite markers were found in all the unigenes, and 17 SNP linked to six salt-stress induced DEGs were revealed. These candidate genes provide novel insights for understanding the salt tolerance mechanism of asparagus bean in the future.

De novo transcriptome assembly: A comprehensive cross-species comparison of short-read RNA-Seq assemblers

BACKGROUND

In recent years, massively parallel complementary DNA sequencing (RNA sequencing [RNA-Seq]) has emerged as a fast, cost-effective, and robust technology to study entire transcriptomes in various manners. In particular, for non-model organisms and in the absence of an appropriate reference genome, RNA-Seq is used to reconstruct the transcriptome de novo. Although the de novo transcriptome assembly of non-model organisms has been on the rise recently and new tools are frequently developing, there is still a knowledge gap about which assembly software should be used to build a comprehensive de novo assembly.

RESULTS

Here, we present a large-scale comparative study in which 10 de novo assembly tools are applied to 9 RNA-Seq data sets spanning different kingdoms of life. Overall, we built >200 single assemblies and evaluated their performance on a combination of 20 biological-based and reference-free metrics. Our study is accompanied by a comprehensive and extensible Electronic Supplement that summarizes all data sets, assembly execution instructions, and evaluation results. Trinity, SPAdes, and Trans-ABySS, followed by Bridger and SOAPdenovo-Trans, generally outperformed the other tools compared. Moreover, we observed species-specific differences in the performance of each assembler. No tool delivered the best results for all data sets.

CONCLUSIONS

We recommend a careful choice and normalization of evaluation metrics to select the best assembling results as a critical step in the reconstruction of a comprehensive de novo transcriptome assembly.

Reference genes selection for Calotropis procera under different salt stress conditions

Calotropis procera is a perennial Asian shrub with significant adaptation to adverse climate conditions and poor soils. Given its increased salt and drought stress tolerance, C. procera stands out as a powerful candidate to provide alternative genetic resources for biotechnological approaches. The qPCR (real-time quantitative polymerase chain reaction), widely recognized among the most accurate methods for quantifying gene expression, demands suitable reference genes (RGs) to avoid over- or underestimations of the relative expression and incorrect interpretation. This study aimed at evaluating the stability of ten RGs for normalization of gene expression of root and leaf of C. procera under different salt stress conditions and different collection times. The selected RGs were used on expression analysis of three target genes. Three independent experiments were carried out in greenhouse with young plants: i) Leaf₁₀₀ = leaf samples collected 30 min, 2 h, 8 h and 45 days after NaCl-stress (100 mM NaCl); ii) Root₅₀ and iii) Root₂₀₀ = root samples collected 30 min, 2 h, 8 h and 1day after NaCl-stress (50 and 200 mM NaCl, respectively). Stability rank among the three algorithms used showed high agreement for the four most stable RGs. The four most stable RGs showed high congruence among all combination of collection time, for each software studied, with minor disagreements. CYP23 was the best RG (rank of top four) for all experimental conditions (Leaf₁₀₀, Root_50, and Root₂₀₀). Using appropriated RGs, we validated the relative expression level of three differentially expressed target genes (NAC78, CNBL4, and ND1) in Leaf₁₀₀ and Root₂₀₀ samples. This study provides the first selection of stable reference genes for C. procera under salinity. Our results emphasize the need for caution when evaluating the stability RGs under different amplitude of variable factors.

Transcriptional Profile Corroborates that bml Mutant Plays likely Role in Premature Leaf Senescence of Rice (Oryza sativa L.)

Leaf senescence is the last period of leaf growth and a dynamic procedure associated with its death. The adaptability of the plants to changing environments occurs thanks to leaf senescence. Hence, transcriptional profiling is important to figure out the exact mechanisms of inducing leaf senescence in a particular crop, such as rice. From this perspective, leaf samples of two different rice genotypes, the brown midrib leaf (bml) mutant and its wild type (WT) were sampled for transcriptional profiling to identify differentially-expressed genes (DEGs). We identified 2670 DEGs, among which 1657 genes were up- and 1013 genes were down-regulated. These DEGs were enriched in binding and catalytic activity, followed by the single organism process and metabolic process through gene ontology (GO), while the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the DEGs were related to the plant hormone signal transduction and photosynthetic pathway enrichment. The expression pattern and the clustering of DEGs revealed that the WRKY and NAC family, as well as zinc finger transcription factors, had greater effects on early-senescence of leaf compared to other transcription factors. These findings will help to elucidate the precise functional role of bml rice mutant in the early-leaf senescence.

Transcriptome-based identification and validation of reference genes for plant-bacteria interaction studies using Nicotiana benthamiana

RT-qPCR is a widely used technique for the analysis of gene expression. Accurate estimation of transcript abundance relies strongly on a normalization that requires the use of reference genes that are stably expressed in the conditions analyzed. Initially, they were adopted from those used in Northern blot experiments, but an increasing number of publications highlight the need to find and validate alternative reference genes for the particular system under study. The development of high-throughput sequencing techniques has facilitated the identification of such stably expressed genes. Nicotiana benthamiana has been extensively used as a model in the plant research field. In spite of this, there is scarce information regarding suitable RT-qPCR reference genes for this species. Employing RNA-seq data previously generated from tomato plants, combined with newly generated data from N. benthamiana leaves infiltrated with Pseudomonas fluorescens, we identified and tested a set of 9 candidate reference genes. Using three different algorithms, we found that NbUbe35, NbNQO and NbErpA exhibit less variable gene expression in our pathosystem than previously used genes. Furthermore, the combined use of the first two is sufficient for robust gene expression analysis. We encourage employing these novel reference genes in future RT-qPCR experiments involving N. benthamiana and Pseudomonas spp.

Transcriptome Profiling and Characterization of Drought-Tolerant Potato Plant (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is the third most important food crop, and breeding drought-tolerant varieties is vital research goal. However, detailed molecular mechanisms in response to drought stress in potatoes are not well known. In this study, we developed EMS-mutagenized potatoes that showed significant tolerance to drought stress compared to the wild-type (WT) 'Desiree' cultivar. In addition, changes to transcripts as a result of drought stress in WT and drought-tolerant (DR) plants were investigated by de novo assembly using the Illumina platform. One-week-old WT and DR plants were treated with -1.8 Mpa polyethylene glycol-8000, and total RNA was prepared from plants harvested at 0, 6, 12, 24, and 48 h for subsequent RNA sequencing. In total, 61,100 transcripts and 5,118 differentially expressed genes (DEGs) displaying up- or down-regulation were identified in pairwise comparisons of WT and DR plants following drought conditions. Transcriptome profiling showed the number of DEGs with up-regulation and down-regulation at 909, 977, 1181, 1225 and 826 between WT and DR plants at 0, 6, 12, 24, and 48 h, respectively. Results of KEGG enrichment showed that the drought tolerance mechanism of the DR plant can mainly be explained by two aspects, the 'photosynthetic-antenna protein' and 'protein processing of the endoplasmic reticulum'. We also divided eight expression patterns in four pairwise comparisons of DR plants (DR0 vs DR6, DR12, DR24, DR48) under PEG treatment. Our comprehensive transcriptome data will further enhance our understanding of the mechanisms regulating drought tolerance in tetraploid potato cultivars.

A differential k-mer analysis pipeline for comparing RNA-Seq transcriptome and meta-transcriptome datasets without a reference

Next-generation DNA sequencing technologies, such as RNA-Seq, currently dominate genome-wide gene expression studies. A standard approach to analyse this data requires mapping sequence reads to a reference and counting the number of reads which map to each gene. However, for many transcriptome studies, a suitable reference genome is unavailable, especially for meta-transcriptome studies which assay gene expression from mixed populations of organisms. Where a reference is unavailable, it is possible to generate a reference by the de novo assembly of the sequence reads. However, the high cost of generating high-coverage data for de novo assembly hinders this approach and more importantly the accurate assembly of such data is challenging, especially for meta-transcriptome data, and resulting assemblies frequently suffer from collapsed regions or chimeric sequences. As an alternative to the standard reference mapping approach, we have developed a k-mer-based analysis pipeline (DiffKAP) to identify differentially expressed reads between RNA-Seq datasets without the requirement for a reference. We compared the DiffKAP approach with the traditional Tophat/Cuffdiff method using RNA-Seq data from soybean, which has a suitable reference genome. We subsequently examined differential gene expression for a coral meta-transcriptome where no reference is available, and validated the results using qRT-PCR. We conclude that DiffKAP is an accurate method to study differential gene expression in complex meta-transcriptomes without the requirement of a reference genome.

Transcriptome and Expression Profiling Analysis of Recalcitrant Tea (Camellia sinensis L.) Seeds Sensitive to Dehydration

The tea plant (Camellia sinensis (L.) O. Kuntze) is an economically important woody perennial nonalcoholic health beverage crop. Tea seeds are categorized as recalcitrant and are sensitive to dehydration treatment. However, the molecular basis of this phenomenon has not been investigated. Thus, we analyzed the genome-wide expression profiles of three dehydration stages using RNA-Seq and digital gene expression (DGE) technologies. We performed de novo assembly and obtained a total of 91,925 nonredundant unigenes, of which 58,472 were extensively annotated. By a hierarchical clustering of differentially expressed genes (DEGs), we found that 8929 DEGs were downregulated and 5875 DEGs were upregulated during dehydration treatment. A series of genes related to ABA biosynthesis and signal transduction, transcription factor, antioxidant enzyme, LEA protein, and proline metabolism that have been reported to function in dehydration process were found to be downregulated. Additionally, the expression profiles of 12 selected genes related to tea seed dehydration treatment were confirmed by qRT-PCR analysis. To our knowledge, this is the first genome-wide study elucidating the possible molecular mechanisms of sensitivity of recalcitrant tea seeds to dehydration. The results obtained in this study contribute to the preservation of tea seeds as genetic resources and can also be used to explore the mechanism of dehydration sensitivity of other recalcitrant seeds.

Identification and characterization of novel cecropins from the Oxysternon conspicillatum neotropic dung beetle

Dung beetles are exposed to a complex microbiological ecosystem during their life cycle. Characterization of novel host-defense peptides (HDP) is essential to understanding the host innate immune response in insects. It constitutes a promising alternative to look for new therapeutic agents against pathogenic microbes. We identified four new HDP, Oxysterlins 1, 2, 3, and 4 from the transcriptome of the Oxysternon conspicillatum dung beetle. These HDP display a highly conserved signal peptide and a mature peptide, characterized by an overall positive charge (cationic) (pI: 10.23–11.49), a hydrophobic ratio (ΦH: 35–41), and amphipathicity. Oxysterlins 1, 2, and 3 have a linear α-helix structure, whilst Oxysterlin 4 has a mixture of both α-helix and β-sheet structures without disulfide bonds through bioinformatics prediction and circular dichroism. Oxysterlins are part of the cecropin family group in an exclusive clade related to beetle cecropins. They have predominant antimicrobial activity against Gram-negative bacteria, including multidrug resistant strains (3.12–50 μg/mL) measured by plate microdilution. Their kinetics, in a time-killing curve showed concentration-dependent bactericidal activity. Furthermore, these HDP have low toxicity against human erythrocytes (62.5–500 μg/mL) and Vero cells (250–500 μg/mL). This article describes new HDP of the cecropin family from the Oxysternon conspicillatum dung beetle, with antimicrobial activity against multidrug resistant bacteria and low toxicity.

Complementary RNA-Sequencing Based Transcriptomics and iTRAQ Proteomics Reveal the Mechanism of the Alleviation of Quinclorac Stress by Salicylic Acid in Oryza sativa ssp. japonica

To uncover the alleviation mechanism of quinclorac stress by salicylic acid (SA), leaf samples of Oryza sativa ssp. Japonica under quinclorac stress with and without SA pre-treatment were analyzed for transcriptional and proteomic profiling to determine the differentially expressed genes (DEGs) and proteins (DEPs), respectively. Results showed that quinclorac stress altered the expression of 2207 DEGs (1427 up-regulated, 780 down-regulated) and 147 DEPs (98 down-regulated, 49 up-regulated). These genes and proteins were enriched in glutathione (GSH) metabolism, porphyrin and chlorophyll metabolism, the biosynthesis of secondary metabolites, glyoxylate and dicarboxylate metabolism, and so on. It also influenced apetala2- ethylene-responsive element binding protein (AP2-EREBP) family, myeloblastosis (MYB) family and WRKY family transcription factors. After SA pre-treatment, 697 genes and 124 proteins were differentially expressed. Pathway analysis showed similar enrichments in GSH, glyoxylate and dicarboxylate metabolism. Transcription factors were distributed in basic helix-loop-helix (bHLH), MYB, Tify and WRKY families. Quantitative real-time PCR results revealed that quinclorac stress induced the expression of glutathion reductase (GR) genes (OsGR2, OsGR3), which was further pronounced by SA pre-treatment. Quinclorac stress further mediated the accumulation of acetaldehyde in rice, while SA enhanced the expression of OsALDH2B5 and OsALDH7 to accelerate the metabolism of herbicide quinclorac for the protection of rice. Correlation analysis between transcriptome and proteomics demonstrated that, under quinclorac stress, correlated proteins/genes were mainly involved in the inhibition of intermediate steps in the biosynthesis of chlorophyll. Other interesting proteins/genes and pathways regulated by herbicide quinclorac and modulated by SA pre-treatment were also discussed, based on the transcriptome and proteomics results.