Search

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.

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.

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.

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.

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 (COL R) 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 COL R Ab. Our study investigated two Italian XDR isogenic colistin-susceptible/resistant (COL S/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 COL R Ab isolates was conducted. COL R vs. COL S sRNAome, analyzed integrating the DE sRNAs obtained from the two different libraries, revealed some statistically significant DE sRNAs in COL R Ab. In detail, we found: (i) two different under-expressed cis -acting sRNAs (Ab sRNA and Ab sRNA) mapping in antisense orientation the 16S rRNA gene A1S_r01, (ii) one under-expressed cis -acting sRNA (Ab sRNA) targeting the A1S_2505 gene (hypothetical protein), (iii) one under-expressed microRNA-size small RNA fragment (Ab sRNA) and its pre-micro Ab sRNA targeting the A1S_0501 gene (hypothetical protein), (iv) as well as an over-expressed microRNA-size small RNA fragment (Ab sRNA) and its pre-micro Ab sRNA 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 COL R 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.

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.

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. - 2019 Elsevier Inc. Scopus

[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., This work has been carried out with financial support from the European Union, through the Horizon H2020 research and innovation programme under grant agreement 634429 (ParaFishControl). MCP was contracted under CSIC PIE project no. 201740E013. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

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.

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.

Chinese hawthorn (Crataegus pinnatifida) is an important ornamental and economic horticultural plant. However, the lack of molecular markers has limited the development and utilization of hawthorn germplasm resources. Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) allow precise and effective cultivar characterization and are routinely used for genetic diversity analysis. Thus, we first reported the development of polymorphic EST-SSR markers in C. pinnatifida with perfect repeats using Illumina RNA-Seq technique. In total, we investigated 14,364 unigenes, from which 5091 EST-SSR loci were mined. Di-nucleotides (2012, 39.52%) were the most abundant SSRs, followed by mono- (1989, 39.07%), and tri-nucleotides (1024, 20.11%). On the basis of these EST-SSRs, a total of 300 primer pairs were designed and used for polymorphism analysis in 70 accessions collected from different geographical regions of China. Of 239 (79.67%) pairs of primer-generated amplification products, 163 (54.33%) pairs of primers showed polymorphism. Finally, 33 primers with high polymorphism were selected for genetic diversity analysis and tested on 70 individuals with low-cost fluorescence-labeled M13 primers using capillary electrophoresis genotyping platform. A total of 108 alleles were amplified by 33 SSR markers, with the number of alleles (Na) ranging from 2 to 14 per locus (mean: 4.939), and the effective number of alleles (Ne) ranging from 1.258 to 3.214 (mean: 2.221). The mean values of gene diversity (He), observed heterozygosity (Ho), and polymorphism information content (PIC) were 0.524 (range 0.205&ndash, 0.689), 0.709 (range 0.132&ndash, 1.000), and 0.450 (range 0.184&ndash, 0.642), respectively. Furthermore, the dendrogram constructed based on the EST-SSR separated the cultivars into two main clusters. In sum, our study was the first comprehensive study on the development and analysis of a large set of SSR markers in hawthorn. The results suggested that the use of NGS techniques for SSR development represented a powerful tool for genetic studies. Additionally, fluorescence-labeled M13 markers proved to be a valuable method for genotyping. All of these EST-SSR markers have agronomic potential and constitute a scientific basis for future studies on the identification, classification, and innovation of hawthorn germplasms.

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.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12870-017-1064-9) contains supplementary material, which is available to authorized users.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12870-015-0696-x) contains supplementary material, which is available to authorized users.

Background The polyadenylation of RNA is critical for gene functioning, but the conserved sequence motifs (often called signal or signature motifs), motif locations and abundances, and base composition patterns around mRNA polyadenylation [poly(A)] sites are still uncharacterized in most species. The evolutionary tendency for poly(A) site selection is still largely unknown. Results We analyzed the poly(A) site regions of 31 species or phyla. Different groups of species showed different poly(A) signal motifs: UUACUU at the poly(A) site in the parasite Trypanosoma cruzi; UGUAAC (approximately 13 bases upstream of the site) in the alga Chlamydomonas reinhardtii; UGUUUG (or UGUUUGUU) at mainly the fourth base downstream of the poly(A) site in the parasite Blastocystis hominis; and AAUAAA at approximately 16 bases and approximately 19 bases upstream of the poly(A) site in animals and plants, respectively. Polyadenylation signal motifs are usually several hundred times more abundant around poly(A) sites than in whole genomes. These predominant motifs usually had very specific locations, whether upstream of, at, or downstream of poly(A) sites, depending on the species or phylum. The poly(A) site was usually an adenosine (A) in all analyzed species except for B. hominis, and there was weak A predominance in C. reinhardtii. Fungi, animals, plants, and the protist Phytophthora infestans shared a general base abundance pattern (or base composition pattern) of “U-rich—A-rich—U-rich—Poly(A) site—U-rich regions”, or U-A-U-A-U for short, with some variation for each kingdom or subkingdom. Conclusion This study identified the poly(A) signal motifs, motif locations, and base composition patterns around mRNA poly(A) sites in protists, fungi, plants, and animals and provided insight into poly(A) site evolution.