Applications of next-generation sequencing in plant biology

Morphometrics of Daucus (Apiaceae): a counterpart to a phylogenomic study

PREMISE OF STUDY

Molecular phylogenetics of genome-scale data sets (phylogenomics) often produces phylogenetic trees with unprecedented resolution. A companion phylogenomics analysis of Daucus using 94 conserved nuclear orthologs supported many of the traditional species but showed unexpected results that require morphological analyses to help interpret them in a practical taxonomic context.

METHODS

We evaluated character state distributions, stepwise discriminant analyses, canonical variate analyses, and hierarchical cluster analyses from 40 morphological characters from 81 accessions of 14 taxa of Daucus and eight species in related genera in an experimental plot.

KEY RESULTS

Most characters showed tremendous variation with character state overlap across many taxa. Multivariate analyses separated the outgroup taxa easily from the Daucus ingroup. Concordant with molecular analyses, most species form phenetic groups, except the same taxa that are problematical in the molecular results: (1) the subspecies of D. carota, (2) D. sahariensis and D. syrticus, and (3) D. broteri and D. guttatus.

CONCLUSIONS

Phenetic analyses, in combination with molecular data, support many Daucus species, but mostly by overlapping ranges of size and meristic variation. The subspecies of D. carota are poorly separated morphologically, are paraphyletic, and all could be recognized at the subspecies rank under D. carota. Daucus sahariensis and D. syrticus are so similar morphologically that they could be placed in synonymy. Combined molecular and morphological data support three species in accessions previously identified as D. broteri and D. guttatus. Molecular and morphological results support the new combination Daucus carota subsp. capillifolius.

A molecular identification protocol for roots of boreal forest tree species

PREMISE OF THE STUDY

Roots play a key role in many ecological processes, yet our ability to identify species from bulk root samples is limited. Molecular tools may be used to identify species from root samples, but they have not yet been developed for most systems. Here we present a PCR-based method previously used to identify roots of grassland species, modified for use in boreal forests. •

METHODS

We used repeatable interspecific size differences in fluorescent amplified fragment length polymorphisms of three noncoding chloroplast DNA regions to identify seven woody species common to boreal forests in Alberta, Canada. •

RESULTS

Abies balsamea, Alnus crispa, Betula papyrifera, Pinus contorta, and Populus tremuloides were identifiable to species, while Picea glauca and Picea mariana were identifiable to genus. In mixtures of known composition of foliar DNA, species were identified with 98% accuracy using one region. Mixed root samples of unknown composition were identified with 100% accuracy; four species were identified using one region, while three species were identified using two regions. •

DISCUSSION

This methodology is accurate, efficient, and inexpensive, and thus a valuable approach for ecological studies of roots. Furthermore, this method has now been validated for both grassland and boreal forest systems, and thus may also have applications in any plant community.

Phylogenetic properties of 50 nuclear loci in Medicago (Leguminosae) generated using multiplexed sequence capture and next-generation sequencing

Next-generation sequencing technology has increased the capacity to generate molecular data for plant biological research, including phylogenetics, and can potentially contribute to resolving complex phylogenetic problems. The evolutionary history of Medicago L. (Leguminosae: Trifoliae) remains unresolved due to incongruence between published phylogenies. Identification of the processes causing this genealogical incongruence is essential for the inference of a correct species phylogeny of the genus and requires that more molecular data, preferably from low-copy nuclear genes, are obtained across different species. Here we report the development of 50 novel LCN markers in Medicago and assess the phylogenetic properties of each marker. We used the genomic resources available for Medicago truncatula Gaertn., hybridisation-based gene enrichment (sequence capture) techniques and Next-Generation Sequencing to generate sequences. This alternative proves to be a cost-effective approach to amplicon sequencing in phylogenetic studies at the genus or tribe level and allows for an increase in number and size of targeted loci. Substitution rate estimates for each of the 50 loci are provided, and an overview of the variation in substitution rates among a large number of low-copy nuclear genes in plants is presented for the first time. Aligned sequences of major species lineages of Medicago and its sister genus are made available and can be used in further probe development for sequence-capture of the same markers.

Phylogenomics of the carrot genus (Daucus, Apiaceae)

•

PREMISE OF THE STUDY

We explored the utility of multiple nuclear orthologs for the taxonomic resolution of wild and cultivated carrot, Daucus species.•

METHODS

We studied the phylogeny of 92 accessions of 13 species and two subspecies of Daucus and 15 accessions of related genera (107 accessions total) with DNA sequences of 94 nuclear orthologs. Reiterative analyses examined data of both alleles using ambiguity codes or a single allele with the highest coverage, trimmed vs. untrimmed homopolymers; pure exonic vs. pure intronic data; the use of all 94 markers vs. a reduced subset of markers; and analysis of a concatenated data set vs. a coalescent (species tree) approach.•

KEY RESULTS

Our maximum parsimony and maximum likelihood trees were highly resolved, with 100% bootstrap support for most of the external and many of the internal clades. They resolved multiple accessions of many different species as monophyletic with strong support, but failed to support other species. The single allele analysis gave slightly better topological resolution; trimming homopolymers failed to increase taxonomic resolution; the exonic data had a smaller proportion of parsimony-informative characters. Similar results demonstrating the same dominant topology can be obtained with many fewer markers. A Bayesian concordance analysis provided an overall similar phylogeny, but the coalescent analysis provided drastic changes in topology to all the above.•

CONCLUSIONS

Our research highlights some difficult species groups in Daucus and misidentifications in germplasm collections. It highlights a useful subset of markers and approaches for future studies of dominant topologies in Daucus.

Genomics: a potential panacea for the perennial problem

Perennial crops represent important fresh and processed food sources worldwide, but advancements in breeding perennials are often impeded due to their very nature. The perennial crops we rely on most for food take several years to reach production maturity and require large spaces to grow, which make breeding new cultivars costly compared with most annual crops. Because breeding perennials is inefficient and expensive, they are often grown in monocultures consisting of small numbers of elite cultivars that are vegetatively propagated for decades or even centuries. This practice puts many perennial crops at risk for calamity since they remain stationary in the face of evolving pest and disease pressures. Although there is tremendous genetic diversity available to them, perennial crop breeders often struggle to generate commercially successful cultivars in a timely and cost-effective manner because of the high costs of breeding. Moreover, consumers often expect the same cultivars to be available indefinitely, and there is often little or no incentive for growers and retailers to take the risk of adopting new cultivars. While genomics studies linking DNA variants to commercially important traits have been performed in diverse perennial crops, the translation of these studies into accelerated breeding of improved cultivars has been limited. Here we explain the "perennial problem" in detail and demonstrate how modern genomics tools can significantly improve the cost effectiveness of breeding perennial crops and thereby prevent crucial food sources from succumbing to the perils of perpetual propagation.

Deep RNA-Seq to unlock the gene bank of floral development in Sinapis arvensis

Sinapis arvensis is a weed with strong biological activity. Despite being a problematic annual weed that contaminates agricultural crop yield, it is a valuable alien germplasm resource. It can be utilized for broadening the genetic background of Brassica crops with desirable agricultural traits like resistance to blackleg (Leptosphaeria maculans), stem rot (Sclerotinia sclerotium) and pod shatter (caused by FRUITFULL gene). However, few genetic studies of S. arvensis were reported because of the lack of genomic resources. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive dataset for S. arvensis for the first time. We used Illumina paired-end sequencing technology to sequence the S. arvensis flower transcriptome and generated 40,981,443 reads that were assembled into 131,278 transcripts. We de novo assembled 96,562 high quality unigenes with an average length of 832 bp. A total of 33,662 full-length ORF complete sequences were identified, and 41,415 unigenes were mapped onto 128 pathways using the KEGG Pathway database. The annotated unigenes were compared against Brassica rapa, B. oleracea, B. napus and Arabidopsis thaliana. Among these unigenes, 76,324 were identified as putative homologs of annotated sequences in the public protein databases, of which 1194 were associated with plant hormone signal transduction and 113 were related to gibberellin homeostasis/signaling. Unigenes that did not match any of those sequence datasets were considered to be unique to S. arvensis. Furthermore, 21,321 simple sequence repeats were found. Our study will enhance the currently available resources for Brassicaceae and will provide a platform for future genomic studies for genetic improvement of Brassica crops.

Future of portable devices for plant pathogen diagnosis

The demand for rapid and accurate diagnosis of plant diseases has risen in the last decade. On-site diagnosis of single or multiple pathogens using portable devices is the first step in this endeavour. Despite extensive attempts to develop portable devices for pathogen detection, current technologies are still restricted to detecting known pathogens with limited detection accuracy. Developing new detection techniques for rapid and accurate detection of multiple plant pathogens and their associated variants is essential. Recent single DNA sequencing technologies are a promising new avenue for developing future portable devices for plant pathogen detection. In this review, we detail the current progress in portable devices and technologies used for detecting plant pathogens, the current position of emerging sequencing technologies for analysis of plant genomics, and the future of portable devices for rapid pathogen diagnosis.

The role of auxin transporters in monocots development

Auxin is a key regulator of plant growth and development, orchestrating cell division, elongation and differentiation, embryonic development, root and stem tropisms, apical dominance, and transition to flowering. Auxin levels are higher in undifferentiated cell populations and decrease following organ initiation and tissue differentiation. This differential auxin distribution is achieved by polar auxin transport (PAT) mediated by auxin transport proteins. There are four major families of auxin transporters in plants: PIN-FORMED (PIN), ATP-binding cassette family B (ABCB), AUXIN1/LIKE-AUX1s, and PIN-LIKES. These families include proteins located at the plasma membrane or at the endoplasmic reticulum (ER), which participate in auxin influx, efflux or both, from the apoplast into the cell or from the cytosol into the ER compartment. Auxin transporters have been largely studied in the dicotyledon model species Arabidopsis, but there is increasing evidence of their role in auxin regulated development in monocotyledon species. In monocots, families of auxin transporters are enlarged and often include duplicated genes and proteins with high sequence similarity. Some of these proteins underwent sub- and neo-functionalization with substantial modification to their structure and expression in organs such as adventitious roots, panicles, tassels, and ears. Most of the present information on monocot auxin transporters function derives from studies conducted in rice, maize, sorghum, and Brachypodium, using pharmacological applications (PAT inhibitors) or down-/up-regulation (over-expression and RNA interference) of candidate genes. Gene expression studies and comparison of predicted protein structures have also increased our knowledge of the role of PAT in monocots. However, knockout mutants and functional characterization of single genes are still scarce and the future availability of such resources will prove crucial to elucidate the role of auxin transporters in monocots development.

Next generation restoration genetics: applications and opportunities

Restoration ecology is a young scientific discipline underpinning improvements in the rapid global expansion of ecological restoration. The application of molecular tools over the past 20 years has made an important contribution to understanding genetic factors influencing ecological restoration success. Here we illustrate how recent advances in next generation sequencing (NGS) methods are revolutionising the practical contribution of genetics to restoration. Novel applications include a dramatically enhanced capacity to measure adaptive variation for optimal seed sourcing, high-throughput assessment and monitoring of natural and restored biological communities aboveground and belowground, and gene expression analysis as a measure of genetic resilience of restored populations. Challenges remain in data generation, handling and analysis, and how best to apply NGS for practical outcomes in restoration.

De novo assembly and transcriptome analysis of the rubber tree (Hevea brasiliensis) and SNP markers development for rubber biosynthesis pathways

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next-generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra-individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.

Whole genome shotgun sequences for microsatellite discovery and application in cultivated and wild Macadamia (Proteaceae)

PREMISE OF THE STUDY

Next-generation sequencing (NGS) data are widely used for single-nucleotide polymorphism discovery and genetic marker development in species with limited available genome information. We developed microsatellite primers for the Proteaceae nut crop species Macadamia integrifolia and assessed cross-species transferability in all congeners to investigate genetic identification of cultivars and gene flow. •

METHODS AND RESULTS

Primers were designed from both raw and assembled Illumina NGS paired-end reads. The final 12 microsatellite markers selected were polymorphic among wild individuals of all four Macadamia species-M. integrifolia, M. tetraphylla, M. ternifolia, and M. jansenii-and in commercial macadamia cultivars including hybrids. •

CONCLUSIONS

We demonstrate the utility of raw and assembled Illumina NGS reads from total genomic DNA for the rapid development of microsatellites in Macadamia. These primers will facilitate future studies of population structure, hybridization, parentage, and cultivar identification in cultivated and wild Macadamia populations.

High-throughput transcriptome sequencing and preliminary functional analysis in four Neotropical tree species

Background

The Amazonian rainforest is predicted to suffer from ongoing environmental changes. Despite the need to evaluate the impact of such changes on tree genetic diversity, we almost entirely lack genomic resources.

Results

In this study, we analysed the transcriptome of four tropical tree species (Carapa guianensis, Eperua falcata, Symphonia globulifera and Virola michelii) with contrasting ecological features, belonging to four widespread botanical families (respectively Meliaceae, Fabaceae, Clusiaceae and Myristicaceae). We sequenced cDNA libraries from three organs (leaves, stems, and roots) using 454 pyrosequencing. We have developed an R and bioperl-based bioinformatic procedure for de novo assembly, gene functional annotation and marker discovery. Mismatch identification takes into account single-base quality values as well as the likelihood of false variants as a function of contig depth and number of sequenced chromosomes. Between 17103 (for Symphonia globulifera) and 23390 (for Eperua falcata) contigs were assembled. Organs varied in the numbers of unigenes they apparently express, with higher number in roots. Patterns of gene expression were similar across species, with metabolism of aromatic compounds standing out as an overrepresented gene function. Transcripts corresponding to several gene functions were found to be over- or underrepresented in each organ. We identified between 4434 (for Symphonia globulifera) and 9076 (for Virola surinamensis) well-supported mismatches. The resulting overall mismatch density was comprised between 0.89 (S. globulifera) and 1.05 (V. surinamensis) mismatches/100 bp in variation-containing contigs.

Conclusion

The relative representation of gene functions in the four transcriptomes suggests that secondary metabolism may be particularly important in tropical trees. The differential representation of transcripts among tissues suggests differential gene expression, which opens the way to functional studies in these non-model, ecologically important species. We found substantial amounts of mismatches in the four species. These newly identified putative variants are a first step towards acquiring much needed genomic resources for tropical tree species.

Electronic supplementary material

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

Peroxidase gene discovery from the horseradish transcriptome

BACKGROUND

Horseradish peroxidases (HRPs) from Armoracia rusticana have long been utilized as reporters in various diagnostic assays and histochemical stainings. Regardless of their increasing importance in the field of life sciences and suggested uses in medical applications, chemical synthesis and other industrial applications, the HRP isoenzymes, their substrate specificities and enzymatic properties are poorly characterized. Due to lacking sequence information of natural isoenzymes and the low levels of HRP expression in heterologous hosts, commercially available HRP is still extracted as a mixture of isoenzymes from the roots of A. rusticana.

RESULTS

In this study, a normalized, size-selected A. rusticana transcriptome library was sequenced using 454 Titanium technology. The resulting reads were assembled into 14871 isotigs with an average length of 1133 bp. Sequence databases, ORF finding and ORF characterization were utilized to identify peroxidase genes from the 14871 isotigs generated by de novo assembly. The sequences were manually reviewed and verified with Sanger sequencing of PCR amplified genomic fragments, resulting in the discovery of 28 secretory peroxidases, 23 of them previously unknown. A total of 22 isoenzymes including allelic variants were successfully expressed in Pichia pastoris and showed peroxidase activity with at least one of the substrates tested, thus enabling their development into commercial pure isoenzymes.

CONCLUSIONS

This study demonstrates that transcriptome sequencing combined with sequence motif search is a powerful concept for the discovery and quick supply of new enzymes and isoenzymes from any plant or other eukaryotic organisms. Identification and manual verification of the sequences of 28 HRP isoenzymes do not only contribute a set of peroxidases for industrial, biological and biomedical applications, but also provide valuable information on the reliability of the approach in identifying and characterizing a large group of isoenzymes.

A target enrichment method for gathering phylogenetic information from hundreds of loci: An example from the Compositae

PREMISE OF THE STUDY

The Compositae (Asteraceae) are a large and diverse family of plants, and the most comprehensive phylogeny to date is a meta-tree based on 10 chloroplast loci that has several major unresolved nodes. We describe the development of an approach that enables the rapid sequencing of large numbers of orthologous nuclear loci to facilitate efficient phylogenomic analyses. •

METHODS AND RESULTS

We designed a set of sequence capture probes that target conserved orthologous sequences in the Compositae. We also developed a bioinformatic and phylogenetic workflow for processing and analyzing the resulting data. Application of our approach to 15 species from across the Compositae resulted in the production of phylogenetically informative sequence data from 763 loci and the successful reconstruction of known phylogenetic relationships across the family. •

CONCLUSIONS

These methods should be of great use to members of the broader Compositae community, and the general approach should also be of use to researchers studying other families.

TR-DB: an open-access database of compounds affecting the ethylene-induced triple response in Arabidopsis

Small molecules which act as hormone agonists or antagonists represent useful tools in fundamental research and are widely applied in agriculture to control hormone effects. High-throughput screening of large chemical compound libraries has yielded new findings in plant biology, with possible future applications in agriculture and horticulture. To further understand ethylene biosynthesis/signaling and its crosstalk with other hormones, we screened a 12,000 compound chemical library based on an ethylene-related bioassay of dark-grown Arabidopsis thaliana (L.) Heynh. seedlings. From the initial screening, 1313 (∼11%) biologically active small molecules altering the phenotype triggered by the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC), were identified. Selection and sorting in classes were based on the angle of curvature of the apical hook, the length and width of the hypocotyl and the root. A MySQL-database was constructed (https://chaos.ugent.be/WE15/) including basic chemical information on the compounds, images illustrating the phenotypes, phenotype descriptions and classification. The research perspectives for different classes of hit compounds will be evaluated, and some general screening tips for customized high-throughput screening and pitfalls will be discussed.

Plant taxonomy: a historical perspective, current challenges, and perspectives

Taxonomy is the science that explores, describes, names, and classifies all organisms. In this introductory chapter, we highlight the major historical steps in the elaboration of this science that provides baseline data for all fields of biology and plays a vital role for society but is also an independent, complex, and sound hypothesis-driven scientific discipline.In a first part, we underline that plant taxonomy is one of the earliest scientific disciplines that emerged thousands of years ago, even before the important contributions of Greeks and Romans (e.g., Theophrastus, Pliny the Elder, and Dioscorides). In the fifteenth to sixteenth centuries, plant taxonomy benefited from the Great Navigations, the invention of the printing press, the creation of botanic gardens, and the use of the drying technique to preserve plant specimens. In parallel with the growing body of morpho-anatomical data, subsequent major steps in the history of plant taxonomy include the emergence of the concept of natural classification, the adoption of the binomial naming system (with the major role of Linnaeus) and other universal rules for the naming of plants, the formulation of the principle of subordination of characters, and the advent of the evolutionary thought. More recently, the cladistic theory (initiated by Hennig) and the rapid advances in DNA technologies allowed to infer phylogenies and to propose true natural, genealogy-based classifications.In a second part, we put the emphasis on the challenges that plant taxonomy faces nowadays. The still very incomplete taxonomic knowledge of the worldwide flora (the so-called taxonomic impediment) is seriously hampering conservation efforts that are especially crucial as biodiversity enters its sixth extinction crisis. It appears mainly due to insufficient funding, lack of taxonomic expertise, and lack of communication and coordination. We then review recent initiatives to overcome these limitations and to anticipate how taxonomy should and could evolve. In particular, the use of molecular data has been era-splitting for taxonomy and may allow an accelerated pace of species discovery. We examine both strengths and limitations of such techniques in comparison to morphology-based investigations, we give broad recommendations on the use of molecular tools for plant taxonomy, and we highlight the need for an integrative taxonomy based on evidence from multiple sources.

Microsatellite markers from the Ion Torrent: a multi-species contrast to 454 shotgun sequencing

The development and screening of microsatellite markers have been accelerated by next-generation sequencing (NGS) technology and in particular GS-FLX pyro-sequencing (454). More recent platforms such as the PGM semiconductor sequencer (Ion Torrent) offer potential benefits such as dramatic reductions in cost, but to date have not been well utilized. Here, we critically compare the advantages and disadvantages of microsatellite development using PGM semiconductor sequencing and GS-FLX pyro-sequencing for two gymnosperm (a conifer and a cycad) and one angiosperm species. We show that these NGS platforms differ in the quantity of returned sequence data, unique microsatellite data and primer design opportunities, mostly consistent with the differences in read length. The strength of the PGM lies in the large amount of data generated at a comparatively lower cost and time. The strength of GS-FLX lies in the return of longer average length sequences and therefore greater flexibility in producing markers with variable product length, due to longer flanking regions, which is ideal for capillary multiplexing. These differences need to be considered when choosing a NGS method for microsatellite discovery. However, the ongoing improvement in read lengths of the NGS platforms will reduce the disadvantage of the current short read lengths, particularly for the PGM platform, allowing greater flexibility in primer design coupled with the power of a larger number of sequences.

Students' perspective on genomics: from sample to sequence using the case study of blueberry

Advances in genomic sequencing technologies in the past decade have revolutionized the field of genomics, resulting in faster and less expensive sequencing. Holding back the potential for innovation, however, is a widespread lack of understanding of genomics and sequencing by the general public. In an attempt to remedy this problem, this paper presents an introduction to the fields of genomics, bioinformatics, and proteomics using the blueberry genome as a model case study of the plant genomics field. The blueberry (Vaccinium sect. Cyanococcus) is often cited as a “super food” in the media due to its nutritional benefits and global economic importance. There have been a number of related genomic publications in the past 20 years; however, a completed genome and a full analysis into the health-related pathways are still needed. As exemplified by this blueberry case study, there are opportunities for future genomic research into numerous beneficial plant species. The solid background presented in this paper provides future researchers the foundation to explore these uncharted areas.

Peculiar evolutionary history of miR390-guided TAS3-like genes in land plants

PCR-based approach was used as a phylogenetic profiling tool to probe genomic DNA samples from representatives of evolutionary distant moss taxa, namely, classes Bryopsida, Tetraphidopsida, Polytrichopsida, Andreaeopsida, and Sphagnopsida. We found relatives of all Physcomitrella patens miR390 and TAS3-like loci in these plant taxa excluding Sphagnopsida. Importantly, cloning and sequencing of Marchantia polymorpha genomic DNA showed miR390 and TAS3-like sequences which were also found among genomic reads of M. polymorpha at NCBI database. Our data suggest that the ancient plant miR390-dependent TAS molecular machinery firstly evolved to target AP2-like mRNAs in Marchantiophyta and only then both ARF- and AP2-specific mRNAs in mosses. The presented analysis shows that moss TAS3 families may undergone losses of tasiAP2 sites during evolution toward ferns and seed plants. These data confirm that miR390-guided genes coding for ARF- and AP2-specific ta-siRNAs have been gradually changed during land plant evolution.

Construction of a reference genetic linkage map for carnation (Dianthus caryophyllus L.)

Background

Genetic linkage maps are important tools for many genetic applications including mapping of quantitative trait loci (QTLs), identifying DNA markers for fingerprinting, and map-based gene cloning. Carnation (Dianthus caryophyllus L.) is an important ornamental flower worldwide. We previously reported a random amplified polymorphic DNA (RAPD)-based genetic linkage map derived from Dianthus capitatus ssp. andrezejowskianus and a simple sequence repeat (SSR)-based genetic linkage map constructed using data from intraspecific F₂ populations; however, the number of markers was insufficient, and so the number of linkage groups (LGs) did not coincide with the number of chromosomes (x = 15). Therefore, we aimed to produce a high-density genetic map to improve its usefulness for breeding purposes and genetic research.

Results

We improved the SSR-based genetic linkage map using SSR markers derived from a genomic library, expression sequence tags, and RNA-seq data. Linkage analysis revealed that 412 SSR loci (including 234 newly developed SSR loci) could be mapped to 17 linkage groups (LGs) covering 969.6 cM. Comparison of five minor LGs covering less than 50 cM with LGs in our previous RAPD-based genetic map suggested that four LGs could be integrated into two LGs by anchoring common SSR loci. Consequently, the number of LGs corresponded to the number of chromosomes (x = 15). We added 192 new SSRs, eight RAPD, and two sequence-tagged site loci to refine the RAPD-based genetic linkage map, which comprised 15 LGs consisting of 348 loci covering 978.3 cM. The two maps had 125 SSR loci in common, and most of the positions of markers were conserved between them. We identified 635 loci in carnation using the two linkage maps. We also mapped QTLs for two traits (bacterial wilt resistance and anthocyanin pigmentation in the flower) and a phenotypic locus for flower-type by analyzing previously reported genotype and phenotype data.

Conclusions

The improved genetic linkage maps and SSR markers developed in this study will serve as reference genetic linkage maps for members of the genus Dianthus, including carnation, and will be useful for mapping QTLs associated with various traits, and for improving carnation breeding programs.

A more complete picture of metal hyperaccumulation through next-generation sequencing technologies

The mechanistic understanding of metal hyperaccumulation has benefitted immensely from the use of molecular genetics tools developed for Arabidopsis thaliana. The revolution in DNA sequencing will enable even greater strides in the near future, this time not restricted to the family Brassicaceae. Reference genomes are within reach for many ecologically interesting species including heterozygous outbreeders. They will allow deep RNA-seq transcriptome studies and the re-sequencing of contrasting individuals to unravel the genetic basis of phenotypic variation. Cell-type specific transcriptome analyses, which will be essential for the dissection of metal translocation pathways in hyperaccumulators, can be achieved through the combination of RNA-seq and translatome approaches. Affordable high-resolution genotyping of many individuals enables the elucidation of quantitative trait loci in intra- and interspecific crosses as well as through genome-wide association mapping across large panels of accessions. Furthermore, genome-wide scans have the power to detect loci under recent selection. Together these approaches will lead to a detailed understanding of the evolutionary path towards the emergence of hyperaccumulation traits.

SNP discovery by illumina-based transcriptome sequencing of the olive and the genetic characterization of Turkish olive genotypes revealed by AFLP, SSR and SNP markers

Background

The olive tree (Olea europaea L.) is a diploid (2n = 2x = 46) outcrossing species mainly grown in the Mediterranean area, where it is the most important oil-producing crop. Because of its economic, cultural and ecological importance, various DNA markers have been used in the olive to characterize and elucidate homonyms, synonyms and unknown accessions. However, a comprehensive characterization and a full sequence of its transcriptome are unavailable, leading to the importance of an efficient large-scale single nucleotide polymorphism (SNP) discovery in olive. The objectives of this study were (1) to discover olive SNPs using next-generation sequencing and to identify SNP primers for cultivar identification and (2) to characterize 96 olive genotypes originating from different regions of Turkey.

Methodology/Principal Findings

Next-generation sequencing technology was used with five distinct olive genotypes and generated cDNA, producing 126,542,413 reads using an Illumina Genome Analyzer IIx. Following quality and size trimming, the high-quality reads were assembled into 22,052 contigs with an average length of 1,321 bases and 45 singletons. The SNPs were filtered and 2,987 high-quality putative SNP primers were identified. The assembled sequences and singletons were subjected to BLAST similarity searches and annotated with a Gene Ontology identifier. To identify the 96 olive genotypes, these SNP primers were applied to the genotypes in combination with amplified fragment length polymorphism (AFLP) and simple sequence repeats (SSR) markers.

Conclusions/Significance

This study marks the highest number of SNP markers discovered to date from olive genotypes using transcriptome sequencing. The developed SNP markers will provide a useful source for molecular genetic studies, such as genetic diversity and characterization, high density quantitative trait locus (QTL) analysis, association mapping and map-based gene cloning in the olive. High levels of genetic variation among Turkish olive genotypes revealed by SNPs, AFLPs and SSRs allowed us to characterize the Turkish olive genotype.

Augmenting transcriptome assembly by combining de novo and genome-guided tools

Researchers interested in studying and constructing transcriptomes, especially for non-model species, face the conundrum of choosing from a number of available de novo and genome-guided assemblers. None of the popular assembly tools in use today achieve requisite sensitivity, specificity or recovery of full-length transcripts on their own. Here, we present a comprehensive comparative study of the performance of various assemblers. Additionally, we present an approach to combinatorially augment transciptome assembly by using both de novo and genome-guided tools. In our study, we obtained the best recovery and most full-length transcripts with Trinity and TopHat1-Cufflinks, respectively. The sensitivity of the assembly and isoform recovery was superior, without compromising much on the specificity, when transcripts from Trinity were augmented with those from TopHat1-Cufflinks.

Transcriptome characteristics and six alternative expressed genes positively correlated with the phase transition of annual cambial activities in Chinese Fir (Cunninghamia lanceolata (Lamb.) Hook)

Background

The molecular mechanisms that govern cambial activity in angiosperms are well established, but little is known about these molecular mechanisms in gymnosperms. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook), a diploid (2n  = 2x  = 22) gymnosperm, is one of the most important industrial and commercial timber species in China. Here, we performed transcriptome sequencing to identify the repertoire of genes expressed in cambium tissue of Chinese fir.

Methodology/Principal Findings

Based on previous studies, the four stage-specific cambial tissues of Chinese fir were defined using transmission electron microscopy (TEM). In total, 20 million sequencing reads (3.6 Gb) were obtained using Illumina sequencing from Chinese fir cambium tissue collected at active growth stage, with a mean length of 131 bp and a N50 of 90 bp. SOAPdenovo software was used to assemble 62,895 unigenes. These unigenes were further functionally annotated by comparing their sequences to public protein databases. Expression analysis revealed that the altered expression of six homologous genes (ClWOX1, ClWOX4, ClCLV1-like, ClCLV-like, ClCLE12, and ClPIN1-like) correlated positively with changes in cambial activities; moreover, these six genes might be directly involved in cambial function in Chinese fir. Further, the full-length cDNAs and DNAs for ClWOX1 and ClWOX4 were cloned and analyzed.

Conclusions

In this study, a large number of tissue/stage-specific unigene sequences were generated from the active growth stage of Chinese fir cambium. Transcriptome sequencing of Chinese fir not only provides extensive genetic resources for understanding the molecular mechanisms underlying cambial activities in Chinese fir, but also is expected to be an important foundation for future genetic studies of Chinese fir. This study indicates that ClWOX1 and ClWOX4 could be possible reverse genetic target genes for revealing the molecular mechanisms of cambial activities in Chinese fir.

Plant tolerance to high temperature in a changing environment: scientific fundamentals and production of heat stress-tolerant crops

Global warming is predicted to have a general negative effect on plant growth due to the damaging effect of high temperatures on plant development. The increasing threat of climatological extremes including very high temperatures might lead to catastrophic loss of crop productivity and result in wide spread famine. In this review, we assess the impact of global climate change on the agricultural crop production. There is a differential effect of climate change both in terms of geographic location and the crops that will likely show the most extreme reductions in yield as a result of expected extreme fluctuations in temperature and global warming in general. High temperature stress has a wide range of effects on plants in terms of physiology, biochemistry and gene regulation pathways. However, strategies exist to crop improvement for heat stress tolerance. In this review, we present recent advances of research on all these levels of investigation and focus on potential leads that may help to understand more fully the mechanisms that make plants tolerant or susceptible to heat stress. Finally, we review possible procedures and methods which could lead to the generation of new varieties with sustainable yield production, in a world likely to be challenged both by increasing population, higher average temperatures and larger temperature fluctuations.

Single-copy nuclear gene primers for Streptanthus and other Brassicaceae from genomic scans, published data, and ESTs

PREMISE OF THE STUDY

We report 11 primer sets for nine single-copy nuclear genes in Streptanthus and other Thelypodieae (Brassicaceae) and their utility at tribal-level and species-level phylogenetics in this poorly resolved group. •

METHODS AND RESULTS

We selected regions based on a cross-referenced matrix of previous studies and public Brassica expressed sequence tags. To design primers, we used alignments of low-depth-coverage Illumina sequencing of genomic DNA for two species of Brassica mapped onto Arabidopsis thaliana. We report several primer combinations for five regions that consistently amplified a single band and yielded high-quality sequences for at least 70% of the species assayed, and for four additional regions whose utility might be clade specific. •

CONCLUSIONS

Our primers will be useful in improving resolution at shallow depths across the Thelypodieae, and likely in other Brassicaceae.

Comparative genome analysis and phylogenetic relationship of order Liliales insight from the complete plastid genome sequences of two Lilies (Lilium longiflorum and Alstroemeria aurea)

Monocots are one of the most diverse, successful and economically important clades of angiosperms. We attempt to analyse the complete plastid genome sequences of two lilies and their lengths were 152,793bp in Lilium longiflorum (Liliaceae) and 155,510bp in Alstroemeria aurea (Alstroemeriaceae). Phylogenetic analyses were performed for 28 taxa including major lineages of monocots using the sequences of 79 plastid genes for clarifying the phylogenetic relationship of the order Liliales. The sister relationship of Liliales and Asparagales-commelinids was improved with high resolution. Comparative analyses of inter-familial and inter-specific sequence variation were also carried out among three families of Liliaceae, Smilacaceae, and Alstroemeriaceae, and between two Lilium species of L. longflorum and L. superbum. Gene content and order were conserved in the order Liliales except infA loss in Smilax and Alstroemeria. IR boundaries were similar in IRa, however, IRb showed different extension patterns as JLB of Smilax and JSB in Alstroemeria. Ka/Ks ratio was high in matK among the pair-wise comparison of three families and the most variable genes were psaJ, ycf1, rpl32, rpl22, matK, and ccsA among the three families and rps15, rpoA, matK, and ndhF between Lilium.

Evolutionary consequences, constraints and potential of polyploidy in plants

Polyploidy, the possession of more than 2 complete genomes, is a major force in plant evolution known to affect the genetic and genomic constitution and the phenotype of an organism, which will have consequences for its ecology and geography as well as for lineage diversification and speciation. In this review, we discuss phylogenetic patterns in the incidence of polyploidy including possible underlying causes, the role of polyploidy for diversification, the effects of polyploidy on geographical and ecological patterns, and putative underlying mechanisms as well as chromosome evolution and evolution of repetitive DNA following polyploidization. Spurred by technological advances, a lot has been learned about these aspects both in model and increasingly also in nonmodel species. Despite this enormous progress, long-standing questions about polyploidy still cannot be unambiguously answered, due to frequently idiosyncratic outcomes and insufficient integration of different organizational levels (from genes to ecology), but likely this will change in the near future. See also the sister article focusing on animals by Choleva and Janko in this themed issue.

Phylogenomic analyses of nuclear genes reveal the evolutionary relationships within the BEP clade and the evidence of positive selection in Poaceae

BEP clade of the grass family (Poaceae) is composed of three subfamilies, i.e. Bambusoideae, Ehrhartoideae, and Pooideae. Controversies on the phylogenetic relationships among three subfamilies still persist in spite of great efforts. However, previous evidence was mainly provided from plastid genes with only a few nuclear genes utilized. Given different evolutionary histories recorded by plastid and nuclear genes, it is indispensable to uncover their relationships based on nuclear genes. Here, eleven species with whole-sequenced genome and six species with transcriptomic data were included in this study. A total of 121 one-to-one orthologous groups (OGs) were identified and phylogenetic trees were reconstructed by different tree-building methods. Genes which might have undergone positive selection and played important roles in adaptive evolution were also investigated from 314 and 173 one-to-one OGs in two bamboo species and 14 grass species, respectively. Our results support the ((B, P) E) topology with high supporting values. Besides, our findings also indicate that 24 and nine orthologs with statistically significant evidence of positive selection are mainly involved in abiotic and biotic stress response, reproduction and development, plant metabolism and enzyme etc. from two bamboo species and 14 grass species, respectively. In summary, this study demonstrates the power of phylogenomic approach to shed lights on the evolutionary relationships within the BEP clade, and offers valuable insights into adaptive evolution of the grass family.

Comparative transcriptome analysis to investigate the high starch accumulation of duckweed (Landoltia punctata) under nutrient starvation

BACKGROUND

Duckweed can thrive on anthropogenic wastewater and produce tremendous biomass production. Due to its relatively high starch and low lignin percentage, duckweed is a good candidate for bioethanol fermentation. Previous studies have observed that water devoid of nutrients is good for starch accumulation, but its molecular mechanism remains unrevealed.

RESULTS

This study globally analyzed the response to nutrient starvation in order to investigate the starch accumulation in duckweed (Landoltia punctata). L. punctata was transferred from nutrient-rich solution to distilled water and sampled at different time points. Physiological measurements demonstrated that the activity of ADP-glucose pyrophosphorylase, the key enzyme of starch synthesis, as well as the starch percentage in duckweed, increased continuously under nutrient starvation. Samples collected at 0 h, 2 h and 24 h time points respectively were used for comparative gene expression analysis using RNA-Seq. A comprehensive transcriptome, comprising of 74,797 contigs, was constructed by a de novo assembly of the RNA-Seq reads. Gene expression profiling results showed that the expression of some transcripts encoding key enzymes involved in starch biosynthesis was up-regulated, while the expression of transcripts encoding enzymes involved in starch consumption were down-regulated, the expression of some photosynthesis-related transcripts were down-regulated during the first 24 h, and the expression of some transporter transcripts were up-regulated within the first 2 h. Very interestingly, most transcripts encoding key enzymes involved in flavonoid biosynthesis were highly expressed regardless of starvation, while transcripts encoding laccase, the last rate-limiting enzyme of lignifications, exhibited very low expression abundance in all three samples.

CONCLUSION

Our study provides a comprehensive expression profiling of L. punctata under nutrient starvation, which indicates that nutrient starvation down-regulated the global metabolic status, redirects metabolic flux of fixed CO2 into starch synthesis branch resulting in starch accumulation in L. punctata.

Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

The paper reviews the current state of low and single copy nuclear markers that have been applied successfully in plant phylogenetics to date, and discusses case studies highlighting the potential of massively parallel high throughput or next-generation sequencing (NGS) approaches for molecular phylogenetic and evolutionary investigations. The current state, prospects and challenges of specific single- or low-copy plant nuclear markers as well as phylogenomic case studies are presented and evaluated.

Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest

Plants produce a vast array of specialized metabolites, many of which are used as pharmaceuticals, flavors, fragrances, and other high-value fine chemicals. However, most of these compounds occur in non-model plants for which genomic sequence information is not yet available. The production of a large amount of nucleotide sequence data using next-generation technologies is now relatively fast and cost-effective, especially when using the latest Roche-454 and Illumina sequencers with enhanced base-calling accuracy. To investigate specialized metabolite biosynthesis in non-model plants we have established a data-mining framework, employing next-generation sequencing and computational algorithms, to construct and analyze the transcriptomes of 75 non-model plants that produce compounds of interest for biotechnological applications. After sequence assembly an extensive annotation approach was applied to assign functional information to over 800,000 putative transcripts. The annotation is based on direct searches against public databases, including RefSeq and InterPro. Gene Ontology (GO), Enzyme Commission (EC) annotations and associated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps are also collected. As a proof-of-concept, the selection of biosynthetic gene candidates associated with six specialized metabolic pathways is described. A web-based BLAST server has been established to allow public access to assembled transcriptome databases for all 75 plant species of the PhytoMetaSyn Project (www.phytometasyn.ca).

Short read sequencing in studies of natural variation and adaptation

Today's high throughput sequencing approaches, coupled with equally revolutionary advances in bioinformatics, allow us to describe and analyze genomes in unprecedented detail. Short Read Sequencing (SRS) approaches have been especially instrumental in bringing genomic analysis to a wide range of questions and species in plant biology. We can now connect genotypes and phenotypes with greater efficiency, and investigate the molecular basis of natural variation and adaptation in a genomic framework. New and creative applications of SRS and other genomic approaches are not only reshaping how we study natural variation, but also our overall understanding of gene and genome evolution. Here we discuss examples of the application of SRS technologies to the characterization of genetic diversity, genome evolution and adaptation in plants.

Screening of genomic libraries

Microsatellites, or simple sequence repeats (SSRs), have proven to be an important molecular marker in plant genetics and breeding research. The main strategies to obtain these markers can be through genomic DNA and from expressed sequence tags (ESTs) from mRNA/cDNA libraries. Genetic studies using microsatellite markers have increased rapidly because they can be highly polymorphic, codominant markers and they show heterozygous conserved sequences. Here, we describe a methodology to obtain microsatellite using the enrichment library of DNA genomic sequences. This method is highly efficient to development microsatellite markers especially in plants that do not have available ESTs or genome databases. This methodology has been used to enrich SSR marker libraries in Citrus spp., an important tool to genotype germplasm, to select zygotic hybrids, and to saturate genetic maps in breeding programs.

Revisiting the phylogeny of papilionoid legumes: New insights from comprehensively sampled early-branching lineages

PREMISE OF STUDY

Phylogenetic relationships of the papilionoid legumes (Papilionoideae) reveal that the early branches are more highly diverse in floral morphology than are other clades of Papilionoideae. This study attempts for the first time to comprehensively sample the early-branching clades of this economically and ecologically important legume subfamily and thus to resolve relationships among them. •

METHODS

Parsimony and Bayesian phylogenetic analyses of the plastid matK and trnL intron sequences included 29 genera not yet sampled in matK phylogenies of the Papilionoideae, 11 of which were sampled for DNA sequence data for the first time. •

KEY RESULTS

The comprehensively sampled matK phylogeny better resolved the deep-branching relationships and increased support for many clades within Papilionoideae. The potentially earliest-branching papilionoid clade does not include any genus traditionally assigned to tribe Swartzieae. Dipterygeae is monophyletic with the inclusion of Monopteryx. The genera Aldina and Amphimas represent two of the nine main but as yet unresolved lineages comprising the large 50-kb inversion clade within papilionoids. The quinolizidine-alkaloid-accumulating genistoid clade is expanded to include a strongly supported subclade containing Ormosia and the previously unplaced Clathrotropis s.s., Panurea, and Spirotropis. Camoensia is the first-branching genus of the core genistoids. •

CONCLUSIONS

The well-resolved phylogeny of the earliest-branching papilionoids generated in this study will greatly facilitate the efforts to redefine and stabilize the classification of this legume subfamily. Many key floral traits did not often predict phylogenetic relationships, so comparative studies on floral evolution and plant-animal interactions, for example, should also benefit from this study.

Next generation sequencing and de novo transcriptome analysis of Costus pictus D. Don, a non-model plant with potent anti-diabetic properties

Background

Phyto-remedies for diabetic control are popular among patients with Type II Diabetes mellitus (DM), in addition to other diabetic control measures. A number of plant species are known to possess diabetic control properties. Costus pictus D. Don is popularly known as “Insulin Plant” in Southern India whose leaves have been reported to increase insulin pools in blood plasma. Next Generation Sequencing is employed as a powerful tool for identifying molecular signatures in the transcriptome related to physiological functions of plant tissues. We sequenced the leaf transcriptome of C. pictus using Illumina reversible dye terminator sequencing technology and used combination of bioinformatics tools for identifying transcripts related to anti-diabetic properties of C. pictus.

Results

A total of 55,006 transcripts were identified, of which 69.15% transcripts could be annotated. We identified transcripts related to pathways of bixin biosynthesis and geraniol and geranial biosynthesis as major transcripts from the class of isoprenoid secondary metabolites and validated the presence of putative norbixin methyltransferase, a precursor of Bixin. The transcripts encoding these terpenoids are known to be Peroxisome Proliferator-Activated Receptor (PPAR) agonists and anti-glycation agents. Sequential extraction and High Performance Liquid Chromatography (HPLC) confirmed the presence of bixin in C. pictus methanolic extracts. Another significant transcript identified in relation to anti-diabetic, anti-obesity and immuno-modulation is of Abscisic Acid biosynthetic pathway. We also report many other transcripts for the biosynthesis of antitumor, anti-oxidant and antimicrobial metabolites of C. pictus leaves.

Conclusion

Solid molecular signatures (transcripts related to bixin, abscisic acid, and geranial and geraniol biosynthesis) for the anti-diabetic properties of C. pictus leaves and vital clues related to the other phytochemical functions like antitumor, anti-oxidant, immuno-modulatory, anti-microbial and anti-malarial properties through the secondary metabolite pathway annotations are reported. The data provided will be of immense help to researchers working in the treatment of DM using herbal therapies.

The lemur revolution starts now: the genomic coming of age for a non-model organism

Morris Goodman was a revolutionary. Together with a mere handful of like-minded scientists, Morris established himself as a leader in the molecular phylogenetic revolution of the 1960s. The effects of this revolution are most evident in this journal, which he founded in 1992. Happily for lemur biologists, one of Morris Goodman's primary interests was in reconstructing the phylogeny of the primates, including the tooth-combed Lorisifomes of Africa and Asia, and the Lemuriformes of Madagascar (collectively referred to as the suborder Strepsirrhini). This paper traces the development of molecular phylogenetic and evolutionary genetic trends and methods over the 50-year expanse of Morris Goodman's career, particularly as they apply to our understanding of lemuriform phylogeny, biogeography, and biology. Notably, this perspective reveals that the lemuriform genome is sufficiently rich in phylogenetic signal such that the very earliest molecular phylogenetic studies - many of which were conducted by Goodman himself - have been validated by contemporary studies that have exploited advanced computational methods applied to phylogenomic scale data; studies that were beyond imagining in the earliest days of phylogeny reconstruction. Nonetheless, the frontier still beckons. New technologies for gathering and analyzing genomic data will allow investigators to build upon what can now be considered a nearly-known phylogeny of the Lemuriformes in order to ask innovative questions about the evolutionary mechanisms that generate and maintain the extraordinary breadth and depth of biological diversity within this remarkable clade of primates.