Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

De novo transcriptome assembly of Dalbergia sissoo Roxb. (Fabaceae) under Botryodiplodia theobromae-induced dieback disease

Dalbergia sissoo Roxb. (Shisham) is a timber-producing species of economic, cultural, and medicinal importance in the Indian subcontinent. In the past few decades, Shisham's dieback disease caused by the fungus Botryodiplodia theobromae has become an evolving issue in the subcontinent endangering its survival. To gain insights into this issue, a standard transcriptome assembly was deployed to assess the response of D. sissoo at the transcriptomic level under the stress of B. theobromae infection. For RNA isolation, the control and infected leaf tissue samples were taken from 1-year-old greenhouse-grown D. sissoo plants after 20 days of stem-base spore inoculation. cDNA synthesis was performed from these freshly isolated RNA samples that were then sent for sequencing. About 18.14 Gb (Giga base) of data was generated using the BGISEQ-500 sequencing platform. In terms of Unigenes, 513,821 were identified after a combined assembly of all samples and then filtering the abundance. The total length of Unigenes, their average length, N50, and GC-content were 310,523,693 bp, 604 bp, 1,101 bp, and 39.95% respectively. The Unigenes were annotated using 7 functional databases i.e., 200,355 (NR: 38.99%), 164,973 (NT: 32.11%), 123,733 (Swissprot: 24.08%), 142,580 (KOG: 27.75%), 139,588 (KEGG: 27.17%), 99,752 (GO: 19.41%), and 137,281 (InterPro: 26.72%). Furthermore, the Transdecoder detected 115,762 CDS. In terms of SSR (Simple Sequence Repeat) markers, 62,863 of them were distributed on 51,508 Unigenes and on the predicted 4673 TF (Transcription Factor) coding Unigenes. A total of 16,018 up- and 19,530 down-regulated Differentially Expressed Genes (DEGs) were also identified. Moreover, the Plant Resistance Genes (PRGs) had a count of 9230. We are hopeful that in the future, these identified Unigenes, SSR markers, DEGs and PRGs will provide the prerequisites for managing Shisham dieback disease, its breeding, and in tree improvement programs.

Resistant and Susceptible Pinus thunbergii ParL. Show Highly Divergent Patterns of Differentially Expressed Genes during the Process of Infection by Bursaphelenchus xylophilus

Pine wilt disease (PWD) is a devastating disease that threatens pine forests worldwide, and breeding resistant pines is an important management strategy used to reduce its impact. A batch of resistant seeds of P. thunbergii was introduced from Japan. Based on the resistant materials, we obtained somatic plants through somatic embryogenesis. In this study, we performed transcriptome analysis to further understand the defense response of resistant somatic plants of P. thunbergii to PWD. The results showed that, after pine wood nematode (PWN) infection, resistant P. thunbergii stimulated more differential expression genes (DEGs) and involved more regulatory pathways than did susceptible P. thunbergii. For the first time, the alpha-linolenic acid metabolism and linoleic acid metabolism were intensively observed in pines resisting PWN infection. The related genes disease resistance protein RPS2 (SUMM2) and pathogenesis-related genes (PR1), as well as reactive oxygen species (ROS)-related genes were significantly up-expressed in order to contribute to protection against PWN inoculation in P. thunbergii. In addition, the diterpenoid biosynthesis pathway was significantly enriched only in resistant P. thunbergii. These findings provided valuable genetic information for future breeding of resistant conifers, and could contribute to the development of new diagnostic tools for early screening of resistant pine seedlings based on specific PWN-tolerance-related markers.

Comprehensive Organ-Specific Profiling of Douglas Fir (Pseudotsuga menziesii) Proteome

The Douglas fir (Pseudotsuga menziesii) is a conifer native to North America that has become increasingly popular in plantations in France due to its many advantages as timber: rapid growth, quality wood, and good adaptation to climate change. Tree genetic improvement programs require knowledge of a species' genetic structure and history and the development of genetic markers. The very slow progress in this field, for Douglas fir as well as the entire genus Pinus, can be explained using the very large size of their genomes, as well as by the presence of numerous highly repeated sequences. Proteomics, therefore, provides a powerful way to access genomic information of otherwise challenging species. Here, we present the first Douglas fir proteomes acquired using nLC-MS/MS from 12 different plant organs or tissues. We identified 3975 different proteins and quantified 3462 of them, then examined the distribution of specific proteins across plant organs/tissues and their implications in various molecular processes. As the first large proteomic study of a resinous tree species with organ-specific profiling, this short note provides an important foundation for future genomic annotations of conifers and other trees.

Development of EST-SSRs based on the transcriptome of Castanopsis carlesii and cross-species transferability in other Castanopsis species

Castanopsis carlesii (Hemsl.) Hay. is a widely distributed and dominant tree species native to subtropical China with significant ecological and economic value. Due to serious human-related disturbance, its wild resources have been increasingly reduced, and whether may result in the loss of genetic diversity. However, no population genetics studies of natural C. carlesii have been reported to date. Microsatellite markers have been a useful tool in population genetics. Therefore, we developed EST-SSR markers based on the transcriptome sequencing of C. carlesii leaves. A total of 149,380,224 clean reads were obtained, and 63,012 nonredundant unigenes with a mean length of 1,034 bp were assembled and annotated based on sequence similarity searches in the Nr, Nt, KO, SwissProt, PFAM, KOG, and GO databases. The results showed that only 5,559 (8.82%) unigenes were annotated in all seven databases, but 46,338 (73.53%) could be annotated in at least one database. A total of 31,459 potential EST-SSRs were identified in 18,690 unigenes, with an average frequency of one SSR approximately 2 kb. Among the 100 EST-SSR primer pairs designed, 49 primer pairs successfully produced the expected product by amplification, with a success rate of 49%, but only 20 primer pairs showed abundant polymorphisms. Polymorphisms were verified using 25 samples from C. carlesii in Qimen, Anhui. A total of 119 alleles were detected, with a mean number of alleles (Na) of 5.95 per locus and a mean polymorphism information content (PIC) of 0.6125. All the 20 newly developed EST-SSR markers were verified in other Castanopsis species (C. sclerophylla, C. lamontii, C. fargesii, C. eyrei and C. jucunda). Sixteen primer pairs showed successful amplification in all five Castanopsis species (80%), and the transferability ratios ranged from 90% to 100%. These developed EST-SSR markers can be applied to population genetic and germplasm evaluations of C. carlesii and related species.

Cyanidin-3-O-glucoside Contributes to Leaf Color Change by Regulating Two bHLH Transcription Factors in Phoebe bournei

Anthocyanins produce different-colored pigments in plant organs, which provide ornamental value. Thus, this study was conducted to understand the mechanism of anthocyanin synthesis in ornamental plants. Phoebe bournei, a Chinese specialty tree, has high ornamental and economic value due to its rich leaf color and diverse metabolic products. Here, the metabolic data and gene expression of red P. bournei leaves at the three developmental stages were evaluated to elucidate the color-production mechanism in the red-leaved P. bournei species. First, metabolomic analysis identified 34 anthocyanin metabolites showing high levels of cyanidin-3-O-glucoside (cya-3-O-glu) in the S1 stage, which may suggest that it is a characteristic metabolite associated with the red coloration of the leaves. Second, transcriptome analysis showed that 94 structural genes were involved in anthocyanin biosynthesis, especially flavanone 3'-hydroxy-lase (PbF3'H), and were significantly correlated with the cya-3-O-glu level. Third, K-means clustering analysis and phylogenetic analyses identified PbbHLH1 and PbbHLH2, which shared the same expression pattern as most structural genes, indicating that these two PbbHLH genes may be regulators of anthocyanin biosynthesis in P. bournei. Finally, overexpression of PbbHLH1 and PbbHLH2 in Nicotiana tabacum leaves triggered anthocyanin accumulation. These findings provide a basis for cultivating P. bournei varieties that have high ornamental value.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

Transcriptome Analysis and Identification of a Female-Specific SSR Marker in Pistacia chinensis Based on Illumina Paired-End RNA Sequencing

Pistacia chinensis Bunge (P. chinensis), a dioecious plant species, has been widely found in China. The female P. chinensis plants are more important than male plants in agricultural production, as their seeds can serve as an ideal feedstock for biodiesel. However, the sex of P. chinensis plants is hard to distinguish during the seedling stage due to the scarcity of available transcriptomic and genomic information. In this work, Illumina paired-end RNA sequencing assay was conducted to unravel the transcriptomic profiles of female and male P. chinensis flower buds. In total, 50,925,088 and 51,470,578 clean reads were obtained from the female and male cDNA libraries, respectively. After quality checks and de novo assembly, a total of 83,370 unigenes with a mean length of 1.3 kb were screened. Overall, 64,539 unigenes (77.48%) could be matched in at least one of the NR, NT, Swiss-Prot, COG, KEGG, and GO databases, 71 of which were putatively related to the floral development of P. chinensis. Additionally, 21,662 simple sequence repeat (SSR) motifs were identified in 17,028 unigenes of P. chinensis, and the mononucleotide motif was the most dominant type of repeats (52.59%) in P. chinensis, followed by dinucleotide (22.29%), trinucleotide (20.15%). The most abundant repeats were AG/CT (13.97%), followed by AAC/GTT (6.75%) and AT/TA (6.10%). Based on these SSR, 983 EST-SSR primers were designed, 151 of which were randomly chosen for validation. Of these validated EST-SSR markers, 25 SSR markers were found to be polymorphic between male and female plants. One SSR marker, namelyPCSSR55, displayed excellent specificity in female plants, which could clearly distinguish between male and female P. chinensis. Altogether, our findings not only reveal that the EST-SSR marker is extremely effective in distinguishing between male and female P. chinensis but also provide a solid framework for sex determination of plant seedlings.

Transcriptome analysis provides new insights into cold adaptation of corsac fox ( Vulpes Corsac )

Vulpesare widely distributed throughout the world and have undergone drastic physiological and phenotypic changes in response to their environment. However, little is known about the underlying genetic causes of these traits, especially Vulpes corsac. In this study, RNA‐Seq was used to obtain a comprehensive dataset for multiple pooled tissues of corsac fox, and selection analysis of orthologous genes was performed to identify the genes that may be influenced by the low‐temperature environment. More than 6.32 Gb clean reads were obtained and assembled into a total of 173,353 unigenes with an average length of 557 bp for corsac fox. Selective pressure analysis showed that 16 positively selected genes (PSGs) were identified in corsac fox, red fox, and arctic fox. Enrichment analysis of PSGs showed that the LRP11 gene was enriched in several pathways related to the low‐temperature response and might play a key role in response to environmental stimuli of foxes. In addition, several positively selected genes were related to DNA damage repair (ELP2 and CHAF1A), innate immunity (ARRDC4 and S100A12), and the respiratory chain (NDUFA5), and these positively selected genes might play a role in adaptation to harsh wild fox environments. The results of common orthologous gene analysis showed that gene flow or convergent evolution might be an important factor in promoting regional differentiation of foxes. Our study provides a valuable transcriptomic resource for the evolutionary history of the corsac fox and the adaptations to the extreme environments.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

De novo assembly, transcriptome characterization and marker discovery in Indian major carp, Labeo rohita through pyrosequencing

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

Genome-Wide Transcriptome Profiling, Characterization, and Functional Identification of NAC Transcription Factors in Sorghum under Salt Stress

Salinity stress has become a significant concern to global food security. Revealing the mechanisms that enable plants to survive under salinity has immense significance. Sorghum has increasingly attracted researchers interested in understanding the survival and adaptation strategies to high salinity. However, systematic analysis of the DEGs (differentially expressed genes) and their relative expression has not been reported in sorghum under salt stress. The de novo transcriptomic analysis of sorghum under different salinity levels from 60 to 120 mM NaCl was generated using Illumina HiSeq. Approximately 323.49 million high-quality reads, with an average contig length of 1145 bp, were assembled de novo. On average, 62% of unigenes were functionally annotated to known proteins. These DEGs were mainly involved in several important metabolic processes, such as carbohydrate and lipid metabolism, cell wall biogenesis, photosynthesis, and hormone signaling. SSG 59-3 alleviated the adverse effects of salinity by suppressing oxidative stress (H₂O₂) and stimulating enzymatic and non-enzymatic antioxidant activities (SOD, APX, CAT, APX, POX, GR, GSH, ASC, proline, and GB), as well as protecting cell membrane integrity (MDA and electrolyte leakage). Significant up-regulation of transcripts encoding the NAC, MYB, and WRYK families, NHX transporters, the aquaporin protein family, photosynthetic genes, antioxidants, and compatible osmolyte proteins were observed. The tolerant line (SSG 59-3) engaged highly efficient machinery in response to elevated salinity, especially during the transport and influx of K⁺ ions, signal transduction, and osmotic homeostasis. Our data provide insights into the evolution of the NAC TFs gene family and further support the hypothesis that these genes are essential for plant responses to salinity. The findings may provide a molecular foundation for further exploring the potential functions of NAC TFs in developing salt-resistant sorghum lines.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

A robust KASP marker for selection of four pairs of linked leaf rust and stripe rust resistance genes introgressed on chromosome arm 5DS from different wheat genomes

Stripe rust and leaf rust are among the most devastating diseases of wheat, limiting its production globally. Wheat wild relatives harbour genetic diversity for new genes and alleles for all major wheat diseases. However, the use of this genetic variation from wild progenitor and non-progenitor species has been limited in the breeding programs. Reasons include limited recombination of donor and recipient genomes and the lack of tertiary gene pool markers. Here, we describe the development of a SNP based marker from the flow-sorted and sequenced Aegilops umbellulata chromosome 5U which can be used for marker assisted selection of four pair of alien leaf rust and stripe rust resistance genes. Lr57-Yr40_CAPS16 marker was reported earlier to be linked with alien leaf and stripe rust resistance genes introgressed on wheat chromosome 5DS. Due to its dominant nature and laborious to work with, a new SNP-based KASP marker, XTa5DS-2754099_kasp23, was developed from the same CAPS marker contig. XTa5DS-2754099_kasp23 was tested in Aegilops umbellulata, Ae. geniculata, Ae. peregrina and Ae. caudata derived alien introgression lines, which harbour four pairs of linked leaf and stripe rust genes; Lr76-Yr70, Lr57-Yr40, LrP- YrP, LrAc-YrAc, respectively. This KASP marker was found to be effective for the selection of the aforesaid four pairs of leaf rust and stripe rust resistance genes. Further, we tested and validated XTa5DS-2754099_kasp23 on commercial varieties and advanced breeding lines from four countries (India, Egypt, Australia and UK) including hexaploid and durum wheat. Our results provide evidence that KASP marker, XTa5DS-2754099_kasp23 can be used in marker-assisted selection of the four pairs of rust resistance alien genes in wheat breeding programmes.

de novo transcriptomic profiling of differentially expressed genes in grass halophyte Urochondra setulosa under high salinity

Soil salinity is one of the major limiting factors for crop productivity across the world. Halophytes have recently been a source of attraction for exploring the survival and tolerance mechanisms at extreme saline conditions. Urochondra setulosa is one of the obligate grass halophyte that can survive in up to 1000 mM NaCl. The de novo transcriptome of Urochondra leaves at different salt concentrations of 300-500 mM NaCl was generated on Illumina HiSeq. Approximately 352.78 million high quality reads with an average contig length of 1259 bp were assembled de novo. A total of 120,231 unigenes were identified. On an average, 65% unigenes were functionally annotated to known proteins. Approximately 35% unigenes were specific to Urochondra. Differential expression revealed significant enrichment (P < 0.05) of transcription factors, transporters and metabolites suggesting the transcriptional regulation of ion homeostasis and signalling at high salt concentrations in this grass. Also, about 143 unigenes were biologically related to salt stress responsive genes. Randomly selected genes of important pathways were validated for functional characterization. This study provides useful information to understand the gene regulation at extremely saline levels. The study offers the first comprehensive evaluation of Urochondra setulosa leaf transcriptome. Examining non-model organisms that can survive in harsh environment can provide novel insights into the stress coping mechanisms which can be useful to develop improved agricultural crops.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Raw transcriptomics data to gene specific SSRs: a validated free bioinformatics workflow for biologists

Recent advances in next-generation sequencing technologies have paved the path for a considerable amount of sequencing data at a relatively low cost. This has revolutionized the genomics and transcriptomics studies. However, different challenges are now created in handling such data with available bioinformatics platforms both in assembly and downstream analysis performed in order to infer correct biological meaning. Though there are a handful of commercial software and tools for some of the procedures, cost of such tools has made them prohibitive for most research laboratories. While individual open-source or free software tools are available for most of the bioinformatics applications, those components usually operate standalone and are not combined for a user-friendly workflow. Therefore, beginners in bioinformatics might find analysis procedures starting from raw sequence data too complicated and time-consuming with the associated learning-curve. Here, we outline a procedure for de novo transcriptome assembly and Simple Sequence Repeats (SSR) primer design solely based on tools that are available online for free use. For validation of the developed workflow, we used Illumina HiSeq reads of different tissue samples of Santalum album (sandalwood), generated from a previous transcriptomics project. A portion of the designed primers were tested in the lab with relevant samples and all of them successfully amplified the targeted regions. The presented bioinformatics workflow can accurately assemble quality transcriptomes and develop gene specific SSRs. Beginner biologists and researchers in bioinformatics can easily utilize this workflow for research purposes.

Transcriptome sequencing reveals signatures of positive selection in the Spot-Tailed Earless Lizard

The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.

SNP marker development in Pinus sylvestris L. in stress-responsive genes characterized from Pinus cembra L. transcriptomes

Massively parallel sequencing of cDNA is an efficient route for generating sequence collections that represent expressed genes under different environmental control. The analysis of their sequence helps in developing molecular markers, such as SNPs, which represent a useful tool in detecting adaptive signals in populations. In this study novel PCR markers, based on stress responsive genes, were designed from the transcriptome of the haploxylon Swiss stone pine (Pinus cembra L.) and tested for SNPs in the diploxylon Scots pine (Pinus sylvestris L.). 84 primers were tested on P. sylvestris DNA samples originating from three different types of habitat. After sequencing and BLAST search of the amplified products, parts of 19 different candidate genes were analysed by considering the polymorphic sites, insertions/deletions as well as synonymous and non-synonymous SNPs. In a total of 3735 sites no indels, eight synonymous and 11 non-synonymous SNPs were found. By providing de novo molecular markers developed in P. cembra and tested for transferability in Scots pine, our results give support for the use of de novo markers targeting conserved regions across different pines. The SNPs detected may have important applications in further studies of adaptive genetic variation, providing tools to study relevant genes important in the long-term adaptation of pine species.

Combining independent de novo assemblies to optimize leaf transcriptome of Persian walnut

Transcriptome resources can facilitate to increase yield and quality of walnuts. Finding the best transcriptome assembly has not been the subject of walnuts research as yet. This research generated 240,179,782 reads from 11 walnut leaves according to cDNA libraries. The reads provided a complete de novo transcriptome assembly. Fifteen different transcriptome assemblies were constructed from five different well-known assemblers used in scientific literature with different k-mer lengths (Bridger, BinPacker, SOAPdenovo-Trans, Trinity and SPAdes) as well as two merging approaches (EvidentialGene and Transfuse). Based on the four quality metrics of assembly, the results indicated an efficiency in the process of merging the assemblies after being generated by de novo assemblers. Finally, EvidentialGene was recognized as the best assembler for the de novo assembly of the leaf transcriptome in walnut. Among a total number of 183,191 transcripts which were generated by EvidentialGene, there were 109,413 transcripts capable of protein potential (59.72%) and 104,926 were recognized as ORFs (57.27%). In addition, 79,185 transcripts were predicted to exist with at least one hit to the Pfam database. A number of 3,931 transcription factors were identified by BLAST searching against PlnTFDB. Furthermore, 6,591 of the predicted peptide sequences contained signaling peptides, while 92,704 contained transmembrane domains. Comparison of the assembled transcripts with transcripts of the walnut and published genome assembly for the 'Chandler' cultivar using the BLAST algorithm led to identify a total number of 27,304 and 19,178 homologue transcripts, respectively. De novo transcriptomes in walnut leaves can be developed for the future studies in functional genomics and genetic studies of walnuts.

Draft genome of multiple resistance donor plant Sinapis alba: An insight into SSRs, annotations and phylogenetics

BACKGROUND

Sinapis alba is a wild member of the Brassicaceae family reported to possess genetic resistance against major biotic and abiotic stresses of oilseed brassicas. However, the resistance nature of S. alba was not exploited generously due to the unavailability of usable genome sequences in public databases. Therefore, the present study was conducted to assemble the first draft genome from raw whole genome shotgun sequences with annotation and develop simple sequence repeat markers for molecular genetics and marker-assisted breeding.

RESULTS

The raw genome sequences had 96x coverage on the Illumina platform with 170 Gbp data. The developed assembly by SOAPdenovo2 has ~459 Mbp genome size covered in 403,423 contigs with an average size of 1138.04 bp. The assembly was BLASTX with Arabidopsis thaliana which showed 32.9% positive hits between both plants. The top hit species distribution analysis showed the highest similarity with A. thaliana. A total of 809,597 GO level annotations were recorded after BLASTX results, and 34,012 sequences were annotated with different enzyme codes grouped under seven classes. The gene prediction tool AUGUSTUS identified 113,107 probable genes with an average size of 684 bp. The biochemical pathway annotation assigned 16,119 potential genes to 152 KEGG maps and 1751 enzyme codes. The development of potential SSRs from the de-novo assembly yielded 70731 unique primer pairs. Out of 159 randomly selected SSR markers for validation, 149 successfully amplified in S. alba. However, 10 SSR markers did not amplify during the validation experiment.

CONCLUSION

The annotated genome assembly with a large number of SSRs was developed in the present study. To the best of our knowledge, this is the first report of S. alba genome assembly development, annotation, and SSRs mining to date. The data presented here will be a very important resource for future crop improvement programs, especially for resistant breeding.

Functional genomics by integrated analysis of transcriptome of sweet potato (Ipomoea batatas (L.) Lam.) during root formation

BACKGROUND

Sweet potato is easily propagated by cuttings. But the molecular biological mechanism of adventitious root formation are not yet clear.

OBJECTIVE

To understand the molecular mechanisms of adventitious root formation from stem cuttings in sweet potato.

METHODS

RNA-seq analysis was performed using un-rooted stem (0 day) and rooted stem (3 days). Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, comparison with Arabidopsis transcription factors (TFs) of DEGs were conducted to investigate the characteristics of genes and TFs involved in root formation. In addition, qRT-PCR analysis using roots at 0, 3, 6, 9, and 12 days after planting was performed to confirm RNA-seq reliability and related genes expression.

RESULTS

42,459 representative transcripts and 2092 DEGs were obtained through the RNA-seq analysis. The DEGs indicated the GO terms related to the single-organism metabolic process and cell periphery, and involved in the biosynthesis of secondary metabolites, and phenylpropanoid biosynthesis in KEGG pathways. The comparison with Arabidopsis thaliana TF database showed that 3 TFs (WRKY, NAC, bHLH) involved in root formation of sweet potato. qRT-PCR analysis, which was conducted to confirm the reliability of RNA-seq analysis, indicated that some metabolisms including oxidative stress and wounding, transport, hormone may be involved in adventitious root formation.

CONCLUSIONS

The detected genes related to secondary metabolism, some hormone (auxin, gibberellin), transports, etc. and 3 TFs (WRKY, NAC, bHLH) may have functions in adventitious roots formation. This results provide valuable resources for future research on the adventitious root formation of sweet potato.

Transcriptome sequencing and whole genome expression profiling of hexaploid sweetpotato under salt stress

BACKGROUND

Purple-fleshed sweetpotato (PFSP) is one of the most important crops in the word which helps to bridge the food gap and contribute to solve the malnutrition problem especially in developing countries. Salt stress is seriously limiting its production and distribution. Due to lacking of reference genome, transcriptome sequencing is offering a rapid approach for crop improvement with promising agronomic traits and stress adaptability.

RESULTS

Five cDNA libraries were prepared from the third true leaf of hexaploid sweetpotato at seedlings stage (Xuzi-8 cultivar) treated with 200 mM NaCl for 0, 1, 6, 12, 48 h. Using second and third generation technology, Illumina sequencing generated 170,344,392 clean high-quality long reads that were assembled into 15,998 unigenes with an average length 2178 base pair and 96.55% of these unigenes were functionally annotated in the NR protein database. A number of 537 unigenes failed to hit any homologs which may be considered as novel genes. The current results indicated that sweetpotato plants behavior during the first hour of salt stress was different than the other three time points. Furthermore, expression profiling analysis identified 4, 479, 281, 508 significantly expressed unigenes in salt stress treated samples at the different time points including 1, 6, 12, 48 h, respectively as compared to control. In addition, there were 4, 1202, 764 and 2195 transcription factors differentially regulated DEGs by salt stress at different time points including 1, 6, 12, 48 h of salt stress. Validation experiment was done using 6 randomly selected unigenes and the results was in agree with the DEG results. Protein kinases include many genes which were found to play a vital role in phosphorylation process and act as a signal transductor/ receptor proteins in membranes. These findings suggest that salt stress tolerance in hexaploid sweetpotato plants may be mainly affected by TFs, PKs, Protein Detox and hormones related genes which contribute to enhance salt tolerance.

CONCLUSION

These transcriptome sequencing data of hexaploid sweetpotato under salt stress conditions can provide a valuable resource for sweetpotato breeding research and focus on novel insights into hexaploid sweetpotato responses to salt stress. In addition, it offers new candidate genes or markers that can be used as a guide to the future studies attempting to breed salt tolerance sweetpotato cultivars.

Transcriptome landscape of Rafflesia cantleyi floral buds reveals insights into the roles of transcription factors and phytohormones in flower development

Rafflesia possesses unique biological features and known primarily for producing the world’s largest and existing as a single flower. However, to date, little is known about key regulators participating in Rafflesia flower development. In order to further understand the molecular mechanism that regulates Rafflesia cantleyi flower development, RNA-seq data from three developmental stages of floral bud, representing the floral organ primordia initiation, floral organ differentiation, and floral bud outgrowth, were analysed. A total of 89,890 transcripts were assembled of which up to 35% could be annotated based on homology search. Advanced transcriptome analysis using K-mean clustering on the differentially expressed genes (DEGs) was able to identify 12 expression clusters that reflect major trends and key transitional states, which correlate to specific developmental stages. Through this, comparative gene expression analysis of different floral bud stages identified various transcription factors related to flower development. The members of WRKY, NAC, bHLH, and MYB families are the most represented among the DEGs, suggesting their important function in flower development. Furthermore, pathway enrichment analysis also revealed DEGs that are involved in various phytohormone signal transduction events such as auxin and auxin transport, cytokinin and gibberellin biosynthesis. Results of this study imply that transcription factors and phytohormone signalling pathways play major role in Rafflesia floral bud development. This study provides an invaluable resource for molecular studies of the flower development process in Rafflesia and other plant species.

Improving lipid production by strain development in microalgae: Strategies, challenges and perspectives

Over the past decade, the number of original articles and reviews presenting microalgae as a promising feedstock for biodiesel has increased tremendously. Many improvements of microalgae have been achieved through selection and strain development for industrial applications. However, the large-scale production of lipids for commercialization is not yet realistic because the production is still much more expensive than that of agricultural products. This review summarizes recent research on the induction of lipid biosynthesis in microalgae and the various strategies of genetic and metabolic engineering for enhancing lipid production. Strain engineering targets are proposed based on these strategies. To address current limitations of strain engineering for lipid production, this review provides insights on recent engineering strategies based on molecular tools and methods, and also discusses further perspectives.

Grouper: graph-based clustering and annotation for improved de novo transcriptome analysis

Motivation

De novo transcriptome analysis using RNA-seq offers a promising means to study gene expression in non-model organisms. Yet, the difficulty of transcriptome assembly means that the contigs provided by the assembler often represent a fractured and incomplete view of the transcriptome, complicating downstream analysis. We introduce Grouper, a new method for clustering contigs from de novo assemblies that are likely to belong to the same transcripts and genes; these groups can subsequently be analyzed more robustly. When provided with access to the genome of a related organism, Grouper can transfer annotations to the de novo assembly, further improving the clustering.

Results

On de novo assemblies from four different species, we show that Grouper is able to accurately cluster a larger number of contigs than the existing state-of-the-art method. The Grouper pipeline is able to map greater than 10% more reads against the contigs, leading to accurate downstream differential expression analyses. The labeling module, in the presence of a closely related annotated genome, can efficiently transfer annotations to the contigs and use this information to further improve clustering. Overall, Grouper provides a complete and efficient pipeline for processing de novo transcriptomic assemblies.

Availability and implementation

The Grouper software is freely available at https://github.com/COMBINE-lab/grouper under the 2-clause BSD license.

Supplementary information

Supplementary data are available at Bioinformatics online.

Interspecific hybridisation and LTR-retrotransposon mobilisation-related structural variation in plants: A case study

The dynamics of long-terminal-repeat retrotransposons in two poplar species (Populus deltoides and P. nigra) and in an interspecific hybrid, recently synthesized, were investigated by analyzing the genomic abundance and transcription levels of a collection of 828 full-length retroelements identified in the genome sequence of P. trichocarpa, all occurring also in the genomes of P. deltoides and P. nigra. Overall, genomic abundance and transcription levels of many retrotransposons in the hybrid resulted higher or lower than expected by calculating the mean of the parental values. A bioinformatics procedure was established to ascertain the occurrence of the same retrotransposon loci in the three genotypes. The results indicated that retrotransposon abundance variations between the hybrid and the mean value of the parents were due to i) co-segregation of retrotransposon high- or low-abundant haplotypes; ii) new retroelement insertions; iii) retrotransposon loss. Concerning retrotransposon expression, this was generally low, with only 14/828 elements over- or under-expressed in the hybrid than expected by calculating the mean of the parents. It is concluded that interspecific hybridisation between the two poplar species determine quantitative variation and differential expression of some retrotransposons, with possible consequences for the genetic differentiation of the hybrid.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Transcriptome analysis and codominant markers development in caper, a drought tolerant orphan crop with medicinal value

Caper (Capparis spinosa L.) is a xerophytic shrub cultivated for its flower buds and fruits, used as food and for their medicinal properties. Breeding programs and even proper taxonomic classification of the genus Capparis has been hampered so far by the lack of reliable genetic information and molecular markers. Here, we present the first genomic resource for C. spinosa, generated by transcriptomic approach and de novo assembly. The sequencing effort produced nearly 80 million clean reads assembled into 124,723 unitranscripts. Careful annotation and comparison with public databases revealed homologs to genes with a key role in important metabolic pathways linked to abiotic stress tolerance and bio-compounds production, such purine, thiamine and phenylpropanoid biosynthesis, α-linolenic acid and lipid metabolism. Additionally, a panel of genes involved in stomatal development/distribution and encoding for Stress Associated Proteins (SAPs) was also identified. We also used the transcriptomic data to uncover novel molecular markers for caper. Out of 50 SSRs tested, 14 proved polymorphic and represent the first set of SSR markers for the genus Capparis. This transcriptome will be an important contribution to future studies and breeding programs for this orphan crop, aiding to the development of improved varieties to sustain agriculture in arid conditions.

Molecular cloning, expression profiling, and functional analysis of a broad-complex isoform 2/3 (Br-Z2/Z3) transcription factor in the diamondback moth, Plutella xylostella (L.)

The diamondback moth, Plutella xylostella (L.) (Lepidoptera: Plutellidae), is a widespread and destructive pest of cruciferous crops. New strategies for controlling it are needed because it is rapidly developing resistance to conventional pesticides. In insects, transcription factors (TFs) including broad-complex (Br-C) are thought to be useful for insecticide development because they are able to regulate the transcription of functional genes involved in responses to external stimuli including insecticides. In the present study, we cloned and sequenced the open reading frames (ORFs) of three BTB-ZF encoding genes from the diamondback moth deposited in the National Center for Biotechnology Information (NCBI) database under accessions MG753773, MG288674, and MG753772. The lengths of these ORFs were 1,680, 1,428, and 1,647 bp, respectively. The phylogenetic analysis based on the predicted amino acid sequences of ZF domains showed that MG753773 and MG288674 belonged to Z2/Z3 and Z7 of Br-C while MG753772 belonged to Ttk types. In the agreement, the highest expression level of MG753773 occurred during the prepupal stage, MG288674 and MG753772 were expressed during all stages and peaked in the adult and egg stages, respectively. RNA interference silencing of MG753773 in the late third instar larvae significantly decreased survival and pupation of the insects. With precocene II, transcription of MG753773 increased (4×) in the fourth instar larva 24 hr later; 48 hr later the rate of prepupation and pupation was significantly higher. These findings will contribute to the development of new regulators of the growth and development for diamondback moth control.

De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through Illumina sequencing

Curcuma alismatifolia widely used as an ornamental plant in Thailand and Cambodia. This species of herbaceous perennial from the Zingiberaceae family, includes cultivars with a wide range of colours and long postharvest life, and is used as an ornamental cut flower, as a potted plant, and in exterior landscapes. For further genetic improvement, however, little genomic information and no specific molecular markers are available. The present study used Illumina sequencing and de novo transcriptome assembly of two C. alismatifolia cvs, 'Chiang Mai Pink' and 'UB Snow 701', to develop simple sequence repeat markers for genetic diversity studies. After de novo assembly, 62,105 unigenes were generated and 48,813 (78.60%) showed significant similarities versus six functional protein databases. In addition, 9,351 expressed sequence tag-simple sequence repeats (EST-SSRs) were identified with a distribution frequency of 12.5% total unigenes. Out of 8,955 designed EST-SSR primers, 150 primers were selected for the development of potential molecular markers. Among these markers, 17 EST-SSR markers presented a moderate level of genetic diversity among three C. alismatifolia cultivars, one hybrid, three Curcuma, and two Zingiber species. Three different genetic groups within these species were revealed using EST-SSR markers, indicating that the markers developed in this study can be effectively applied to the population genetic analysis of Curcuma and Zingiber species. This report describes the first analysis of transcriptome data of an important ornamental ginger cultivars, also provides a valuable resource for gene discovery and marker development in the genus Curcuma.

De novo transcriptome sequencing and SSR markers development for Cedrela balansae C.DC., a native tree species of northwest Argentina

The endangered Cedrela balansae C.DC. (Meliaceae) is a high-value timber species with great potential for forest plantations that inhabits the tropical forests in Northwestern Argentina.Research on this species is scarce because of the limited genetic and genomic information available. Here, we explored the transcriptome of C. balansae using 454 GS FLX Titanium next-generation sequencing (NGS) technology. Following de novo assembling, we identified 27,111 non-redundant unigenes longer than 200 bp, and considered these transcripts for further downstream analysis. The functional annotation was performed searching the 27,111 unigenes against the NR-Protein and the Interproscan databases. This analysis revealed 26,977 genes with homology in at least one of the Database analyzed. Furthermore, 7,774 unigenes in 142 different active biological pathways in C. balansae were identified with the KEGG database. Moreover, after in silico analyses, we detected 2,663 simple sequence repeats (SSRs) markers. A subset of 70 SSRs related to important "stress tolerance" traits based on functional annotation evidence, were selected for wet PCR-validation in C. balansae and other Cedrela species inhabiting in northwest and northeast of Argentina (C. fissilis, C. saltensis and C. angustifolia). Successful transferability was between 77% and 93% and thanks to this study, 32 polymorphic functional SSRs for all analyzed Cedrela species are now available. The gene catalog and molecular markers obtained here represent a starting point for further research, which will assist genetic breeding programs in the Cedrela genus and will contribute to identifying key populations for its preservation.

De novo sequencing and transcriptome assembly of Arisaema heterophyllum Blume and identification of genes involved in isoflavonoid biosynthesis

Arisaema heterophyllum Blume (AhBl) is one of the valued medicinal plants. However, its genetic information is limited, which impedes further studies of this valuable resource. To investigate the genes involved in the isoflavonoid biosynthesis, we deeply performed transcriptome sequencing for AhBl. An average of 10.98 Gb clean reads were obtained based on root, tuber and leaf tissues, and 109,937 unigenes were yielded after de novo assembly. In total, 72,287 of those unigenes were annotated in at least one public database. The numbers of expressed unigenes in each tissue were 35,686, 43,363 and 47,783, respectively. The overall expression levels of transcripts in leaf were higher than those in root and tuber. Differentially expressed genes analysis indicated that a total of 12,448 shared unigenes were detected in all three tissues, 10,215 of which were higher expressed in tuber than that in root and leaf. Besides, 87 candidate unigenes that encode for enzymes involved in biosynthesis of isoflavonoid were identified and analyzed, and some key enzyme genes were experimentally validated by quantitative Real-Time PCR (qRT-PCR). This study provides a unique dataset for the systematic analysis of AhBl functional genes and expression characteristics, and facilitates the future study of the pharmacological mechanism of AhBl.

Generation and classification of transcriptomes in two Croomia species and molecular evolution of CYC/TB1 genes in Stemonaceae

The genus Croomia (Stemonaceae) is an excellent model for studying the evolution of the Eastern Asia (EA)-Eastern North America (ENA) floristic disjunction and the genetic mechanisms of floral zygomorphy formation. In addition to the presence of both actinomorphic and zygomorphic flowers within the genus, species are disjunctively distributed between EA and ENA. However, due to the limited availability of genomic resources, few studies of Croomia have examined these questions. In this study, we sequenced the floral and leaf transcriptomes of the zygomorphic flowered C roomia heterosepala and the actinomorphic flowered Croomia japonica, and used comparative genomic approaches to investigate the transcriptome evolution of the two closely related species. The sequencing and de novo assembly of transcriptomes from flowers of C. heterosepala (ChFlower), flowers of C. japonica (CjFlower), and leaves of C. japonica (CjLeaf) yielded 57,193, 62,131 and 64,448 unigenes, respectively. In addition, estimation of Ka/Ks ratios for 11,566 potential orthologous groups between ChFlower and CjFlower revealed that only six pairs had Ka/Ks ratios significantly greater than 1 and are likely under positive selection. A total of 429 single copy nuclear genes (SCNGs) and 21,460 expression sequence tags-simple sequence repeats (EST-SSRs) were identified in this study. Specifically, we identified seven CYC/TB1-like genes from Stemonaceae. Phylogenetic and molecular evolution analyses indicated that these CYC/TB1-like genes formed a monophyletic clade (SteTBL1) and were subject to strong purifying selection. The shifts of floral symmetry in Stemonaceae do not appear to be correlated with TBL copy number.

Good Citizenship Made Easy: A Step-by-Step Guide to Submitting RNA-Seq Data to NCBI

The analysis of transcriptome data from non-model organisms contributes to our understanding of diverse aspects of evolutionary biology, including developmental processes, speciation, adaptation, and extinction. Underlying this diversity is one shared feature, the generation of enormous amounts of sequence data. Data availability requirements in most journals oblige researchers to make their raw transcriptome data publicly available, and the databases housed at the National Center for Biotechnology Information (NCBI) are a popular choice for data deposition. Unfortunately, the successful submission of raw sequences to the Sequence Read Archive (SRA) and transcriptome assemblies to the Transcriptome Shotgun Assembly (TSA) can be challenging for novice users, significantly delaying data availability and publication. Here we present two comprehensive protocols for submitting RNA-Seq data to NCBI databases, accompanied by an easy-to-use website that facilitates the timely submission of data by researchers of any experience level. © 2018 by John Wiley & Sons, Inc.

Integrated tool for microsatellite isolation and validation from the reference genome and their application in the study of breeding turnover in an endangered avian population

Accurate individual identification is required to estimate survival rates in avian populations. For endangered species, non-invasive methods of obtaining individual identification, such as using molted feathers as a source of DNA for microsatellite markers, are preferred because of less disturbance, easy sample preparation and high efficiency. With the availability of many avian genomes, a few pipelines isolating genome-wide microsatellites have been published, but it is still a challenge to isolate microsatellites from the reference genome efficiently. Here, we have developed an integrated tool comprising a bioinformatic pipeline and experimental procedures for microsatellite isolation and validation based on the reference genome. We have identified over 95 000 microsatellite loci and established a system comprising 10 highly polymorphic markers (PIC value: 0.49-0.93, mean: 0.79) for an endangered species, saker falcon (Falco cherrug). These markers (except 1) were successfully amplified in 126 molted feathers, exhibiting high amplification success rates (83.9-99.7%), high quality index (0.90-0.97) and low allelic dropout rates (1-9.5%). To further assess the efficiency of this marker system in a population study, we identified individual sakers using these molted feathers (adult) and 146 plucked feathers (offspring). The use of parent and offspring samples enabled us to infer the genotype of missing samples (N = 28), and all adult genotypes were used to ascertain that breeding turnover is a useful proxy for survival estimation in sakers. Our study presents a cost-effective tool for microsatellite isolation based on publicly available reference genomes and demonstrates the power of this tool in estimating key parameters of avian population dynamics.

Transcriptome sequencing of Pinus kesiya var. langbianensis and comparative analysis in the Pinus phylogeny

Background

Pines are widely distributed in the Northern Hemisphere and have a long evolutionary history. The availability of transcriptome data has facilitated comparative transcriptomics for studying the evolutionary patterns associated with the different geographical distributions of species in the Pinus phylogeny.

Results

The transcriptome of Pinus kesiya var. langbianensis was sequenced using the Illumina HiSeq 2000 platform, and a total of 68,881 unigenes were assembled by Trinity. Transcriptome sequences of another 12 conifer species were downloaded from public databases. All of the pairwise orthologues were identified by comparative transcriptome analysis in 13 conifer species, from which the rate of diversification was calculated and a phylogenetic tree inferred. All of the fast-evolving positive selection sequences were identified, and some salt-, drought-, and abscisic acid-resistance genes were discovered.

Conclusions

mRNA sequences of P. kesiya var. langbianensis were obtained by transcriptome sequencing, and a large number of simple sequence repeat and short nucleotide polymorphism loci were detected. These data can be used in molecular marker-assisted selected in pine breeding. Divergence times were estimated in the 13 conifer species using comparative transcriptomic analysis. A number of positive selection genes were found to be related to environmental factors. Salt- and abscisic acid-related genes exhibited different selection patterns between coastal and inland Pinus. Our findings help elucidate speciation patterns in the Pinus lineage.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5127-6) contains supplementary material, which is available to authorized users.

Evolutionary insights from comparative transcriptome and transcriptome-wide coalescence analyses in Tetrastigma hemsleyanum

BACKGROUND

Tetrastigma hemsleyanum is of great medicinal importance and used as a model system to address the evolutionary history of warm-temperate evergreen (WTE) forest biomes in East Asia over Neogene time scales. However, further studies on the neutral and adaptive divergence processes of T. hemsleyanum are currently impeded by a lack of genomic resources. In this study, we de novo assembled and annotated a reference transcriptome for two cpDNA lineages (Central-South-East vs. Southwest) of T. hemsleyanum. We further used comparative genomic and multilocus coalescent approaches to investigate the tempo and mode of lineage diversification in T. hemsleyanum.

RESULTS

A total of 52,838 and 65,197 unigenes with an N50 of 1,667 and 1,841 bp for Central-South-East (CSE) and Southwest (SW) lineages, respectively, were recovered, and 6,692 putative orthologs were identified between the two lineages. Estimation of K_a/K_s ratios for these orthologs revealed that ten genes had K_a/K_s values significantly greater than 0.5 (P < 0.05), whereas 2,099 (K_a/K_s < 0.5, P < 0.05) were inferred to be under purifying selection. Based on three bioinformatic strategies, we identified a total of 1,018 single-copy nuclear genes (SCNGs) from the orthologs. We successfully designed eight nuclear gene primer pairs with high intraspecific variation (e.g. h_T = 0.923, π_T = 1.68×10^-3), when surveyed across a subset of T. hemsleyanum individuals. Concordant with the previous cpDNA data, the haplotype networks constructed for most nuclear gene loci clearly identified the two lineages. A multilocus coalescence analysis suggested that the separation between the two lineages appears to have occurred during the mid-Pliocene. Despite their ancient divergence, both lineages experienced expansion at rather localized scales and have continued to exchange genes at a low rate.

CONCLUSIONS

This study demonstrated the utility of transcriptome sequencing as a basis for SCNG development in non-model species and the advantages of integrating multiple nuclear loci for phylogeographic and phylogenetic studies.

Transcriptomic Profiling and Gene Disruption Revealed that Two Genes Related to PUFAs/DHA Biosynthesis May be Essential for Cell Growth of Aurantiochytrium sp

Aurantiochytrium sp. PKU#SW7 is a thraustochytrid strain that was found to exhibit high potential for docosahexaenoic acid (DHA, C22:6n-3) production. In this work, the transcriptome of Aurantiochytrium sp. PKU#SW7 was analyzed for the study of genes involved in basic metabolic functions and especially in the mechanisms of DHA biosynthesis. Sequence annotation and functional analysis revealed that the strain contains components of fatty acid synthesis (FAS) and polyketide synthase (PKS) pathways. Fatty acid desaturases and elongases were identified as components of FAS pathway, whilst key components of PKS pathway were also found in the cDNA library. The relative contribution of the two pathways to the synthesis of DHA was unknown, as both pathways appeared to be lacking full complement of genes for standalone synthesis of DHA. Further analysis of two putative genes encoding the very-long-chain (3R)-3-hydroxyacyl-CoA dehydratase and dehydrase/isomerase involved in FAS and PKS pathways, respectively, revealed that under various salinity conditions, their relative expression levels changed corresponding to the variation of DHA content in Aurantiochytrium sp. Independent knock outs of these genes in Aurantiochytrium sp. resulted in poor cell growth, probably due to little or no intracellular DHA accumulation. Hence, it can be speculated that both genes are engaged in DHA biosynthesis and DHA in Aurantiochytrium sp. could be produced by jointed actions of both FAS and PKS systems.

Comparative Transcriptomics Among Four White Pine Species

Conifers are the dominant plant species throughout the high latitude boreal forests as well as some lower latitude temperate forests of North America, Europe, and Asia. As such, they play an integral economic and ecological role across much of the world. This study focused on the characterization of needle transcriptomes from four ecologically important and understudied North American white pines within the Pinus subgenus Strobus. The populations of many Strobus species are challenged by native and introduced pathogens, native insects, and abiotic factors. RNA from the needles of western white pine (Pinus monticola), limber pine (Pinus flexilis), whitebark pine (Pinus albicaulis), and sugar pine (Pinus lambertiana) was sampled, Illumina short read sequenced, and de novo assembled. The assembled transcripts and their subsequent structural and functional annotations were processed through custom pipelines to contend with the challenges of non-model organism transcriptome validation. Orthologous gene family analysis of over 58,000 translated transcripts, implemented through Tribe-MCL, estimated the shared and unique gene space among the four species. This revealed 2025 conserved gene families, of which 408 were aligned to estimate levels of divergence and reveal patterns of selection. Specific candidate genes previously associated with drought tolerance and white pine blister rust resistance in conifers were investigated.

A transcriptome screen for positive selection in domesticated breadfruit and its wild relatives (Artocarpus spp.)

PREMISE OF THE STUDY

Underutilized crops, such as breadfruit (Artocarpus altilis, Moraceae) have the potential to improve global food security. Humans have artificially selected many cultivars of breadfruit since its domestication began approximately 3500 years ago. The goal of this research was to identify transcriptomic signals of positive selection and to develop genomic resources that may facilitate the development of improved breadfruit cultivars in the future.

METHODS

A reference transcriptome of breadfruit was assembled de novo and annotated. Twenty-four transcriptomes of breadfruit and its wild relatives were generated and analyzed to reveal signals of positive selection that may have resulted from local adaptation or natural selection. Emphasis was placed on MADS-box genes, which are important because they often regulate fruiting timing and structures, and on carotenoid biosynthesis genes, which can impact the nutritional quality of the fruit.

KEY RESULTS

Over 1000 genes showed signals of positive selection, and these genes were enriched for localization to plastids. Nucleotide sites and individuals under positive selection were discovered in MADS-box genes and carotenoid biosynthesis genes, with several sites located in cofactor or DNA-binding domains. A McDonald-Kreitman test comparing wild to cultivated samples revealed selection in one of the carotenoid biosynthesis genes, abscisic acid 8'-hydroxylase 3.

CONCLUSIONS

This research highlights some of the many genes that may have been intentionally or unintentionally selected for during the human-mediated dispersal of breadfruit and stresses the importance of conserving a varied germplasm collection. It has revealed candidate genes for further study and produced new genomic resources for breadfruit.

Generation and characterization of the blood transcriptome of Macaca thibetana and comparative analysis with M. mulatta

Transcriptome profiles provide a large transcript sequence data set for genomic study, particularly in organisms that have no accurate genome data published. The Tibetan macaque (Macaca thibetana) is commonly considered to be an endemic species to China and an important animal in biomedical research in the present day. In the present study, we report the de novo assembly and characterization of the blood transcriptome of the Tibetan macaque from three individuals, and we also sequenced the blood transcriptome of the rhesus macaque (Macaca mulatta) for comparison. Using RNA-seq technology, 138 million sequencing reads of the M. thibetana transcriptome were generated. The assembly included 327 871 transcripts with an N50 of 1571 bp. According to the sequence similarity search, 80 317 (24.5%) transcripts were annotated in the nr protein database. All transcripts from M. thibetana and M. mulatta were functionally classified and compared using GO and KEGG analyses. The two transcriptomes were different in the GO term of nutrient reservoir activity, and in the KEGG subcategories of signaling molecules and interaction, infectious diseases, cell growth and death, and immune system. The transcriptomes in this study would provide a valuable resource for future functional and comparative genomic studies, and even for biological studies of this non-human primate.

Development and Application of Genomic Resources in an Endangered Palaeoendemic Tree, Parrotia subaequalis (Hamamelidaceae) From Eastern China

Parrotia subaequalis is an endangered palaeoendemic tree from disjunct montane sites in eastern China. Due to the lack of effective genomic resources, the genetic diversity and population structure of this endangered species are not clearly understood. In this study, we conducted paired-end shotgun sequencing (2 × 125 bp) of genomic DNA for two individuals of P. subaequalis on the Illumina HiSeq platform. Based on the resulting sequences, we have successfully assembled the complete chloroplast genome of P. subaequalis, as well as identified the polymorphic chloroplast microsatellites (cpSSRs), nuclear microsatellites (nSSRs) and mutational hotspots of chloroplast. Ten polymorphic cpSSR loci and 12 polymorphic nSSR loci were used to genotype 96 individuals of P. subaequalis from six populations to estimate genetic diversity and population structure. Our results revealed that P. subaequalis exhibited abundant genetic diversity (e.g., cpSSRs: Hcp = 0.862; nSSRs: HT = 0.559) and high genetic differentiation (e.g., cpSSRs: RST = 0.652; nSSRs: RST = 0.331), and characterized by a low pollen-to-seed migration ratio (r ≈ 1.78). These genetic patterns are attributable to its long evolutionary histories and low levels of contemporary inter-population gene flow by pollen and seed. In addition, lack of isolation-by-distance pattern and strong population genetic structuring in both marker systems, suggests that long-term isolation and/or habitat fragmentation as well as genetic drift may have also contributed to the geographic differentiation of P. subaequalis. Therefore, long-term habitat protection is the most important methods to prevent further loss of genetic variation and a decrease in effective population size. Furthermore, both cpSSRs and nSSRs revealed that P. subaequalis populations consisted of three genetic clusters, which should be considered as separated conservation units.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Transcriptional profiling of cork oak phellogenic cells isolated by laser microdissection

The phenylpropanoid pathway impacts the cork quality development. In cork of bad quality, the flavonoid route is favored, whereas in good quality, cork lignin and suberin production prevails. Cork oaks develop a thick cork tissue as a protective shield that results of the continuous activity of a secondary meristem, the cork cambium, or phellogen. Most studies applied to developmental processes do not consider the cell types from which the samples were extracted. Here, laser microdissection (LM) coupled with transcript profiling using RNA sequencing (454 pyrosequencing) was applied to phellogen cells of trees producing low- and good quality cork. Functional annotation and functional enrichment analyses showed that stress-related genes are enriched in samples extracted from trees producing good quality cork (GQC). This process is under tight transcriptional (transcription factors, kinases) regulation and also hormonal control involving ABA, ethylene, and auxins. The phellogen cells collected from trees producing bad quality cork (BQC) show a consistent up-regulation of genes belonging to the flavonoid pathway as a response to stress. They also display a different modulation of cell wall genes resulting into a thinner cork layer, i.e., less meristematic activity. Based on the analysis of the phenylpropanoid pathway regulating genes, in GQC, the synthesis of lignin and suberin is promoted, whereas in BQC, the same pathway favors the biosynthesis of free phenolic compounds. This study provided new insights of how cell-specific gene expression can determine tissue and organ morphology and physiology and identified robust candidate genes that can be used in breeding programs aiming at improving cork quality.

Genomic abundance and transcriptional activity of diverse gypsy and copia long terminal repeat retrotransposons in three wild sunflower species

BACKGROUND

Long terminal repeat (LTR) retrotransposons are highly abundant in plant genomes and require transcriptional activity for their proliferative mode of replication. These sequences exist in plant genomes as diverse sublineages within the main element superfamilies (i.e., gypsy and copia). While transcriptional activity of these elements is increasingly recognized as a regular attribute of plant transcriptomes, it is currently unknown the extent to which different sublineages of these elements are transcriptionally active both within and across species. In the current report, we utilize next generation sequencing methods to examine genomic copy number abundance of diverse LTR retrotransposon sublineages and their corresponding levels of transcriptional activity in three diploid wild sunflower species, Helianthus agrestis, H. carnosus and H. porteri.

RESULTS

The diploid sunflower species under investigation differ in genome size 2.75-fold, with 2C values of 22.93 for H. agrestis, 12.31 for H. carnosus and 8.33 for H. porteri. The same diverse gypsy and copia sublineages of LTR retrotransposons were identified across species, but with gypsy sequences consistently more abundant than copia and with global gypsy sequence abundance positively correlated with nuclear genome size. Transcriptional activity was detected for multiple copia and gypsy sequences, with significantly higher activity levels detected for copia versus gypsy. Interestingly, of 11 elements identified as transcriptionally active, 5 exhibited detectable expression in all three species and 3 exhibited detectable expression in two species.

CONCLUSIONS

Combined analyses of LTR retrotransposon genomic abundance and transcriptional activity across three sunflower species provides novel insights into genome size evolution and transposable element dynamics in this group. Despite considerable variation in nuclear genome size among species, relatively conserved patterns of LTR retrotransposon transcriptional activity were observed, with a highly overlapping set of copia and gypsy sequences observed to be transcriptionally active across species. A higher proportion of copia versus gypsy elements were found to be transcriptionally active and these sequences also were expressed at higher levels.