Characterization and development of transcriptome-derived novel EST-SSR markers to assess genetic diversity in Chaetomium globosum.

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. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03794-7.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features.

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.

Combined PacBio Iso-Seq and Illumina RNA-Seq Analysis of the Tuta absoluta (Meyrick) Transcriptome and Cytochrome P450 Genes.

Tuta absoluta (Meyrick) is a devastating invasive pest worldwide. The abamectin and chlorantraniliprole complex have become an alternative option for chemical control because they can enhance insecticidal activity and delay increased drug resistance. Notably, pests are inevitably resistant to various types of insecticides, and compound insecticides are no exception. To identify potential genes involved in the detoxification of abamectin and chlorantraniliprole complex in T. absoluta, PacBio SMRT-seq transcriptome sequencing and Illumina RNA-seq analysis of abamectin and chlorantraniliprole complex-treated T. absoluta were performed. We obtained 80,492 non-redundant transcripts, 62,762 (77.97%) transcripts that were successfully annotated, and 15,524 differentially expressed transcripts (DETs). GO annotation results showed that most of these DETs were involved in the biological processes of life-sustaining activities, such as cellular, metabolic, and single-organism processes. The KEGG pathway enrichment results showed that the pathways related to glutathione metabolism, fatty acid and amino acid synthesis, and metabolism were related to the response to abamectin and chlorantraniliprole complex in T. absoluta. Among these, 21 P450s were differentially expressed (11 upregulated and 10 downregulated). The qRT-PCR results for the eight upregulated P450 genes after abamectin and chlorantraniliprole complex treatment were consistent with the RNA-Seq data. Our findings provide new full-length transcriptional data and information for further studies on detoxification-related genes in T. absoluta.

Gene-centric coverage of the human liver transcriptome: QPCR, Illumina, and Oxford Nanopore RNA-Seq.

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).

Salt Stress Tolerance in Casuarina glauca: Insights from the Branchlets Transcriptome.

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.

Liposomal Formulation of Botanical Extracts may Enhance Yield Triggering PR Genes and Phenylpropanoid Pathway in Barley (Hordeum vulgare).

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.

De Novo Transcriptomic Characterization Enables Novel Microsatellite Identification and Marker Development in Betta splendens.

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.

Comparative physiological and full-length transcriptome analyses reveal the molecular mechanism of melatonin-mediated salt tolerance in okra (Abelmoschus esculentus L.).

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.

Innate immune-gene expression during experimental amyloodiniosis in European seabass (Dicentrarchus labrax).

The ectoparasite protozoan Amyloodinium ocellatum (AO) is the causative agent of amyloodiniosis in European seabass (ESB, Dicentrarchus labrax). There is a lack of information about basic molecular immune response mechanisms of ESB during AO infestation. Therefore, to compare gene expression between experimental AO-infested ESB tissues and uninfested ESB tissues (gills and head kidney) RNA-seq was adopted. The RNA-seq revealed multiple differentially expressed genes (DEG), namely 679 upregulated genes and 360 downregulated genes in the gills, and 206 upregulated genes and 170 downregulated genes in head kidney. In gills, genes related to the immune system (perforin, CC1) and protein binding were upregulated. Several genes involved in IFN related pathways were upregulated in the head kidney. Subsequently, to validate the DEG from amyloodiniosis, 26 ESB (mean weight 14 g) per tank in triplicate were bath challenged for 2 h with AO (3.5 × 106/tank; 70 dinospores/mL) under controlled conditions (26-28 °C and 34‰ salinity). As a control group (non-infested), 26 ESB per tank in triplicate were also used. Changes in the expression of innate immune genes in gills and head kidney at 2, 3, 5, 7 and 23 dpi were analysed using real-time PCR. The results indicated that the expression of cytokines (CC1, IL-8) and antimicrobial peptide (Hep) were strongly stimulated and reached a peak at 5 dpi in the early infestation stage, followed by a gradual reduction in the recovery stage (23 dpi). Noticeably, the immunoglobulin (IgM) expression was higher at 23 dpi compared to 7 dpi. Furthermore, in-situ hybridization showed positive signals of CC1 mRNA in AO infested gills compared to the control group. Altogether, chemokines were involved in the immune process under AO infestation and this evidence allows a better understanding of the immune response in European seabass during amyloodiniosis.

SMRT and Illumina RNA sequencing reveal novel insights into the heat stress response and crosstalk with leaf senescence in tall fescue.

BACKGROUND: As a cool-season grass species, tall fescue (Festuca arundinacea) is challenged by increasing temperatures. Heat acclimation or activation of leaf senescence, are two main strategies when tall fescue is exposed to heat stress (HS). However, lacking a genome sequence, the complexity of hexaploidy nature, and the short read of second-generation sequencing hinder a comprehensive understanding of the mechanism. This study aims to characterize the molecular mechanism of heat adaptation and heat-induced senescence at transcriptional and post-transcriptional levels. RESULTS: Transcriptome of heat-treated (1 h and 72 h) and senescent leaves of tall fescue were generated by combining single-molecular real-time and Illumina sequencing. In total, 4076; 6917, and 11,918 differentially expressed genes (DEGs) were induced by short- and long-term heat stress (HS), and senescence, respectively. Venn and bioinformatics analyses of DEGs showed that short-term HS strongly activated heat shock proteins (Hsps) and heat shock factors (Hsfs), as well as specifically activated FK506-binding proteins (FKBPs), calcium signaling genes, glutathione S-transferase genes, photosynthesis-related genes, and phytohormone signaling genes. By contrast, long-term HS shared most of DEGs with senescence, including the up-regulated chlorophyll catabolic genes, phytohormone synthesis/degradation genes, stress-related genes, and NACs, and the down-regulated photosynthesis-related genes, FKBPs, and catalases. Subsequently, transient overexpression in tobacco showed that FaHsfA2a (up-regulated specifically by short-term HS) reduced cell membrane damages caused by HS, but FaNAC029 and FaNAM-B1 (up-regulated by long-term HS and senescence) increased the damages. Besides, alternative splicing was widely observed in HS and senescence responsive genes, including Hsps, Hsfs, and phytohormone signaling/synthesis genes. CONCLUSIONS: The short-term HS can stimulate gene responses and improve thermotolerance, but long-term HS is a damage and may accelerate leaf senescence. These results contribute to our understanding of the molecular mechanism underlying heat adaptation and heat-induced senescence.

SMRT sequencing of the full-length transcriptome of the Rhynchophorus ferrugineus (Coleoptera: Curculionidae).

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.

Low Long Terminal Repeat (LTR)-Retrotransposon Expression in Leaves of the Marine Phanerogam Posidonia Oceanica L.

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.

Transcriptome assembly and abiotic related gene expression analysis of date palm reveal candidate genes involved in response to cadmium stress.

Date palm in Tunisia is of major economic importance but are also factors of social, environmental and economic stability. An annotated assembly of the transcriptome of cultivar Deglet Nour was reported. RNA was isolated from plant Cd-contaminated leaves, and 37,049 unique Illumina RNA-seq reads were used in a transcriptome assembly. The draft transcriptome assembly consists of 6789 contigs and 17.285 singletons with a means length of 858 bp and 1.042 bp, respectively. The final assembly was functionally annotated using Blast2GO software, allowing the identification of putative genes controlling important agronomic traits. The annotated transcriptome data sets were used to query all known Kyoto Encyclopedia of Genes and Genomes pathways. The most represented molecular functions and biological processes were nucleotide binding and transcription, transport and response to stress and abiotic and biotic stimuli. A prediction of the genes interaction network was proposed by selecting corresponding functionally similar genes from Arabidopsis datasets, downloaded by GeneMANIA version 2.1. Several Cd-responsive genes expression was monitored in in vitro isolated explant of Cd stressed Deglet Nour. Some chelators encoding genes were upregulated confirming in silico findings. Genes encoding HMs transporters in date palm showed expression enhancement more pronounced after 20 days of exposure. P. dactylifera transcriptome provides a valuable resource for future functional analysis of candidate genes involved in metal stress response.

The complexity of alternative splicing and landscape of tissue-specific expression in lotus (Nelumbo nucifera) unveiled by Illumina- and single-molecule real-time-based RNA-sequencing.

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.

Acting locally - affecting globally: RNA sequencing of gilthead sea bream with a mild Sparicotyle chrysophrii infection reveals effects on apoptosis, immune and hypoxia related genes.

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.

COLR Acinetobacter baumannii sRNA Signatures: Computational Comparative Identification and Biological Targets.

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 (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 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.

Validation of QTL mapping and transcriptome profiling for identification of candidate genes associated with nitrogen stress tolerance in sorghum.

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.

De novo assembly and characterization of leaf and floral transcriptomes of the hybridizing bromeliad species (Pitcairnia spp.) adapted to Neotropical Inselbergs.

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.

A nonlethal sampling method to obtain, generate and assemble whole blood transcriptomes from small, wild mammals.

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.

De Novo Transcriptome Assembly and Development of Novel Microsatellite Markers for the Traditional Chinese Medicinal Herb, Veratrilla baillonii Franch (Gentianaceae).

Veratrilla baillonii Franch is an important Chinese medicinal herb for treating liver-related diseases, which has been over-collected in the recent decades. However, the effective conservation and related population genetic study has been hindered because of the lack of genome sequences and genetic markers in the natural population. We have conducted RNA-seq on V. baillonii. We performed de novo assembly of these data to characterize the V. baillonii transcriptome, resulting in 133,019 contigs with size >200 bp. These contigs were annotated using the NCBI nonredundant database and Gene Ontology (GO) terms. From these contigs, we developed novel microsatellite simple sequence repeat (SSR) markers, identifying a total of 40,885 SSRs. SSRs with repeat motifs of 1-4 bp (mono-, di-, tri-, and tetranucleotides) accounted for 99.8% of all SSRs, with mononucleotide repeats most common, followed by dinucleotide (16.2%) and trinucleotide repeats (14.7%). We selected 151 SSRs for experimental validation, of which 74 were confirmed by polymerase chain reaction. Fourteen SSRs were determined to be polymorphic by screening 40 individuals from six distant populations. The number of alleles per locus ranged from two to four, and the expected heterozygosity varied from 0.2637 to 0.8571, suggesting that these SSR markers are highly polymorphic and effective for further genetic analysis in the nature population. In addition, we explored the genetic structure of V. baillonii using five SSRs in four geographic populations and found that the identified genotypes were clustered into two phylogenetic clades: the Mekong River clade and Jinsha River clade. This result indicates that these two regions may harbor highly divergent genetic lineages and enriched genetic diversity. The de novo transcriptome sequences and new SSR markers discovered by this study provide an initial step for understanding the population genetics of V. baillonii, and a valuable resource for effective conservation management.