Search

Chaetomium globosum Kunze, an internationally recognized biocontrol fungus. It mycoparasitizes various plant pathogens and produce antifungal metabolites to suppress the growth of pathogenic fungi. Lack of detailed genome level diversity studies has delimited the development and utilization of potential C. globosum strains. The present study was taken to reveal the distribution, identification, and characterization of expressed sequence tag-simple sequence repeats (EST-SSRs) in C. globosum. RNA-Seq experiment was performed for C. globosum potential isolate Cg2 (AY429049) using Illumina HiSeq 2500. Reference-guided de novo assembly yielded 45,582 transcripts containing 27,957 unigenes. We generated a new set of 8485 EST-SSR markers distributed in 5908 unigene sequences with one SSR locus distribution density per 6.1 kb. Six distinct classes of SSR repeat motifs were identified. The most abundant were mononucleotide repeats (51.67%), followed by tri-nucleotides (36.61%). Out of 5034 EST-SSR primers, 50 primer pairs were selected and validated for the polymorphic study of 15 C. globosum isolates. Twenty-two SSR markers showed average genetic polymorphism among C. globosum isolates. The number of alleles (Na) per marker ranges from 2 to 4, with a total of 74 alleles detected for 22 markers with a mean polymorphism information content (PIC) value of 0.4. UPGMA hierarchical clustering analysis generated three main clusters of C. globosum isolates and exhibited a lower similarity index range from 0.59 to 0.85. Thus, the newly developed EST-SSR markers could replace traditional methods for determining diversity. The study will also enhance the genomic research in C. globosum to explore its biocontrol potential against phytopathogens. [ABSTRACT FROM AUTHOR]

Liriodendron chinense is a very popular tree used in landscaping because of its beautiful flowers and unique leaf shape. Although the L. chinense genome was released in 2018 based on PacBio long reads and genetic linkage map and provides a valuable resource for gene function discovery and transcriptome analyses of L. chinense. Yet, that version of the L. chinense genome relied on ab initio annotation and a small amount of transcriptome data, and lacks 5′ UTRs and 3′ UTRs and alternative splicing (AS) annotation; hence, it is imperative that an improved annotation be sought. Herein, we re-annotated the structure and function of genes across the entire L. chinense nuclear genome based on PacBio long reads and Illumina short reads retrieved from near-round organization types. This updated annotation of the L. chinense genome, Lchi2.1, includes a total of 42,831 gene models with a high proportion (92.8%) of complete BUSCOs. Among the Lchi2.1 annotation, 21,324 genes were modified, 15,770 novel genes were added, 14,851 genes were augmented with information on the 5′ UTRs and 3′ UTRs, and 6751 genes were found as alternatively spliced isoforms at the whole-genome scale. Additionally, we re-analyzed the transcriptome data across various stages of leaf shape development for L. chinense, founding many genes enriched to post-embryonic plant morphogenesis, plant organ formation, and anatomical structure formation involved in morphogenesis, which suggests these genes could be related to leaf morphological processes. Overall, the updated annotation of the L. chinense genome presented in this study will contribute to gene functional studies in L. chinense and its comparative genomic analysis with related plant species. [ABSTRACT FROM AUTHOR]

The impediment of inadequate shrimp gonad development has adversely affected their production in the aquaculture sector. One promising strategy to enhance the quality and timing of offspring for commercial purposes is the regulation of their reproductive activity. Nevertheless, the availability of information pertaining to ovarian development in Fenneropenaeus penicillatus remains limited. To address this knowledge gap, our study employed a multi-faceted approach, integrating cutting-edge technologies such as Iso-seq, RNA-seq, RT-qPCR, histological and morphological observations. The Iso-seq generating 81,964,598 subreads with an average length of 1600 bp. This effort yielded 926,713 CCS reads and 557,043 FLNC reads, culminating in a robust set of 16,587 non-redundant full-length transcripts. This study utilized CPC, CPAT, PLEK, and pfam, to predict the numbers of long non-coding RNAs (lncRNAs). The outcomes were as follows: 5204, 4478, 7055, and 8175 lncRNAs, respectively. Our findings unveiled the intricate molecular underpinnings of ovarian development, encompassing pivotal pathways such as oocyte meiosis, insulin signaling pathway, oxytocin signaling pathway, progesterone-mediated oocyte maturation, estrogen signaling pathway, and GnRH signaling pathway. Notably, we identified specific genes—including Vtg, VgR, Sxl4, UCH, Hsp70, PP1, CDK2, Rheb, MKK3, and CDC42—that exhibited dynamic expression patterns during ovarian development. This study has produced a comprehensive dataset of full-length transcripts that holds significant potential as a genomic reference for future research on F. penicillatus. Furthermore, these insights may substantially contribute to the advancement of artificial breeding techniques and the ongoing development of sustainable aquaculture practices for F. penicillatus.

Simple Summary: Batocera horsfieldi (Hope) is an important wood-boring pest in China. This pest primarily infests tree trunks by feeding on the woody tissue, creating a network of interconnected tunnels within the trunk. These tunnels become blocked with insect feces and wood debris, causing damage, decay, and even the death of the host plant's tissues. In this study, single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies were employed to conduct full-length transcriptome sequencing of male and female adults of B. horsfieldi. A total of 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. Clustering and redundancy removal of the full-length non-chimera reads resulted in 39,912 consensus reads. Additionally, functional annotation was performed on a total of 84,650 transcripts in seven different databases. This study provides an important foundation for future exploration of gene regulation in the interaction between B. horsfieldi and host plants using RNA interference (RNAi), and it offers a scientific basis for the prevention and control of B. horsfieldi. Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi. [ABSTRACT FROM AUTHOR]

Thuja koraiensis Nakai (Cupressaceae) is an endangered and ecologically important conifer endemic to Korea. Previously, we established a protocol for micropropagation in T. koraiensis, which involved somatic embryogenesis from embryogenic callus of T. koraiensis. However, the molecular mechanisms underlying somatic embryogenesis remain unclear. Herein, we performed transcriptomic analysis to identify somatic embryogenesis-related genes of T. koraiensis via Illumina RNA sequencing. We conducted de novo transcriptome assembly using a Trinity assembler, which produced 274,077 transcript contigs clustered into 205,843 transcripts (unigenes), with an average length of 825 base pairs. Of all the unigenes, 14.69%, 18.62%, and 7.4% had homologs in the Gene Ontology, NCBI Non-redundant Protein, and NCBI Nucleotide databases, respectively. Among these mRNA sequences, expression of putative embryogenesis-associated transcription factors, namely BABYBOOM (BBM), WUSCHEL-RELATED HOMEOBOX (WOX), and SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK), was analyzed during somatic embryogenesis. RT-PCR analysis revealed that TkBBM, TkWOX, and TkSERK were highly expressed in embryogenic callus and seedling roots, whereas their expression was reduced in seedling leaves. Our findings provide new insights into the roles of BBM, WOX, and SERK in somatic embryogenesis. Our results may serve as a reference for comparative transcriptome analysis in related species and further aid functional genomics studies in conifers. [ABSTRACT FROM AUTHOR]

This work aimed to study the plant conditioning effect and mode of action of a plant-based biostimulant used in organic farming. This new generation plant biostimulant, named ELICE16INDURES®, is rich in plant bio-active ingredients containing eleven supercritical botanical extracts encapsulated in nano-scale liposomes. The dose–response (10 to 240 g ha−1) relationship was tested in a field population of autumn barley (Hordeum vulgare) test crop, and underlying molecular mechanisms were studied. Applying nanotechnology, cell-identical nanoparticles may help the better uptake and delivery of active ingredients increasing resilience, vitality, and crop yield. The amount of harvested crops showed a significant increase of 27.5% and 39.9% interconnected to higher normalized difference vegetation index (NDVI) of 20% and 25% after the treatment of low and high dosages (20 and 240 g ha−1), respectively. Illumina NextSeq 550 sequencing, gene expression profiling, and KEGG-pathway analysis of outstanding dosages indicated the upregulation of pathogenesis-related (PR) and other genes—associated with induced resistance—which showed dose dependency as well. [ABSTRACT FROM AUTHOR]

Climate change and the accelerated rate of population growth are imposing a progressive degradation of natural ecosystems worldwide. In this context, the use of pioneer trees represents a powerful approach to reverse the situation. Among others, N2-fixing actinorhizal trees constitute important elements of plant communities and have been successfully used in land reclamation at a global scale. In this study, we have analyzed the transcriptome of the photosynthetic organs of Casuarina glauca (branchlets) to unravel the molecular mechanisms underlying salt stress tolerance. For that, C. glauca plants supplied either with chemical nitrogen (KNO3+) or nodulated by Frankia (NOD+) were exposed to a gradient of salt concentrations (200, 400, and 600 mM NaCl) and RNA-Seq was performed. An average of ca. 25 million clean reads was obtained for each group of plants, corresponding to 86,202 unigenes. The patterns of differentially expressed genes (DEGs) clearly separate two groups: (i) control- and 200 mM NaCl-treated plants, and (ii) 400 and 600 mM NaCl-treated plants. Additionally, although the number of total transcripts was relatively high in both plant groups, the percentage of significant DEGs was very low, ranging from 6 (200 mM NaCl/NOD+) to 314 (600 mM NaCl/KNO3+), mostly involving down-regulation. The vast majority of up-regulated genes was related to regulatory processes, reinforcing the hypothesis that some ecotypes of C. glauca have a strong stress-responsive system with an extensive set of constitutive defense mechanisms, complemented by a tight mechanism of transcriptional and post-transcriptional regulation. The results suggest that the robustness of the stress response system in C. glauca is regulated by a limited number of genes that tightly regulate detoxification and protein/enzyme stability, highlighting the complexity of the molecular interactions leading to salinity tolerance in this species. [ABSTRACT FROM AUTHOR]

It has been shown that the best coverage of the HepG2 cell line transcriptome encoded by genes of a single chromosome, chromosome 18, is achieved by a combination of two sequencing platforms, Illumina RNA-Seq and Oxford Nanopore Technologies (ONT), using cut-off levels of FPKM > 0 and TPM > 0, respectively. In this study, we investigated the extent to which the combination of these transcriptomic analysis methods makes it possible to achieve a high coverage of the transcriptome encoded by the genes of other human chromosomes. A comparative analysis of transcriptome coverage for various types of biological material was carried out, and the HepG2 cell line transcriptome was compared with the transcriptome of liver tissue cells. In addition, the contribution of variability in the coverage of expressed genes in human transcriptomes to the creation of a draft human transcriptome was evaluated. For human liver tissues, ONT makes an extremely insignificant contribution to the overall coverage of the transcriptome. Thus, to ensure maximum coverage of the liver tissue transcriptome, it is sufficient to apply only one technology: Illumina RNA-Seq (FPKM > 0).

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Abstract Background Melatonin, a multifunctional signal molecule, has been reported to play crucial roles in growth and development and stress responses in various plant species. Okra (Abelmoschus esculentus L.) is a food crop with extremely high values of nutrition and healthcare. Recent reports have revealed the protective role of melatonin in alleviating salt stress. However, little is known about its regulatory mechanisms in response to salt stress in okra. Results In this study, we explored whether exogenous melatonin pretreatment could alleviate salt stress (300 mM NaCl) of okra plants. Results showed that exogenous application of melatonin (50 μM) significantly enhanced plant tolerance to salt stress, as demonstrated by the plant resistant phenotype, as well as by the higher levels of the net photosynthetic rate, chlorophyll fluorescence and chlorophyll content in comparison with nontreated salt-stressed plants. Additionally, melatonin pretreatment remarkably decreased the levels of lipid peroxidation and H2O2 content and scavenged O2 •- in melatonin-pretreated plants, which may be attributed to the higher levels of enzyme activities including POD and GR. Moreover, a combination of third- (PacBio) and second-generation (Illumina) sequencing technologies was applied to sequence full-length transcriptomes of okra. A total of 121,360 unigenes was obtained, and the size of transcript lengths ranged from 500 to 6000 bp. Illumina RNA-seq analysis showed that: Comparing with control, 1776, 1063 and 1074 differential expression genes (DEGs) were identified from the three treatments (NaCl, MT50 and MT + NaCl, respectively). These genes were enriched in more than 10 GO terms and 34 KEGG pathways. Nitrogen metabolism, sulfur metabolism, and alanine, aspartate and glutamate metabolism were significantly enriched in all three treatments. Many transcription factors including MYB, WRKY, NAC etc., were also identified as DEGs. Conclusions Our preliminary results suggested that melatonin pretreatment enhanced salt tolerance of okra plants for the first time. These data provide the first set of full-length isoforms in okra and more comprehensive insights into the molecular mechanism of melatonin responses to salt stress.

Thuja koraiensis Nakai (Cupressaceae) is an endangered and ecologically important conifer endemic to Korea. Previously, we established a protocol for micropropagation in T. koraiensis, which involved somatic embryogenesis from embryogenic callus of T. koraiensis. However, the molecular mechanisms underlying somatic embryogenesis remain unclear. Herein, we performed transcriptomic analysis to identify somatic embryogenesis-related genes of T. koraiensis via Illumina RNA sequencing. We conducted de novo transcriptome assembly using a Trinity assembler, which produced 274,077 transcript contigs clustered into 205,843 transcripts (unigenes), with an average length of 825 base pairs. Of all the unigenes, 14.69%, 18.62%, and 7.4% had homologs in the Gene Ontology, NCBI Non-redundant Protein, and NCBI Nucleotide databases, respectively. Among these mRNA sequences, expression of putative embryogenesis-associated transcription factors, namely BABYBOOM (BBM), WUSCHEL-RELATED HOMEOBOX (WOX), and SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK), was analyzed during somatic embryogenesis. RT-PCR analysis revealed that TkBBM, TkWOX, and TkSERK were highly expressed in embryogenic callus and seedling roots, whereas their expression was reduced in seedling leaves. Our findings provide new insights into the roles of BBM, WOX, and SERK in somatic embryogenesis. Our results may serve as a reference for comparative transcriptome analysis in related species and further aid functional genomics studies in conifers.

Climate change and the accelerated rate of population growth are imposing a progressive degradation of natural ecosystems worldwide. In this context, the use of pioneer trees represents a powerful approach to reverse the situation. Among others, N2-fixing actinorhizal trees constitute important elements of plant communities and have been successfully used in land reclamation at a global scale. In this study, we have analyzed the transcriptome of the photosynthetic organs of Casuarina glauca (branchlets) to unravel the molecular mechanisms underlying salt stress tolerance. For that, C. glauca plants supplied either with chemical nitrogen (KNO3+) or nodulated by Frankia (NOD+) were exposed to a gradient of salt concentrations (200, 400, and 600 mM NaCl) and RNA-Seq was performed. An average of ca. 25 million clean reads was obtained for each group of plants, corresponding to 86,202 unigenes. The patterns of differentially expressed genes (DEGs) clearly separate two groups: (i) control- and 200 mM NaCl-treated plants, and (ii) 400 and 600 mM NaCl-treated plants. Additionally, although the number of total transcripts was relatively high in both plant groups, the percentage of significant DEGs was very low, ranging from 6 (200 mM NaCl/NOD+) to 314 (600 mM NaCl/KNO3+), mostly involving down-regulation. The vast majority of up-regulated genes was related to regulatory processes, reinforcing the hypothesis that some ecotypes of C. glauca have a strong stress-responsive system with an extensive set of constitutive defense mechanisms, complemented by a tight mechanism of transcriptional and post-transcriptional regulation. The results suggest that the robustness of the stress response system in C. glauca is regulated by a limited number of genes that tightly regulate detoxification and protein/enzyme stability, highlighting the complexity of the molecular interactions leading to salinity tolerance in this species.

This work aimed to study the plant conditioning effect and mode of action of a plant-based biostimulant used in organic farming. This new generation plant biostimulant, named ELICE16INDURES®, is rich in plant bio-active ingredients containing eleven supercritical botanical extracts encapsulated in nano-scale liposomes. The dose–response (10 to 240 g ha−1) relationship was tested in a field population of autumn barley (Hordeum vulgare) test crop, and underlying molecular mechanisms were studied. Applying nanotechnology, cell-identical nanoparticles may help the better uptake and delivery of active ingredients increasing resilience, vitality, and crop yield. The amount of harvested crops showed a significant increase of 27.5% and 39.9% interconnected to higher normalized difference vegetation index (NDVI) of 20% and 25% after the treatment of low and high dosages (20 and 240 g ha−1), respectively. Illumina NextSeq 550 sequencing, gene expression profiling, and KEGG-pathway analysis of outstanding dosages indicated the upregulation of pathogenesis-related (PR) and other genes—associated with induced resistance—which showed dose dependency as well.

Abstract Background Monogenean flatworms are the main fish ectoparasites inflicting serious economic losses in aquaculture. The polyopisthocotylean Sparicotyle chrysophrii parasitizes the gills of gilthead sea bream (GSB, Sparus aurata) causing anaemia, lamellae fusion and sloughing of epithelial cells, with the consequent hypoxia, emaciation, lethargy and mortality. Currently no preventive or curative measures against this disease exist and therefore information on the host-parasite interaction is crucial to find mitigation solutions for sparicotylosis. The knowledge about gene regulation in monogenean-host models mostly comes from freshwater monopysthocotyleans and almost nothing is known about polyopisthocotyleans. The current study aims to decipher the host response at local (gills) and systemic (spleen, liver) levels in farmed GSB with a mild natural S. chrysophrii infection by transcriptomic analysis. Results Using Illumina RNA sequencing and transcriptomic analysis, a total of 2581 differentially expressed transcripts were identified in infected fish when compared to uninfected controls. Gill tissues in contact with the parasite (P gills) displayed regulation of fewer genes (700) than gill portions not in contact with the parasite (NP gills) (1235), most likely due to a local silencing effect of the parasite. The systemic reaction in the spleen was much higher than that at the parasite attachment site (local) (1240), and higher than in liver (334). NP gills displayed a strong enrichment of genes mainly related to immune response and apoptosis. Processes such as apoptosis, inflammation and cell proliferation dominated gills, whereas inhibition of apoptosis, autophagy, platelet activation, signalling and aggregation, and inflammasome were observed in spleen. Proteasome markers were increased in all tissues, whereas hypoxia-related genes were down-regulated in gills and spleen. Conclusions Contrasting forces seem to be acting at local and systemic levels. The splenic down-regulation could be part of a hypometabolic response, to counteract the hypoxia induced by the parasite damage to the gills and to concentrate the energy on defence and repair responses. Alternatively, it can be also interpreted as the often observed action of helminths to modify host immunity in its own interest. These results provide the first toolkit for future studies towards understanding and management of this parasitosis.

Salix cupularis, an afforestation tree for soil and water conservation, is a deciduous willow species commonly cultivated in Sichuan Province, China. The research on the molecular mechanism of drought resistance of S. cupularis will contribute to the subsequent study of soil and water conservation, biodiversity, and ecological protection. Here, the results of an RNA-Seq analysis of transcriptomic data using Illumina deep sequencing technology to compare control (T1) and drought-stressed (T2) seedlings of S. cupularis were reported. Polyethylene glycol was used to induce drought stress in S. cupularis seedlings; 15-day-old seedlings were used for RNA-Seq analysis. De novo gene assembly was used to generate the comprehensive transcriptome, which was comprised of 65,228 unigenes, among which 16,192 were> 1 kb and accounted for 24.82% of the total unigene library. Analysis of the gene structure based on the unigene database, including open reading frame (ORF) predictions, simple sequence repeat (SSR) analysis, and inter-sample single nucleotide polymorphism (SNP) analysis, resulted in the identification of 758,779 SNP and 6779 SSR. After comparing the T1 and T2 gene profiles, 4289 differentially expressed genes (DEGs) were identified, of which 2340 were upregulated and 1949 were downregulated. The gene ontology (GO) analysis classified genes into the cellular component, molecular function, and biological process categories. The clusters of orthologous genes (COG) database was constructed based on the evolutionary phylogeny of bacteria, algae, and eukaryotes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs between T1 and T2 revealed that the top 20 enriched KEGG metabolic pathways were primarily in porphyrin and chlorophyll metabolism, photosynthesis, phenylpropanoid biosynthesis of phenylalanine metabolism, pentose and glucoronate interconversions, and glycan degradation, among other metabolic pathways. Collectively, these findings provide a theoretical basis for future studies on the genetic resistance of S. cupularis, and it also lay a foundation for the flow-up research of soil and water conservation, biodiversity, and ecological protection.

Background: Melatonin, a multifunctional signal molecule, has been reported to play crucial roles in growth and development and stress responses in various plant species. Okra (Abelmoschus esculentus L.) is a food crop with extremely high values of nutrition and healthcare. Recent reports have revealed the protective role of melatonin in alleviating salt stress. However, little is known about its regulatory mechanisms in response to salt stress in okra. Results: In this study, we explored whether exogenous melatonin pretreatment could alleviate salt stress (300 mM NaCl) of okra plants. Results showed that exogenous application of melatonin (50 μM) significantly enhanced plant tolerance to salt stress, as demonstrated by the plant resistant phenotype, as well as by the higher levels of the net photosynthetic rate, chlorophyll fluorescence and chlorophyll content in comparison with nontreated salt-stressed plants. Additionally, melatonin pretreatment remarkably decreased the levels of lipid peroxidation and H2O2 content and scavenged O2•- in melatonin-pretreated plants, which may be attributed to the higher levels of enzyme activities including POD and GR. Moreover, a combination of third- (PacBio) and second-generation (Illumina) sequencing technologies was applied to sequence full-length transcriptomes of okra. A total of 121,360 unigenes was obtained, and the size of transcript lengths ranged from 500 to 6000 bp. Illumina RNA-seq analysis showed that: Comparing with control, 1776, 1063 and 1074 differential expression genes (DEGs) were identified from the three treatments (NaCl, MT50 and MT + NaCl, respectively). These genes were enriched in more than 10 GO terms and 34 KEGG pathways. Nitrogen metabolism, sulfur metabolism, and alanine, aspartate and glutamate metabolism were significantly enriched in all three treatments. Many transcription factors including MYB, WRKY, NAC etc., were also identified as DEGs. Conclusions: Our preliminary results suggested that melatonin pretreatment enhanced salt tolerance of okra plants for the first time. These data provide the first set of full-length isoforms in okra and more comprehensive insights into the molecular mechanism of melatonin responses to salt stress. [ABSTRACT FROM AUTHOR]

Background Red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most destructive insects for palm trees in the world. However, its genome resources are still in the blank stage, which limits the study of molecular and growth development analysis. Methods In this study, we used PacBio Iso-Seq and Illumina RNA-seq to first generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th larva, female and male) to increase our understanding of the life cycle and molecular characteristics of R. ferrugineus. Results A total of 63,801 nonredundant full-length transcripts were generated with an average length of 2,964 bp from three developmental stages, including the 7th instar larva, pupa, female adult and male adult. These transcripts showed a high annotation rate in seven public databases, with 54,999 (86.20%) successfully annotated. Meanwhile, 2,184 alternative splicing (AS) events, 2,084 transcription factors (TFs), 66,230 simple sequence repeats (SSR) and 9,618 Long noncoding RNAs (lncRNAs) were identified. In summary, our results provide a new source of full-length transcriptional data and information for the further study of gene expression and genetics in R. ferrugineus.

The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens.

Multidrug-Resistant (MDR) and Extensively Drug Resistant (XDR) Acinetobacter baumannii (Ab) represent a serious cause of healthcare-associated infections worldwide. Currently, the available treatment options are very restricted and colistin-based therapies are last-line treatments of these infections, even though colistin resistant (COLR) Ab have rarely been isolated yet. In bacteria, small non-coding RNAs (sRNAs) have been implicated in regulatory pathways of different biological functions, however, no knowledge exists about the sRNA role on the biological adaptation in COLRAb. Our study investigated two Italian XDR isogenic colistin-susceptible/resistant (COLS/R) Ab strain-pairs to discover new sRNA signatures. Comparative sRNA transcriptome (sRNAome) analyses were carried out by Illumina RNA-seq using both a Tru-Seq and a Short Insert library, whilst Ab ATCC 17978 and ACICU Reference Genome assembly, mapping, annotation and statistically significant differential expression (q-value ≤ 0.01) of the raw reads were performed by the Rockhopper tool. A computational filtering, sorting only similarly statistically significant differentially expressed (DE) sRNAs mapping on the same gene in both COLRAb isolates was conducted. COLR vs. COLS sRNAome, analyzed integrating the DE sRNAs obtained from the two different libraries, revealed some statistically significant DE sRNAs in COLRAb. In detail, we found: (i) two different under-expressed cis-acting sRNAs (AbsRNA1 and AbsRNA2) mapping in antisense orientation the 16S rRNA gene A1S_r01, (ii) one under-expressed cis-acting sRNA (AbsRNA3) targeting the A1S_2505 gene (hypothetical protein), (iii) one under-expressed microRNA-size small RNA fragment (AbsRNA4) and its pre-microAbsRNA4 targeting the A1S_0501 gene (hypothetical protein), (iv) as well as an over-expressed microRNA-size small RNA fragment (AbsRNA5) and its pre-microAbsRNA5 targeting the A1S_3097 gene (signal peptide). Custom TaqMan® probe-based real-time qPCRs validated the expression pattern of the selected sRNA candidates shown by RNA-seq. Furthermore, analysis on sRNA ΔA1S_r01, ΔA1S_2505 as well as the over-expressed A1S_3097 mutants revealed no effects on colistin resistance. Our study, for the first time, found the sRNAome signatures of clinical COLRAb with a computational prediction of their targets related to protein synthesis, host-microbe interaction and other different biological functions, including biofilm production, cell-cycle control, virulence, and antibiotic-resistance.

Multidrug-Resistant (MDR) and Extensively Drug Resistant (XDR) Acinetobacter baumannii (Ab) represent a serious cause of healthcare-associated infections worldwide. Currently, the available treatment options are very restricted and colistin-based therapies are last-line treatments of these infections, even though colistin resistant (COLR) Ab have rarely been isolated yet. In bacteria, small non-coding RNAs (sRNAs) have been implicated in regulatory pathways of different biological functions, however, no knowledge exists about the sRNA role on the biological adaptation in COLR Ab. Our study investigated two Italian XDR isogenic colistin-susceptible/resistant (COLS/R) Ab strain-pairs to discover new sRNA signatures. Comparative sRNA transcriptome (sRNAome) analyses were carried out by Illumina RNA-seq using both a Tru-Seq and a Short Insert library, whilst Ab ATCC 17978 and ACICU Reference Genome assembly, mapping, annotation and statistically significant differential expression (q -value ≤ 0.01) of the raw reads were performed by the Rockhopper tool. A computational filtering, sorting only similarly statistically significant differentially expressed (DE) sRNAs mapping on the same gene in both COLR Ab isolates was conducted. COLR vs. COLS sRNAome, analyzed integrating the DE sRNAs obtained from the two different libraries, revealed some statistically significant DE sRNAs in COLR Ab. In detail, we found: (i) two different under-expressed cis -acting sRNAs (Ab sRNA1 and Ab sRNA2) mapping in antisense orientation the 16S rRNA gene A1S_r01, (ii) one under-expressed cis -acting sRNA (Ab sRNA3) targeting the A1S_2505 gene (hypothetical protein), (iii) one under-expressed microRNA-size small RNA fragment (Ab sRNA4) and its pre-micro Ab sRNA4 targeting the A1S_0501 gene (hypothetical protein), (iv) as well as an over-expressed microRNA-size small RNA fragment (Ab sRNA5) and its pre-micro Ab sRNA5 targeting the A1S_3097 gene (signal peptide). Custom TaqMan® probe-based real-time qPCRs validated the expression pattern of the selected sRNA candidates shown by RNA-seq. Furthermore, analysis on sRNA ΔA1S_r01, ΔA1S_2505 as well as the over-expressed A1S_3097 mutants revealed no effects on colistin resistance. Our study, for the first time, found the sRNAome signatures of clinical COLR Ab with a computational prediction of their targets related to protein synthesis, host-microbe interaction and other different biological functions, including biofilm production, cell-cycle control, virulence, and antibiotic-resistance. [ABSTRACT FROM AUTHOR]

Alternative splicing (AS) plays a critical role in regulating different physiological and developmental processes in eukaryotes, by dramatically increasing the diversity of the transcriptome and the proteome. However, the saturation and complexity of AS remain unclear in lotus due to its limitation of rare obtainment of full-length multiple-splice isoforms. In this study, we apply a hybrid assembly strategy by combining single-molecule real-time sequencing and Illumina RNA-seq to get a comprehensive insight into the lotus transcriptomic landscape. We identified 211,802 high-quality full-length non-chimeric reads, with 192,690 non-redundant isoforms, and updated the lotus reference gene model. Moreover, our analysis identified a total of 104,288 AS events from 16,543 genes, with alternative 3ʹ splice-site being the predominant model, following by intron retention. By exploring tissue datasets, 370 tissue-specific AS events were identified among 12 tissues. Both the tissue-specific genes and isoforms might play important roles in tissue or organ development, and are suitable for 'ABCE' model partly in floral tissues. A large number of AS events and isoform variants identified in our study enhance the understanding of transcriptional diversity in lotus, and provide valuable resource for further functional genomic studies. [ABSTRACT FROM AUTHOR]

Seagrasses as Posidonia oceanica reproduce mostly by vegetative propagation, which can reduce genetic variability within populations. Since, in clonally propagated species, insurgence of genetic variability can be determined by the activity of transposable elements, we have estimated the activity of such repeat elements by measuring their expression level in the leaves of plants from a Mediterranean site, for which Illumina complementary DNA (cDNA) sequence reads (produced from RNAs isolated by leaves of plants from deep and shallow meadows) were publicly available. Firstly, we produced a collection of retrotransposon-related sequences and then mapped Illumina cDNA reads onto these sequences. With this approach, it was evident that Posidonia retrotransposons are, in general, barely expressed; only nine elements resulted transcribed at levels comparable with those of reference genes encoding tubulins and actins. Differences in transcript abundance were observed according to the superfamily and the lineage to which the retrotransposons belonged. Only small differences were observed between retrotransposon expression levels in leaves of shallow and deep Posidonia meadow stands, whereas one TAR/Tork element resulted differentially expressed in deep plants exposed to heat. It can be concluded that, in P. oceanica, the contribution of retrotransposon activity to genetic variability is reduced, although the nine specific active elements could actually produce new structural variations.

Spotted scat (Scatophagus argus) is an economically important farmed fish, particularly in East and Southeast Asia. Because there has been little research on reproductive development and regulation in this species, the lack of a mature artificial reproduction technology remains a barrier for the sustainable development of the aquaculture industry. More genetic and genomic background knowledge is urgently needed for an in-depth understanding of the molecular mechanism of reproductive process and identification of functional genes related to sexual differentiation, gonad maturation and gametogenesis. For these reasons, we performed transcriptomic analysis on spotted scat using a multiple tissue sample mixing strategy. The Illumina RNA sequencing generated 118 510 486 raw reads. After trimming, de novo assembly was performed and yielded 99 888 unigenes with an average length of 905.75 bp. A total of 45 015 unigenes were successfully annotated to the Nr, Swiss-Prot, KOG and KEGG databases. Additionally, 23 783 and 27 183 annotated unigenes were assigned to 56 Gene Ontology (GO) functional groups and 228 KEGG pathways, respectively. Subsequently, 2 474 transcripts associated with reproduction were selected using GO term and KEGG pathway assignments, and a number of reproduction-related genes involved in sex differentiation, gonad development and gametogenesis were identified. Furthermore, 22 279 simple sequence repeat (SSR) loci were discovered and characterized. The comprehensive transcript dataset described here greatly increases the genetic information available for spotted scat and contributes valuable sequence resources for functional gene mining and analysis. Candidate transcripts involved in reproduction would make good starting points for future studies on reproductive mechanisms, and the putative sex differentiation-related genes will be helpful for sex-determining gene identification and sex-specific marker isolation. Lastly, the SSRs can serve as marker resources for future research into genetics, marker-assisted selection (MAS) and conservation biology. [ABSTRACT FROM AUTHOR]

Background: Quantitative trait loci (QTLs) detected in one mapping population may not be detected in other mapping populations at all the time. Therefore, before being used for marker assisted breeding, QTLs need to be validated in different environments and/or genetic backgrounds to rule out statistical anomalies. In this regard, we mapped the QTLs controlling various agronomic traits in a recombinant inbred line (RIL) population in response to Nitrogen (N) stress and validated these with the reported QTLs in our earlier study to find the stable and consistent QTLs across populations. Also, with Illumina RNA-sequencing we checked the differential expression of gene (DEG) transcripts between parents and pools of RILs with high and low nitrogen use efficiency (NUE) and overlaid these DEGs on to the common validated QTLs to find candidate genes associated with N-stress tolerance in sorghum. Results: An F7 RIL population derived from a cross between CK60 (N-stress sensitive) and San Chi San (N-stress tolerant) inbred sorghum lines was used to map QTLs for 11 agronomic traits tested under different N-levels. Composite interval mapping analysis detected a total of 32 QTLs for 11 agronomic traits. Validation of these QTLs revealed that of the detected, nine QTLs from this population were consistent with the reported QTLs in earlier study using CK60/China17 RIL population. The validated QTLs were located on chromosomes 1, 6, 7, 8, and 9. In addition, root transcriptomic profiling detected 55 and 20 differentially expressed gene (DEG) transcripts between parents and pools of RILs with high and low NUE respectively. Also, overlay of these DEG transcripts on to the validated QTLs found candidate genes transcripts for NUE and also showed the expected differential expression. For example, DEG transcripts encoding Lysine histidine transporter 1 (LHT1) had abundant expression in San Chi San and the tolerant RIL pool, whereas DEG transcripts encoding seed storage albumin, transcription factor IIIC (TFIIIC) and dwarfing gene (DW2) encoding multidrug resistance-associated protein-9 homolog showed abundant expression in CK60 parent, similar to earlier study. Conclusions: The validated QTLs among different mapping populations would be the most reliable and stable QTLs across germplasm. The DEG transcripts found in the validated QTL regions will serve as future candidate genes for enhancing NUE in sorghum using molecular approaches. [ABSTRACT FROM AUTHOR]

We present the leaf and floral transcriptomes of two hybridizing bromeliad species that differ in their major pollinator systems. Here we identified candidate genes responsible for pollinator attraction and reproductive isolation in these two species. We searched for candidate genes involved in floral traits, such as colour. Approximately 34 Gbp of cDNA sequence data were produced from both tissues and species, resulting in a total of 424 506 914 raw reads. The de novo-assembled transcriptomes consisted of a total of 263 955 contigs, further clustered into 110 977 unigenes. Over 58% of the unigenes were functionally annotated and assigned to one or more Gene Ontology terms. The transcriptomes revealed 144 unique transcripts that encode key enzymes in the flavonoid and anthocyanin biosynthesis pathways. The domain/family annotation and phylogenetic analysis allowed us to infer, by homology, potential functions of the genes encoding MYB, HD- ZIP and bZIP- HY5 transcription factors, as well as WD40 protein, which may be involved in anthocyanin and flavonoid regulation in these species. These candidate genes are associated with natural regulation in flower colour in other plant species and will facilitate future studies aimed at elucidating the molecular basis of adaptive differentiation and the evolution of mechanisms of pollinator-mediated reproductive isolation in these two bromeliads. In addition, we identified a total of 49 439 microsatellite loci. These resources will assist future research into adaptation and speciation events in bromeliad species, thus providing a starting point for investigation of the molecular mechanisms of the traits responsible for their reproductive isolation. [ABSTRACT FROM AUTHOR]

Background: Sorghum is an important C4 crop which relies on applied Nitrogen fertilizers (N) for optimal yields, of which substantial amounts are lost into the atmosphere. Understanding the genetic variation of sorghum in response to limited nitrogen supply is important for elucidating the underlying genetic mechanisms of nitrogen utilization. Results: A bi-parental mapping population consisting of 131 recombinant inbred lines (RILs) was used to map quantitative trait loci (QTLs) influencing different agronomic traits evaluated under normal N (100 kg.ha-1 fertilizer) and low N (0 kg.ha-1 fertilizer) conditions. A linkage map spanning 1614 cM was developed using 642 polymorphic single nucleotide polymorphisms (SNPs) detected in the population using Genotyping-By-Sequencing (GBS) technology. Composite interval mapping detected a total of 38 QTLs for 11 agronomic traits tested under different nitrogen levels. The phenotypic variation explained by individual QTL ranged from 6.2 to 50.8 %. Illumina RNA sequencing data generated on seedling root tissues revealed 726 differentially expressed gene (DEG) transcripts between parents, of which 108 were mapped close to the QTL regions. Conclusions: Co-localized regions affecting multiple traits were detected on chromosomes 1, 5, 6, 7 and 9. These potentially pleiotropic regions were coincident with the genomic regions of cloned QTLs, including genes associated with flowering time, Ma3 on chromosome 1 and Ma1 on chromosome 6, gene associated with plant height, Dw2 on chromosome 6. In these regions, RNA sequencing data showed differential expression of transcripts related to nitrogen metabolism (Ferredoxin-nitrate reductase), glycolysis (Phosphofructo-2-kinase), seed storage proteins, plant hormone metabolism and membrane transport. The differentially expressed transcripts underlying the pleiotropic QTL regions could be potential targets for improving sorghum performance under limited N fertilizer through marker assisted selection. [ABSTRACT FROM AUTHOR]

The acquisition of tissue samples from wild populations is a constant challenge in conservation biology, especially for endangered species and protected species where nonlethal sampling is the only option. Whole blood has been suggested as a nonlethal sample type that contains a high percentage of bodywide and genomewide transcripts and therefore can be used to assess the transcriptional status of an individual, and to infer a high percentage of the genome. However, only limited quantities of blood can be nonlethally sampled from small species and it is not known if enough genetic material is contained in only a few drops of blood, which represents the upper limit of sample collection for some small species. In this study, we developed a nonlethal sampling method, the laboratory protocols and a bioinformatic pipeline to sequence and assemble the whole blood transcriptome, using Illumina RNA-Seq, from wild greater mouse-eared bats (Myotis myotis). For optimal results, both ribosomal and globin RNAs must be removed before library construction. Treatment of DNase is recommended but not required enabling the use of smaller amounts of starting RNA. A large proportion of protein-coding genes (61%) in the genome were expressed in the blood transcriptome, comparable to brain (65%), kidney (63%) and liver (58%) transcriptomes, and up to 99% of the mitogenome (excluding D-loop) was recovered in the RNA-Seq data. In conclusion, this nonlethal blood sampling method provides an opportunity for a genomewide transcriptomic study of small, endangered or critically protected species, without sacrificing any individuals. [ABSTRACT FROM AUTHOR]

It has been shown that the best coverage of the HepG2 cell line transcriptome encoded by genes of a single chromosome, chromosome 18, is achieved by a combination of two sequencing platforms, Illumina RNA-Seq and Oxford Nanopore Technologies (ONT), using cut-off levels of FPKM > 0 and TPM > 0, respectively. In this study, we investigated the extent to which the combination of these transcriptomic analysis methods makes it possible to achieve a high coverage of the transcriptome encoded by the genes of other human chromosomes. A comparative analysis of transcriptome coverage for various types of biological material was carried out, and the HepG2 cell line transcriptome was compared with the transcriptome of liver tissue cells. In addition, the contribution of variability in the coverage of expressed genes in human transcriptomes to the creation of a draft human transcriptome was evaluated. For human liver tissues, ONT makes an extremely insignificant contribution to the overall coverage of the transcriptome. Thus, to ensure maximum coverage of the liver tissue transcriptome, it is sufficient to apply only one technology: Illumina RNA-Seq (FPKM > 0)., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ilgisonis, Ponomarenko, Tarbeeva, Lisitsa, Zgoda, Radko and Archakov.)

NAC [no apical meristem (NAM), Arabidopsis thaliana transcription activation factor [ATAF1/2] and cup-shaped cotyledon (CUC2)] proteins belong to one of the largest plant-specific transcription factor (TF) families and play important roles in plant development processes, response to biotic and abiotic cues and hormone signalling. Our genome-wide analysis identified 110 StNAC genes in potato encoding for 136 proteins, including 14 membrane-bound TFs. The physical map positions of StNAC genes on 12 potato chromosomes were non-random, and 40 genes were found to be distributed in 16 clusters. The StNAC proteins were phylogenetically clustered into 12 subgroups. Phylogenetic analysis of StNACs along with their Arabidopsis and rice counterparts divided these proteins into 18 subgroups. Our comparative analysis has also identified 36 putative TNAC proteins, which appear to be restricted to Solanaceae family. In silico expression analysis, using Illumina RNA-seq transcriptome data, revealed tissue-specific, biotic, abiotic stress and hormone-responsive expression profile of StNAC genes. Several StNAC genes, including StNAC072 and StNAC101that are orthologs of known stress-responsive Arabidopsis RESPONSIVE TO DEHYDRATION 26 (RD26) were identified as highly abiotic stress responsive. Quantitative real-time polymerase chain reaction analysis largely corroborated the expression profile of StNAC genes as revealed by the RNA-seq data. Taken together, this analysis indicates towards putative functions of several StNAC TFs, which will provide blue-print for their functional characterization and utilization in potato improvement. [ABSTRACT FROM PUBLISHER]

The soft-shelled turtle Pelodiscus sinensis is a high-profile turtle species because of its nutritional and medicinal value in Asian countries. However, little is known about the genes that are involved in formation of their nutritional quality traits, especially the molecular mechanisms responsible for unsaturated fatty acid and collagen biosynthesis. In the present study, the transcriptomes from six tissues from Pelodiscus sinensis were sequenced using an Illumina paired-end sequencing platform. We obtained more than 47 million sequencing reads and 73 954 unigenes with an average size of 754 bp by de-novo assembly. In total, 55.19% of the unigenes (40 814) had significant similarity with proteins in the National Center of Biotechnology Information (NCBI) non-redundant protein database and Swiss-Prot database ( E-value <10). Of these annotated unigenes, 9 156 and 11 947 unigenes were assigned to 52 gene ontology categories (GO) and 25 clusters of orthologous groups (COG), respectively. In total, 26 496 (35.83%) unigenes were assigned to 242 pathways using the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, we found a number of highly expressed genes involved in the regulation of P. sinensis unsaturated fatty acid biosynthesis and collagen formation, including desaturases, growth factors, transcription factors, and extracellular matrix components. Our data represent the most comprehensive sequence resource available for the Chinese soft-shelled turtle and could provide a basis for new research on this turtle as well as the molecular genetics and functional genomics of other terrapins. To our knowledge, we report for the first time, the large-scale RNA sequencing (RNA-Seq) of terrapin animals and would enrich the knowledge of turtles for future research. [ABSTRACT FROM AUTHOR]

The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens. [ABSTRACT FROM AUTHOR]