Western white pine SNP discovery and high-throughput genotyping for breeding and conservation applications

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

SNP Detection in Pinus pinaster Transcriptome and Association with Resistance to Pinewood Nematode

Pinewood nematode (PWN, Bursaphelenchus xylophilus) is the causal agent of pine wilt disease (PWD), which severely affects Pinus pinaster stands in southwestern Europe. Despite the high susceptibility of P. pinaster, individuals of selected half-sib families have shown genetic variability in survival after PWN inoculation, indicating that breeding for resistance can be a valuable strategy to control PWD. In this work, RNA-seq data from susceptible and resistant plants inoculated with PWN were used for SNP discovery and analysis. A total of 186,506 SNPs were identified, of which 31 were highly differentiated between resistant and susceptible plants, including SNPs in genes involved in cell wall lignification, a process previously linked to PWN resistance. Fifteen of these SNPs were selected for validation through Sanger sequencing and 14 were validated. To evaluate SNP-phenotype associations, 40 half-sib plants were genotyped for six validated SNPs. Associations with phenotype after PWN inoculation were found for two SNPs in two different genes (MEE12 and PCMP-E91), as well as two haplotypes of HIPP41, although significance was not maintained following Bonferroni correction. SNPs here detected may be useful for the development of molecular markers for PWD resistance and should be further investigated in future association studies.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

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.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

In-vitro anti-fungal assay and association analysis reveal a role for the Pinus monticola PR10 gene (PmPR10-3.1) in quantitative disease resistance to white pine blister rust

Pathogenesis-related (PR) proteins play important roles in plant defense response. However, functional investigation of PR10 genes is still limited and their physiological roles have not been conclusively characterized in biological processes of conifer trees. Here, we identified multiple novel members in the western white pine (Pinus monticola) PmPR10 family by bioinformatic mining available transcriptomic data. Phylogenetic analysis of protein sequences revealed four PR10 and two PR10-like clusters with a high synteny across different species of five-needle pines. Of 10 PmPR10 genes, PmPR10-3.1 was selected and expressed in Escherichia coli. The purified recombinant protein exhibited inhibitory effects on spore hyphal growth of fungal pathogens Cronartium ribicola, Phoma exigua, and Phoma argillacea by in-vitro anti-fungal analysis. Genetic variation analysis detected a total of 21 single nucleotide polymorphisms (SNPs) within PmPR10-3.1 in a collection of P. monticola seed families. A nonsynonymous SNP (t178g) showed significant association with relative levels of quantitative disease resistance (QDR), explaining about 8.7% of phenotypic variation as the peak value across all SNPs. Our results provide valuable insight into the genetic architecture underlying P. monticola QDR and imply that PmPR10-3.1 may function as an important component in conifer basal immunity for non-specific resistance to a wide spectrum of pathogens.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

A 3K Axiom SNP array from a transcriptome-wide SNP resource sheds new light on the genetic diversity and structure of the iconic subtropical conifer tree Araucaria angustifolia (Bert.) Kuntze

High-throughput SNP genotyping has become a precondition to move to higher precision and wider genome coverage genetic analysis of natural and breeding populations of non-model species. We developed a 44,318 annotated SNP catalog for Araucaria angustifolia, a grandiose subtropical conifer tree, one of the only two native Brazilian gymnosperms, critically endangered due to its valuable wood and seeds. Following transcriptome assembly and annotation, SNPs were discovered from RNA-seq and pooled RAD-seq data. From the SNP catalog, an Axiom® SNP array with 3,038 validated SNPs was developed and used to provide a comprehensive look at the genetic diversity and structure of 15 populations across the natural range of the species. RNA-seq was a far superior source of SNPs when compared to RAD-seq in terms of conversion rate to polymorphic markers on the array, likely due to the more efficient complexity reduction of the huge conifer genome. By matching microsatellite and SNP data on the same set of A. angustifolia individuals, we show that SNPs reflect more precisely the actual genome-wide patterns of genetic diversity and structure, challenging previous microsatellite-based assessments. Moreover, SNPs corroborated the known major north-south genetic cline, but allowed a more accurate attribution to regional versus among-population differentiation, indicating the potential to select ancestry-informative markers. The availability of a public, user-friendly 3K SNP array for A. angustifolia and a catalog of 44,318 SNPs predicted to provide ~29,000 informative SNPs across ~20,000 loci across the genome, will allow tackling still unsettled questions on its evolutionary history, toward a more comprehensive picture of the origin, past dynamics and future trend of the species' genetic resources. Additionally, but not less importantly, the SNP array described, unlocks the potential to adopt genomic prediction methods to accelerate the still very timid efforts of systematic tree breeding of A. angustifolia.

Characterization of Cronartium ribicola dsRNAs reveals novel members of the family Totiviridae and viral association with fungal virulence

Background

Mycoviruses were recently discovered in the white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.). Detection and characterization of their double stranded RNA (dsRNA) would facilitate understanding of pathogen virulence and disease pathogenesis in WPBR systems.

Methods

Full-length cDNAs were cloned from the dsRNAs purified from viral-infected C. ribicola, and their cDNA sequences were determined by DNA sequencing. Evolutionary relationships of the dsRNAs with related mycoviruses were determined by phylogenetic analysis. Dynamic distributions of the viral RNAs within samples of their fungal host C. ribicola were investigated by measurement of viral genome prevalence and viral gene expression.

Results

In this study we identified and characterized five novel dsRNAs from C. ribicola, designated as Cronartium ribicola totivirus 1–5 (CrTV1 to CrTV5). These dsRNA sequences encode capsid protein and RNA-dependent RNA polymerase with significant homologies to dsRNA viruses of the family Totiviridae. Phylogenetic analysis showed that the CrTVs were grouped into two distinct clades. CrTV2 through CrTV5 clustered within the genus Totivirus. CrTV1 along with a few un-assigned dsRNAs constituted a distinct phyletic clade that is genetically distant from presently known genera in the Totiviridae family, indicating that CrTV1 represents a novel genus in the Totiviridae family. The CrTVs were prevalent in fungal samples obtained from infected western white pine, whitebark pine, and limber pines. Viral RNAs were generally expressed at higher levels during in planta mycelium growth than in aeciospores and urediniospores. CrTV4 was significantly associated with C. ribicola virulent pathotype and specific C. ribicola host tree species, suggesting dsRNAs as potential tools for dissection of pathogenic mechanisms of C. ribicola and diagnosis of C. ribicola pathotypes.

Conclusion

Phylogenetic and expression analyses of viruses in the WPBR pathogen, C. ribicola, have enchanced our understanding of virus diversity in the family Totiviridae, and provided a potential strategy to utilize pathotype-associated mycoviruses to control fungal forest diseases.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

Core set construction and association analysis of Pinus massoniana from Guangdong province in southern China using SLAF-seq

Germplasm resource collection and utilization are important in forestry species breeding. High-through sequencing technologies have been playing increasing roles in forestry breeding. In this study, specific-locus amplified fragment sequencing (SLAF-seq) was employed to analyze 149 masson pine (Pinus massoniana) accessions collected from Guangdong in China. A large number of 471,660 SNPs in the total collection were identified from 599,164 polymorphic SLAF tags. Population structure analysis showed that 149 masson pines could not be obviously divided into subpopulations. Two core sets, containing 29 masson pine accessions for increasing resin and wood yield respectively, were obtained from the total collection. Phenotypic analyses of five traits showed abundant variations, 25 suggestive and 9 significant SNPs were associated with the resin-yielding capacity (RYC') and volume of wood (VW) using EMMAX and FaST-LMM; 22 suggestive and 11 significant SNPs were associated with RYC' and VW using mrMLM and FASTmrMLM. Moreover, a large number of associated SNPs were detected in trait HT, DBH, RW and RYC using mrMLM, FASTmrMLM, FASTmrEMMA and ISIS EM-BLASSO. The core germplasm sets would be a valuable resource for masson pine improvement and breeding. In addition, the associated SNP markers would be meaningful for masson pine resource selection.

Multiomics analysis reveals a distinct mechanism of oleaginousness in the emerging model alga Chromochloris zofingiensis

Chromochloris zofingiensis, featured due to its capability to simultaneously synthesize triacylglycerol (TAG) and astaxanthin, is emerging as a leading candidate alga for production uses. To better understand the oleaginous mechanism of this alga, we conducted a multiomics analysis by systematically integrating time-resolved transcriptomes, lipidomes and metabolomes in response to nitrogen deprivation. The data analysis unraveled the distinct mechanism of TAG accumulation, which involved coordinated stimulation of multiple biological processes including supply of energy and reductants, carbon reallocation from protein and starch, and 'pushing' and 'pulling' carbon to TAG synthesis. Unlike the model alga Chlamydomonas, de novo fatty acid synthesis in C. zofingiensis was promoted, together with enhanced turnover of both glycolipids and phospholipids, supporting the drastic need of acyls for TAG assembly. Moreover, genomewide analysis identified many key functional enzymes and transcription factors that had engineering potential for TAG modulation. Two genes encoding glycerol-3-phosphate acyltransferase (GPAT), the first committed enzyme for TAG assembly, were found in the C. zofingiensis genome; in vivo functional characterization revealed that extrachloroplastic GPAT instead of chloroplastic GPAT played a central role in TAG synthesis. These findings illuminate distinct oleaginousness mechanisms in C. zofingiensis and pave the way towards rational manipulation of this alga to becone an emerging model for trait improvements.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.

Profiling methyl jasmonate-responsive transcriptome for understanding induced systemic resistance in whitebark pine (Pinus albicaulis)

RNA-seq analysis on whitebark pine needles demonstrated that methyl jasmonate (MeJA)-triggered transcriptome re-programming substantially overlapped with defense responses against insects and fungal pathogens in Pinus species, increasing current knowledge regarding induced systemic resistance (ISR) to pathogens and pests in whitebark pine. Many whitebark pine populations are in steep decline due to high susceptibility to mountain pine beetle and the non-native white pine blister rust (WPBR). Resistance, including induced systemic resistance (ISR), is not well characterized in whitebark pine, narrowing the current options for increasing the success of restoration and breeding programs. Exogenous jasmonates are known to trigger ISR by activating the plant's immune system through regulation of gene expression to produce chemical defense compounds. This study reports profiles of whitebark pine needle transcriptomes, following methyl jasmonate (MeJA) treatment using RNA-seq. A MeJA-responsive transcriptome was de novo assembled and transcriptome profiling identified a set of differentially expressed genes (DEGs), revealing 1422 up- and 999 down-regulated transcripts with at least twofold change (FDR corrected p < 0.05) in needle tissues in response to MeJA application. GO analysis revealed that these DEGs have putative functions in plant defense signalling, transcription regulation, biosyntheses of secondary metabolites, and other biological processes. Lineage-specific expression of defense-related genes was characterized through comparison with MeJA signalling in model plants. In particular, MeJA-triggered transcriptome re-programming substantially overlapped with defense responses against WPBR and insects in related Pinus species, suggesting that MeJA may be used to improve whitebark pine resistance to pathogens/pests. Our study provides new insights into molecular mechanisms and metabolic pathways involved in whitebark pine ISR. DEGs identified in this study can be used as candidates to facilitate identification of genomic variation contributing to host resistance and aid in breeding selection of elite genotypes with better adaptive fitness to environmental stressors in this endangered tree species.

Overexpression of the Prunus sogdiana NBS-LRR Subgroup Gene PsoRPM2 Promotes Resistance to the Root-Knot Nematode Meloidogyne incognita in Tobacco

Root-knot nematodes (RKNs), particularly Meloidogyne incognita, are the most devastating soil-borne pathogens that significantly affect the production of Prunus spp. fruit. RKN infection is difficult to control and consequently causes massive yield losses each year. However, several germplasms of wild Prunus spp. have been shown to display resistance to M. incognita. Consequently, both the isolation of novel plant resistance (R) genes and the characterization of their resistance mechanisms are important strategies for future disease control. R proteins require the co-chaperone protein HSP90-SGT1-RAR1 to achieve correct folding, maturation, and stabilization. Here, we used homologous cloning to isolate the R gene PsoRPM2 from the RKN-resistant species Prunus sogdiana. PsoRPM2 was found to encode a TIR-NB-LRR-type protein and react with significantly elevated PsoRPM2 expression levels in response to RKN infection. Transient expression assays indicated PsoRPM2 to be located in both the cytoplasm and the nucleus. Four transgenic tobacco lines that heterologously expressed PsoRPM2 showed enhanced resistance to M. incognita. Yeast two-hybrid analysis and bimolecular fluorescence complementation analysis demonstrated that both PsoRAR1 and PsoRPM2 interacted with PsoHSP90-1 and PsoSGT1, but not with one another. These results indicate that the observed PsoRPM2-mediated RKN resistance requires both PsoHSP90-1 and PsoSGT1, further suggesting that PsoRAR1 plays a functionally redundant role in the HSP90-SGT1-RAR1 co-chaperone.

Saturated genic SNP mapping identified functional candidates and selection tools for the Pinus monticola Cr2 locus controlling resistance to white pine blister rust

Molecular breeding incorporates efficient tools to increase rust resistance in five-needle pines. Susceptibility of native five-needle pines to white pine blister rust (WPBR), caused by the non-native invasive fungus Cronartium ribicola (J.C. Fisch.), has significantly reduced wild populations of these conifers in North America. Major resistance (R) genes against specific avirulent pathotypes have been found in several five-needle pine species. In this study, we screened genic SNP markers by comparative transcriptome and genetic association analyses and constructed saturated linkage maps for the western white pine (Pinus monticola) R locus (Cr2). Phenotypic segregation was measured by a hypersensitive reaction (HR)-like response on the needles and disease symptoms of cankered stems post inoculation by the C. ribicola avcr2 race. SNP genotypes were determined by HRM- and TaqMan-based SNP genotyping. Saturated maps of the Cr2-linkage group (LG) were constructed in three seed families using a total of 34 SNP markers within 21 unique genes. Cr2 was consistently flanked by contig_2142 (encoding a ruvb-like protein) and contig_3772 (encoding a delta-fatty acid desaturase) across the three seed families. Cr2 was anchored to the Pinus consensus LG-1, which differs from LGs where other R loci of Pinus species were mapped. GO annotation identified a set of NBS-LRR and other resistance-related genes as R candidates in the Cr2 region. Association of one nonsynonymous SNP locus of an NBS-LRR gene with Cr2-mediated phenotypes provides a valuable tool for marker-assisted selection (MAS), which will shorten the breeding cycle of resistance screening and aid in the restoration of WPBR-disturbed forest ecosystems.

Discovery of Putative Herbicide Resistance Genes and Its Regulatory Network in Chickpea Using Transcriptome Sequencing

Background: Chickpea (Cicer arietinum L.) contributes 75% of total pulse production. Being cheaper than animal protein, makes it important in dietary requirement of developing countries. Weed not only competes with chickpea resulting into drastic yield reduction but also creates problem of harboring fungi, bacterial diseases and insect pests. Chemical approach having new herbicide discovery has constraint of limited lead molecule options, statutory regulations and environmental clearance. Through genetic approach, transgenic herbicide tolerant crop has given successful result but led to serious concern over ecological safety thus non-transgenic approach like marker assisted selection is desirable. Since large variability in tolerance limit of herbicide already exists in chickpea varieties, thus the genes offering herbicide tolerance can be introgressed in variety improvement programme. Transcriptome studies can discover such associated key genes with herbicide tolerance in chickpea. Results: This is first transcriptomic studies of chickpea or even any legume crop using two herbicide susceptible and tolerant genotypes exposed to imidazoline (Imazethapyr). Approximately 90 million paired-end reads generated from four samples were processed and assembled into 30,803 contigs using reference based assembly. We report 6,310 differentially expressed genes (DEGs), of which 3,037 were regulated by 980 miRNAs, 1,528 transcription factors associated with 897 DEGs, 47 Hub proteins, 3,540 putative Simple Sequence Repeat-Functional Domain Marker (SSR-FDM), 13,778 genic Single Nucleotide Polymorphism (SNP) putative markers and 1,174 Indels. Randomly selected 20 DEGs were validated using qPCR. Pathway analysis suggested that xenobiotic degradation related gene, glutathione S-transferase (GST) were only up-regulated in presence of herbicide. Down-regulation of DNA replication genes and up-regulation of abscisic acid pathway genes were observed. Study further reveals the role of cytochrome P450, xyloglucan endotransglucosylase/hydrolase, glutamate dehydrogenase, methyl crotonoyl carboxylase and of thaumatin-like genes in herbicide resistance. Conclusion: Reported DEGs can be used as genomic resource for future discovery of candidate genes associated with herbicide tolerance. Reported markers can be used for future association studies in order to develop marker assisted selection (MAS) for refinement. In endeavor of chickpea variety development programme, these findings can be of immense use in improving productivity of chickpea germplasm.

Genetic Diversity and Population Structure of Whitebark Pine (Pinus albicaulis Engelm.) in Western North America

Whitebark pine (WBP, Pinus albicaulis Engelm.) is an endangered conifer species due to heavy mortality from white pine blister rust (WPBR, caused by Cronartium ribicola) and mountain pine beetle (Dendroctonus ponderosae). Information about genetic diversity and population structure is of fundamental importance for its conservation and restoration. However, current knowledge on the genetic constitution and genomic variation is still limited for WBP. In this study, an integrated genomics approach was applied to characterize seed collections from WBP breeding programs in western North America. RNA-seq analysis was used for de novo assembly of the WBP needle transcriptome, which contains 97,447 protein-coding transcripts. Within the transcriptome, single nucleotide polymorphisms (SNPs) were discovered, and more than 22,000 of them were non-synonymous SNPs (ns-SNPs). Following the annotation of genes with ns-SNPs, 216 ns-SNPs within candidate genes with putative functions in disease resistance and plant defense were selected to design SNP arrays for high-throughput genotyping. Among these SNP loci, 71 were highly polymorphic, with sufficient variation to identify a unique genotype for each of the 371 individuals originating from British Columbia (Canada), Oregon and Washington (USA). A clear genetic differentiation was evident among seed families. Analyses of genetic spatial patterns revealed varying degrees of diversity and the existence of several genetic subgroups in the WBP breeding populations. Genetic components were associated with geographic variables and phenotypic rating of WPBR disease severity across landscapes, which may facilitate further identification of WBP genotypes and gene alleles contributing to local adaptation and quantitative resistance to WPBR. The WBP genomic resources developed here provide an invaluable tool for further studies and for exploitation and utilization of the genetic diversity preserved within this endangered conifer and other five-needle pines.

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.

Characterization of Five Novel Mitoviruses in the White Pine Blister Rust Fungus Cronartium ribicola

The white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.) is an exotic invasive forest pathogen causing severe stem canker disease of native white pine trees (subgenus Strobus) in North America. The present study reports discovery of five novel mitoviruses in C. ribicola by deep RNA sequencing. The complete genome of each mitovirus was determined by rapid amplification of cDNA ends (RACE) and reverse transcriptase-polymerase chain reaction (RT-PCR). A single open reading frame (ORF) encoding a putative RNA-dependent RNA polymerase (RdRp) was detected in each of the viral genomes using mitochondrial genetic codes. Phylogenetic analysis indicated that the C. ribicola mitoviruses (CrMV1 to CrMV5) are new putative species of the genus Mitovirus. qRT-PCR and RNA-Seq analyses revealed that viral RNAs were significantly increased in fungal mycelia in cankered pine stems compared to expression during two different stages of spore development, suggesting that viral genome replication and transcription benefit from active growth of the host fungus. CrMVs were widespread with relatively high levels of minor allele frequency (MAF) in western North America. As the first report of mitoviruses in the Class Pucciniomycetes, this work allows further investigation of the dynamics of a viral community in the WPBR pathosystem, including potential impacts that may affect pathogenicity and virulence of the host fungus.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors

Background

The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq.

Results

cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses.

Conclusion

This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Electronic supplementary material

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