Comparative transcriptomics of a complex of four European pine species

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

Full-Length Transcriptome Sequencing-Based Analysis of Pinus sylvestris var. mongolica in Response to Sirex noctilio Venom

Simple Summary

Sirex noctilio, as a devastating international forestry quarantine pest whose venom can cause a series of physiological changes in the host plants, such as needle wilting, yellowing, decreased transpiration rate and increased respiration rate, etc. In this study, a full-length reference transcript of Pinus sylvestris var. mongolica was constructed by combining second- and third-generation transcriptome sequencing technologies. We also identified the specific expression genes and transcription factors of P. sylvestris var. mongolica under S. noctilio venom and wounding stress. S. noctilio venom mainly induced the expression of genes related to ROS, GAPDH and GPX, and mechanical damage mainly induced the photosynthesis−related genes. The results provide a better understanding of the molecular regulation of pine trees in response to S. noctilio venom.

Abstract

Sirex noctilio is a major international quarantine pest that recently emerged in northeast China to specifically invade conifers. During female oviposition, venom is injected into the host together with its symbiotic fungus to alter the normal Pinus physiology and weaken or even kill the tree. In China, the Mongolian pine (Pinus sylvestris var. mongolica), an important wind-proof and sand-fixing species, is the unique host of S. noctilio. To explore the interplay between S. noctilio venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with S. noctilio venom. The analysis was performed at 12 h, 24 h and 72 h. PacBio ISO-seq was used and integrated with RNA-seq to construct an accurate full-length transcriptomic database. We obtained 52,963 high-precision unigenes, consisting of 48,654 (91.86%) unigenes that were BLASTed to known sequences in the public database and 4309 unigenes without any annotation information, which were presumed to be new genes. The number of differentially expressed genes (DEGs) increased with the treatment time, and the DEGs were most abundant at 72 h. A total of 706 inoculation-specific DEGs (475 upregulated and 231 downregulated) and 387 wounding-specific DEGs (183 upregulated and 204 downregulated) were identified compared with the control. Under venom stress, we identified 6 DEGs associated with reactive oxygen species (ROS) and 20 resistance genes in Mongolian pine. Overall, 52 transcription factors (TFs) were found under venom stress, 45 of which belonged to the AP2/ERF TF family and were upregulated. A total of 13 genes related to the photosystem, 3 genes related photo-regulation, and 9 TFs were identified under wounding stress. In conclusion, several novel putative genes were found in Mongolian pine by PacBio ISO seq. Meanwhile, we also identified various genes that were resistant to S. noctilio venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.

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.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Molecular Insight into Genetic Structure and Diversity of Putative Hybrid Swarms of Pinus sylvestris × P. mugo in Slovakia

The genetic structures of the four putative hybrid swarms of Pinus sylvestris × P. mugo in Slovakia were analyzed in terms of individual admixture proportions calculated via inter-primer binding site (iPBS) marker loci. This work aimed to reevaluate the hybrid swarms’ differentiation status as postulated in the previous studies at both population and genomic levels. The study confirmed intermediate divergence of each of the swarms examined. Based on 80-loci hybrid index scores, we have revealed the presence of introgressive and intermediate hybrids with frequencies corresponding to differentiation estimates. Surprisingly, irrespective of individual phenotypes, the most frequent intermediates were found in Sucha Hora (29.5%) and Obsivanka (28.6%) populations, which resemble rather pure P. mugo and were previously considered as bimodal hybrid zones with a negligible amount of hybrids. The remaining hybrid zone population Zuberec seems to be highly introgressed to P. sylvestris, while Tisovnica is clearly inclined to P. mugo. The results and different methodologies are discussed.

Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome. Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.

Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches

Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Transcript Dynamics in Wounded and Inoculated Scots Pine

Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.

Genetic Consequences of Hybridization in Relict Isolated Trees Pinus sylvestris and the Pinus mugo Complex

Scots pine (Pinus sylvestris L.) and the taxa from the P. mugo complex can hybridize in the contact zones and produce fertile hybrids. A unique example of an early Holocene relict population of P. sylvestris and P. uliginosa (a taxon from the P. mugo complex) growing on the tops of Jurassic sandstone rocks is located in Błędne Skały (Sudetes). Phenotypically, there are trees resembling P. sylvestris, P. uliginosa and intermediate forms between them. We expected that some of P. sylvestris and/or P. uliginosa-like trees could be in fact cryptic hybrids resembling one of the parental phenotypes. To address this question, we examined randomly sampled individuals, using a set of plastid (cpDNA), nuclear (nDNA) and mitochondrial (mtDNA) markers as well as biometric characteristics of needles and cones. The results were compared to the same measurements of allopatric reference populations of the P. sylvestris and the P. mugo complex (Pinus mugo s.s, P. uncinata and P. uliginosa). We detected cpDNA barcodes of the P. mugo complex in most individuals with the P. sylvestris phenotype, while we did not detect cpDNA diagnostic of P. sylvestris within P. uliginosa-like trees. These results indicate the presence of cryptic hybrids of the P. sylvestris phenotype. We found only three typical P. sylvestris individuals that were clustered with the species reference populations based on needle and cone characteristics. Most trees showed intermediate characteristics between P. sylvestris and P. uliginosa-like trees, indicating intensive and probably long-lasting hybridization of the taxa at this area and subsequent gene erosion of parental species.

Development of a single nucleotide polymorphism array for population genomic studies in four European pine species

Pines are some of the most ecologically and economically important tree species in the world, and many have enormous natural distributions or have been extensively planted. However, a lack of rapid genotyping capability is hampering progress in understanding the molecular basis of genetic variation in these species. Here, we deliver an efficient tool for genotyping thousands of single nucleotide polymorphism (SNP) markers across the genome that can be applied to genetic studies in pines. Polymorphisms from resequenced candidate genes and transcriptome sequences of P. sylvestris, P. mugo, P. uncinata, P. uliginosa and P. radiata were used to design a 49,829 SNP array (Axiom_PineGAP, Thermo Fisher). Over a third (34.68%) of the unigenes identified from the P. sylvestris transcriptome were represented on the array, which was used to screen samples of four pine species. The conversion rate for the array on all samples was 42% (N = 20,795 SNPs) and was similar for SNPs sourced from resequenced candidate gene and transcriptome sequences. The broad representation of gene ontology terms by unigenes containing converted SNPs reflected their coverage across the full transcriptome. Over a quarter of successfully converted SNPs were polymorphic among all species, and the data were successful in discriminating among the species and some individual populations. The SNP array provides a valuable new tool to advance genetic studies in these species and demonstrates the effectiveness of the technology for rapid genotyping in species with large and complex genomes.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

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.

Improved recovery of ancient DNA from subfossil wood - application to the world's oldest Late Glacial pine forest

Ancient DNA from historical and subfossil wood has a great potential to provide new insights into the history of tree populations. However, its extraction and analysis have not become routine, mainly because contamination of the wood with modern plant material can complicate the verification of genetic information. Here, we used sapwood tissue from 22 subfossil pines that were growing c. 13 000 yr bp in Zurich, Switzerland. We developed and evaluated protocols to eliminate surface contamination, and we tested ancient DNA authenticity based on plastid DNA metabarcoding and the assessment of post-mortem DNA damage. A novel approach using laser irradiation coupled with bleaching and surface removal was most efficient in eliminating contaminating DNA. DNA metabarcoding confirmed which ancient DNA samples repeatedly amplified pine DNA and were free of exogenous plant taxa. Pine DNA sequences of these samples showed a high degree of cytosine to thymine mismatches, typical of post-mortem damage. Stringent decontamination of wood surfaces combined with DNA metabarcoding and assessment of post-mortem DNA damage allowed us to authenticate ancient DNA retrieved from the oldest Late Glacial pine forest. These techniques can be applied to any subfossil wood and are likely to improve the accessibility of relict wood for genome-scale ancient DNA studies.

Chemical characterization and identification of Pinaceae pollen by infrared microspectroscopy

FTIR microspectroscopy, in combination with spectral averaging procedure, enables precise analysis of pollen grains for chemical characterization and identification studies of fresh and fossilised pollen in botany, ecology and palaeosciences. Infrared microspectroscopy (µFTIR) of Pinaceae pollen can provide valuable information on plant phenology, ecophysiology and paleoecology, but measurements are challenging, resulting in unreproducible spectra. The comparative analysis of µFTIR spectra belonging to morphologically different Pinaceae pollen, namely bisaccate Pinus and monosaccate Tsuga pollen, was conducted. The study shows that the main cause of spectral variability is non-radial symmetry of bisaccate pollen grains, while additional variation is caused by Mie scattering. Averaging over relatively small number of single pollen grain spectra (approx. 5-10) results with reproducible data on pollen chemical composition. The practical applicability of the µFTIR spectral averaging method has been demonstrated by the partial least-squares regression-based differentiation of the two closely related Pinus species with morphologically indistinguishable pollen: Pinus mugo (mountain pine) and Pinus sylvestris (Scots pine). The study has demonstrated that the µFTIR approach can be used for identification, differentiation and chemical characterization of pollen with complex morphology. The methodology enables analysis of fresh pollen, as well as fossil pollen from sediment core samples, and can be used in botany, ecology and paleoecology for study of biotic and abiotic effects on plants.

Early phenology and growth trait variation in closely related European pine species

Closely related taxa occupying different environments are valuable systems for studying evolution. In this study, we examined differences in early phenology (bud set, bud burst) and early growth in a common garden trial of closely related pine species: Pinus sylvestris, P. mugo, and P. uncinata. Seeds for the trial were sourced from populations across the ranges of each species in Europe. Over first 4 years of development, clear differences were observed between species, while the most significant intraspecific differentiation was observed among plants from P. sylvestris populations from continental European locations. Trait differences within P. sylvestris were highly correlated with altitude and latitude of the site of origin. Meanwhile, P. mugo populations from the Carpathians had the earliest bud set and bud flush compared to other populations of the species. Overall, populations from the P. mugo complex from heterogeneous mountain environments and P. sylvestris from the Scottish Highlands showed the highest within-population variation for the focal traits. Although the three species have been shown to be genetically highly similar, this study reveals large differences in key adaptive traits both among and within species.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping

Background

Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola.

Results

Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed.

Conclusions

We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Electronic supplementary material

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