Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach

Fire-mediated germination syndromes in Leucadendron (Proteaceae) and their functional correlates

A mechanistic understanding of fire-driven seedling recruitment is essential for effective conservation management of fire-prone vegetation, such as South African fynbos, especially with rare and threatened taxa. The genus Leucadendron (Proteaceae) is an ideal candidate for comparative germination studies, comprising 85 species with a mixture of contrasting life-history traits (killed by fire vs able to resprout; serotinous vs geosporous) and seed morphologies (nutlets vs winged achenes). Individual and combined effects of heat and smoke on seed germination of 40 species were quantified in the laboratory, and Bayesian inference applied to distinguish biologically meaningful treatment effects from non-zero, but biologically trivial, effects. Three germination syndromes were identified based on whether germination was dependent on, enhanced by, or independent of direct fire cues (heat and smoke). Seed storage location was the most reliable predictor of germination syndromes, with soil-stored seeds c. 80% more likely to respond to direct fire cues (primarily smoke) than canopy-stored seeds. Notable exceptions were L. linifolium, with an absolute requirement for smoke to germinate (the third serotinous species so reported), and two other serotinous species with smoke-enhanced germination. Nutlet-bearing species, whether serotinous or geosporous, were c. 70% more likely to respond to fire cues than winged seeds, but there was no evidence for an effect of phylogeny or persistence strategy on germination. This comprehensive account of seed germination characteristics and identification of germination syndromes and their predictors, supports propagation, conservation and restoration initiatives in this iconic fynbos genus and other fire-prone shrubs with canopy or soil-stored seeds.

Identifying genetic markers for a range of phylogenetic utility-From species to family level

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies.

Unraveling the Developmental and Genetic Mechanisms Underpinning Floral Architecture in Proteaceae

Proteaceae are a basal eudicot family with a highly conserved floral groundplan but which displays considerable variation in other aspects of floral and inflorescence morphology. Their morphological diversity and phylogenetic position make them good candidates for understanding the evolution of floral architecture, in particular the question of the homology of the undifferentiated perianth with the differentiated perianth of core eudicots, and the mechanisms underlying the repeated evolution of zygomorphy. In this paper, we combine a morphological approach to explore floral ontogenesis and a transcriptomic approach to access the genes involved in floral organ identity and development, focusing on Grevillea juniperina, a species from subfamily Grevilleoideae. We present developmental data for Grevillea juniperina and three additional species that differ in their floral symmetry using stereomicroscopy, SEM and High Resolution X-Ray Computed Tomography. We find that the adnation of stamens to tepals takes place at early developmental stages, and that the establishment of bilateral symmetry coincides with the asymmetrical growth of the single carpel. To set a framework for understanding the genetic basis of floral development in Proteaceae, we generated and annotated de novo a reference leaf/flower transcriptome from Grevillea juniperina. We found Grevillea homologs of all lineages of MADS-box genes involved in floral organ identity. Using Arabidopsis thaliana gene expression data as a reference, we found homologs of other genes involved in floral development in the transcriptome of G. juniperina. We also found at least 21 class I and class II TCP genes, a gene family involved in the regulation of growth processes, including floral symmetry. The expression patterns of a set of floral genes obtained from the transcriptome were characterized during floral development to assess their organ specificity and asymmetry of expression.

A New Perspective on Polyploid Fragaria (Strawberry) Genome Composition Based on Large-Scale, Multi-Locus Phylogenetic Analysis

The subgenomic compositions of the octoploid (2n = 8× = 56) strawberry (Fragaria) species, including the economically important cultivated species Fragaria x ananassa, have been a topic of long-standing interest. Phylogenomic approaches utilizing next-generation sequencing technologies offer a new window into species relationships and the subgenomic compositions of polyploids. We have conducted a large-scale phylogenetic analysis of Fragaria (strawberry) species using the Fluidigm Access Array system and 454 sequencing platform. About 24 single-copy or low-copy nuclear genes distributed across the genome were amplified and sequenced from 96 genomic DNA samples representing 16 Fragaria species from diploid (2×) to decaploid (10×), including the most extensive sampling of octoploid taxa yet reported. Individual gene trees were constructed by different tree-building methods. Mosaic genomic structures of diploid Fragaria species consisting of sequences at different phylogenetic positions were observed. Our findings support the presence in octoploid species of genetic signatures from at least five diploid ancestors (F. vesca, F. iinumae, F. bucharica, F. viridis, and at least one additional allele contributor of unknown identity), and questions the extent to which distinct subgenomes are preserved over evolutionary time in the allopolyploid Fragaria species. In addition, our data support divergence between the two wild octoploid species, F. virginiana and F. chiloensis.

Life-History Traits Evolved Jointly with Climatic Niche and Disturbance Regime in the Genus Leucadendron (Proteaceae)

Organisms have evolved a diversity of life-history strategies to cope with variation in their environment. Persistence as adults and/or seeds across recruitment events allows species to dampen the effects of environmental fluctuations. The evolution of life cycles with overlapping generations should thus permit the colonization of environments with uncertain recruitment. We tested this hypothesis in Leucadendron (Proteaceae), a genus with high functional diversity native to fire-prone habitats in the South African fynbos. We analyzed the joint evolution of life-history traits (adult survival and seed-bank strategies) and ecological niches (climate and fire regime), using comparative methods and accounting for various sources of uncertainty. In the fynbos, species with canopy seed banks that are unable to survive fire as adults display nonoverlapping generations. In contrast, resprouters with an underground seed bank may be less threatened by extreme climatic events and fire intervals, given their iteroparity and long-lasting seed bank. Life cycles with nonoverlapping generations indeed jointly evolved with niches with less exposure to frost but not with those with less exposure to drought. Canopy seed banks jointly evolved with niches with more predictable fire return, compared to underground seed banks. The evolution of extraordinary functional diversity among fynbos plants thus reflects, at least in part, the diversity of both climates and fire regimes in this region.

Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae)

PREMISE OF THE STUDY

We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline.

METHODS AND RESULTS

Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a "pseudoreference" closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae.

CONCLUSIONS

This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.

Anchored phylogenomics improves the resolution of evolutionary relationships in the rapid radiation of Protea L

PREMISE OF THE STUDY

Estimating phylogenetic relationships in relatively recent evolutionary radiations is challenging, especially if short branches associated with recent divergence result in multiple gene tree histories. We combine anchored enrichment next-generation sequencing with species tree analyses to produce a robust estimate of phylogenetic relationships in the genus Protea (Proteaceae), an iconic radiation in South Africa.

METHODS

We sampled multiple individuals within 59 out of 112 species of Protea and 6 outgroup species for a total of 163 individuals, and obtained sequences for 498 low-copy, orthologous nuclear loci using anchored phylogenomics. We compare several approaches for building species trees, and explore gene tree-species tree discrepancies to determine whether poor phylogenetic resolution reflects a lack of informative sites, incomplete lineage sorting, or hybridization.

KEY RESULTS

Phylogenetic estimates from species tree approaches are similar to one another and recover previously well-supported clades within Protea, in addition to providing well-supported phylogenetic hypotheses for previously poorly resolved intrageneric relationships. Individual gene trees are markedly different from one another and from species trees. Nonetheless, analyses indicate that differences among gene trees occur primarily concerning clades supported by short branches.

CONCLUSIONS

Species tree methods using hundreds of nuclear loci provided strong support for many previously unresolved relationships in the radiation of the genus Protea. In cases where support for particular relationships remains low, these appear to arise from few informative sites and lack of information rather than strongly supported disagreement among gene trees.

A genome-scale mining strategy for recovering novel rapidly-evolving nuclear single-copy genes for addressing shallow-scale phylogenetics in Hydrangea

Background

Identifying orthologous molecular markers that potentially resolve relationships at and below species level has been a major challenge in molecular phylogenetics over the past decade. Non-coding regions of nuclear low- or single-copy markers are a vast and promising source of data providing information for shallow-scale phylogenetics. Taking advantage of public transcriptome data from the One Thousand Plant Project (1KP), we developed a genome-scale mining strategy for recovering potentially orthologous single-copy markers to address low-scale phylogenetics. Our marker design targeted the amplification of intron-rich nuclear single-copy regions from genomic DNA. As a case study we used Hydrangea section Cornidia, one of the most recently diverged lineages within Hydrangeaceae (Cornales), for comparing the performance of three of these nuclear markers to other “fast” evolving plastid markers.

Results

Our data mining and filtering process retrieved 73 putative nuclear single-copy genes which are potentially useful for resolving phylogenetic relationships at a range of divergence depths within Cornales. The three assessed nuclear markers showed considerably more phylogenetic signal for shallow evolutionary depths than conventional plastid markers. Phylogenetic signal in plastid markers increased less markedly towards deeper evolutionary divergences. Potential phylogenetic noise introduced by nuclear markers was lower than their respective phylogenetic signal across all evolutionary depths. In contrast, plastid markers showed higher probabilities for introducing phylogenetic noise than signal at the deepest evolutionary divergences within the tribe Hydrangeeae (Hydrangeaceae).

Conclusions

While nuclear single-copy markers are highly informative for shallow evolutionary depths without introducing phylogenetic noise, plastid markers might be more appropriate for resolving deeper-level divergences such as the backbone relationships of the Hydrangeaceae family and deeper, at which non-coding parts of nuclear markers could potentially introduce noise due to elevated rates of evolution. The herein developed and demonstrated transcriptome based mining strategy has a great potential for the design of novel and highly informative nuclear markers for a range of plant groups and evolutionary scales.

Electronic supplementary material

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

Efficient Detection of Novel Nuclear Markers for Brassicaceae by Transcriptome Sequencing

The lack of DNA sequence information for most non-model organisms impairs the design of primers that are universally applicable for the study of molecular polymorphisms in nuclear markers. Next-generation sequencing (NGS) techniques nowadays provide a powerful approach to overcome this limitation. We present a flexible and inexpensive method to identify large numbers of nuclear primer pairs that amplify in most Brassicaceae species. We first obtained and mapped NGS transcriptome sequencing reads from two of the distantly related Brassicaceae species, Cardamine hirsuta and Arabis alpina, onto the Arabidopsis thaliana reference genome, and then identified short conserved sequence motifs among the three species bioinformatically. From these, primer pairs to amplify coding regions (nuclear protein coding loci, NPCL) and exon-primed intron-crossing sequences (EPIC) were developed. We identified 2,334 universally applicable primer pairs, targeting 1,164 genes, which provide a large pool of markers as readily usable genomic resource that will help addressing novel questions in the Brassicaceae family. Testing a subset of the newly designed nuclear primer pairs revealed that a great majority yielded a single amplicon in all of the 30 investigated Brassicaceae taxa. Sequence analysis and phylogenetic reconstruction with a subset of these markers on different levels of phylogenetic divergence in the mustard family were compared with previous studies. The results corroborate the usefulness of the newly developed primer pairs, e.g., for phylogenetic analyses or population genetic studies. Thus, our method provides a cost-effective approach for designing nuclear loci across a broad range of taxa and is compatible with current NGS technologies.

MarkerMiner 1.0: A new application for phylogenetic marker development using angiosperm transcriptomes

PREMISE OF THE STUDY

Targeted sequencing using next-generation sequencing (NGS) platforms offers enormous potential for plant systematics by enabling economical acquisition of multilocus data sets that can resolve difficult phylogenetic problems. However, because discovery of single-copy nuclear (SCN) loci from NGS data requires both bioinformatics skills and access to high-performance computing resources, the application of NGS data has been limited.

METHODS AND RESULTS

We developed MarkerMiner 1.0, a fully automated, open-access bioinformatic workflow and application for discovery of SCN loci in angiosperms. Our new tool identified as many as 1993 SCN loci from transcriptomic data sampled as part of four independent test cases representing marker development projects at different phylogenetic scales.

CONCLUSIONS

MarkerMiner is an easy-to-use and effective tool for discovery of putative SCN loci. It can be run locally or via the Web, and its tabular and alignment outputs facilitate efficient downstream assessments of phylogenetic utility, locus selection, intron-exon boundary prediction, and primer or probe development.

Conserved genetic regions across angiosperms as tools to develop single-copy nuclear markers in gymnosperms: an example using cycads

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single-copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis-Populus-Vitis-Oryza Shared Single-Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single-copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.