Comparative analyses of two Geraniaceae transcriptomes using next-generation sequencing

Well-resolved phylogeny supports repeated evolution of keel flowers as a synergistic contributor to papilionoid legume diversification

The butterfly-shaped keel flower is a highly successful floral form in angiosperms. These flowers steer the mechanical interaction with bees and thus are hypothesized to accelerate pollinator-driven diversification. The exceptionally labile evolution of keel flowers in Papilionoideae (Fabaceae) provides a suitable system to test this hypothesis. Using 1456 low-copy nuclear loci, we confidently resolve the early divergence history of Papilionoideae. Constrained by this backbone phylogeny, we generated a time tree for 3326 Fabales to evaluate the tempo and mode of diversification within a state-dependent evolutionary framework. The first keel flowers emerged c. 59.0 million years ago in Papilionoideae, predating the earliest fossil by 3-4 million years. The Miocene diversification of Papilionoideae coincided with the rapid evolution of keel flowers. At least six independent origins and 32 losses of keel flowers were identified in Papilionoideae, Cercidoideae, and Polygalaceae. However, the state-dependent diversification model was not favored. Lack of radiation associated with keel flowers suggests that diversification within Papilionoideae was not solely driven by pollinator-mediated selection, but instead an outcome of the synergistic effects of multiple innovations, including nitrogen fixation and chemical defense, as well as dispersal into subtropical and temperate regions.

Plastid-nuclear coevolution of ribosomal protein genes in papilionoid legumes

In plants, cellular function is orchestrated by three distinct genomes located within the nucleus, mitochondrion, and plastid. These genomes are interdependent, requiring tightly coordinated maintenance and expression. Plastids host several multisubunit protein complexes encoded by both the plastid and nuclear genomes. To investigate plastid-nuclear coevolution, this study focused on plastid ribosomal protein genes that are encoded by both plastid and nuclear genomes from 50 taxa across 15 of the 22 early branching major clades of the legume subfamily Papilionoideae. Comparative analysis of substitution rates was conducted across five gene sets: nuclear-encoded plastid-targeted ribosomal protein genes (NuCpRP), nuclear-encoded cytosol-targeted ribosomal genes (NuCyRP), other nuclear-encoded plastid-targeted genes that are not involved in ribosomes (NuCpOT), plastid-encoded ribosomal protein genes (CpRP) and plastid-encoded photosynthesis genes (CpPS).1 Elevated nonsynonymous substitution rates (dN) and ratios of nonsynonymous to synonymous substitution rates (dN/dS; ω) were observed in both CpRP and NuCpRP compared to the other gene sets. Significant differences in dN for CpRP and NuCpRP were found between the papilionoid 50-kb inversion clade and other legumes. Using coevolution statistics and evolutionary rate covariation, strong signals of cytonuclear coevolution were identified, where nonsynonymous substitutions in CpRP and NuCpRP genes co-occur along the same branches of the Papilionoideae phylogeny. Increased ω in a few CpRP genes was due to intensified positive selection whereas most of the CpRP and NuCpRP increased ω was caused by relaxed purifying selection. This pattern not only underscores the role of cytonuclear incompatibility in driving speciation but also highlights its constraints on the genetic enhancement of papilionoid crop species.

Contrasting Metabolisms in Green and White Leaf Sectors of Variegated Pelargonium zonale—An Integrative Transcriptomic and Metabolomic Study

The photosynthetically active green leaf (GL) and non-active white leaf (WL) tissues of variegated Pelargonium zonale provide an excellent model system for studying processes associated with photosynthesis and sink-source interactions, enabling the same microenvironmental conditions. By combining differential transcriptomics and metabolomics, we identified the main differences between these two metabolically contrasting tissues. Genes related to photosynthesis and associated pigments, the Calvin–Benson cycle, fermentation, and glycolysis were strongly repressed in WL. On the other hand, genes related to nitrogen and protein metabolism, defence, cytoskeletal components (motor proteins), cell division, DNA replication, repair and recombination, chromatin remodelling, and histone modifications were upregulated in WL. A content of soluble sugars, TCA intermediates, ascorbate, and hydroxybenzoic acids was lower, while the concentration of free amino acids (AAs), hydroxycinnamic acids, and several quercetin and kaempferol glycosides was higher in WL than in GL. Therefore, WL presents a carbon sink and depends on photosynthetic and energy-generating processes in GL. Furthermore, the upregulated nitrogen metabolism in WL compensates for the insufficient energy from carbon metabolism by providing alternative respiratory substrates. At the same time, WL serves as nitrogen storage. Overall, our study provides a new genetic data resource for the use of this excellent model system and for ornamental pelargonium breeding and contributes to uncovering molecular mechanisms underlying variegation and its adaptive ecological value.

Rampant Nuclear Transfer and Substitutions of Plastid Genes in Passiflora

Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora.

Differentially Expressed Genes Shared by Two Distinct Cytoplasmic Male Sterility (CMS) Types of Silene vulgaris Suggest the Importance of Oxidative Stress in Pollen Abortion

Cytoplasmic male sterility (CMS), encoded by the interacting mitochondrial and nuclear genes, causes pollen abortion or non-viability. CMS is widely used in agriculture and extensively studied in crops. Much less is known about CMS in wild species. We performed a comparative transcriptomic analysis of male sterile and fertile individuals of Silene vulgaris, a model plant for the study of gynodioecy, to reveal the genes responsible for pollen abortion in this species. We used RNA-seq datasets previously employed for the analysis of mitochondrial and plastid transcriptomes of female and hermaphrodite flower buds, making it possible to compare the transcriptomes derived from three genomes in the same RNA specimen. We assembled de novo transcriptomes for two haplotypes of S. vulgaris and identified differentially expressed genes between the females and hermaphrodites, associated with stress response or pollen development. The gene for alternative oxidase was downregulated in females. The genetic pathways controlling CMS in S. vulgaris are similar to those in crops. The high number of the differentially expressed nuclear genes contrasts with the uniformity of organellar transcriptomes across genders, which suggests these pathways are evolutionarily conserved and that selective mechanisms may shield organellar transcription against changes in the cytoplasmic transcriptome.

Large-scale phylogenomics reveals ancient introgression in Asian Hepatica and new insights into the origin of the insular endemic Hepatica maxima

Hepatica maxima is native to Ulleungdo, which is one of the oceanic islands in Korea, and it likely originated via anagenetic speciation from the Korean mainland species H. asiatica. However, the relationships among the Asian lineages remain unresolved. Phylogenomics based on plant genomes can provide new insights into the evolutionary history of plants. We first generated plastid, mitochondrial and transcriptome sequences of the insular endemic species H. maxima. Using the genomic data for H. maxima, we obtained a phylogenomic dataset consisting of 76 plastid, 37 mitochondrial and 413 nuclear genes from Asian Hepatica and two outgroups. Coalescent- and concatenation-based methods revealed cytonuclear and organellar discordance in the lineage. The presence of gynodioecy with cytoplasmic male sterility in Asian Hepatica suggests that the discordance is correlated with potential disruption of linkage disequilibrium between the organellar genomes. Species network analyses revealed a deep history of hybridization and introgression in Asian Hepatica. We discovered that ancient and recent introgression events occurred throughout the evolutionary history of the insular endemic species H. maxima. The introgression may serve as an important source of genetic variation to facilitate adaptation to the Ulleungdo environment.

The Phytochemical Composition, Biological Effects and Biotechnological Approaches to the Production of High-Value Essential Oil from Geranium

Geraniaceae family plants are highly reputed aromatic and medicinal perennial branched herbs. The high economic value of these plants is due to their secondary metabolites, especially essential oil of foliage, which is a complex mixture of volatile phytochemicals, such as terpenes, esters, aldehydes, alcohols, ketones and phenols. The main phytoconstituents of the essential oil belong to the terpenoid group of metabolites, such as monoterpenes, sesquiterpenes, diterpenes and their esters. Of these, geraniol, linalool, citronellol and their esters (50–70%) generally constitute a major portion of essential oil, responsible for its fragrance. Essential oil is biosynthesized in specialized tissues known as glandular trichomes present in leaves, green branches and fresh flowers. Geraniaceae family plants have been highly useful in the perfumery, cosmetics, aromatherapy, pharmaceuticals and food industries. Several pharmacological properties such as antifungal, anti-inflammatory, anti-cancerous, anti-depressant, antibacterial, antioxidant, antiseptic, anti-dysentery, and antidiabetic properties are attributed to the presence of geranium oil. Further, it improves blood circulation, treats congestion, cleans the lymphatic system, strengthens the immune system, and is effective in combating nervousness, constipation, insomnia, anxiety and high blood pressure. The chapter discusses the phytochemical composition, pharmacological properties, genomics of essential oil biosynthetic pathway, enhancement of essential oil yield, and several biotechnological approaches to enhance the quantity as well as quality of essential oil in geranium.

RNA-seq highlights parallel and contrasting patterns in the evolution of the nuclear genome of fully mycoheterotrophic plants

BACKGROUND

While photosynthesis is the most notable trait of plants, several lineages of plants (so-called full heterotrophs) have adapted to obtain organic compounds from other sources. The switch to heterotrophy leads to profound changes at the morphological, physiological and genomic levels.

RESULTS

Here, we characterize the transcriptomes of three species representing two lineages of mycoheterotrophic plants: orchids (Epipogium aphyllum and Epipogium roseum) and Ericaceae (Hypopitys monotropa). Comparative analysis is used to highlight the parallelism between distantly related fully heterotrophic plants. In both lineages, we observed genome-wide elimination of nuclear genes that encode proteins related to photosynthesis, while systems associated with protein import to plastids as well as plastid transcription and translation remain active. Genes encoding components of plastid ribosomes that have been lost from the plastid genomes have not been transferred to the nuclear genomes; instead, some of the encoded proteins have been substituted by homologs. The nuclear genes of both Epipogium species accumulated nucleotide substitutions twice as rapidly as their photosynthetic relatives; in contrast, no increase in the substitution rate was observed in H. monotropa.

CONCLUSIONS

Full heterotrophy leads to profound changes in nuclear gene content. The observed increase in the rate of nucleotide substitutions is lineage specific, rather than a universal phenomenon among non-photosynthetic plants.

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.

Sequencing and De Novo Assembly of the Toxicodendron radicans (Poison Ivy) Transcriptome

Contact with poison ivy plants is widely dreaded because they produce a natural product called urushiol that is responsible for allergenic contact delayed-dermatitis symptoms lasting for weeks. For this reason, the catchphrase most associated with poison ivy is "leaves of three, let it be", which serves the purpose of both identification and an appeal for avoidance. Ironically, despite this notoriety, there is a dearth of specific knowledge about nearly all other aspects of poison ivy physiology and ecology. As a means of gaining a more molecular-oriented understanding of poison ivy physiology and ecology, Next Generation DNA sequencing technology was used to develop poison ivy root and leaf RNA-seq transcriptome resources. De novo assembled transcriptomes were analyzed to generate a core set of high quality expressed transcripts present in poison ivy tissue. The predicted protein sequences were evaluated for similarity to SwissProt homologs and InterProScan domains, as well as assigned both GO terms and KEGG annotations. Over 23,000 simple sequence repeats were identified in the transcriptome, and corresponding oligo nucleotide primer pairs were designed. A pan-transcriptome analysis of existing Anacardiaceae transcriptomes revealed conserved and unique transcripts among these species.

Transcriptomic analysis of salt stress responsive genes in Rhazya stricta

Rhazya stricta is an evergreen shrub that is widely distributed across Western and South Asia, and like many other members of the Apocynaceae produces monoterpene indole alkaloids that have anti-cancer properties. This species is adapted to very harsh desert conditions making it an excellent system for studying tolerance to high temperatures and salinity. RNA-Seq analysis was performed on R. stricta exposed to severe salt stress (500 mM NaCl) across four time intervals (0, 2, 12 and 24 h) to examine mechanisms of salt tolerance. A large number of transcripts including genes encoding tetrapyrroles and pentatricopeptide repeat (PPR) proteins were regulated only after 12 h of stress of seedlings grown in controlled greenhouse conditions. Mechanisms of salt tolerance in R. stricta may involve the upregulation of genes encoding chaperone protein Dnaj6, UDP-glucosyl transferase 85a2, protein transparent testa 12 and respiratory burst oxidase homolog protein b. Many of the highly-expressed genes act on protecting protein folding during salt stress and the production of flavonoids, key secondary metabolites in stress tolerance. Other regulated genes encode enzymes in the porphyrin and chlorophyll metabolic pathway with important roles during plant growth, photosynthesis, hormone signaling and abiotic responses. Heme biosynthesis in R. stricta leaves might add to the level of salt stress tolerance by maintaining appropriate levels of photosynthesis and normal plant growth as well as by the participation in reactive oxygen species (ROS) production under stress. We speculate that the high expression levels of PPR genes may be dependent on expression levels of their targeted editing genes. Although the results of PPR gene family indicated regulation of a large number of transcripts under salt stress, PPR actions were independent of the salt stress because their RNA editing patterns were unchanged.

De novo transcriptome analysis of rose-scented geranium provides insights into the metabolic specificity of terpene and tartaric acid biosynthesis

Background

Rose-scented geranium (Pelargonium sp.) is a perennial herb that produces a high value essential oil of fragrant significance due to the characteristic compositional blend of rose-oxide and acyclic monoterpenoids in foliage. Recently, the plant has also been shown to produce tartaric acid in leaf tissues. Rose-scented geranium represents top-tier cash crop in terms of economic returns and significance of the plant and plant products. However, there has hardly been any study on its metabolism and functional genomics, nor any genomic expression dataset resource is available in public domain. Therefore, to begin the gains in molecular understanding of specialized metabolic pathways of the plant, de novo sequencing of rose-scented geranium leaf transcriptome, transcript assembly, annotation, expression profiling as well as their validation were carried out.

Results

De novo transcriptome analysis resulted a total of 78,943 unique contigs (average length: 623 bp, and N50 length: 752 bp) from 15.44 million high quality raw reads. In silico functional annotation led to the identification of several putative genes representing terpene, ascorbic acid and tartaric acid biosynthetic pathways, hormone metabolism, and transcription factors. Additionally, a total of 6,040 simple sequence repeat (SSR) motifs were identified in 6.8% of the expressed transcripts. The highest frequency of SSR was of tri-nucleotides (50%). Further, transcriptome assembly was validated for randomly selected putative genes by standard PCR-based approach. In silico expression profile of assembled contigs were validated by real-time PCR analysis of selected transcripts.

Conclusion

Being the first report on transcriptome analysis of rose-scented geranium the data sets and the leads and directions reflected in this investigation will serve as a foundation for pursuing and understanding molecular aspects of its biology, and specialized metabolic pathways, metabolic engineering, genetic diversity as well as molecular breeding.

Electronic supplementary material

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

The Evolution of the FT/TFL1 Genes in Amaranthaceae and Their Expression Patterns in the Course of Vegetative Growth and Flowering in Chenopodium rubrum

The FT/TFL1 gene family controls important aspects of plant development: MFT-like genes affect germination, TFL1-like genes act as floral inhibitors, and FT-like genes are floral activators. Gene duplications produced paralogs with modified functions required by the specific lifestyles of various angiosperm species. We constructed the transcriptome of the weedy annual plant Chenopodium rubrum and used it for the comprehensive search for the FT/TFL1 genes. We analyzed their phylogenetic relationships across Amaranthaceae and all angiosperms. We discovered a very ancient phylogenetic clade of FT genes represented by the CrFTL3 gene of C. rubrum Another paralog CrFTL2 showed an unusual structural rearrangement which might have contributed to the functional shift. We examined the transcription patterns of the FT/TFL1 genes during the vegetative growth and floral transition in C. rubrum to get clues about their possible functions. All the genes except for the constitutively expressed CrFTL2 gene, and the CrFTL3 gene, which was transcribed only in seeds, exhibited organ-specific expression influenced by the specific light regime. The CrFTL1 gene was confirmed as a single floral activator from the FT/TFL1 family in C. rubrum Its floral promoting activity may be counteracted by CrTFL1 C. rubrum emerges as an easily manipulated model for the study of floral induction in weedy fast-cycling plants lacking a juvenile phase.

The nuclear genome of Rhazya stricta and the evolution of alkaloid diversity in a medically relevant clade of Apocynaceae

Alkaloid accumulation in plants is activated in response to stress, is limited in distribution and specific alkaloid repertoires are variable across taxa. Rauvolfioideae (Apocynaceae, Gentianales) represents a major center of structural expansion in the monoterpenoid indole alkaloids (MIAs) yielding thousands of unique molecules including highly valuable chemotherapeutics. The paucity of genome-level data for Apocynaceae precludes a deeper understanding of MIA pathway evolution hindering the elucidation of remaining pathway enzymes and the improvement of MIA availability in planta or in vitro. We sequenced the nuclear genome of Rhazya stricta (Apocynaceae, Rauvolfioideae) and present this high quality assembly in comparison with that of coffee (Rubiaceae, Coffea canephora, Gentianales) and others to investigate the evolution of genome-scale features. The annotated Rhazya genome was used to develop the community resource, RhaCyc, a metabolic pathway database. Gene family trees were constructed to identify homologs of MIA pathway genes and to examine their evolutionary history. We found that, unlike Coffea, the Rhazya lineage has experienced many structural rearrangements. Gene tree analyses suggest recent, lineage-specific expansion and diversification among homologs encoding MIA pathway genes in Gentianales and provide candidate sequences with the potential to close gaps in characterized pathways and support prospecting for new MIA production avenues.

A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives

We produced a unique large data set of reference transcriptomes to obtain new knowledge about the evolution of plant genomes and crop domestication. For this purpose, we validated a RNA-Seq data assembly protocol to perform comparative population genomics. For the validation, we assessed and compared the quality of de novo Illumina short-read assemblies using data from two crops for which an annotated reference genome was available, namely grapevine and sorghum. We used the same protocol for the release of 26 new transcriptomes of crop plants and wild relatives, including still understudied crops such as yam, pearl millet and fonio. The species list has a wide taxonomic representation with the inclusion of 15 monocots and 11 eudicots. All contigs were annotated using BLAST, prot4EST and Blast2GO. A strong originality of the data set is that each crop is associated with close relative species, which will permit whole-genome comparative evolutionary studies between crops and their wild-related species. This large resource will thus serve research communities working on both crops and model organisms. All the data are available at http://arcad-bioinformatics.southgreen.fr/.

Here We Are, But Where Do We Go? A Systematic Review of Crustacean Transcriptomic Studies from 2014-2015

Despite their economic, ecological, and experimental importance, genomic resources remain scarce for crustaceans. In lieu of genomes, many researchers have taken advantage of technological advancements to instead sequence and assemble crustacean transcriptomes de novo However, there is little consensus on what standard operating procedures are, or should be, for the field. Here, we systematically reviewed 53 studies published during 2014-2015 that utilized transcriptomic resources from this taxonomic group in an effort to identify commonalities as well as potential weaknesses that have applicability beyond just crustaceans. In general, these studies utilized RNA-Seq data, both novel and publicly available, to characterize transcriptomes and/or identify differentially expressed genes (DEGs) between treatments. Although the software suite Trinity was popular in assembly pipelines and other programs were also commonly employed, many studies failed to report crucial details regarding bioinformatic methodologies, including read mappers and the utilized parameters in identifying and characterizing DEGs. Annotation percentages for assembled transcriptomic contigs were low, averaging 32% overall. While other metrics, such as numbers of contigs and DEGs reported, correlated with the number of sequence reads utilized per sample, these did reach apparent saturation with increasing sequencing depth. Most disturbingly, a number of studies (55%) reported DEGs based on non-replicated experimental designs and single biological replicates for each treatment. Given this, we suggest future RNA-Seq experiments targeting transcriptome characterization conduct deeper (i.e., 50-100 M reads) sequencing while those examining differential expression instead focus more on increased biological replicates at shallower (i.e., ∼10-20 M reads/sample) sequencing depths. Moreover, the community must avoid submitting for review, or accepting for publication, non-replicated differential expression studies. Finally, mining the ever growing publicly available transcriptomic data from crustaceans will allow future studies to focus on hypothesis-driven research instead of continuing to simply characterize transcriptomes. As an example of this, we utilized neurotoxin sequences from the recently described remipede venom gland transcriptome in conjunction with publicly available crustacean transcriptomic data to derive preliminary results and hypotheses regarding the evolution of venom in crustaceans.

Plastid-Nuclear Interaction and Accelerated Coevolution in Plastid Ribosomal Genes in Geraniaceae

Plastids and mitochondria have many protein complexes that include subunits encoded by organelle and nuclear genomes. In animal cells, compensatory evolution between mitochondrial and nuclear-encoded subunits was identified and the high mitochondrial mutation rates were hypothesized to drive compensatory evolution in nuclear genomes. In plant cells, compensatory evolution between plastid and nucleus has rarely been investigated in a phylogenetic framework. To investigate plastid–nuclear coevolution, we focused on plastid ribosomal protein genes that are encoded by plastid and nuclear genomes from 27 Geraniales species. Substitution rates were compared for five sets of genes representing plastid- and nuclear-encoded ribosomal subunit proteins targeted to the cytosol or the plastid as well as nonribosomal protein controls. We found that nonsynonymous substitution rates (dN) and the ratios of nonsynonymous to synonymous substitution rates (ω) were accelerated in both plastid- (CpRP) and nuclear-encoded subunits (NuCpRP) of the plastid ribosome relative to control sequences. Our analyses revealed strong signals of cytonuclear coevolution between plastid- and nuclear-encoded subunits, in which nonsynonymous substitutions in CpRP and NuCpRP tend to occur along the same branches in the Geraniaceae phylogeny. This coevolution pattern cannot be explained by physical interaction between amino acid residues. The forces driving accelerated coevolution varied with cellular compartment of the sequence. Increased ω in CpRP was mainly due to intensified positive selection whereas increased ω in NuCpRP was caused by relaxed purifying selection. In addition, the many indels identified in plastid rRNA genes in Geraniaceae may have contributed to changes in plastid subunits.

Divergence of RNA polymerase α subunits in angiosperm plastid genomes is mediated by genomic rearrangement

Genes for the plastid-encoded RNA polymerase (PEP) persist in the plastid genomes of all photosynthetic angiosperms. However, three unrelated lineages (Annonaceae, Passifloraceae and Geraniaceae) have been identified with unusually divergent open reading frames (ORFs) in the conserved region of rpoA, the gene encoding the PEP α subunit. We used sequence-based approaches to evaluate whether these genes retain function. Both gene sequences and complete plastid genome sequences were assembled and analyzed from each of the three angiosperm families. Multiple lines of evidence indicated that the rpoA sequences are likely functional despite retaining as low as 30% nucleotide sequence identity with rpoA genes from outgroups in the same angiosperm order. The ratio of non-synonymous to synonymous substitutions indicated that these genes are under purifying selection, and bioinformatic prediction of conserved domains indicated that functional domains are preserved. One of the lineages (Pelargonium, Geraniaceae) contains species with multiple rpoA-like ORFs that show evidence of ongoing inter-paralog gene conversion. The plastid genomes containing these divergent rpoA genes have experienced extensive structural rearrangement, including large expansions of the inverted repeat. We propose that illegitimate recombination, not positive selection, has driven the divergence of rpoA.

Coevolution between Nuclear-Encoded DNA Replication, Recombination, and Repair Genes and Plastid Genome Complexity

Disruption of DNA replication, recombination, and repair (DNA-RRR) systems has been hypothesized to cause highly elevated nucleotide substitution rates and genome rearrangements in the plastids of angiosperms, but this theory remains untested. To investigate nuclear–plastid genome (plastome) coevolution in Geraniaceae, four different measures of plastome complexity (rearrangements, repeats, nucleotide insertions/deletions, and substitution rates) were evaluated along with substitution rates of 12 nuclear-encoded, plastid-targeted DNA-RRR genes from 27 Geraniales species. Significant correlations were detected for nonsynonymous (dN) but not synonymous (dS) substitution rates for three DNA-RRR genes (uvrB/C, why1, and gyrA) supporting a role for these genes in accelerated plastid genome evolution in Geraniaceae. Furthermore, correlation between dN of uvrB/C and plastome complexity suggests the presence of nucleotide excision repair system in plastids. Significant correlations were also detected between plastome complexity and 13 of the 90 nuclear-encoded organelle-targeted genes investigated. Comparisons revealed significant acceleration of dN in plastid-targeted genes of Geraniales relative to Brassicales suggesting this correlation may be an artifact of elevated rates in this gene set in Geraniaceae. Correlation between dN of plastid-targeted DNA-RRR genes and plastome complexity supports the hypothesis that the aberrant patterns in angiosperm plastome evolution could be caused by dysfunction in DNA-RRR systems.

De Novo Transcriptome Sequencing of Oryza officinalis Wall ex Watt to Identify Disease-Resistance Genes

Oryza officinalis Wall ex Watt is one of the most important wild relatives of cultivated rice and exhibits high resistance to many diseases. It has been used as a source of genes for introgression into cultivated rice. However, there are limited genomic resources and little genetic information publicly reported for this species. To better understand the pathways and factors involved in disease resistance and accelerating the process of rice breeding, we carried out a de novo transcriptome sequencing of O. officinalis. In this research, 137,229 contigs were obtained ranging from 200 to 19,214 bp with an N50 of 2331 bp through de novo assembly of leaves, stems and roots in O. officinalis using an Illumina HiSeq 2000 platform. Based on sequence similarity searches against a non-redundant protein database, a total of 88,249 contigs were annotated with gene descriptions and 75,589 transcripts were further assigned to GO terms. Candidate genes for plant–pathogen interaction and plant hormones regulation pathways involved in disease-resistance were identified. Further analyses of gene expression profiles showed that the majority of genes related to disease resistance were all expressed in the three tissues. In addition, there are two kinds of rice bacterial blight-resistant genes in O. officinalis, including two Xa1 genes and three Xa26 genes. All 2 Xa1 genes showed the highest expression level in stem, whereas one of Xa26 was expressed dominantly in leaf and other 2 Xa26 genes displayed low expression level in all three tissues. This transcriptomic database provides an opportunity for identifying the genes involved in disease-resistance and will provide a basis for studying functional genomics of O. officinalis and genetic improvement of cultivated rice in the future.

A transcriptome approach towards understanding the development of ripening capacity in 'Bartlett' pears (Pyrus communis L.)

BACKGROUND

The capacity of European pear fruit (Pyrus communis L.) to ripen after harvest develops during the final stages of growth on the tree. The objective of this study was to characterize changes in 'Bartlett' pear fruit physico-chemical properties and transcription profiles during fruit maturation leading to attainment of ripening capacity.

RESULTS

The softening response of pear fruit held for 14 days at 20 °C after harvest depended on their maturity. We identified four maturity stages: S1-failed to soften and S2- displayed partial softening (with or without ET-ethylene treatment); S3 - able to soften following ET; and S4 - able to soften without ET. Illumina sequencing and Trinity assembly generated 68,010 unigenes (mean length of 911 bp), of which 32.8 % were annotated to the RefSeq plant database. Higher numbers of differentially expressed transcripts were recorded in the S3-S4 and S1-S2 transitions (2805 and 2505 unigenes, respectively) than in the S2-S3 transition (2037 unigenes). High expression of genes putatively encoding pectin degradation enzymes in the S1-S2 transition suggests pectic oligomers may be involved as early signals triggering the transition to responsiveness to ethylene in pear fruit. Moreover, the co-expression of these genes with Exps (Expansins) suggests their collaboration in modifying cell wall polysaccharide networks that are required for fruit growth. K-means cluster analysis revealed that auxin signaling associated transcripts were enriched in cluster K6 that showed the highest gene expression at S3. AP2/EREBP (APETALA 2/ethylene response element binding protein) and bHLH (basic helix-loop-helix) transcripts were enriched in all three transition S1-S2, S2-S3, and S3-S4. Several members of Aux/IAA (Auxin/indole-3-acetic acid), ARF (Auxin response factors), and WRKY appeared to play an important role in orchestrating the S2-S3 transition.

CONCLUSIONS

We identified maturity stages associated with the development of ripening capacity in 'Bartlett' pear, and described the transcription profile of fruit at these stages. Our findings suggest that auxin is essential in regulating the transition of pear fruit from being ethylene-unresponsive (S2) to ethylene-responsive (S3), resulting in fruit softening. The transcriptome will be helpful for future studies about specific developmental pathways regulating the transition to ripening.

Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers

The exchange of genetic material between cellular organelles through intracellular gene transfer (IGT) or between species by horizontal gene transfer (HGT) has played an important role in plant mitochondrial genome evolution. The mitochondrial genomes of Geraniaceae display a number of unusual phenomena including highly accelerated rates of synonymous substitutions, extensive gene loss and reduction in RNA editing. Mitochondrial DNA sequences assembled for 17 species of Geranium revealed substantial reduction in gene and intron content relative to the ancestor of the Geranium lineage. Comparative analyses of nuclear transcriptome data suggest that a number of these sequences have been functionally relocated to the nucleus via IGT. Evidence for rampant HGT was detected in several Geranium species containing foreign organellar DNA from diverse eudicots, including many transfers from parasitic plants. One lineage has experienced multiple, independent HGT episodes, many of which occurred within the past 5.5 Myr. Both duplicative and recapture HGT were documented in Geranium lineages. The mitochondrial genome of Geranium brycei contains at least four independent HGT tracts that are absent in its nearest relative. Furthermore, G. brycei mitochondria carry two copies of the cox1 gene that differ in intron content, providing insight into contrasting hypotheses on cox1 intron evolution.

Exploring triacylglycerol biosynthetic pathway in developing seeds of Chia (Salvia hispanica L.): a transcriptomic approach

Chia (Salvia hispanica L.), a member of the mint family (Lamiaceae), is a rediscovered crop with great importance in health and nutrition and is also the highest known terrestrial plant source of heart-healthy omega-3 fatty acid, alpha linolenic acid (ALA). At present, there is no public genomic information or database available for this crop, hindering research on its genetic improvement through genomics-assisted breeding programs. The first comprehensive analysis of the global transcriptome profile of developing Salvia hispanica L. seeds, with special reference to lipid biosynthesis is presented in this study. RNA from five different stages of seed development was extracted and sequenced separately using the Illumina GAIIx platform. De novo assembly of processed reads in the pooled transcriptome using Trinity yielded 76,014 transcripts. The total transcript length was 66,944,462 bases (66.9 Mb), with an average length of approximately 880 bases. In the molecular functions category of Gene Ontology (GO) terms, ATP binding and nucleotide binding were found to be the most abundant and in the biological processes category, the metabolic process and the regulation of transcription-DNA-dependent and oxidation-reduction process were abundant. From the EuKaryotic Orthologous Groups of proteins (KOG) classification, the major category was “Metabolism” (31.97%), of which the most prominent class was ‘carbohydrate metabolism and transport’ (5.81% of total KOG classifications) followed by ‘secondary metabolite biosynthesis transport and catabolism’ (5.34%) and ‘lipid metabolism’ (4.57%). A majority of the candidate genes involved in lipid biosynthesis and oil accumulation were identified. Furthermore, 5596 simple sequence repeats (SSRs) were identified. The transcriptome data was further validated through confirmative PCR and qRT-PCR for select lipid genes. Our study provides insight into the complex transcriptome and will contribute to further genome-wide research and understanding of chia. The identified novel UniGenes will facilitate gene discovery and creation of genomic resource for this crop.

RNA-Seq analysis of the wild barley (H. spontaneum) leaf transcriptome under salt stress

Wild salt-tolerant barley (Hordeum spontaneum) is the ancestor of cultivated barley (Hordeum vulgare or H. vulgare). Although the cultivated barley genome is well studied, little is known about genome structure and function of its wild ancestor. In the present study, RNA-Seq analysis was performed on young leaves of wild barley treated with salt (500mM NaCl) at four different time intervals. Transcriptome sequencing yielded 103 to 115 million reads for all replicates of each treatment, corresponding to over 10 billion nucleotides per sample. Of the total reads, between 74.8 and 80.3% could be mapped and 77.4 to 81.7% of the transcripts were found in the H. vulgare unigene database (unigene-mapped). The unmapped wild barley reads for all treatments and replicates were assembled de novo and the resulting contigs were used as a new reference genome. This resulted in 94.3 to 95.3% of the unmapped reads mapping to the new reference. The number of differentially expressed transcripts was 9277, 3861 of which were unigene-mapped. The annotated unigene- and de novo-mapped transcripts (5100) were utilized to generate expression clusters across time of salt stress treatment. Two-dimensional hierarchical clustering classified differential expression profiles into nine expression clusters, four of which were selected for further analysis. Differentially expressed transcripts were assigned to the main functional categories. The most important groups were "response to external stimulus" and "electron-carrier activity". Highly expressed transcripts are involved in several biological processes, including electron transport and exchanger mechanisms, flavonoid biosynthesis, reactive oxygen species (ROS) scavenging, ethylene production, signaling network and protein refolding. The comparisons demonstrated that mRNA-Seq is an efficient method for the analysis of differentially expressed genes and biological processes under salt stress.

NDH expression marks major transitions in plant evolution and reveals coordinate intracellular gene loss

BACKGROUND

Key innovations have facilitated novel niche utilization, such as the movement of the algal predecessors of land plants into terrestrial habitats where drastic fluctuations in light intensity, ultraviolet radiation and water limitation required a number of adaptations. The NDH (NADH dehydrogenase-like) complex of Viridiplantae plastids participates in adapting the photosynthetic response to environmental stress, suggesting its involvement in the transition to terrestrial habitats. Although relatively rare, the loss or pseudogenization of plastid NDH genes is widely distributed across diverse lineages of photoautotrophic seed plants and mutants/transgenics lacking NDH function demonstrate little difference from wild type under non-stressed conditions. This study analyzes large transcriptomic and genomic datasets to evaluate the persistence and loss of NDH expression across plants.

RESULTS

Nuclear expression profiles showed accretion of the NDH gene complement at key transitions in land plant evolution, such as the transition to land and at the base of the angiosperm lineage. While detection of transcripts for a selection of non-NDH, photosynthesis related proteins was independent of the state of NDH, coordinate, lineage-specific loss of plastid NDH genes and expression of nuclear-encoded NDH subunits was documented in Pinaceae, gnetophytes, Orchidaceae and Geraniales confirming the independent and complete loss of NDH in these diverse seed plant taxa.

CONCLUSION

The broad phylogenetic distribution of NDH loss and the subtle phenotypes of mutants suggest that the NDH complex is of limited biological significance in contemporary plants. While NDH activity appears dispensable under favorable conditions, there were likely sufficiently frequent episodes of abiotic stress affecting terrestrial habitats to allow the retention of NDH activity. These findings reveal genetic factors influencing plant/environment interactions in a changing climate through 450 million years of land plant evolution.

Coordinated rates of evolution between interacting plastid and nuclear genes in Geraniaceae

Although gene coevolution has been widely observed within individuals and between different organisms, rarely has this phenomenon been investigated within a phylogenetic framework. The Geraniaceae is an attractive system in which to study plastid-nuclear genome coevolution due to the highly elevated evolutionary rates in plastid genomes. In plants, the plastid-encoded RNA polymerase (PEP) is a protein complex composed of subunits encoded by both plastid (rpoA, rpoB, rpoC1, and rpoC2) and nuclear genes (sig1-6). We used transcriptome and genomic data for 27 species of Geraniales in a systematic evaluation of coevolution between genes encoding subunits of the PEP holoenzyme. We detected strong correlations of dN (nonsynonymous substitutions) but not dS (synonymous substitutions) within rpoB/sig1 and rpoC2/sig2, but not for other plastid/nuclear gene pairs, and identified the correlation of dN/dS ratio between rpoB/C1/C2 and sig1/5/6, rpoC1/C2 and sig2, and rpoB/C2 and sig3 genes. Correlated rates between interacting plastid and nuclear sequences across the Geraniales could result from plastid-nuclear genome coevolution. Analyses of coevolved amino acid positions suggest that structurally mediated coevolution is not the major driver of plastid-nuclear coevolution. The detection of strong correlation of evolutionary rates between SIG and RNAP genes suggests a plausible explanation for plastome-genome incompatibility in Geraniaceae.

Optimal assembly strategies of transcriptome related to ploidies of eukaryotic organisms

BACKGROUND

Several de novo transcriptome assemblers have been developed recently to assemble the short reads generated from the next-generation sequencing platforms and different strategies were employed for assembling transcriptomes of various eukaryotes without genome sequences. Though there are some comparisons among these de novo assembly tools for assembling transcriptomes of different eukaryotic organisms, there is no report about the relationship between assembly strategies and ploidies of the organisms.

RESULTS

When we de novo assembled transcriptomes of sweet potato (hexaploid), Trametes gallica (a diploid fungus), Oryza meyeriana (a diploid wild rice), five assemblers, including Edena, Oases, Soaptrans, IDBA-tran and Trinity, were used in different strategies (Single-Assembler Single-Parameter, SASP; Single-Assembler Multiple-Parameters, SAMP; Combined De novo Transcriptome Assembly, CDTA, that is multiple assembler multiple parameter). It was found that CDTA strategy has the best performance compared with other two strategies for assembling transcriptome of the hexaploid sweet potato, whereas SAMP strategy with assembler Oases is better than other strategies for assembling transcriptomes of diploid fungus and the wild rice transcriptomes.

CONCLUSION

Based on the results from ours and others, it is suggested that CDTA strategy is better used for transcriptome assembly of polyploidy organisms and SAMP strategy of Oases is outperformed for those diploid organisms without genome sequences.

Evolutionary and biotechnology implications of plastid genome variation in the inverted-repeat-lacking clade of legumes

Land plant plastid genomes (plastomes) provide a tractable model for evolutionary study in that they are relatively compact and gene dense. Among the groups that display an appropriate level of variation for structural features, the inverted-repeat-lacking clade (IRLC) of papilionoid legumes presents the potential to advance general understanding of the mechanisms of genomic evolution. Here, are presented six complete plastome sequences from economically important species of the IRLC, a lineage previously represented by only five completed plastomes. A number of characters are compared across the IRLC including gene retention and divergence, synteny, repeat structure and functional gene transfer to the nucleus. The loss of clpP intron 2 was identified in one newly sequenced member of IRLC, Glycyrrhiza glabra. Using deeply sequenced nuclear transcriptomes from two species helped clarify the nature of the functional transfer of accD to the nucleus in Trifolium, which likely occurred in the lineage leading to subgenus Trifolium. Legumes are second only to cereal crops in agricultural importance based on area harvested and total production. Genetic improvement via plastid transformation of IRLC crop species is an appealing proposition. Comparative analyses of intergenic spacer regions emphasize the need for complete genome sequences for developing transformation vectors for plastid genetic engineering of legume crops.

Complete sequences of organelle genomes from the medicinal plant Rhazya stricta (Apocynaceae) and contrasting patterns of mitochondrial genome evolution across asterids

Background

Rhazya stricta is native to arid regions in South Asia and the Middle East and is used extensively in folk medicine to treat a wide range of diseases. In addition to generating genomic resources for this medicinally important plant, analyses of the complete plastid and mitochondrial genomes and a nuclear transcriptome from Rhazya provide insights into inter-compartmental transfers between genomes and the patterns of evolution among eight asterid mitochondrial genomes.

Results

The 154,841 bp plastid genome is highly conserved with gene content and order identical to the ancestral organization of angiosperms. The 548,608 bp mitochondrial genome exhibits a number of phenomena including the presence of recombinogenic repeats that generate a multipartite organization, transferred DNA from the plastid and nuclear genomes, and bidirectional DNA transfers between the mitochondrion and the nucleus. The mitochondrial genes sdh3 and rps14 have been transferred to the nucleus and have acquired targeting presequences. In the case of rps14, two copies are present in the nucleus; only one has a mitochondrial targeting presequence and may be functional. Phylogenetic analyses of both nuclear and mitochondrial copies of rps14 across angiosperms suggests Rhazya has experienced a single transfer of this gene to the nucleus, followed by a duplication event. Furthermore, the phylogenetic distribution of gene losses and the high level of sequence divergence in targeting presequences suggest multiple, independent transfers of both sdh3 and rps14 across asterids. Comparative analyses of mitochondrial genomes of eight sequenced asterids indicates a complicated evolutionary history in this large angiosperm clade with considerable diversity in genome organization and size, repeat, gene and intron content, and amount of foreign DNA from the plastid and nuclear genomes.

Conclusions

Organelle genomes of Rhazya stricta provide valuable information for improving the understanding of mitochondrial genome evolution among angiosperms. The genomic data have enabled a rigorous examination of the gene transfer events. Rhazya is unique among the eight sequenced asterids in the types of events that have shaped the evolution of its mitochondrial genome. Furthermore, the organelle genomes of R. stricta provide valuable genomic resources for utilizing this important medicinal plant in biotechnology applications.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-405) contains supplementary material, which is available to authorized users.