SNP discovery with EST and NextGen sequencing in switchgrass (Panicum virgatum L.)

Transcriptome divergence during leaf development in two contrasting switchgrass (Panicum virgatum L.) cultivars

The genetics and responses to biotic stressors of tetraploid switchgrass (Panicum virgatum L.) lowland cultivar 'Kanlow' and upland cultivar Summer are distinct and can be exploited for trait improvement. In general, there is a paucity of data on the basal differences in transcription across tissue developmental times for switchgrass cultivars. Here, the changes in basal and temporal expression of genes related to leaf functions were evaluated for greenhouse grown 'Kanlow', and 'Summer' plants. Three biological replicates of the 4th leaf pooled from 15 plants per replicate were harvested at regular intervals beginning from leaf emergence through senescence. Increases and decreases in leaf chlorophyll and N content were similar for both cultivars. Likewise, multidimensional scaling (MDS) analysis indicated both cultivar-independent and cultivar-specific gene expression. Cultivar-independent genes and gene-networks included those associated with leaf function, such as growth/senescence, carbon/nitrogen assimilation, photosynthesis, chlorophyll biosynthesis, and chlorophyll degradation. However, many genes encoding nucleotide-binding leucine rich repeat (NB-LRRs) proteins and wall-bound kinases associated with detecting and responding to environmental signals were differentially expressed. Several of these belonged to unique cultivar-specific gene co-expression networks. Analysis of genomic resequencing data provided several examples of NB-LRRs genes that were not expressed and/or apparently absent in the genomes of Summer plants. It is plausible that cultivar (ecotype)-specific genes and gene-networks could be one of the drivers for the documented differences in responses to leaf-borne pathogens between these two cultivars. Incorporating broad resistance to plant pathogens in elite switchgrass germplasm could improve sustainability of biomass production under low-input conditions.

Enriching Genomic Resources and Marker Development from Transcript Sequences of Jatropha curcas for Microgravity Studies

Jatropha (Jatropha curcas L.) is an economically important species with a great potential for biodiesel production. To enrich the jatropha genomic databases and resources for microgravity studies, we sequenced and annotated the transcriptome of jatropha and developed SSR and SNP markers from the transcriptome sequences. In total 1,714,433 raw reads with an average length of 441.2 nucleotides were generated. De novo assembling and clustering resulted in 115,611 uniquely assembled sequences (UASs) including 21,418 full-length cDNAs and 23,264 new jatropha transcript sequences. The whole set of UASs were fully annotated, out of which 59,903 (51.81%) were assigned with gene ontology (GO) term, 12,584 (10.88%) had orthologs in Eukaryotic Orthologous Groups (KOG), and 8,822 (7.63%) were mapped to 317 pathways in six different categories in Kyoto Encyclopedia of Genes and Genome (KEGG) database, and it contained 3,588 putative transcription factors. From the UASs, 9,798 SSRs were discovered with AG/CT as the most frequent (45.8%) SSR motif type. Further 38,693 SNPs were detected and 7,584 remained after filtering. This UAS set has enriched the current jatropha genomic databases and provided a large number of genetic markers, which can facilitate jatropha genetic improvement and many other genetic and biological studies.

In situ embryo rescue for generation of wide intra- and interspecific hybrids of Panicum virgatum L

Wide crosses have been used for decades as a method for transferring novel genetic material and traits in plant breeding. Historically, many products of wide crosses require tedious and inefficient surgical embryo rescue prior to embryo abortion to recover single plantlets. We have utilized transgenic switchgrass (Panicum virgatum L. cv Alamo) as a pollen donor in conjunction with antibiotic or herbicide selection for recovery of intra-and interspecific F1 crosses by using developing ovules from the female parent and selecting for embryogenic cultures derived from the in situ immature embryo. Using this approach, several intravarietial crosses were generated between transgenic Alamo and the switchgrass varieties Kanlow, Blackwell and Cave-in-Rock as well as an interspecific cross with Atlantic coastal panicgrass. This procedure selected F1 embryogenic callus produced from the developing embryo contained within isolated immature ovules. Several clonal plants were successfully regenerated from each cross. Southern blot, PCR, phenotypic analyses and genomic analysis confirmed F1 hybrids. Using genotyping-by-sequencing shows the hybridization of the recovered plants by determining the ratio of transgressive markers to total compared markers between parents and their potential offspring. The ratio of transgressive markers to total compared markers was significantly lower between parents and their predicted offspring than between parents and offspring unrelated to them. This approach provides the possibility to move useful transgenes into varieties that are recalcitrant to direct transformation which can be optionally segregated thus useful to create new hybrids, as well as recovery of wide crosses that are either difficult or impossible using traditional techniques.

Genetic variation, population structure and linkage disequilibrium in Switchgrass with ISSR, SCoT and EST-SSR markers

Background

To evaluate genetic variation, population structure, and the extent of linkage disequilibrium (LD), 134 switchgrass (Panicum virgatum L.) samples were analyzed with 51 markers, including 16 ISSRs, 20 SCoTs, and 15 EST-SSRs.

Results

In this study, a high level of genetic variation was observed in the switchgrass samples and they had an average Nei’s gene diversity index (H) of 0.311. A total of 793 bands were obtained, of which 708 (89.28 %) were polymorphic. Using a parameter marker index (MI), the efficiency of the three types of markers (ISSR, SCoT, and EST-SSR) in the study were compared and we found that SCoT had a higher marker efficiency than the other two markers. The 134 switchgrass samples could be divided into two sub-populations based on STRUCTURE, UPGMA clustering, and principal coordinate analyses (PCA), and upland and lowland ecotypes could be separated by UPGMA clustering and PCA analyses. Linkage disequilibrium analysis revealed an average r² of 0.035 across all 51 markers, indicating a trend of higher LD in sub-population 2 than that in sub-population 1 (P < 0.01).

Conclusions

The population structure revealed in this study will guide the design of future association studies using these switchgrass samples.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

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

SWEEP: A Tool for Filtering High-Quality SNPs in Polyploid Crops

High-throughput next-generation sequence-based genotyping and single nucleotide polymorphism (SNP) detection opens the door for emerging genomics-based breeding strategies such as genome-wide association analysis and genomic selection. In polyploids, SNP detection is confounded by a highly similar homeologous sequence where a polymorphism between subgenomes must be differentiated from a SNP. We have developed and implemented a novel tool called SWEEP: Sliding Window Extraction of Explicit Polymorphisms. SWEEP uses subgenome polymorphism haplotypes as contrast to identify true SNPs between genotypes. The tool is a single command script that calls a series of modules based on user-defined options and takes sorted/indexed bam files or vcf files as input. Filtering options are highly flexible and include filtering based on sequence depth, alternate allele ratio, and SNP quality on top of the SWEEP filtering procedure. Using real and simulated data we show that SWEEP outperforms current SNP filtering methods for polyploids. SWEEP can be used for high-quality SNP discovery in polyploid crops.

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Understanding the relationship between genotype and phenotype is a major biological question and being able to predict phenotypes based on molecular genotypes is integral to molecular breeding. Whole-genome duplications have shaped the history of all flowering plants and present challenges to elucidating the relationship between genotype and phenotype, especially in neopolyploid species. Although single nucleotide polymorphisms (SNPs) have become popular tools for genetic mapping, discovery and application of SNPs in polyploids has been difficult. Here, we summarize common experimental approaches to SNP calling, highlighting recent polyploid successes. To examine the impact of software choice on these analyses, we called SNPs among five peanut genotypes using different alignment programs (BWA-mem and Bowtie 2) and variant callers (SAMtools, GATK, and Freebayes). Alignments produced by Bowtie 2 and BWA-mem and analyzed in SAMtools shared 24.5% concordant SNPs, and SAMtools, GATK, and Freebayes shared 1.4% concordant SNPs. A subsequent analysis of simulated Brassica napus chromosome 1A and 1C genotypes demonstrated that, of the three software programs, SAMtools performed with the highest sensitivity and specificity on Bowtie 2 alignments. These results, however, are likely to vary among species, and we therefore propose a series of best practices for SNP calling in polyploids.

Contrasting metabolism in perenniating structures of upland and lowland switchgrass plants late in the growing season

Background

Switchgrass (Panicum virgatum L.) is being developed as a bioenergy crop for many temperate regions of the world. One way to increase biomass yields is to move southern adapted lowland cultivars to more northern latitudes. However, many southerly adapted switchgrass germplasm can suffer significant winter kill in northerly climes.

Materials and Methods

Here, we have applied next-generation sequencing in combination with biochemical analyses to query the metabolism of crowns and rhizomes obtained from two contrasting switchgrass cultivars. Crowns and rhizomes from field-grown lowland (cv Kanlow) and upland (cv Summer) switchgrass cultivars were collected from three randomly selected post-flowering plants. Summer plants were senescing, whereas Kanlow plants were not at this harvest date.

Results

Principal component analysis (PCA) differentiated between both the Summer and Kanlow transcriptomes and metabolomes. Significant differences in transcript abundances were detected for 8,050 genes, including transcription factors such as WRKYs and those associated with phenylpropanoid biosynthesis. Gene-set enrichment analyses showed that a number of pathways were differentially up-regulated in the two populations. For both populations, protein levels and enzyme activities agreed well with transcript abundances for genes involved in the phenylpropanoid pathway that were up-regulated in Kanlow crowns and rhizomes. The combination of these datasets suggests that dormancy-related mechanisms had been triggered in the crowns and rhizomes of the Summer plants, whereas the crowns and rhizomes of Kanlow plants had yet to enter dormancy.

Conclusions

Delayed establishment of dormancy at more northerly latitudes could be one factor that reduces winter-survival in the high-yielding Kanlow plants. Understanding the cellular signatures that accompany the transition to dormancy can be used in the future to select plants with improved winter hardiness.

Transcriptome sequencing for high throughput SNP development and genetic mapping in Pea

BACKGROUND

Pea has a complex genome of 4.3 Gb for which only limited genomic resources are available to date. Although SNP markers are now highly valuable for research and modern breeding, only a few are described and used in pea for genetic diversity and linkage analysis.

RESULTS

We developed a large resource by cDNA sequencing of 8 genotypes representative of modern breeding material using the Roche 454 technology, combining both long reads (400 bp) and high coverage (3.8 million reads, reaching a total of 1,369 megabases). Sequencing data were assembled and generated a 68 K unigene set, from which 41 K were annotated from their best blast hit against the model species Medicago truncatula. Annotated contigs showed an even distribution along M. truncatula pseudochromosomes, suggesting a good representation of the pea genome. 10 K pea contigs were found to be polymorphic among the genetic material surveyed, corresponding to 35 K SNPs.We validated a subset of 1538 SNPs through the GoldenGate assay, proving their ability to structure a diversity panel of breeding germplasm. Among them, 1340 were genetically mapped and used to build a new consensus map comprising a total of 2070 markers. Based on blast analysis, we could establish 1252 bridges between our pea consensus map and the pseudochromosomes of M. truncatula, which provides new insight on synteny between the two species.

CONCLUSIONS

Our approach created significant new resources in pea, i.e. the most comprehensive genetic map to date tightly linked to the model species M. truncatula and a large SNP resource for both academic research and breeding.

Comprehensive analysis of expressed sequence tags from cultivated and wild radish (Raphanus spp.)

BACKGROUND

Radish (Raphanus sativus L., 2n = 2× = 18) is an economically important vegetable crop worldwide. A large collection of radish expressed sequence tags (ESTs) has been generated but remains largely uncharacterized.

RESULTS

In this study, approximately 315,000 ESTs derived from 22 Raphanus cDNA libraries from 18 different genotypes were analyzed, for the purpose of gene and marker discovery and to evaluate large-scale genome duplication and phylogenetic relationships among Raphanus spp. The ESTs were assembled into 85,083 unigenes, of which 90%, 65%, 89% and 89% had homologous sequences in the GenBank nr, SwissProt, TrEMBL and Arabidopsis protein databases, respectively. A total of 66,194 (78%) could be assigned at least one gene ontology (GO) term. Comparative analysis identified 5,595 gene families unique to radish that were significantly enriched with genes related to small molecule metabolism, as well as 12,899 specific to the Brassicaceae that were enriched with genes related to seed oil body biogenesis and responses to phytohormones. The analysis further indicated that the divergence of radish and Brassica rapa occurred approximately 8.9-14.9 million years ago (MYA), following a whole-genome duplication event (12.8-21.4 MYA) in their common ancestor. An additional whole-genome duplication event in radish occurred at 5.1-8.4 MYA, after its divergence from B. rapa. A total of 13,570 simple sequence repeats (SSRs) and 28,758 high-quality single nucleotide polymorphisms (SNPs) were also identified. Using a subset of SNPs, the phylogenetic relationships of eight different accessions of Raphanus was inferred.

CONCLUSION

Comprehensive analysis of radish ESTs provided new insights into radish genome evolution and the phylogenetic relationships of different radish accessions. Moreover, the radish EST sequences and the associated SSR and SNP markers described in this study represent a valuable resource for radish functional genomics studies and breeding.

Transcriptome analysis of heat stress response in switchgrass (Panicum virgatum L.)

BACKGROUND

Global warming predictions indicate that temperatures will increase by another 2-6°C by the end of this century. High temperature is a major abiotic stress limiting plant growth and productivity in many areas of the world. Switchgrass (Panicum virgatum L.) is a model herbaceous bioenergy crop, due to its rapid growth rate, reliable biomass yield, minimal requirements of water and nutrients, adaptability to grow on marginal lands and widespread distribution throughout North America. The effect of high temperature on switchgrass physiology, cell wall composition and biomass yields has been reported. However, there is void in the knowledge of the molecular responses to heat stress in switchgrass.

RESULTS

We conducted long-term heat stress treatment (38°/30°C, day/night, for 50 days) in the switchgrass cultivar Alamo. A significant decrease in the plant height and total biomass was evident in the heat stressed plants compared to controls. Total RNA from control and heat stress samples were used for transcriptome analysis with switchgrass Affymetrix genechips. Following normalization and pre-processing, 5365 probesets were identified as differentially expressed using a 2-fold cutoff. Of these, 2233 probesets (2000 switchgrass unigenes) were up-regulated, and 3132 probesets (2809 unigenes) were down-regulated. Differential expression of 42 randomly selected genes from this list was validated using RT-PCR. Rice orthologs were retrieved for 78.7% of the heat stress responsive switchgrass probesets. Gene ontology (GOs) enrichment analysis using AgriGO program showed that genes related to ATPase regulator, chaperone binding, and protein folding was significantly up-regulated. GOs associated with protein modification, transcription, phosphorus and nitrogen metabolic processes, were significantly down-regulated by heat stress.

CONCLUSIONS

Plausible connections were identified between the identified GOs, physiological responses and heat response phenotype observed in switchgrass plants. Comparative transcriptome analysis in response to heat stress among four monocots - switchgrass, rice, wheat and maize identified 16 common genes, most of which were associated with protein refolding processes. These core genes will be valuable biomarkers for identifying heat sensitive plant germplasm since they are responsive to both short duration as well as chronic heat stress treatments, and are also expressed in different plant growth stages and tissue types.

Towards uncovering the roles of switchgrass peroxidases in plant processes

Herbaceous perennial plants selected as potential biofuel feedstocks had been understudied at the genomic and functional genomic levels. Recent investments, primarily by the U.S. Department of Energy, have led to the development of a number of molecular resources for bioenergy grasses, such as the partially annotated genome for switchgrass (Panicum virgatum L.), and some related diploid species. In its current version, the switchgrass genome contains 65,878 gene models arising from the A and B genomes of this tetraploid grass. The availability of these gene sequences provides a framework to exploit transcriptomic data obtained from next-generation sequencing platforms to address questions of biological importance. One such question pertains to discovery of genes and proteins important for biotic and abiotic stress responses, and how these components might affect biomass quality and stress response in plants engineered for a specific end purpose. It can be expected that production of switchgrass on marginal lands will expose plants to diverse stresses, including herbivory by insects. Class III plant peroxidases have been implicated in many developmental responses such as lignification and in the adaptive responses of plants to insect feeding. Here, we have analyzed the class III peroxidases encoded by the switchgrass genome, and have mined available transcriptomic datasets to develop a first understanding of the expression profiles of the class III peroxidases in different plant tissues. Lastly, we have identified switchgrass peroxidases that appear to be orthologs of enzymes shown to play key roles in lignification and plant defense responses to hemipterans.