De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq

Evolution of glucuronoxylan side chain variability in vascular plants and the compensatory adaptations of cell wall-degrading hydrolases

Polysaccharide structural complexity not only influences cell wall strength and extensibility but also hinders pathogenic and biotechnological attempts to saccharify the wall. In certain species and tissues, glucuronic acid side groups on xylan exhibit arabinopyranose or galactose decorations whose genetic and evolutionary basis is completely unknown, impeding efforts to understand their function and engineer wall digestibility. Genetics and polysaccharide profiling were used to identify the responsible loci in Arabidopsis and Eucalyptus from proposed candidates, while phylogenies uncovered a shared evolutionary origin. GH30-family endo-glucuronoxylanase activities were analysed by electrophoresis, and their differing specificities were rationalised by phylogeny and structural analysis. The newly identified xylan arabinopyranosyltransferases comprise an overlooked subfamily in the GT47-A family of Golgi glycosyltransferases, previously assumed to comprise mainly xyloglucan galactosyltransferases, highlighting an unanticipated adaptation of both donor and acceptor specificities. Further neofunctionalisation has produced a Myrtaceae-specific xylan galactosyltransferase. Simultaneously, GH30 endo-glucuronoxylanases have convergently adapted to overcome these decorations, suggesting a role for these structures in defence. The differential expression of glucuronoxylan-modifying genes across Eucalyptus tissues, however, hints at further functions. Our results demonstrate the rapid adaptability of biosynthetic and degradative carbohydrate-active enzyme activities, providing insight into plant-pathogen interactions and facilitating plant cell wall biotechnological utilisation.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

De novo sequencing, assembly, and characterization of Asparagus racemosus transcriptome and analysis of expression profile of genes involved in the flavonoid biosynthesis pathway

Asparagus racemosus is known for its diverse content of secondary metabolites, i.e., saponins, alkaloids, and a wide range of flavonoids. Flavonoids, including phenols and polyphenols, have a significant role in plant physiology and are synthesized in several tissues. Despite the diverse role of flavonoids, genetic information is limited for flavonoid biosynthesis pathways in A. racemosus. The current study explores full-scale functional genomics information of A. racemosus by de novo transcriptome sequencing using Illumina paired-end sequencing technology to elucidate the genes involved in flavonoid biosynthesis pathways. The de novo assembly of high-quality paired-end reads resulted in ∼2.3 million high-quality reads with a pooled transcript of 45,647 comprising ∼76 Mb transcriptome with a mean length (bp) of 1,674 and N50 of 1,868bp. Furthermore, the coding sequence (CDS) prediction analysis from 45,647 pooled transcripts resulted in 45,444 CDS with a total length and mean length of 76,398,686 and 1,674, respectively. The Gene Ontology (GO) analysis resulted in a high number of CDSs assigned to 25,342 GO terms, which grouped the predicted CDS into three main domains, i.e., Biological Process (19,550), Molecular Function (19,873), and Cellular Component (14,577). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to categorize 6,353 CDS into 25 distinct biological pathway categories, in which the majority of mapped CDS were shown to be related to translation (645), followed by signal transduction (532), carbohydrate metabolism (524), folding, sorting, and degradation (522). Among these, only ∼64 and 14 CDSs were found to be involved in the phenylpropanoid and flavonoid biosynthesis pathways, respectively. Quantitative Real-time PCR was used to check the expression profile of fourteen potential flavonoid biosynthesis pathway genes. The qRT-PCR analysis result matches the transcriptome sequence data validating the Illumina sequence results. Moreover, a large number of genes associated with the flavonoids biosynthesis pathway were found to be upregulated under the induction of methyl jasmonate. The present-day study on transcriptome sequence data of A. racemosus can be utilized for characterizing genes involved in flavonoid biosynthesis pathways and for functional genomics analysis in A. racemosus using the reverse genetics approach (CRISPR/Cas9 technology).

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses

Metal tolerance proteins (MTPs) as Me²⁺/H⁺(K⁺) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd²⁺ and Cu²⁺ exposure might be responsible for metal translocation from roots to leaves.

Haplogenome assembly reveals structural variation in Eucalyptus interspecific hybrids

BACKGROUND

De novo phased (haplo)genome assembly using long-read DNA sequencing data has improved the detection and characterization of structural variants (SVs) in plant and animal genomes. Able to span across haplotypes, long reads allow phased, haplogenome assembly in highly outbred organisms such as forest trees. Eucalyptus tree species and interspecific hybrids are the most widely planted hardwood trees with F1 hybrids of Eucalyptus grandis and E. urophylla forming the bulk of fast-growing pulpwood plantations in subtropical regions. The extent of structural variation and its effect on interspecific hybridization is unknown in these trees. As a first step towards elucidating the extent of structural variation between the genomes of E. grandis and E. urophylla, we sequenced and assembled the haplogenomes contained in an F1 hybrid of the two species.

FINDINGS

Using Nanopore sequencing and a trio-binning approach, we assembled the separate haplogenomes (566.7 Mb and 544.5 Mb) to 98.0% BUSCO completion. High-density SNP genetic linkage maps of both parents allowed scaffolding of 88.0% of the haplogenome contigs into 11 pseudo-chromosomes (scaffold N50 of 43.8 Mb and 42.5 Mb for the E. grandis and E. urophylla haplogenomes, respectively). We identify 48,729 SVs between the two haplogenomes providing the first detailed insight into genome structural rearrangement in these species. The two haplogenomes have similar gene content, 35,572 and 33,915 functionally annotated genes, of which 34.7% are contained in genome rearrangements.

CONCLUSIONS

Knowledge of SV and haplotype diversity in the two species will form the basis for understanding the genetic basis of hybrid superiority in these trees.

Co-overexpression of AtSAT1 and EcPAPR improves seed nutritional value in maize

Maize seeds synthesize insufficient levels of the essential amino acid methionine (Met) to support animal and livestock growth. Serine acetyltransferase1 (SAT1) and 3'-phosphoadenosine-5'-phosphosulfate reductase (PAPR) are key control points for sulfur assimilation into Cys and Met biosynthesis. Two high-MET maize lines pRbcS:AtSAT1 and pRbcS:EcPAPR were obtained through metabolic engineering recently, and their total Met was increased by 1.4- and 1.57-fold, respectively, compared to the wild type. The highest Met maize line, pRbcS:AtSAT1-pRbcS:EcPAPR, was created by stacking the two transgenes, causing total Met to increase 2.24-fold. However, the pRbcS:AtSAT1-pRbcS:EcPAPR plants displayed progressively severe defects in plant growth, including early senescence, stunting, and dwarfing, indicating that excessive sulfur assimilation has an adverse effect on plant development. To explore the mechanism of correlation between Met biosynthesis in maize leaves and storage proteins in developing endosperm, the transcriptomes of the sixth leaf at stage V9 and 18 DAP endosperm of pRbcS:AtSAT1, pRbcS:AtSAT1-pRbcS:EcPAPR, and the null segregants were quantified and analyzed. In pRbcS:AtSAT1-pRbcS:EcPAPR, 3274 genes in leaves (1505 up- and 1769 downregulated) and 679 genes in the endosperm (327 up- and 352 downregulated) were differentially expressed. Gene ontology (GO) and KEGG (Kyoto encyclopedia of genes and genomes) analyses revealed that many genes were associated with Met homeostasis, including transcription factors and genes involved in cysteine and Met metabolism, glutathione metabolism, plant hormone signal transduction, and oxidation-reduction. The data from gene network analysis demonstrated that two genes, serine/threonine-protein kinase (CCR3) and heat shock 70 kDa protein (HSP), were localized in the core of the leaves and endosperm regulation networks, respectively. The results of this study provide insights into the diverse mechanisms that underlie the ideal establishment of enhanced Met levels in maize seeds.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Genome sequencing and identification of cellulase genes in Bacillus paralicheniformis strains from the Red Sea

Background

Cellulolytic microorganisms are considered a key player in the degradation of plant biomass in various environments. These microorganisms can be isolated from various environments, such as soils, the insect gut, the mammalian rumen and oceans. The Red Sea exhibits a unique environment in terms of presenting a high seawater temperature, high salinity, low nutrient levels and high biodiversity. However, there is little information regarding cellulase genes in the Red Sea environment. This study aimed to examine whether the Red Sea can be a resource for the bioprospecting of microbial cellulases by isolating cellulase-producing microorganisms from the Red Sea environment and characterizing cellulase genes.

Results

Three bacterial strains were successfully isolated from the plankton fraction and the surface of seagrass. The isolated strains were identified as Bacillus paralicheniformis and showed strong cellulase activity. These results suggested that these three isolates secreted active cellulases. By whole genome sequencing, we found 10 cellulase genes from the three isolates. We compared the expression of these cellulase genes under cellulase-inducing and non-inducing conditions and found that most of the cellulase genes were generally upregulated during cellulolysis in the isolates. Our operon structure analysis also showed that cellulase genes form operons with genes involved in various kinds of cellular reactions, such as protein metabolism, which suggests the existence of crosstalk between cellulolysis and other metabolic pathways in the bacterial isolates. These results suggest that multiple cellulases are playing important roles in cellulolysis.

Conclusions

Our study reports the isolation and characterization of cellulase-producing bacteria from the Red Sea. Our whole-genome sequencing classified our three isolates as Bacillus paralicheniformis, and we revealed the presence of ten cellulase orthologues in each of three isolates’ genomes. Our comparative expression analysis also identified that most of the cellulase genes were upregulated under the inducing conditions in general. Although cellulases have been roughly classified into three enzyme groups of beta-glucosidase, endo-β-1,4-glucanase and exoglucanase, these findings suggest the importance to consider microbial cellulolysis as a more complex reaction with various kinds of cellulase enzymes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02316-w.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Fritillaria hupehensis

One of the most commonly utilized medicinal plants in China is Fritillaria hupehensis (Hsiao et K.C. Hsia). However, due to a lack of genomic resources, little is known about the biosynthesis of relevant compounds, particularly the flavonoid biosynthesis pathway. A PacBio RS II sequencing generated a total of 342,044 reads from the bulb, leaf, root, and stem, of which 316,438 were full-length (FL) non-redundant reads with an average length of 1365 bp and a N50 of 1888 bp. There were also 38,607 long non-coding RNAs and 7914 simple sequence repeats detected. To improve our understanding of processes implicated in regulating secondary metabolite biosynthesis in F. hupehensis tissues, we evaluated potential metabolic pathways. Overall, this study provides a repertoire of FL transcripts in F. hupehensis for the first time, and it will be a valuable resource for marker-assisted breeding and research into bioactive compounds for medicinal and pharmacological applications.

Genomic evidence of introgression and adaptation in a model subtropical tree species, Eucalyptus grandis

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses.

Metabolome and Transcriptome Analysis of Hexaploid Solidago canadensis Roots Reveals its Invasive Capacity Related to Polyploidy

Polyploid plants are more often invasive species than their diploid counterparts. As the invasiveness of a species is often linked to its production of allelopathic compounds, we hypothesize that differences in invasive ability between cytotypes may be due to their different ability to synthesize allelopathic metabolites. We test this using two cytotypes of Solidago canadensis as the model and use integrated metabolome and transcriptome data to resolve the question. Metabolome analysis identified 122 metabolites about flavonoids, phenylpropanoids and terpenoids, of which 57 were differentially accumulated between the two cytotypes. Transcriptome analysis showed that many differentially expressed genes (DEGs) were enriched in ‘biosynthesis of secondary metabolites’, ‘plant hormone signal transduction’, and ‘MAPK signaling’, covering most steps of plant allelopathic metabolite synthesis. Importantly, the differentially accumulated flavonoids, phenylpropanoids and terpenoids were closely correlated with related DEGs. Furthermore, 30 miRNAs were found to be negatively associated with putative targets, and they were thought to be involved in target gene expression regulation. These miRNAs probably play a vital role in the regulation of metabolite synthesis in hexaploid S. canadensis. The two cytotypes of S. canadensis differ in the allelopathic metabolite synthesis and this difference is associated with regulation of expression of a range of genes. These results suggest that changes in gene expression may underlying the increased invasive potential of the polyploidy.

Organellar carbon metabolism is coordinated with distinct developmental phases of secondary xylem

Subcellular compartmentation of plant biosynthetic pathways in the mitochondria and plastids requires coordinated regulation of nuclear encoded genes, and the role of these genes has been largely ignored by wood researchers. In this study, we constructed a targeted systems genetics coexpression network of xylogenesis in Eucalyptus using plastid and mitochondrial carbon metabolic genes and compared the resulting clusters to the aspen xylem developmental series. The constructed network clusters reveal the organization of transcriptional modules regulating subcellular metabolic functions in plastids and mitochondria. Overlapping genes between the plastid and mitochondrial networks implicate the common transcriptional regulation of carbon metabolism during xylem secondary growth. We show that the central processes of organellar carbon metabolism are distinctly coordinated across the developmental stages of wood formation and are specifically associated with primary growth and secondary cell wall deposition. We also demonstrate that, during xylogenesis, plastid-targeted carbon metabolism is partially regulated by the central clock for carbon allocation towards primary and secondary xylem growth, and we discuss these networks in the context of previously established associations with wood-related complex traits. This study provides a new resolution into the integration and transcriptional regulation of plastid- and mitochondrial-localized carbon metabolism during xylogenesis.

Snakin-1 affects reactive oxygen species and ascorbic acid levels and hormone balance in potato

Snakin-1 is a member of the Solanum tuberosum Snakin/GASA family. We previously demonstrated that Snakin-1 is involved in plant defense to pathogens as well as in plant growth and development, but its mechanism of action has not been completely elucidated yet. Here, we showed that leaves of Snakin-1 silenced potato transgenic plants exhibited increased levels of reactive oxygen species and significantly reduced content of ascorbic acid. Furthermore, Snakin-1 silencing enhanced salicylic acid content in accordance with an increased expression of SA-inducible PRs genes. Interestingly, gibberellic acid levels were also enhanced and transcriptome analysis revealed that a large number of genes related to sterol biosynthesis were downregulated in these silenced lines. Moreover, we demonstrated that Snakin-1 directly interacts with StDIM/DWF1, an enzyme involved in plant sterols biosynthesis. Additionally, the analysis of the expression pattern of PStSN1::GUS in potato showed that Snakin-1 is present mainly in young tissues associated with active growth and cell division zones. Our comprehensive analysis of Snakin-1 silenced lines demonstrated for the first time in potato that Snakin-1 plays a role in redox balance and participates in a complex crosstalk among different hormones.

The plastid and mitochondrial genomes of Eucalyptus grandis

BACKGROUND

Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes.

RESULTS

In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis.

CONCLUSIONS

The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

Identification and expression profiling of genes governing lignin biosynthesis in Casuarina equisetifolia L

Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.

Characterization of Brassinazole resistant (BZR) gene family and stress induced expression in Eucalyptus grandis

Brassinosteroids (BRs) are a group of plant hormones which play a pivotal role in modulating cell elongation, stress responses, vascular differentiation and senescence. In response to BRs, BRASSINAZOLE-RESISTANT (BZR) transcription factors (TFs) accumulate in the nucleus, where they modulate thousands of target genes and coordinate many biological processes, especially in regulating defense against biotic and abiotic stresses. In this study, 6 BZR TFs of Eucalyptus grandis (EgrBZR) from a genome-wide survey were characterized by sequence analysis and expression profiling against several abiotic stresses. The results showed that BZR gene family in Eucalyptus was slightly smaller compared to Populus and Arabidopsis, but all phylogenetic groups were represented. Various systematic in silico analysis of these TFs validated the basic properties of BZRs, whereas comparative studies showed a high degree of similarity with recognized BZRs of other plant species. In the organ-specific expression analyses, 4 EgrBZRs were expressed in vascular tissue indicating their possible functions in wood formation. Meanwhile, almost all EgrBZR genes showed differential transcript abundance levels in response to exogenously applied BR, MeJA, and SA, and salt and cold stresses. Besides, protein interaction analysis showed that all EgrBZR genes were associated with BR signaling directly or indirectly. These TFs were proposed as transcriptional activators or repressors of abiotic stress response and growth and development pathways of E. grandis by participating in BR signaling processes. These findings would be helpful in resolving the regulatory mechanism of EgrBZRs in stress resistance conditions but require further functional study of these potential TFs in Eucalyptus.

EUCANEXT: an integrated database for the exploration of genomic and transcriptomic data from Eucalyptus species

Tree species of the genus Eucalyptus are the most valuable and widely planted hardwoods in the world. Given the economic importance of Eucalyptus trees, much effort has been made towards the generation of specimens with superior forestry properties that can deliver high-quality feedstocks, customized to the industrýs needs for both cellulosic (paper) and lignocellulosic biomass production. In line with these efforts, large sets of molecular data have been generated by several scientific groups, providing invaluable information that can be applied in the development of improved specimens. In order to fully explore the potential of available datasets, the development of a public database that provides integrated access to genomic and transcriptomic data from Eucalyptus is needed. EUCANEXT is a database that analyses and integrates publicly available Eucalyptus molecular data, such as the E. grandis genome assembly and predicted genes, ESTs from several species and digital gene expression from 26 RNA-Seq libraries. The database has been implemented in a Fedora Linux machine running MySQL and Apache, while Perl CGI was used for the web interfaces. EUCANEXT provides a user-friendly web interface for easy access and analysis of publicly available molecular data from Eucalyptus species. This integrated database allows for complex searches by gene name, keyword or sequence similarity and is publicly accessible at http://www.lge.ibi.unicamp.br/eucalyptusdb. Through EUCANEXT, users can perform complex analysis to identify genes related traits of interest using RNA-Seq libraries and tools for differential expression analysis. Moreover, all the bioinformatics pipeline here described, including the database schema and PERL scripts, are readily available and can be applied to any genomic and transcriptomic project, regardless of the organism.

Database URL:http://www.lge.ibi.unicamp.br/eucalyptusdb

De novo sequencing, assembly and characterisation of Aloe vera transcriptome and analysis of expression profiles of genes related to saponin and anthraquinone metabolism

BACKGROUND

Aloe vera is a perennial, succulent, drought-resistant plant that exhibits many pharmacological characteristics such as wound healing ability against skin burns, anti-ulcer, anti-inflammatory, anti-tumor, anti-viral, anti-hypercholesterolemic, anti-hyperglycemic, anti-asthmatic and much more. Despite great medicinal worth, little genomic information is available on Aloe vera. This study is an initiative to explore the full-scale functional genomics of Aloe vera by generating whole transcriptome sequence database, using Illumina HiSeq technology and its progressive annotation specifically with respect to the metabolic specificity of the plant.

RESULTS

Transcriptome sequencing of root and leaf tissue of Aloe vera was performed using Illumina paired-end sequencing technology. De novo assembly of high quality paired-end reads, resulted into 1,61,733 and 2,21,792 transcripts with mean length of 709 and 714 nucleotides for root and leaf respectively. The non-redundant transcripts were clustered using CD-HIT-EST, yielding a total of 1,13,063 and 1,41,310 unigenes for root and leaf respectively. A total of 6114 and 6527 CDS for root and leaf tissue were enriched into 24 different biological pathway categories using KEGG pathway database. DGE profile prepared by calculating FPKM values was analyzed for differential expression of specific gene encoding enzymes involved in secondary metabolite biosynthesis. Sixteen putative genes related to saponin, lignin, anthraquinone, and carotenoid biosynthesis were selected for quantitative expression by real-time PCR. DGE as well as qRT PCR expression analysis represented up-regulation of secondary metabolic genes in root as compared to leaf. Furthermore maximum number of genes was found to be up-regulated after the induction of methyl jasmonate, which stipulates the association of secondary metabolite synthesis with the plant's defense mechanism during stress. Various transcription factors including bHLH, NAC, MYB were identified by searching predicted CDS against PlantTFdb.

CONCLUSIONS

This is the first transcriptome database of Aloe vera and can be potentially utilized to characterize the genes involved in the biosynthesis of important secondary metabolites, metabolic regulation, signal transduction mechanism, understanding function of a particular gene in the biology and physiology of plant of this species as well as other species of Aloe genus.

Genotyping-by-sequencing through transcriptomics: implementation in a range of crop species with varying reproductive habits and ploidy levels

The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker-dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost-effectiveness and sample throughput. Genotyping-by-sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS-transcriptomics (GBS-t) approach developed to be broadly applicable, cost-effective and high-throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS-t was validated as a simple and largely automated, cost-effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS-t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity.

Enriching Genomic Resources and Transcriptional Profile Analysis of Miscanthus sinensis under Drought Stress Based on RNA Sequencing

Miscanthus × giganteus is wildly cultivated as a potential biofuel feedstock around the world; however, the narrow genetic basis and sterile characteristics have become a limitation for its utilization. As a progenitor of M. × giganteus, M. sinensis is widely distributed around East Asia providing well abiotic stress tolerance. To enrich the M. sinensis genomic databases and resources, we sequenced and annotated the transcriptome of M. sinensis by using an Illumina HiSeq 2000 platform. Approximately 316 million high-quality trimmed reads were generated from 349 million raw reads, and a total of 114,747 unigenes were obtained after de novo assembly. Furthermore, 95,897 (83.57%) unigenes were annotated to at least one database including NR, Swiss-Prot, KEGG, COG, GO, and NT, supporting that the sequences obtained were annotated properly. Differentially expressed gene analysis indicates that drought stress 15 days could be a critical period for M. sinensis response to drought stress. The high-throughput transcriptome sequencing of M. sinensis under drought stress has greatly enriched the current genomic available resources. The comparison of DEGs under different periods of drought stress identified a wealth of candidate genes involved in drought tolerance regulatory networks, which will facilitate further genetic improvement and molecular studies of the M. sinensis.

De Novo Transcriptome Sequencing and the Hypothetical Cold Response Mode of Saussurea involucrata in Extreme Cold Environments

Saussurea involucrata grows in high mountain areas covered by snow throughout the year. The temperature of this habitat can change drastically in one day. To gain a better understanding of the cold response signaling pathways and molecular metabolic reactions involved in cold stress tolerance, genome-wide transcriptional analyses were performed using RNA-Seq technologies. A total of 199,758 transcripts were assembled, producing 138,540 unigenes with 46.8 Gb clean data. Overall, 184,416 (92.32%) transcripts were successfully annotated. The 365 transcription factors identified (292 unigenes) belonged to 49 transcription factor families associated with cold stress responses. A total of 343 transcripts on the signal transduction (132 upregulated and 212 downregulated in at least any one of the conditions) were strongly affected by cold temperature, such as the CBL-interacting serine/threonine-protein kinase (CIPKs), receptor-like protein kinases, and protein kinases. The circadian rhythm pathway was activated by cold adaptation, which was necessary to endure the severe temperature changes within a day. There were 346 differentially expressed genes (DEGs) related to transport, of which 138 were upregulated and 22 were downregulated in at least any one of the conditions. Under cold stress conditions, transcriptional regulation, molecular transport, and signal transduction were involved in the adaptation to low temperature in S. involucrata. These findings contribute to our understanding of the adaptation of plants to harsh environments and the survival traits of S. involucrata. In addition, the present study provides insight into the molecular mechanisms of chilling and freezing tolerance.

Root morphogenic pathways in Eucalyptus grandis are modified by the activity of protein arginine methyltransferases

BACKGROUND

Methylation of proteins at arginine residues, catalysed by members of the protein arginine methyltransferase (PRMT) family, is crucial for the regulation of gene transcription and for protein function in eukaryotic organisms. Inhibition of the activity of PRMTs in annual model plants has demonstrated wide-ranging involvement of PRMTs in key plant developmental processes, however, PRMTs have not been characterised or studied in long-lived tree species.

RESULTS

Taking advantage of the recently available genome for Eucalyptus grandis, we demonstrate that most of the major plant PRMTs are conserved in E. grandis as compared to annual plants and that they are expressed in all major plant tissues. Proteomic and transcriptomic analysis in roots suggest that the PRMTs of E. grandis control a number of regulatory proteins and genes related to signalling during cellular/root growth and morphogenesis. We demonstrate here, using chemical inhibition of methylation and transgenic approaches, that plant type I PRMTs are necessary for normal root growth and branching in E. grandis. We further show that EgPRMT1 has a key role in root hair initiation and elongation and is involved in the methylation of β-tubulin, a key protein in cytoskeleton formation.

CONCLUSIONS

Together, our data demonstrate that PRMTs encoded by E. grandis methylate a number of key proteins and alter the transcription of a variety of genes involved in developmental processes. Appropriate levels of expression of type I PRMTs are necessary for the proper growth and development of E. grandis roots.

Transcriptional Responses of Chilean Quinoa (Chenopodium quinoa Willd.) Under Water Deficit Conditions Uncovers ABA-Independent Expression Patterns

HIGHLIGHTS R49 genotype displayed best performance on selected physiological parameters and highest tolerance to drought.R49 drought over-represented transcripts has exhibited 19% of genes (306 contigs) that presented no homology to published databases.Expression pattern for canonical responses to drought such as ABA biosynthesis and other genes induced in response to drought were assessed by qPCR. Global freshwater shortage is one of the biggest challenges of our time, often associated to misuse, increased consumption demands and the effects of climate change, paralleled with the desertification of vast areas. Chenopodium quinoa (Willd.) represents a very promising species, due to both nutritional content and cultivation under water constraint. We characterized drought tolerance of three Chilean genotypes and selected Genotype R49 (Salares ecotype) based upon Relative Water Content (RWC), Electrolyte Leakage (EL) and maximum efficiency of photosystem II (Fv/Fm) after drought treatment, when compared to another two genotypes. Exploratory RNA-Seq of R49 was generated by Illumina paired-ends method comparing drought and control irrigation conditions. We obtained 104.8 million reads, with 54 million reads for control condition and 51 million reads for drought condition. Reads were assembled in 150,952 contigs, were 31,523 contigs have a reading frame of at least 300 nucleotides (100 aminoacids). BLAST2GO annotation showed a 15% of genes without homology to NCBI proteins, but increased to 19% (306 contigs) when focused into drought-induced genes. Expression pattern for canonical drought responses such as ABA biosynthesis and other genes induced were assessed by qPCR, suggesting novelty of R49 drought responses.

The Eucalyptus Tonoplast Intrinsic Protein (TIP) Gene Subfamily: Genomic Organization, Structural Features, and Expression Profiles

Plant aquaporins are water channels implicated in various physiological processes, including growth, development and adaptation to stress. In this study, the Tonoplast Intrinsic Protein (TIP) gene subfamily of Eucalyptus, an economically important woody species, was investigated and characterized. A genome-wide survey of the Eucalyptus grandis genome revealed the presence of eleven putative TIP genes (referred as EgTIP), which were individually assigned by phylogeny to each of the classical TIP1-5 groups. Homology modeling confirmed the presence of the two highly conserved NPA (Asn-Pro-Ala) motifs in the identified EgTIPs. Residue variations in the corresponding selectivity filters, that might reflect differences in EgTIP substrate specificity, were observed. All EgTIP genes, except EgTIP5.1, were transcribed and the majority of them showed organ/tissue-enriched expression. Inspection of the EgTIP promoters revealed the presence of common cis-regulatory elements implicated in abiotic stress and hormone responses pointing to an involvement of the identified genes in abiotic stress responses. In line with these observations, additional gene expression profiling demonstrated increased expression under polyethylene glycol-imposed osmotic stress. Overall, the results obtained suggest that these novel EgTIPs might be functionally implicated in eucalyptus adaptation to stress.

De novo transcriptome sequencing of black pepper (Piper nigrum L.) and an analysis of genes involved in phenylpropanoid metabolism in response to Phytophthora capsici

BACKGROUND

Piper nigrum L., or "black pepper", is an economically important spice crop in tropical regions. Black pepper production is markedly affected by foot rot disease caused by Phytophthora capsici, and genetic improvement of black pepper is essential for combating foot rot diseases. However, little is known about the mechanism of anti- P. capsici in black pepper. The molecular mechanisms underlying foot rot susceptibility were studied by comparing transcriptome analysis between resistant (Piper flaviflorum) and susceptible (Piper nigrum cv. Reyin-1) black pepper species.

RESULTS

116,432 unigenes were acquired from six libraries (three replicates of resistant and susceptible black pepper samples), which were integrated by applying BLAST similarity searches and noted by adopting Kyoto Encyclopaedia of Genes and Gene Ontology (GO) genome orthology identifiers. The reference transcriptome was mapped using two sets of digital gene expression data. Using GO enrichment analysis for the differentially expressed genes, the majority of the genes associated with the phenylpropanoid biosynthesis pathway were identified in P. flaviflorum. In addition, the expression of genes revealed that after susceptible and resistant species were inoculated with P. capsici, the majority of genes incorporated in the phenylpropanoid metabolism pathway were up-regulated in both species. Among various treatments and organs, all the genes were up-regulated to a relatively high degree in resistant species. Phenylalanine ammonia lyase and peroxidase enzyme activity increased in susceptible and resistant species after inoculation with P. capsici, and the resistant species increased faster. The resistant plants retain their vascular structure in lignin revealed by histochemical analysis.

CONCLUSIONS

Our data provide critical information regarding target genes and a technological basis for future studies of black pepper genetic improvements, including transgenic breeding.

Regulation, overexpression, and target gene identification of Potato Homeobox 15 (POTH15) - a class-I KNOX gene in potato

Potato Homeobox 15 (POTH15) is a KNOX-I (Knotted1-like homeobox) family gene in potato that is orthologous to Shoot Meristemless (STM) in Arabidopsis. Despite numerous reports on KNOX genes from different species, studies in potato are limited. Here, we describe photoperiodic regulation of POTH15, its overexpression phenotype, and identification of its potential targets in potato (Solanum tuberosum ssp. andigena). qRT-PCR analysis showed a higher abundance of POTH15 mRNA in shoot tips and stolons under tuber-inducing short-day conditions. POTH15 promoter activity was detected in apical and axillary meristems, stolon tips, tuber eyes, and meristems of tuber sprouts, indicating its role in meristem maintenance and leaf development. POTH15 overexpression altered multiple morphological traits including leaf and stem development, leaflet number, and number of nodes and branches. In particular, the rachis of the leaf was completely reduced and leaves appeared as a bouquet of leaflets. Comparative transcriptomic analysis of 35S::GUS and two POTH15 overexpression lines identified more than 6000 differentially expressed genes, including 2014 common genes between the two overexpression lines. Functional analysis of these genes revealed their involvement in responses to hormones, biotic/abiotic stresses, transcription regulation, and signal transduction. qRT-PCR of selected candidate target genes validated their differential expression in both overexpression lines. Out of 200 randomly chosen POTH15 targets, 173 were found to have at least one tandem TGAC core motif, characteristic of KNOX interaction, within 3.0kb in the upstream sequence of the transcription start site. Overall, this study provides insights to the role of POTH15 in controlling diverse developmental processes in potato.

Transcriptome analysis of immature xylem in the Chinese fir at different developmental phases

Background.Chinese fir [Cunninghamia lanceolata (Lamb.) Hook.] is one of the most important native tree species for timber production in southern China. An understanding of overall fast growing stage, stem growth stage and senescence stage cambium transcriptome variation is lacking. We used transcriptome sequencing to identify the repertoire of genes expressed during development of xylem tissue in Chinese fir, aiming to delineate the molecular mechanisms of wood formation. Results. We carried out transcriptome sequencing at three different cultivation ages (7Y, 15Y and 21Y) generating 68.71 million reads (13.88 Gbp). A total of 140,486 unigenes with a mean size of 568.64 base pairs (bp) were obtained via de novo assembly. Of these, 27,427 unigenes (19.52%) were further annotated by comparison to public protein databases. A total of 5,331 (3.79%) unigenes were mapped into 118 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Differentially expressed genes (DEG) analysis identified 3, 16 and 5,899 DEGs from the comparison of 7Y vs. 15Y, 7Y vs. 21Y and 15Y vs. 21Y, respectively, in the immature xylem tissues, including 2,638 significantly up-regulated and 3,280 significantly down-regulated genes. Besides, five NAC transcription factors, 190 MYB transcription factors, and 34 WRKY transcription factors were identified respectively from Chinese fir transcriptome. Conclusion. Our results revealed the active transcriptional pathways and identified the DEGs at different cultivation phases of Chinese fir wood formation. This transcriptome dataset will aid in understanding and carrying out future studies on the molecular basis of Chinese fir wood formation and contribute to future artificial production and applications.

Tricking Arthrinium malaysianum into Producing Industrially Important Enzymes Under 2-Deoxy D-Glucose Treatment

This study catalogs production of industrially important enzymes and changes in transcript expression caused by 2-deoxy D-glucose (2-DG) treatment in Arthrinium malaysianum cultures. Carbon Catabolite Repression (CCR) induced by 2-DG in this species is cAMP independent unlike many other organisms. Higher levels of secreted endoglucanase (EG), β-glucosidase (BGL), β-xylosidase (BXL), and filter paper activity assay (FPase) enzymes under 2-DG treatment can be exploited for commercial purposes. An integrated RNA sequencing and quantitative proteomic analysis was performed to investigate the cellular response to 2-DG in A. malaysianum. Analysis of RNASeq data under 2-DG treated and control condition reveals that 56% of the unigenes do not have any known similarity to proteins in non-redundant database. Gene Ontology IDs were assigned to 36% of the transcripts (13260) and about 5207 (14%) were mapped to Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). About 1711 genes encoding 2691 transcripts were differentially expressed in treated vs. control samples. Out of the 2691 differentially expressed transcripts, only 582 have any known function. The most up regulated genes belonged to Pentose Phosphate Pathways and carbohydrate degradation class as expected. In addition, genes involved in protein folding, binding, catalytic activity, DNA repair, and secondary metabolites were up-regulated under 2-DG treatment. Whereas genes encoding glycosylation pathways, growth, nutrient reservoir activity was repressed. Gene ontology analysis of the differentially expressed genes indicates metabolic process (35%) is the pre-dominant class followed by carbohydrate degradation (11%), protein folding, and trafficking (6.2%) and transport (5.3%) classes. Unlike other organisms, conventional unfolded protein response (UPR) was not activated in either control or treated conditions. Major enzymes secreted by A. malaysianum are those degrading plant polysaccharides, the most dominant ones being β-glucosidase, as demonstrated by the 2D gel analysis. A set of 7 differentially expressed mRNAs were validated by qPCR. Transmission electron microscopy analyses demonstrated that the 2-DG treated cell walls of hyphae showed significant differences in the cell-wall thickness. Overall 2-DG treatment in A. malaysianum induced secretion of large amount of commercially viable enzymes compared to other known species.

RNA-Seq Library Construction Methods for Transcriptome Analysis

Next-generation sequencing (NGS) technologies have revolutionized the study of genomics with an ever-expanding list of applications. RNA-Seq has emerged as a powerful method, applying transcriptome analysis to a wider range of organisms-most significantly, non-model organisms lacking prior genomic sequencing. Whereas an initial concern of NGS datasets was the potential limitation of short read lengths, short read sequences have been successfully employed in creation of de novo transcriptome assemblies that allow for subsequent mapping of reads for expression analysis. Prior genomic sequence knowledge is no longer a requirement for identification of functional transcriptional elements and for global gene expression characterization. Significant cost reductions in generating RNA-Seq data, and improvements in de novo assemblers, has allowed the analysis of transcriptomes in heretofore unsequenced plant species. These protocols describe standard methods for constructing RNA-Seq libraries to be sequenced on Illumina sequencing platforms for comprehensive transcriptome analysis. © 2016 by John Wiley & Sons, Inc.

Transcriptomic analysis of Arabidopsis overexpressing flowering locus T driven by a meristem-specific promoter that induces early flowering

Here we analyzed in leaves the effect of FT overexpression driven by meristem-specific KNAT1 gene homolog of Arabidopsis thaliana (Lincoln et al., 1994; Long et al., 1996) on the transcriptomic response during plant development. Our results demonstrated that meristematic FT overexpression generates a phenotype with an early flowering independent of photoperiod when compared with wild type (WT) plants. Arabidopsis FT-overexpressor lines (AtFTOE) did not show significant differences compared with WT lines neither in leaf number nor in rosette diameter up to day 21, when AtFTOE flowered. After this period AtFTOE plants started flower production and no new rosette leaves were produced. Additionally, WT plants continued on vegetative stage up to day 40, producing 12-14 rosette leaves before flowering. Transcriptomic analysis of rosette leaves studied by sequencing Illumina RNA-seq allowed us to determine the differential expression in mature leaf rosette of 3652 genes, being 626 of them up-regulated and 3026 down-regulated. Overexpressed genes related with flowering showed up-regulated transcription factors such as MADS-box that are known as flowering markers in meristem and which overexpression has been related with meristem identity preservation and the transition from vegetative to floral stage. Genes related with sugar transport have shown a higher demand of monosaccharides derived from the hydrolysis of sucrose to glucose and probably fructose, which can also be influenced by reproductive stage of AtFTOE plants.

Transcriptome sequencing reveals population differentiation in gene expression linked to functional traits and environmental gradients in the South African shrub Protea repens

Understanding the environmental and genetic mechanisms underlying locally adaptive trait variation across the ranges of species is a major focus of evolutionary biology. Combining transcriptome sequencing with common garden experiments on populations spanning geographical and environmental gradients holds promise for identifying such mechanisms. The South African shrub Protea repens displays diverse phenotypes in the wild along drought and temperature gradients. We grew plants from seeds collected at 19 populations spanning this species' range, and sequenced the transcriptomes of these plants to reveal gene pathways associated with adaptive trait variation. We related expression in co-expressed gene networks to trait phenotypes measured in the common garden and to source population climate. We found that expression in gene networks correlated with source-population environment and with plant traits. In particular, the activity of gene networks enriched for growth related pathways correlated strongly with source site minimum winter temperature and with leaf size, stem diameter and height in the garden. Other gene networks with enrichments for photosynthesis related genes showed associations with precipitation. Our results strongly suggest that this species displays population-level differences in gene expression that have been shaped by source population site climate, and that are reflected in trait variation along environmental gradients.

Differential Gene Expression between Leaf and Rhizome in Atractylodes lancea: A Comparative Transcriptome Analysis

The rhizome of Atractylodes lancea is extensively used in the practice of Traditional Chinese Medicine because of its broad pharmacological activities. This study was designed to characterize the transcriptome profiling of the rhizome and leaf of Atractylodes lancea in an attempt to uncover the molecular mechanisms regulating rhizome formation and growth. Over 270 million clean reads were assembled into 92,366 unigenes, 58% of which are homologous with sequences in public protein databases (NR, Swiss-Prot, GO, and KEGG). Analysis of expression levels showed that genes involved in photosynthesis, stress response, and translation were the most abundant transcripts in the leaf, while transcripts involved in stress response, transcription regulation, translation, and metabolism were dominant in the rhizome. Tissue-specific gene analysis identified distinct gene families active in the leaf and rhizome. Differential gene expression analysis revealed a clear difference in gene expression pattern, identifying 1518 up-regulated genes and 3464 down-regulated genes in the rhizome compared with the leaf, including a series of genes related to signal transduction, primary and secondary metabolism. Transcription factor (TF) analysis identified 42 TF families, with 67 and 60 TFs up-regulated in the rhizome and leaf, respectively. A total of 104 unigenes were identified as candidates for regulating rhizome formation and development. These data offer an overview of the gene expression pattern of the rhizome and leaf and provide essential information for future studies on the molecular mechanisms of controlling rhizome formation and growth. The extensive transcriptome data generated in this study will be a valuable resource for further functional genomics studies of A. lancea.

Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes

BACKGROUND

Gastrodia elata Blume (Orchidaceae) is an important Chinese medicine with several functional components. In the life cycle of G. elata, the orchid develops a symbiotic relationship with two compatible mycorrhizal fungi Mycena spp. and Armillaria mellea during seed germination to form vegetative propagation corm and vegetative growth to develop tubers, respectively. Gastrodin (p-hydroxymethylphenol-beta-D-glucoside) is the most important functional component in G. elata, and gastrodin significantly increases from vegetative propagation corms to tubers. To address the gene regulation mechanism in gastrodin biosynthesis in G. elata, a comparative analysis of de novo transcriptome sequencing among the vegetative propagation corms and tubers of G. elata and A. mellea was conducted using deep sequencing.

RESULTS

Transcriptome comparison between the vegetative propagation corms and juvenile tubers of G. elata revealed 703 differentially expressed unigenes, of which 298 and 405 unigenes were, respectively up-regulated (fold-change ≥ 2, q-value < 0.05, the trimmed mean of M-values (TMM)-normalized fragments per kilobase of transcript per Million mapped reads (FPKM) > 10) and down-regulated (fold-change ≤ 0.5, q-value <0.05, TMM-normalized FPKM > 10) in juvenile tubers. After Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 112 up-regulated unigenes with KEGG Ortholog identifiers (KOids) or enzyme commission (EC) numbers were assigned to 159 isogroups involved in seventy-eight different pathways, and 132 down-regulated unigenes with KOids or EC numbers were assigned to 168 isogroups, involved in eighty different pathways. The analysis of the isogroup genes from all pathways revealed that the two unigenes TRINITY_DN54282_c0_g1 (putative monooxygenases) and TRINITY_DN50323_c0_g1 (putative glycosyltransferases) might participate in hydroxylation and glucosylation in the gastrodin biosynthetic pathway.

CONCLUSIONS

The gene expression of the two unique unigenes encoding monooxygenase and glycosyltransferase significantly increases from vegetative propagation corms to tubers, and the molecular basis of gastrodin biosynthesis in the tubers of G. elata is proposed.

Transcriptome Analysis and Discovery of Genes Relevant to Development in Bradysia odoriphaga at Three Developmental Stages

Bradysia odoriphaga (Diptera: Sciaridae) is the most important pest of Chinese chive (Allium tuberosum) in Asia; however, the molecular genetics are poorly understood. To explore the molecular biological mechanism of development, Illumina sequencing and de novo assembly were performed in the third-instar, fourth-instar, and pupal B. odoriphaga. The study resulted in 16.2 Gb of clean data and 47,578 unigenes (≥125 bp) contained in 7,632,430 contigs, 46.21% of which were annotated from non-redundant protein (NR), Gene Ontology (GO), Clusters of Orthologous Groups (COG), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. It was found that 19.67% of unigenes matched the homologous species mainly, including Aedes aegypti, Culex quinquefasciatus, Ceratitis capitata, and Anopheles gambiae. According to differentially expressed gene (DEG) analysis, 143, 490, and 309 DEGs were annotated as involved in the developmental process in the GO database respectively, in the comparisons of third-instar and fourth-instar larvae, third-instar larvae and pupae, and fourth-instar larvae and pupae. Twenty-five genes were closely related to these processes, including developmental process, reproduction process, and reproductive organs development and programmed cell death (PCD). The information of unigenes assembled in B. odoriphaga through transcriptome and DEG analyses could provide a detailed genetic basis and regulated information for elaborating the developmental mechanism from the larval, pre-pupal to pupal stages of B. odoriphaga.

De Novo Assembled Wheat Transcriptomes Delineate Differentially Expressed Host Genes in Response to Leaf Rust Infection

Pathogens like Puccinia triticina, the causal organism for leaf rust, extensively damages wheat production. The interaction at molecular level between wheat and the pathogen is complex and less explored. The pathogen induced response was characterized using mock- or pathogen inoculated near-isogenic wheat lines (with or without seedling leaf rust resistance gene Lr28). Four Serial Analysis of Gene Expression libraries were prepared from mock- and pathogen inoculated plants and were subjected to Sequencing by Oligonucleotide Ligation and Detection, which generated a total of 165,767,777 reads, each 35 bases long. The reads were processed and multiple k-mers were attempted for de novo transcript assembly; 22 k-mers showed the best results. Altogether 21,345 contigs were generated and functionally characterized by gene ontology annotation, mining for transcription factors and resistance genes. Expression analysis among the four libraries showed extensive alterations in the transcriptome in response to pathogen infection, reflecting reorganizations in major biological processes and metabolic pathways. Role of auxin in determining pathogenesis in susceptible and resistant lines were imperative. The qPCR expression study of four LRR-RLK (Leucine-rich repeat receptor-like protein kinases) genes showed higher expression at 24 hrs after inoculation with pathogen. In summary, the conceptual model of induced resistance in wheat contributes insights on defense responses and imparts knowledge of Puccinia triticina-induced defense transcripts in wheat plants.

A De novo Transcriptomic Approach to Identify Flavonoids and Anthocyanins "Switch-Off" in Olive (Olea europaea L.) Drupes at Different Stages of Maturation

Highlights A de novo transcriptome reconstruction of olive drupes was performed in two genotypesGene expression was monitored during drupe development in two olive cultivarsTranscripts involved in flavonoid and anthocyanin pathways were analyzed in Cassanese and Leucocarpa cultivarsBoth cultivar and developmental stage impact gene expression in Olea europaea fruits. During ripening, the fruits of the olive tree (Olea europaea L.) undergo a progressive chromatic change characterized by the formation of a red-brown "spot" which gradually extends on the epidermis and in the innermost part of the mesocarp. This event finds an exception in the Leucocarpa cultivar, in which we observe a destabilized equilibrium between the metabolisms of chlorophyll and other pigments, particularly the anthocyanins whose switch-off during maturation promotes the white coloration of fruits. Despite its importance, genomic information on the olive tree is still lacking. Different RNA-seq libraries were generated from drupes of "Leucocarpa" and "Cassanese" olive genotypes, sampled at 100 and 130 days after flowering (DAF), and were used in order to identify transcripts involved in the main phenotypic changes of fruits during maturation and their corresponding expression patterns. A total of 103,359 transcripts were obtained and 3792 and 3064 were differentially expressed in "Leucocarpa" and "Cassanese" genotypes, respectively, during 100-130 DAF transition. Among them flavonoid and anthocyanin related transcripts such as phenylalanine ammonia lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol 3'-hydrogenase (F3'H), flavonol 3'5 '-hydrogenase (F3'5'H), flavonol synthase (FLS), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), UDP-glucose:anthocianidin: flavonoid glucosyltransferase (UFGT) were identified. These results contribute to reducing the current gap in information regarding metabolic processes, including those linked to fruit pigmentation in the olive.

The transcriptional landscape of insect galls: psyllid (Hemiptera) gall formation in Hawaiian Metrosideros polymorpha (Myrtaceae)

BACKGROUND

Recent studies show that galling Hymenoptera and Diptera are able to synthesize the plant hormone indole-3-acetic acid (auxin) from tryptophan and that plant response to insect-produced auxin is implicated in gall formation. We examined the leaf transcriptome of galled and ungalled leaves of individuals of the Hawaiian endemic plant Metrosideros polymorpha (Myrtaceae) subject to infestation by psyllid (Hemiptera) gall-makers in the genus Trioza (Triozidae).

RESULTS

Transcript libraries were sequenced using Illumina technology and the reads assembled de novo into contigs. Functional identification of contigs followed a two-step procedure, first identifying contigs by comparison to the completely sequenced genome of the related Eucalyptus, followed by identifying the equivalent Arabidopsis gene using a pre-computed mapping between Eucalyptus and Arabidopsis genes. This allowed us to use the rich functional annotation of the Arabidopsis genome to assess the transcriptional landscape of galling in Metrosideros. Comparing galled and ungalled leaves, we find a highly significant enrichment of expressed genes with a gene ontology (GO) annotation to auxin response in the former. One gene consistently expressed in all galled trees examined but not detected in any libraries from ungalled leaves was the Metrosideros version of SMALL AUXIN UPREGULATED (SAUR) 67 which appears to be a marker for leaf-galling in Metrosideros.

CONCLUSIONS

We conclude that an auxin response is involved in galling by Metrosideros psyllids. The possibility should therefore be considered that psyllids (like other insects examined) are able to synthesize auxin.

De novo transcriptome analysis of Medicago falcata reveals novel insights about the mechanisms underlying abiotic stress-responsive pathway

Background

The entire world is facing a deteriorating environment. Understanding the mechanisms underlying plant responses to external abiotic stresses is important for breeding stress-tolerant crops and herbages. Phytohormones play critical regulatory roles in plants in the response to external and internal cues to regulate growth and development. Medicago falcata is one of the stress-tolerant candidate leguminous species and is able to fix atmospheric nitrogen. This ability allows leguminous plants to grow in nitrogen deficient soils.

Methods

We performed Illumina sequencing of cDNA prepared from abiotic stress treated M. falcata. Sequencedreads were assembled to provide a transcriptome resource. Transcripts were annotated using BLASTsearches against the NCBI non-redundant database and gene ontology definitions were assigned. Acomparison among the three abiotic stress treated samples was carried out. The expression of transcriptswas confirmed with qRT-PCR.

Results

We present an abiotic stress-responsive M. falcata transcriptome using next-generation sequencing data from samples grown under standard, dehydration, high salinity, and cold conditions. We combined reads from all samples and de novo assembled 98,515 transcripts to build the M. falcata gene index. A comprehensive analysis of the transcriptome revealed abiotic stress-responsive mechanisms underlying the metabolism and core signalling components of major phytohormones. We identified nod factor signalling pathways during early symbiotic nodulation that are modified by abiotic stresses. Additionally, a global comparison of homology between the M. falcata and M. truncatula transcriptomes, along with five other leguminous species, revealed a high level of global sequence conservation within the family.

Conclusions

M. falcata is shown to be a model candidate for studying abiotic stress-responsive mechanisms in legumes. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to abiotic stresses. Our data provides many gene candidates that might be used for herbage and crop breeding. Additionally, FalcataBase (http://bioinformatics.cau.edu.cn/falcata/) was built for storing these data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2019-x) contains supplementary material, which is available to authorized users.

Large-scale transcriptional profiling of lignified tissues in Tectona grandis

BACKGROUND

Currently, Tectona grandis is one of the most valuable trees in the world and no transcript dataset related to secondary xylem is available. Considering how important the secondary xylem and sapwood transition from young to mature trees is, little is known about the expression differences between those successional processes and which transcription factors could regulate lignin biosynthesis in this tropical tree. Although MYB transcription factors are one of the largest superfamilies in plants related to secondary metabolism, it has not yet been characterized in teak. These results will open new perspectives for studies of diversity, ecology, breeding and genomic programs aiming to understand deeply the biology of this species.

RESULTS

We present a widely expressed gene catalog for T. grandis using Illumina technology and the de novo assembly. A total of 462,260 transcripts were obtained, with 1,502 and 931 genes differentially expressed for stem and branch secondary xylem, respectively, during age transition. Analysis of stem and branch secondary xylem indicates substantial similarity in gene ontologies including carbohydrate enzymes, response to stress, protein binding, and allowed us to find transcription factors and heat-shock proteins differentially expressed. TgMYB1 displays a MYB domain and a predicted coiled-coil (CC) domain, while TgMYB2, TgMYB3 and TgMYB4 showed R2R3-MYB domain and grouped with MYBs from several gymnosperms and flowering plants. TgMYB1, TgMYB4 and TgCES presented higher expression in mature secondary xylem, in contrast with TgMYB2, TgHsp1, TgHsp2, TgHsp3, and TgBi whose expression is higher in young lignified tissues. TgMYB3 is expressed at lower level in secondary xylem.

CONCLUSIONS

Expression patterns of MYB transcription factors and heat-shock proteins in lignified tissues are dissimilar when tree development was evaluated, obtaining more expression of TgMYB1 and TgMYB4 in lignified tissues of 60-year-old trees, and more expression in TgHsp1, TgHsp2, TgHsp3 and TgBi in stem secondary xylem of 12-year-old trees. We are opening a door for further functional characterization by reverse genetics and marker-assisted selection with those genes. Investigation of some of the key regulators of lignin biosynthesis in teak, however, could be a valuable step towards understanding how rigidity of teak wood and extractives content are different from most other woods. The obtained transcriptome data represents new sequences of T. grandis deposited in public databases, representing an unprecedented opportunity to discover several related-genes associated with secondary xylem such as transcription factors and stress-related genes in a tropical tree.

Transcriptomes That Confer to Plant Defense against Powdery Mildew Disease in Lagerstroemia indica

Transcriptome analysis was conducted in two popular Lagerstroemia cultivars: "Natchez" (NAT), a white flower and powdery mildew resistant interspecific hybrid and "Carolina Beauty" (CAB), a red flower and powdery mildew susceptible L. indica cultivar. RNA-seq reads were generated from Erysiphe australiana infected leaves and de novo assembled. A total of 37,035 unigenes from 224,443 assembled contigs in both genotypes were identified. Approximately 85% of these unigenes have known function. Of them, 475 KEGG genes were found significantly different between the two genotypes. Five of the top ten differentially expressed genes (DEGs) involved in the biosynthesis of secondary metabolites (plant defense) and four in flavonoid biosynthesis pathway (antioxidant activities or flower coloration). Furthermore, 5 of the 12 assembled unigenes in benzoxazinoid biosynthesis and 7 of 11 in flavonoid biosynthesis showed higher transcript abundance in NAT. The relative abundance of transcripts for 16 candidate DEGs (9 from CAB and 7 from NAT) detected by qRT-PCR showed general agreement with the abundances of the assembled transcripts in NAT. This study provided the first transcriptome analyses in L. indica. The differential transcript abundance between two genotypes indicates that it is possible to identify candidate genes that are associated with the plant defenses or flower coloration.

The Eucalyptus terpene synthase gene family

Background

Terpenoids are abundant in the foliage of Eucalyptus, providing the characteristic smell as well as being valuable economically and influencing ecological interactions. Quantitative and qualitative inter- and intra- specific variation of terpenes is common in eucalypts.

Results

The genome sequences of Eucalyptus grandis and E. globulus were mined for terpene synthase genes (TPS) and compared to other plant species. We investigated the relative expression of TPS in seven plant tissues and functionally characterized five TPS genes from E. grandis. Compared to other sequenced plant genomes, Eucalyptus grandis has the largest number of putative functional TPS genes of any sequenced plant. We discovered 113 and 106 putative functional TPS genes in E. grandis and E. globulus, respectively. All but one TPS from E. grandis were expressed in at least one of seven plant tissues examined. Genomic clusters of up to 20 genes were identified. Many TPS are expressed in tissues other than leaves which invites a re-evaluation of the function of terpenes in Eucalyptus.

Conclusions

Our data indicate that terpenes in Eucalyptus may play a wider role in biotic and abiotic interactions than previously thought. Tissue specific expression is common and the possibility of stress induction needs further investigation. Phylogenetic comparison of the two investigated Eucalyptus species gives insight about recent evolution of different clades within the TPS gene family. While the majority of TPS genes occur in orthologous pairs some clades show evidence of recent gene duplication, as well as loss of function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1598-x) contains supplementary material, which is available to authorized users.

The Eucalyptus terpene synthase gene family

BACKGROUND

Terpenoids are abundant in the foliage of Eucalyptus, providing the characteristic smell as well as being valuable economically and influencing ecological interactions. Quantitative and qualitative inter- and intra- specific variation of terpenes is common in eucalypts.

RESULTS

The genome sequences of Eucalyptus grandis and E. globulus were mined for terpene synthase genes (TPS) and compared to other plant species. We investigated the relative expression of TPS in seven plant tissues and functionally characterized five TPS genes from E. grandis. Compared to other sequenced plant genomes, Eucalyptus grandis has the largest number of putative functional TPS genes of any sequenced plant. We discovered 113 and 106 putative functional TPS genes in E. grandis and E. globulus, respectively. All but one TPS from E. grandis were expressed in at least one of seven plant tissues examined. Genomic clusters of up to 20 genes were identified. Many TPS are expressed in tissues other than leaves which invites a re-evaluation of the function of terpenes in Eucalyptus.

CONCLUSIONS

Our data indicate that terpenes in Eucalyptus may play a wider role in biotic and abiotic interactions than previously thought. Tissue specific expression is common and the possibility of stress induction needs further investigation. Phylogenetic comparison of the two investigated Eucalyptus species gives insight about recent evolution of different clades within the TPS gene family. While the majority of TPS genes occur in orthologous pairs some clades show evidence of recent gene duplication, as well as loss of function.

RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

BACKGROUND

The evolutionary basis of reproductive success in different environments is of major interest in the study of plant adaptation. Since the reproductive stage is particularly sensitive to drought, genes affecting reproductive success during this stage are key players in the evolution of adaptive mechanisms. We used an ecological genomics approach to investigate the reproductive response of drought-tolerant and sensitive wild barley accessions originating from different habitats in the Levant.

RESULTS

We sequenced mRNA extracted from spikelets at the flowering stage in drought-treated and control plants. The barley genome was used for a reference-guided assembly and differential expression analysis. Our approach enabled to detect biological processes affecting grain production under drought stress. We detected novel candidate genes and differentially expressed alleles associated with drought tolerance. Drought associated genes were shown to be more conserved than non-associated genes, and drought-tolerance genes were found to evolve more rapidly than other drought associated genes.

CONCLUSIONS

We show that reproductive success under drought stress is not a habitat-specific trait but a shared physiological adaptation that appeared to evolve recently in the evolutionary history of wild barley. Exploring the genomic basis of reproductive success under stress in crop wild progenitors is expected to have considerable ecological and economical applications.