Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

The complete chloroplast genome of Malva wigandii (Alef.) M.F. Ray (Malvaceae, Malvoideae).

The complete chloroplast genome sequence of wild sea mallow Malva wigandii (=Lavatera maritima) was determined and characterized in this study. The genome is 158,162 bp long, containing a pair of inverted repeats regions (IRs) of 25,166 bp, which are separated by a large single-copy region of 86,860 bp and a small single-copy region of 20,970 bp. The sea mallow chloroplast genome has 131 known genes, including 85 protein-coding genes, eight ribosomal RNA genes, and 37 tRNA genes. The phylogenomic analysis showed that M. wigandii forms a cluster with Althaea officinalis with a strong bootstrap support and is sister to sequences belonging to the tribe Gossypieae. All of them are grouped in a lineage with other members of the subfamily Malvoideae. This newly sequenced chloroplast genome sequence provides useful genetic information to explore the origin and evolution of the Mediterranean radiation that gave rise to the generic alliance of Malva.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event.

The consequences of the Cretaceous-Paleogene (K-Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous-Paleogene (K-Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events].

Molecular phylogenetics of Lotus (Leguminosae) with emphasis in the tempo and patterns of colonization in the Macaronesian region.

With a wide distribution range including Europe and Asia, Lotus (Leguminosae) represents the largest genus within Loteae. It is particularly diverse in the Mediterreanean region and in the five archipelagos of Macaronesia (Atlantic Ocean). However, little is known about the relationships among the 14 sections currently recognized within Lotus and about the timing and patterns of its colonization in the Macaronesian region. In this investigation, we use four DNA regions (nuclear ribosomal ITS plus three plastid regions) in the most comprehensive sampling of Lotus species to date (some endemic species within the Canary Islands were poorly represented in previous phylogenetic analyses) to infer relationships within this genus and to establish patterns of colonization in Macaronesia. Divergence time estimates and habitat reconstruction analyses indicate that Lotus likely diverged about 7.86 Ma from its sister group, but all colonization events to Macaronesia occurred more recently (ranging from the last 0.23 to 2.70 Ma). The diversification of Lotus in Macaronesia involved between four and six independent colonization events from four sections currently distributed in Africa and Europe. A major aspect shaping the current distribution of taxa involved intra-island colonization of mainly new habitats and inter-island colonization of mostly similar habitats, with Gran Canaria and Tenerife as the major sources of diversification and of further colonization events. Section Pedrosia is the most diverse in terms of colonization events, number of species, and habitat heterogeneity, including a back-colonization event to the continent. Subsections within Pedrosia radiated into diverse habitat types recently (late Pleistocene, ca 0.23-0.29 Ma) and additional molecular markers and sampling would be necessary to understand the most recent dispersal events of this group within the Canary Islands and Cape Verde.

Plastome phylogeography in two African rain forest legume trees reveals that Dahomey Gap populations originate from the Cameroon volcanic line.

Paleo-environmental data show that the distribution of African rain forests was affected by Quaternary climate changes. In particular, the Dahomey Gap (DG) - a 200 km wide savanna corridor currently separating the West African and Central African rain forest blocks and containing relict rain forest fragments - was forested during the mid-Holocene and possibly during previous interglacial periods, whereas it was dominated by open vegetation (savanna) during glacial periods. Genetic signatures of past population fragmentation and demographic changes have been found in some African forest plant species using nuclear markers, but such events appear not to have been synchronous or shared across species. To better understand the colonization history of the DG by rain forest trees through seed dispersal, the plastid genomes of two widespread African forest legume trees, Anthonotha macrophylla and Distemonanthus benthamianus, were sequenced in 47 individuals for each species, providing unprecedented phylogenetic resolution of their maternal lineages (857 and 115 SNPs, respectively). Both species exhibit distinct lineages separating three regions: 1. Upper Guinea (UG, i.e. the West African forest block), 2. the area ranging from the DG to the Cameroon volcanic line (CVL), and 3. Lower Guinea (LG, the western part of the Central African forest block) where three lineages co-occur. In both species, the DG populations (including southern Nigeria west of Cross River) exhibit much lower genetic diversity than UG and LG populations, and their plastid lineages originate from the CVL, confirming the role of the CVL as an ancient forest refuge. Despite the similar phylogeographic structures displayed by A. macrophylla and D. benthamianus, molecular dating indicates very contrasting ages of lineage divergence (UG diverged from LG since c. 7 Ma and 0.7 Ma, respectively) and DG colonization (probably following the Mid Pleistocene Transition and the Last Glacial Maximum, respectively). The stability of forest refuge areas and repeated similar forest shrinking/expanding events during successive glacial periods might explain why similar phylogeographic patterns can be generated over contrasting timescales.

Large-scale genomic sequence data resolve the deepest divergences in the legume phylogeny and support a near-simultaneous evolutionary origin of all six subfamilies.

Phylogenomics is increasingly used to infer deep-branching relationships while revealing the complexity of evolutionary processes such as incomplete lineage sorting, hybridization/introgression and polyploidization. We investigate the deep-branching relationships among subfamilies of the Leguminosae (or Fabaceae), the third largest angiosperm family. Despite their ecological and economic importance, a robust phylogenetic framework for legumes based on genome-scale sequence data is lacking. We generated alignments of 72 chloroplast genes and 7621 homologous nuclear-encoded proteins, for 157 and 76 taxa, respectively. We analysed these with maximum likelihood, Bayesian inference, and a multispecies coalescent summary method, and evaluated support for alternative topologies across gene trees. We resolve the deepest divergences in the legume phylogeny despite lack of phylogenetic signal across all chloroplast genes and the majority of nuclear genes. Strongly supported conflict in the remainder of nuclear genes is suggestive of incomplete lineage sorting. All six subfamilies originated nearly simultaneously, suggesting that the prevailing view of some subfamilies as 'basal' or 'early-diverging' with respect to others should be abandoned, which has important implications for understanding the evolution of legume diversity and traits. Our study highlights the limits of phylogenetic resolution in relation to rapid successive speciation.

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

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

Phylogenomic analyses reveal an exceptionally high number of evolutionary shifts in a florally diverse clade of African legumes.

Detarioideae is well known for its high diversity of floral traits, including flower symmetry, number of organs, and petal size and morphology. This diversity has been characterized and studied at higher taxonomic levels, but limited analyses have been performed among closely related genera with contrasting floral traits due to the lack of fully resolved phylogenetic relationships. Here, we used four representative transcriptomes to develop an exome capture (target enrichment) bait for the entire subfamily and applied it to the Anthonotha clade using a complete data set (61 specimens) representing all extant floral diversity. Our phylogenetic analyses recovered congruent topologies using ML and Bayesian methods. Anthonotha was recovered as monophyletic contrary to the remaining three genera (Englerodendron, Isomacrolobium and Pseudomacrolobium), which form a monophyletic group sister to Anthonotha. We inferred a total of 35 transitions for the seven floral traits (pertaining to flower symmetry, petals, stamens and staminodes) that we analyzed, suggesting that at least 30% of the species in this group display transitions from the ancestral condition reconstructed for the Anthonotha clade. The main transitions were towards a reduction in the number of organs (petals, stamens and staminodes). Despite the high number of transitions, our analyses indicate that the seven characters are evolving independently in these lineages. Petal morphology is the most labile floral trait with a total of seven independent transitions in number and seven independent transitions to modification in petal types. The diverse petal morphology along the dorsoventral axis of symmetry within the flower is not associated with differences at the micromorphology of petal surface, suggesting that in this group all petals within the flower might possess the same petal identity at the molecular level. Our results provide a solid evolutionary framework for further detailed analyses of the molecular basis of petal identity.

Comparative analysis of two sister Erythrophleum species (Leguminosae) reveal contrasting transcriptome-wide responses to early drought stress.

With the ongoing climate change, African rainforests are expected to experience severe drought events in the future. In Africa, the tropical genus Erythrophleum (Fabaceae) includes two forest sister timber tree species displaying contrasting geographical distributions. Erythrophleum ivorense is adapted to wet evergreen Guineo-Congolian forests, whereas E. suaveolens occurs in a wider range of climates, being found in moist dense forests but also in gallery forests under a relatively drier climate. This geographical distribution pattern suggests that the two species might cope differently to drought at the genomic level. Yet, the genetic basis of tolerance response to drought stress in both species is still uncharacterized. To bridge this gap, we performed an RNA-seq approach on seedlings from each species to monitor their transcriptional responses at different levels of drought stress (0, 2 and 6 weeks after stopping watering seedlings). Monitoring of wilting symptoms revealed that E. suaveolens displayed an earlier phenotypic response to drought stress than E. ivorense. At the transcriptomic level, results revealed 2020 (1204 down-regulated/816 up-regulated) and 1495 differentially expressed genes (DEGs) in response to drought stress from a total of 67,432 and 66,605 contigs assembled in E. ivorense and E. suaveolens, respectively. After identifying 30,374 orthologs between species, we found that only 7 of them were DEGs shared between species, while 587 and 458 were differentially expressed only in E. ivorense or E. suaveolens, respectively. GO and KEGG enrichment analysis revealed that the two species differ in terms of significantly regulated pathways as well as the number and expression profile of DEGs (Up/Down) associated with each pathway in the two stress stages. Our results suggested that the two studied species react differently to drought. E. suaveolens seems displaying a prompt response to drought at its early stage strengthened by the down-regulation of many DEGs encoding for signaling and metabolism-related pathways. A considerable up-regulation of these pathways was also found in E. ivorense at the late stage of drought, suggesting this species may be a late responder. Overall, our data may serve as basis for further understanding the genetic control of drought tolerance in tropical trees and favor the selection of crucial genes for genetically enhancing drought resistance.

Genome-Enhanced Detection and Identification (GEDI) of plant pathogens.

Plant diseases caused by fungi and Oomycetes represent worldwide threats to crops and forest ecosystems. Effective prevention and appropriate management of emerging diseases rely on rapid detection and identification of the causal pathogens. The increase in genomic resources makes it possible to generate novel genome-enhanced DNA detection assays that can exploit whole genomes to discover candidate genes for pathogen detection. A pipeline was developed to identify genome regions that discriminate taxa or groups of taxa and can be converted into PCR assays. The modular pipeline is comprised of four components: (1) selection and genome sequencing of phylogenetically related taxa, (2) identification of clusters of orthologous genes, (3) elimination of false positives by filtering, and (4) assay design. This pipeline was applied to some of the most important plant pathogens across three broad taxonomic groups: Phytophthoras (Stramenopiles, Oomycota), Dothideomycetes (Fungi, Ascomycota) and Pucciniales (Fungi, Basidiomycota). Comparison of 73 fungal and Oomycete genomes led the discovery of 5,939 gene clusters that were unique to the targeted taxa and an additional 535 that were common at higher taxonomic levels. Approximately 28% of the 299 tested were converted into qPCR assays that met our set of specificity criteria. This work demonstrates that a genome-wide approach can efficiently identify multiple taxon-specific genome regions that can be converted into highly specific PCR assays. The possibility to easily obtain multiple alternative regions to design highly specific qPCR assays should be of great help in tackling challenging cases for which higher taxon-resolution is needed.

Dissecting the 'bacon and eggs' phenotype: transcriptomics of post-anthesis colour change in Lotus.

Background and Aims: Post-anthesis colour change (PACC) is widely thought to be an adaptation to signal floral suitability to pollinators. Lotus filicaulis and Lotus sessilifolius are insect-pollinated herbaceous legumes with flowers that open yellow, shift to orange and finally red. This study examines the molecular basis for floral colour change in these Lotus species. Methods: Lotus filicaulis was cultivated in a glasshouse from which pollinating insects (bees) were excluded, and the rate of colour change was recorded in both unpollinated and manually pollinated flowers. Unpollinated flowers from both the yellow stage and the red stage were sampled for sequencing. The transcriptomes of L. filicaulis and L. sessilifolius of both colour stages were analysed for differentially expressed genes and enriched ontologies. Key Results: The rate of progression through PACC doubled when L. filicaulis was hand-pollinated. De novo assembly of RNA-Seq reads from non-model Lotus species outperformed heterologous alignment of reads to the L. japonicus genome. Differential expression analysis suggested that the carotenoid biosynthetic pathway is upregulated at anthesis while the flavonoid biosynthetic pathway is upregulated with the onset of PACC in L. filicaulis and L. sessilifolius . Conclusion: Pollination significantly accelerates PACC in L. filicaulis , consistent with the hypothesis that PACC increases pollination efficiency by directing pollinators to unpollinated flowers. RNA-Seq results show the synchronized upregulation of the entire cyanidin biosynthesis pathway in the red stage of PACC in L. filicaulis and L. sessilifolius . The genes implicated offer the basis for further investigations into how gene families, transcription factors and related pathways are likely to be involved in PACC.

A Multi-Species TaqMan PCR Assay for the Identification of Asian Gypsy Moths (Lymantria spp.) and Other Invasive Lymantriines of Biosecurity Concern to North America.

Preventing the introduction and establishment of forest invasive alien species (FIAS) such as the Asian gypsy moth (AGM) is a high-priority goal for countries with extensive forest resources such as Canada. The name AGM designates a group of closely related Lymantria species (Lepidoptera: Erebidae: Lymantriinae) comprising two L. dispar subspecies (L. dispar asiatica, L. dispar japonica) and three closely related Lymantria species (L. umbrosa, L. albescens, L. postalba), all considered potential FIAS in North America. Ships entering Canadian ports are inspected for the presence of suspicious gypsy moth eggs, but those of AGM are impossible to distinguish from eggs of innocuous Lymantria species. To assist regulatory agencies in their identification of these insects, we designed a suite of TaqMan® assays that provide significant improvements over existing molecular assays targeting AGM. The assays presented here can identify all three L. dispar subspecies (including the European gypsy moth, L. dispar dispar), the three other Lymantria species comprising the AGM complex, plus five additional Lymantria species that pose a threat to forests in North America. The suite of assays is built as a "molecular key" (analogous to a taxonomic key) and involves several parallel singleplex and multiplex qPCR reactions. Each reaction uses a combination of primers and probes designed to separate taxa through discriminatory annealing. The success of these assays is based on the presence of single nucleotide polymorphisms (SNPs) in the 5' region of mitochondrial cytochrome c oxidase I (COI) or in its longer, 3' region, as well as on the presence of an indel in the "FS1" nuclear marker, generating North American and Asian alleles, used here to assess Asian introgression into L. dispar dispar. These assays have the advantage of providing rapid and accurate identification of ten Lymantria species and subspecies considered potential FIAS.

Pollinator shifts drive petal epidermal evolution on the Macaronesian Islands bird-flowered species.

Pollinator shifts are considered to drive floral trait evolution, yet little is still known about the modifications of petal epidermal surface at a biogeographic region scale. Here we investigated how independent shifts from insects to passerine birds in the Macaronesian Islands consistently modified this floral trait (i.e. absence of papillate cells). Using current phylogenies and extensive evidence from field observations, we selected a total of 81 plant species and subspecies for petal microscopy and comparative analysis, including 19 of the 23 insular species pollinated by opportunistic passerine birds (Macaronesian bird-flowered element). Species relying on passerine birds as the most effective pollinators (bird-pollinated) independently evolved at least five times and in all instances associated with a loss of papillate cells, whereas species with a mixed pollination system (birds plus insects and/or other vertebrates) evolved at least five times in Macaronesia and papillate cells were lost in only 25% of these transitions. Our findings suggest that petal micromorphology is a labile trait during pollinator shifts and that papillate cells tend to be absent on those species where pollinators have limited mechanical interaction with flowers, including opportunistic passerine birds that forage by hovering or from the ground.

DNA barcodes successfully identified Macaronesian Lotus (Leguminosae) species within early diverged lineages of Cape Verde and mainland Africa.

Plant DNA barcoding currently relies on the application of a two-locus combination, matK + rbcL. Despite the universality of these two gene regions across plants, it is suspected that this combination might not have sufficient variation to discriminate closely related species. In this study, we tested the performance of this two-locus plant barcode along with the additional plastid regions trnH-psbA, rpoC1 and rpoB and the nuclear region internal transcribed spacer (nrITS) in a group of 38 species of Lotus from the Macaronesian region. The group has radiated into the five archipelagos within this region from mid-Miocene to early Pleistocene, and thus provides both early divergent and recent radiations that pose a particularly difficult challenge for barcoding. The group also has 10 species considered under different levels of conservation concern. We found different levels of species discrimination depending on the age of the lineages. We obtained 100 % of the species identification from mainland Africa and Cape Verde when all six regions were combined. These lineages radiated >4.5 Mya; however, in the most recent radiations from the end of the Pliocene to the mid-Pleistocene (3.5-1.5 Mya), only 30 % of the species were identified. Of the regions examined, the intergenic region trnH-psbA was the most variable and had the greatest discriminatory power (18 %) of the plastid regions when analysed alone. The nrITS region was the best region when analysed alone with a discriminatory power of 26 % of the species. Overall, we identified 52 % of the species and 30 % of the endangered or threatened species within this group when all six regions were combined. Our results are consistent with those of other studies that indicate that additional approaches to barcoding will be needed in recently evolved groups, such as the inclusion of faster evolving regions from the nuclear genome.

Single-nucleotide polymorphism discovery in Leptographium longiclavatum, a mountain pine beetle-associated symbiotic fungus, using whole-genome resequencing.

Single-nucleotide polymorphisms (SNPs) are rapidly becoming the standard markers in population genomics studies; however, their use in nonmodel organisms is limited due to the lack of cost-effective approaches to uncover genome-wide variation, and the large number of individuals needed in the screening process to reduce ascertainment bias. To discover SNPs for population genomics studies in the fungal symbionts of the mountain pine beetle (MPB), we developed a road map to discover SNPs and to produce a genotyping platform. We undertook a whole-genome sequencing approach of Leptographium longiclavatum in combination with available genomics resources of another MPB symbiont, Grosmannia clavigera. We sequenced 71 individuals pooled into four groups using the Illumina sequencing technology. We generated between 27 and 30 million reads of 75 bp that resulted in a total of 1, 181 contigs longer than 2 kb and an assembled genome size of 28.9 Mb (N50 = 48 kb, average depth = 125x). A total of 9052 proteins were annotated, and between 9531 and 17,266 SNPs were identified in the four pools. A subset of 206 genes (containing 574 SNPs, 11% false positives) was used to develop a genotyping platform for this species. Using this roadmap, we developed a genotyping assay with a total of 147 SNPs located in 121 genes using the Illumina(®) Sequenom iPLEX Gold. Our preliminary genotyping (success rate = 85%) of 304 individuals from 36 populations supports the utility of this approach for population genomics studies in other MPB fungal symbionts and other fungal nonmodel species.