Post-glacial evolution of Panicum virgatum: centers of diversity and gene pools revealed by SSR markers and cpDNA sequences

Genomic prediction of switchgrass winter survivorship across diverse lowland populations

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

Foxtail mosaic virus-induced gene silencing (VIGS) in switchgrass (Panicum virgatum L.)

BACKGROUND

Although the genome for the allotetraploid bioenergy crop switchgrass (Panicum virgatum) has been established, limitations in mutant resources have hampered in planta gene function studies toward crop optimization. Virus-induced gene silencing (VIGS) is a versatile technique for transient genetic studies. Here we report the implementation of foxtail mosaic virus (FoMV)-mediated gene silencing in switchgrass in above- and below-ground tissues and at different developmental stages.

RESULTS

The study demonstrated that leaf rub-inoculation is a suitable method for systemic gene silencing in switchgrass. For all three visual marker genes, Magnesium chelatase subunit D (ChlD) and I (ChlI) as well as phytoene desaturase (PDS), phenotypic changes were observed in leaves, albeit at different intensities. Gene silencing efficiency was verified by RT-PCR for all tested genes. Notably, systemic gene silencing was also observed in roots, although silencing efficiency was stronger in leaves (~ 63-94%) as compared to roots (~ 48-78%). Plants at a later developmental stage were moderately less amenable to VIGS than younger plants, but also less perturbed by the viral infection.

CONCLUSIONS

Using FoMV-mediated VIGS could be achieved in switchgrass leaves and roots, providing an alternative approach for studying gene functions and physiological traits in this important bioenergy crop.

Genotyping-by-Sequencing and QTL Mapping of Biomass Yield in Two Switchgrass F1 Populations (Lowland x Coastal and Coastal x Upland)

The prevalence of genetic diversity in switchgrass germplasm can be exploited to capture favorable alleles that increase its range of adaptation and biomass yield. The objectives of the study were to analyze the extent of polymorphism and patterns of segregation distortion in two F1 populations and use the linkage maps to locate QTL for biomass yield. We conducted genotyping-by-sequencing on two populations derived from crosses between the allotetraploid lowland genotype AP13 (a selection from "Alamo") and coastal genotype B6 (a selection from PI 422001) with 285 progeny (AB population) and between B6 and the allotetraploid upland VS16 (a selection from "Summer") with 227 progeny (BV population). As predictable from the Euclidean distance between the parents, a higher number of raw variants was discovered in the coastal × upland BV cross (6 M) compared to the lowland × coastal AB cross (2.5 M). The final number of mapped markers was 3,107 on the BV map and 2,410 on the AB map. More segregation distortion of alleles was seen in the AB population, with 75% distorted loci compared to 11% distorted loci in the BV population. The distortion in the AB population was seen across all chromosomes in both the AP13 and B6 maps and likely resulted from zygotic or post-zygotic selection for increased levels of heterozygosity. Our results suggest lower genetic compatibility between the lowland AP13 and the coastal B6 ecotype than between B6 and the upland ecotype VS16. Four biomass QTLs were mapped in the AB population (LG 2N, 6K, 6N, and 8N) and six QTLs in the BV population [LG 1N (2), 8N (2), 9K, and 9N]. The QTL, with the largest and most consistent effect across years, explaining between 8.4 and 11.5% of the variation, was identified on 6N in the AP13 map. The cumulative effect of all the QTLs explained a sizeable portion of the phenotypic variation in both AB and BV populations and the markers associated with them may potentially be used for the marker-assisted improvement of biomass yield. Since switchgrass improvement is based on increasing favorable allele frequencies through recurrent selection, the transmission bias within individuals and loci needs to be considered as this may affect the genetic gain if the favorable alleles are distorted.

Climate and stomatal traits drive covariation in nighttime stomatal conductance and daytime gas exchange rates in a widespread C4 grass

Nighttime stomatal conductance (g_sn ) varies among plant functional types and species, but factors shaping the evolution of g_sn remain unclear. Examinations of intraspecific variation in g_sn as a function of climate and co-varying leaf traits may provide new insight into the evolution of g_sn and its adaptive significance. We grew 11 genotypes of Panicum virgatum (switchgrass) representing differing home-climates in a common garden experiment and measured nighttime and daytime leaf gas exchange, as well as stomatal density (SD) and size during early-, mid-, and late-summer. We used piecewise structural equation modelling to determine direct and indirect relationships between home-climate, gas exchange, and stomatal traits. We found no direct relationship between home-climate and g_sn . However, genotypes from hotter climates possessed higher SD, which resulted in higher g_sn . Across genotypes, higher g_sn was associated with higher daytime stomatal conductance and net photosynthesis. Our results indicate that higher g_sn may arise in genotypes from hotter climates via increased SD. High SD may provide benefits to genotypes from hotter climates through enhanced daytime transpirational cooling or by permitting maximal gas exchange when conditions are suitable. These results highlight the role of climate and trait coordination in shaping genetic differentiation in g_sn .

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass

Long-term climate change and periodic environmental extremes threaten food and fuel security¹ and global crop productivity^2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience⁵, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation⁶-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

Linkage mapping evidence for a syntenic QTL associated with flowering time in perennial C4 rhizomatous grasses Miscanthus and switchgrass

Flowering in perennial species is directed via complex signalling pathways that adjust to developmental regulations and environmental cues. Synchronized flowering in certain environments is a prerequisite to commercial seed production, and so the elucidation of the genetic architecture of flowering time in Miscanthus and switchgrass could aid breeding in these underdeveloped species. In this context, we assessed a mapping population in Miscanthus and two ecologically diverse switchgrass mapping populations over 3 years from planting. Multiple flowering time quantitative trait loci (QTL) were identified in both species. Remarkably, the most significant Miscanthus and switchgrass QTL proved to be syntenic, located on linkage groups 4 and 2, with logarithm of odds scores of 17.05 and 21.8 respectively. These QTL regions contained three flowering time transcription factors: Squamosa Promoter-binding protein-Like, MADS-box SEPELLATA2 and gibberellin-responsive bHLH137. The former is emerging as a key component of the age-related flowering time pathway.

Spatial, climate and ploidy factors drive genomic diversity and resilience in the widespread grass Themeda triandra

Global climate change poses a significant threat to natural communities around the world, with many plant species showing signs of climate stress. Grassland ecosystems are not an exception, with climate change compounding contemporary pressures such as habitat loss and fragmentation. In this study, we assess the climate resilience of Themeda triandra, a foundational species and the most widespread plant in Australia, by assessing the relative contributions of spatial, environmental and ploidy factors to contemporary genomic variation. Reduced-representation genome sequencing on 472 samples from 52 locations was used to test how the distribution of genomic variation, including ploidy polymorphism, supports adaptation to hotter and drier climates. We explicitly quantified isolation by distance (IBD) and isolation by environment (IBE) and predicted genomic vulnerability of populations to future climates based on expected deviation from current genomic composition. We found that a majority (54%) of genomic variation could be attributed to IBD, while an additional 22% (27% when including ploidy information) could be explained by two temperature and two precipitation climate variables demonstrating IBE. Ploidy polymorphisms were common within populations (31/52 populations), indicating that ploidy mixing is characteristic of T. triandra populations. Genomic vulnerabilities were found to be heterogeneously distributed throughout the landscape, and our analysis suggested that ploidy polymorphism, along with other factors linked to polyploidy, reduced vulnerability to future climates by 60% (0.25-0.10). Our data suggests that polyploidy may facilitate adaptation to hotter climates and highlight the importance of incorporating ploidy in adaptive management strategies to promote the resilience of this and other foundation species.

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.

Complex interactions between day length and diurnal patterns of gene expression drive photoperiodic responses in a perennial C4 grass

Photoperiod is a key environmental cue affecting flowering and biomass traits in plants. Key components of the photoperiodic flowering pathway have been identified in many species, but surprisingly few studies have globally examined the diurnal rhythm of gene expression with changes in day length. Using a cost-effective 3'-Tag RNA sequencing strategy, we characterize 9,010 photoperiod responsive genes with strict statistical testing across a diurnal time series in the C₄ perennial grass, Panicum hallii. We show that the vast majority of photoperiod responses are driven by complex interactions between day length and sampling periods. A fine-scale contrast analysis at each sampling time revealed a detailed picture of the temporal reprogramming of cis-regulatory elements and biological processes under short- and long-day conditions. Phase shift analysis reveals quantitative variation among genes with photoperiod-dependent diurnal patterns. In addition, we identify three photoperiod enriched transcription factor families with key genes involved in photoperiod flowering regulatory networks. Finally, coexpression networks analysis of GIGANTEA homolog predicted 1,668 potential coincidence partners, including five well-known GI-interacting proteins. Our results not only provide a resource for understanding the mechanisms of photoperiod regulation in perennial grasses but also lay a foundation to increase biomass yield in biofuel crops.

Genomic Prediction for Winter Survival of Lowland Switchgrass in the Northern USA

The lowland ecotype of switchgrass has generated considerable interest because of its higher biomass yield and late flowering characteristics compared to the upland ecotype. However, lowland ecotypes planted in northern latitudes exhibit very low winter survival. Implementation of genomic selection could potentially enhance switchgrass breeding for winter survival by reducing generation time while eliminating the dependence on weather. The objectives of this study were to assess the potential of genomic selection for winter survival in lowland switchgrass by combining multiple populations in the training set and applying the selected model in two independent testing datasets for validation. Marker data were generated using exome capture sequencing. Validation was conducted using (1) indirect indicators of winter adaptation based on geographic and climatic variables of accessions from different source locations and (2) winter survival estimates of the phenotype. The prediction accuracies were significantly higher when the training dataset comprising all populations was used in fivefold cross validation but its application was not useful in the independent validation dataset. Nevertheless, modeling for population heterogeneity improved the prediction accuracy to some extent but the genetic relationship between the training and validation populations was found to be more influential. The predicted winter survival of lowland switchgrass indicated latitudinal and longitudinal variability, with the northeast USA the region for most cold tolerant lowland populations. Our results suggested that GS could provide valuable opportunities for improving winter survival and accelerate the lowland switchgrass breeding programs toward the development of cold tolerant cultivars suitable for northern latitudes.

Plant compartment and genetic variation drive microbiome composition in switchgrass roots

Switchgrass (Panicum virgatum) is a promising biofuel crop native to the United States with genotypes that are adapted to a wide range of distinct ecosystems. Various plants have been shown to undergo symbioses with plant growth-promoting bacteria and fungi, however, plant-associated microbial communities of switchgrass have not been extensively studied to date. We present 16S ribosomal RNA gene and internal transcribed spacer (ITS) data of rhizosphere and root endosphere compartments of four switchgrass genotypes to test the hypothesis that host selection of its root microbiota prevails after transfer to non-native soil. We show that differences in bacterial, archaeal and fungal community composition and diversity are strongly driven by plant compartment and switchgrass genotypes and ecotypes. Plant-associated microbiota show an enrichment in Alphaproteobacteria and Actinobacteria as well as Sordariales and Pleosporales compared with the surrounding soil. Root associated compartments display low-complexity communities dominated and enriched in Actinobacteria, in particular Streptomyces, in the lowland genotypes, and in Alphaproteobacteria, specifically Sphingobium, in the upland genotypes. Our comprehensive root analysis serves as a snapshot of host-specific bacterial and fungal associations of switchgrass in the field and confirms that host-selected microbiomes persist after transfer to non-native soil.

Quantitative Trait Loci for Freezing Tolerance in a Lowland x Upland Switchgrass Population

Low-temperature related abiotic stress is an important factor affecting winter survival in lowland switchgrass when grown in northern latitudes in the United States. A better understanding of the genetic architecture of freezing tolerance in switchgrass will aid the development of lowland switchgrass cultivars with improved winter survival. The objectives of this study were to conduct a freezing tolerance assessment, generate a genetic map using single nucleotide polymorphism (SNP) markers, and identify QTL (quantitative trait loci) associated with freezing tolerance in a lowland × upland switchgrass population. A pseudo-F2 mapping population was generated from an initial cross between the lowland population Ellsworth and the upland cultivar Summer. The segregating progenies were screened for freezing tolerance in a controlled-environment facility. Two clonal replicates of each genotype were tested at six different treatment temperatures ranging from -15 to -5°C at an interval of 2°C for two time periods. Tiller emergence (days) and tiller number were recorded following the recovery of each genotype with the hypothesis that upland genotype is the source for higher tiller number and early tiller emergence. Survivorship of the pseudo-F2 population ranged from 89% at -5°C to 5% at -15°C with an average LT50 of -9.7°C. Genotype had a significant effect on all traits except tiller number at -15°C. A linkage map was constructed from bi-allelic single nucleotide polymorphism markers generated using exome capture sequencing. The final map consisted of 1618 markers and 2626 cM, with an average inter-marker distance of 1.8 cM. Six significant QTL were identified, one each on chromosomes 1K, 5K, 5N, 6K, 6N, and 9K, for the following traits: tiller number, tiller emergence days and LT50. A comparative genomics study revealed important freezing tolerance genes/proteins, such as COR47, DREB2B, zinc finger-CCCH, WRKY, GIGANTEA, HSP70, and NRT2, among others that reside within the 1.5 LOD confidence interval of the identified QTL.

Breeding progress and preparedness for mass-scale deployment of perennial lignocellulosic biomass crops switchgrass, miscanthus, willow and poplar

Genetic improvement through breeding is one of the key approaches to increasing biomass supply. This paper documents the breeding progress to date for four perennial biomass crops (PBCs) that have high output-input energy ratios: namely Panicum virgatum (switchgrass), species of the genera Miscanthus (miscanthus), Salix (willow) and Populus (poplar). For each crop, we report on the size of germplasm collections, the efforts to date to phenotype and genotype, the diversity available for breeding and on the scale of breeding work as indicated by number of attempted crosses. We also report on the development of faster and more precise breeding using molecular breeding techniques. Poplar is the model tree for genetic studies and is furthest ahead in terms of biological knowledge and genetic resources. Linkage maps, transgenesis and genome editing methods are now being used in commercially focused poplar breeding. These are in development in switchgrass, miscanthus and willow generating large genetic and phenotypic data sets requiring concomitant efforts in informatics to create summaries that can be accessed and used by practical breeders. Cultivars of switchgrass and miscanthus can be seed-based synthetic populations, semihybrids or clones. Willow and poplar cultivars are commercially deployed as clones. At local and regional level, the most advanced cultivars in each crop are at technology readiness levels which could be scaled to planting rates of thousands of hectares per year in about 5 years with existing commercial developers. Investment in further development of better cultivars is subject to current market failure and the long breeding cycles. We conclude that sustained public investment in breeding plays a key role in delivering future mass-scale deployment of PBCs.

Candidate Variants for Additive and Interactive Effects on Bioenergy Traits in Switchgrass (Panicum virgatum L.) Identified by Genome-Wide Association Analyses

Switchgrass ( L.) is a promising herbaceous energy crop, but further gains in biomass yield and quality must be achieved to enable a viable bioenergy industry. Developing DNA markers can contribute to such progress, but depiction of genetic bases should be reliable, involving simple additive marker effects and also interactions with genetic backgrounds (e.g., ecotypes) or synergies with other markers. We analyzed plant height, C content, N content, and mineral concentration in a diverse panel consisting of 512 genotypes of upland and lowland ecotypes. We performed association analyses based on exome capture sequencing and tested 439,170 markers for marginal effects, 83,290 markers for marker × ecotype interactions, and up to 311,445 marker pairs for pairwise interactions. Analyses of pairwise interactions focused on subsets of marker pairs preselected on the basis of marginal marker effects, gene ontology annotation, and pairwise marker associations. Our tests identified 12 significant effects. Homology and gene expression information corroborated seven effects and indicated plausible causal pathways: flowering time and lignin synthesis for plant height; plant growth and senescence for C content and mineral concentration. Four pairwise interactions were detected, including three interactions preselected on the basis of pairwise marker correlations. Furthermore, a marker × ecotype interaction and a pairwise interaction were confirmed in an independent switchgrass panel. Our analyses identified reliable candidate variants for important bioenergy traits. Moreover, they exemplified the importance of interactive effects for depicting genetic bases and illustrated the usefulness of preselecting marker pairs for identifying pairwise marker interactions in association studies.

Natural variation in genes potentially involved in plant architecture and adaptation in switchgrass (Panicum virgatum L.)

Background

Advances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.

Results

A total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima’s test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.

Conclusions

Several non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1193-2) contains supplementary material, which is available to authorized users.

Extensive Genetic Diversity is Present within North American Switchgrass Germplasm

Switchgrass ( is a perennial native North American grass present in two ecotypes: upland, found primarily in the northern range of switchgrass habitats, and lowland, found largely in the southern reaches of switchgrass habitats. Previous studies focused on a diversity panel of primarily northern switchgrass, so to expand our knowledge of genetic diversity in a broader set of North American switchgrass, exome capture sequence data were generated for 632 additional, primarily lowland individuals. In total, over 37 million single nucleotide polymorphisms (SNPs) were identified and a set of 1.9 million high-confidence SNPs were obtained from 1169 individuals from 140 populations (67 upland, 65 lowland, 8 admixed) were used in downstream analyses of genetic diversity and population structure. Seven separate population groups were identified with moderate genetic differentiation [mean fixation index (Fst) estimate of 0.06] between the lowland and the upland populations. Ecotype-specific and population-specific SNPs were identified for use in germplasm evaluations. Relative to rice ( L.), maize ( L.), soybean [ (L.) Merr.], and Gaertn., analyses of nucleotide diversity revealed a high degree of genetic diversity (0.0135) across all individuals, consistent with the outcrossing mode of reproduction and the polyploidy of switchgrass. This study supports the hypothesis that repeated glaciation events, ploidy barriers, and restricted gene flow caused by flowering time differences have resulted in distinct gene pools across ecotypes and geographic regions. These data provide a resource to associate alleles with traits of interest for forage, restoration, and biofuel feedstock efforts in switchgrass.

Quantitative trait loci for cell wall composition traits measured using near-infrared spectroscopy in the model C4 perennial grass Panicum hallii

BACKGROUND

Biofuels derived from lignocellulosic plant material are an important component of current renewable energy strategies. Improvement efforts in biofuel feedstock crops have been primarily focused on increasing biomass yield with less consideration for tissue quality or composition. Four primary components found in the plant cell wall contribute to the overall quality of plant tissue and conversion characteristics, cellulose and hemicellulose polysaccharides are the primary targets for fuel conversion, while lignin and ash provide structure and defense. We explore the genetic architecture of tissue characteristics using a quantitative trait loci (QTL) mapping approach in Panicum hallii, a model lignocellulosic grass system. Diversity in the mapping population was generated by crossing xeric and mesic varietals, comparative to northern upland and southern lowland ecotypes in switchgrass. We use near-infrared spectroscopy with a primary analytical method to create a P. hallii specific calibration model to quickly quantify cell wall components.

RESULTS

Ash, lignin, glucan, and xylan comprise 68% of total dry biomass in P. hallii: comparable to other feedstocks. We identified 14 QTL and one epistatic interaction across these four cell wall traits and found almost half of the QTL to localize to a single linkage group.

CONCLUSIONS

Panicum hallii serves as the genomic model for its close relative and emerging biofuel crop, switchgrass (P. virgatum). We used high throughput phenotyping to map genomic regions that impact natural variation in leaf tissue composition. Understanding the genetic architecture of tissue traits in a tractable model grass system will lead to a better understanding of cell wall structure as well as provide genomic resources for bioenergy crop breeding programs.

Identification of massive molecular markers in Echinochloa phyllopogon using a restriction-site associated DNA approach

Echinochloa phyllopogon proliferation seriously threatens rice production worldwide. We combined a restriction-site associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers for E. phyllopogon. RAD tags were generated from the genomic DNA of two E. phyllopogon plants, and sequenced to produce 5197.7 Mb and 5242.9 Mb high quality sequences, respectively. The GC content of E. phyllopogon was 45.8%, which is high for monocots. In total, 4710 putative SSRs were identified in 4132 contigs, which permitted the design of PCR primers for E. phyllopogon. Most repeat motifs among the SSRs identified were dinucleotide (>82%), and most of these SSRs were four motif-repeats (>75%). The most frequent motif was AT, accounting for 36.3%-37.2%, followed by AG and AC. In total, 78 putative polymorphic SSR loci were found. A total of 49,179 SNPs were discovered between the two samples of E. phyllopogon, 67.1% of which were transversions and 32.9% were transitions. We used eight SSRs to study the genetic diversity of four E. phyllopogon populations collected from rice fields in China and all eight loci tested were polymorphic.

Selection Signatures in Four Lignin Genes from Switchgrass Populations Divergently Selected for In Vitro Dry Matter Digestibility

Switchgrass is undergoing development as a dedicated cellulosic bioenergy crop. Fermentation of lignocellulosic biomass to ethanol in a bioenergy system or to volatile fatty acids in a livestock production system is strongly and negatively influenced by lignification of cell walls. This study detects specific loci that exhibit selection signatures across switchgrass breeding populations that differ in in vitro dry matter digestibility (IVDMD), ethanol yield, and lignin concentration. Allele frequency changes in candidate genes were used to detect loci under selection. Out of the 183 polymorphisms identified in the four candidate genes, twenty-five loci in the intron regions and four loci in coding regions were found to display a selection signature. All loci in the coding regions are synonymous substitutions. Selection in both directions were observed on polymorphisms that appeared to be under selection. Genetic diversity and linkage disequilibrium within the candidate genes were low. The recurrent divergent selection caused excessive moderate allele frequencies in the cycle 3 reduced lignin population as compared to the base population. This study provides valuable insight on genetic changes occurring in short-term selection in the polyploid populations, and discovered potential markers for breeding switchgrass with improved biomass quality.

The Genetic Basis of Upland/Lowland Ecotype Divergence in Switchgrass (Panicum virgatum)

The evolution of locally adapted ecotypes is a common phenomenon that generates diversity within plant species. However, we know surprisingly little about the genetic mechanisms underlying the locally adapted traits involved in ecotype formation. The genetic architecture underlying locally adapted traits dictates how an organism will respond to environmental selection pressures, and has major implications for evolutionary ecology, conservation, and crop breeding. To understand the genetic architecture underlying the divergence of switchgrass (Panicum virgatum) ecotypes, we constructed a genetic mapping population through a four-way outbred cross between two northern upland and two southern lowland accessions. Trait segregation in this mapping population was largely consistent with multiple independent loci controlling the suite of traits that characterizes ecotype divergence. We assembled a joint linkage map using ddRADseq, and mapped quantitative trait loci (QTL) for traits that are divergent between ecotypes, including flowering time, plant size, physiological processes, and disease resistance. Overall, we found that most QTL had small to intermediate effects. While we identified colocalizing QTL for multiple traits, we did not find any large-effect QTL that clearly controlled multiple traits through pleiotropy or tight physical linkage. These results indicate that ecologically important traits in switchgrass have a complex genetic basis, and that similar loci may underlie divergence across the geographic range of the ecotypes.

Genotypic diversity effects on biomass production in native perennial bioenergy cropping systems

The perennial grass species that are being developed as biomass feedstock crops harbor extensive genotypic diversity, but the effects of this diversity on biomass production are not well understood. We investigated the effects of genotypic diversity in switchgrass (Panicum virgatum) and big bluestem (Andropogon gerardii) on perennial biomass cropping systems in two experiments conducted over 2008-2014 at a 5.4-ha fertile field site in northeastern Illinois, USA. We varied levels of switchgrass and big bluestem genotypic diversity using various local and nonlocal cultivars - under low or high species diversity, with or without nitrogen inputs - and quantified establishment, biomass yield, and biomass composition. In one experiment ('agronomic trial'), we compared three switchgrass cultivars in monoculture to a switchgrass cultivar mixture and three different species mixtures, with or without N fertilization. In another experiment ('diversity gradient'), we varied diversity levels in switchgrass and big bluestem (1, 2, 4, or 6 cultivars per plot), with one or two species per plot. In both experiments, cultivar mixtures produced yields equivalent to or greater than the best cultivars. In the agronomic trial, the three switchgrass mixture showed the highest production overall, though not significantly different than best cultivar monoculture. In the diversity gradient, genotypic mixtures had one-third higher biomass production than the average monoculture, and none of the monocultures were significantly higher yielding than the average mixture. Year-to-year variation in yields was lowest in the three-cultivar switchgrass mixtures and Cave-In-Rock (the southern Illinois cultivar) and also reduced in the mixture of switchgrass and big bluestem relative to the species monocultures. The effects of genotypic diversity on biomass composition were modest relative to the differences among species and genotypes. Our findings suggest that local genotypes can be included in biomass cropping systems without compromising yields and that genotypic mixtures could help provide high, stable yields of high-quality biomass feedstocks.

Relative Performance of Non-Local Cultivars and Local, Wild Populations of Switchgrass (Panicum virgatum) in Competition Experiments

The possibility of increased invasiveness in cultivated varieties of native perennial species is a question of interest in biofuel risk assessment. Competitive success is a key factor in the fitness and invasive potential of perennial plants, and thus the large-scale release of high-yielding biomass cultivars warrants empirical comparisons with local conspecifics in the presence of competitors. We evaluated the performance of non-local cultivars and local wild biotypes of the tallgrass species Panicum virgatum L. (switchgrass) in competition experiments during two growing seasons in Ohio and Iowa. At each location, we measured growth and reproductive traits (plant height, tiller number, flowering time, aboveground biomass, and seed production) of four non-locally sourced cultivars and two locally collected wild biotypes. Plants were grown in common garden experiments under three types of competition, referred to as none, moderate (with Schizachyrium scoparium), and high (with Bromus inermis). In both states, the two “lowland” cultivars grew taller, flowered later, and produced between 2x and 7.5x more biomass and between 3x and 34x more seeds per plant than local wild biotypes, while the other two cultivars were comparable to wild biotypes in these traits. Competition did not affect relative differences among biotypes, with the exception of shoot number, which was more similar among biotypes under high competition. Insights into functional differences between cultivars and wild biotypes are crucial for developing biomass crops while mitigating the potential for invasiveness. Here, two of the four cultivars generally performed better than wild biotypes, indicating that these biotypes may pose more of a risk in terms of their ability to establish vigorous feral populations in new regions outside of their area of origin. Our results support an ongoing assessment of switchgrass cultivars developed for large-scale planting for biofuels.

Accuracy of Genomic Prediction in Switchgrass (Panicum virgatum L.) Improved by Accounting for Linkage Disequilibrium

Switchgrass is a relatively high-yielding and environmentally sustainable biomass crop, but further genetic gains in biomass yield must be achieved to make it an economically viable bioenergy feedstock. Genomic selection (GS) is an attractive technology to generate rapid genetic gains in switchgrass, and meet the goals of a substantial displacement of petroleum use with biofuels in the near future. In this study, we empirically assessed prediction procedures for genomic selection in two different populations, consisting of 137 and 110 half-sib families of switchgrass, tested in two locations in the United States for three agronomic traits: dry matter yield, plant height, and heading date. Marker data were produced for the families’ parents by exome capture sequencing, generating up to 141,030 polymorphic markers with available genomic-location and annotation information. We evaluated prediction procedures that varied not only by learning schemes and prediction models, but also by the way the data were preprocessed to account for redundancy in marker information. More complex genomic prediction procedures were generally not significantly more accurate than the simplest procedure, likely due to limited population sizes. Nevertheless, a highly significant gain in prediction accuracy was achieved by transforming the marker data through a marker correlation matrix. Our results suggest that marker-data transformations and, more generally, the account of linkage disequilibrium among markers, offer valuable opportunities for improving prediction procedures in GS. Some of the achieved prediction accuracies should motivate implementation of GS in switchgrass breeding programs.

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

Switchgrass (Panicum virgatum L.) Genotypes Differ between Coastal Sites and Inland Road Corridors in the Northeastern US

Switchgrass (Panicum virgatum L.) is a North American grass that exhibits vast genetic diversity across its geographic range. In the Northeastern US, local switchgrass populations were restricted to a narrow coastal zone before European settlement, but current populations inhabit inland road verges raising questions about their origin and genetics. These questions are important because switchgrass lines with novel traits are being cultivated as a biofuel feedstock, and gene flow could impact the genetic integrity and distribution of local populations. This study was designed to determine if: 1) switchgrass plants collected in the Long Island Sound Coastal Lowland coastal Level IV ecoregion represented local populations, and 2) switchgrass plants collected from road verges in the adjacent inland regions were most closely related to local coastal populations or switchgrass from other geographic regions. The study used 18 microsatellite markers to infer the genetic relationships between 122 collected switchgrass plants and a reference dataset consisting of 28 cultivars representing ecotypes, ploidy levels, and lineages from North America. Results showed that 84% of 88 plants collected in the coastal plants were most closely aligned with the Lowland tetraploid genetic pool. Among this group, 61 coastal plants were similar to, but distinct from, all Lowland tetraploid cultivars in the reference dataset leading to the designation of a genetic sub-population called the Southern New England Lowland Tetraploids. In contrast, 67% of 34 plants collected in road verges in the inland ecoregions were most similar to two Upland octoploid cultivars; only 24% of roadside plants were Lowland tetraploid. These results suggest that cryptic, non-local genotypes exist in road verges and that gene flow from biofuels plantations could contribute to further changes in switchgrass population genetics in the Northeast.

Genome-wide distribution comparative and composition analysis of the SSRs in Poaceae

Background

The Poaceae family is of great importance to human beings since it comprises the cereal grasses which are the main sources for human food and animal feed. With the rapid growth of genomic data from Poaceae members, comparative genomics becomes a convinent method to study genetics of diffierent species. The SSRs (Simple Sequence Repeats) are widely used markers in the studies of Poaceae for their high abundance and stability.

Results

In this study, using the genomic sequences of 9 Poaceae species, we detected 11,993,943 SSR loci and developed 6,799,910 SSR primer pairs. The results show that SSRs are distributed on all the genomic elements in grass. Hexamer is the most frequent motif and AT/TA is the most frequent motif in dimer. The abundance of the SSRs has a positive linear relationship with the recombination rate. SSR sequences in the coding regions involve a higher GC content in the Poaceae than that in the other species. SSRs of 70-80 bp in length showed the highest AT/GC base ratio among all of these loci. The result shows the highest polymorphism rate belongs to the SSRs ranged from 30 bp to 40 bp. Using all the SSR primers of Japonica, nineteen universal primers were selected and located on the genome of the grass family. The information of SSR loci, the SSR primers and the tools of mining and analyzing SSR are provided in the PSSRD (Poaceae SSR Database, http://biodb.sdau.edu.cn/pssrd/).

Conclusions

Our study and the PSSRD database provide a foundation for the comparative study in the Poaceae and it will accelerate the study on markers application, gene mapping and molecular breeding.

Assessment of drought tolerance of 49 switchgrass (Panicum virgatum) genotypes using physiological and morphological parameters

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season C4 grass that is a target lignocellulosic biofuel species. In many regions, drought stress is one of the major limiting factors for switchgrass growth. The objective of this study was to evaluate the drought tolerance of 49 switchgrass genotypes. The relative drought stress tolerance was determined based on a set of parameters including plant height, leaf length, leaf width, leaf sheath length, leaf relative water content (RWC), electrolyte leakage (EL), photosynthetic rate (Pn), stomatal conductance (g s), transpiration rate (Tr), intercellular CO2 concentration (Ci), and water use efficiency (WUE).

RESULTS

SRAP marker analysis determined that the selected 49 switchgrass genotypes represent a diverse genetic pool of switchgrass germplasm. Principal component analysis (PCA) and drought stress indexes (DSI) of each physiological parameter showed significant differences in the drought stress tolerance among the 49 genotypes. Heatmap and PCA data revealed that physiological parameters are more sensitive than morphological parameters in distinguishing the control and drought treatments. Metabolite profiling data found that under drought stress, the five best drought-tolerant genotypes tended to have higher levels of abscisic acid (ABA), spermine, trehalose, and fructose in comparison to the five most drought-sensitive genotypes.

CONCLUSION

Based on PCA ranking value, the genotypes TEM-SEC, TEM-LoDorm, BN-13645-64, Alamo, BN-10860-61, BN-12323-69, TEM-SLC, T-2086, T-2100, T-2101, Caddo, and Blackwell-1 had relatively higher ranking values, indicating that they are more tolerant to drought. In contrast, the genotypes Grif Nebraska 28, Grenville-2, Central Iowa Germplasm, Cave-in-Rock, Dacotah, and Nebraska 28 were found to be relatively sensitive to drought stress. By analyzing physiological response parameters and different metabolic profiles, the methods utilized in this study identified drought-tolerant and drought-sensitive switchgrass genotypes. These results provide a foundation for future research directed at understanding the molecular mechanisms underlying switchgrass tolerance to drought.

The genetics of divergence and reproductive isolation between ecotypes of Panicum hallii

The process of plant speciation often involves the evolution of divergent ecotypes in response to differences in soil water availability between habitats. While the same set of traits is frequently associated with xeric/mesic ecotype divergence, it is unknown whether those traits evolve independently or if they evolve in tandem as a result of genetic colocalization either by pleiotropy or genetic linkage. The self-fertilizing C4 grass species Panicum hallii includes two major ecotypes found in xeric (var. hallii) or mesic (var. filipes) habitats. We constructed the first linkage map for P. hallii by genotyping a reduced representation genomic library of an F2 population derived from an intercross of var. hallii and filipes. We then evaluated the genetic architecture of divergence between these ecotypes through quantitative trait locus (QTL) mapping. Overall, we mapped QTLs for nine morphological traits that are involved in the divergence between the ecotypes. QTLs for five key ecotype-differentiating traits all colocalized to the same region of linkage group five. Leaf physiological traits were less divergent between ecotypes, but we still mapped five physiological QTLs. We also discovered a two-locus Dobzhansky-Muller hybrid incompatibility. Our study suggests that ecotype-differentiating traits may evolve in tandem as a result of genetic colocalization.

Accelerating the switchgrass (Panicum virgatum L.) breeding cycle using genomic selection approaches

Switchgrass (Panicum virgatum L.) is a perennial grass undergoing development as a biofuel feedstock. One of the most important factors hindering breeding efforts in this species is the need for accurate measurement of biomass yield on a per-hectare basis. Genomic selection on simple-to-measure traits that approximate biomass yield has the potential to significantly speed up the breeding cycle. Recent advances in switchgrass genomic and phenotypic resources are now making it possible to evaluate the potential of genomic selection of such traits. We leveraged these resources to study the ability of three widely-used genomic selection models to predict phenotypic values of morphological and biomass quality traits in an association panel consisting of predominantly northern adapted upland germplasm. High prediction accuracies were obtained for most of the traits, with standability having the highest ten-fold cross validation prediction accuracy (0.52). Moreover, the morphological traits generally had higher prediction accuracies than the biomass quality traits. Nevertheless, our results suggest that the quality of current genomic and phenotypic resources available for switchgrass is sufficiently high for genomic selection to significantly impact breeding efforts for biomass yield.

Distribution models for Panicum virgatum (Poaceae) reveal an expanded range in present and future climate regimes in the northeastern United States

PREMISE OF THE STUDY

Expanded area cultivated with the bioenergy crop Panicum virgatum (switchgrass) could alter the genetics of native populations through gene flow, so understanding current and future species distribution is a first step toward estimating ecological impacts. We surveyed switchgrass distribution in the northeastern United States and generated statistical models to address hypotheses about current distribution relative to historical records and responses to climate change.

METHODS

Surveys were conducted on 1600 km of road verges along environmental gradients. Switchgrass abundance became the training data for two multivariate generalized linear models that generated maps representing the probability of switchgrass in road verges. Models were evaluated and the superior model was used with variables from three climate change scenarios for 2050 and 2099.

KEY RESULTS

Switchgrass populations were found in 41% of roadside plots and up to 188 km from the coast. The environmental variables temperature, urban areas, and sandy soils were positively correlated with switchgrass presence, while elevation, soil pH, and distance to the coast were negatively correlated. The model without spatial autocorrelation performed better. Models and maps incorporating climate change predictions showed a sharp northward shift in suitable habitat.

CONCLUSIONS

Switchgrass populations in the northeastern United States occur on inland road verges, supporting the idea that species distribution has expanded relative to historical descriptions of a restricted coastal habitat. The optimal model showed that mean temperature, elevation, and urban development were important in switchgrass distribution today, and climate change will increase suitable habitat for future bioenergy production and wild populations.

Population genomic variation reveals roles of history, adaptation and ploidy in switchgrass

Geographic patterns of genetic variation are shaped by multiple evolutionary processes, including genetic drift, migration and natural selection. Switchgrass (Panicum virgatum L.) has strong genetic and adaptive differentiation despite life history characteristics that promote high levels of gene flow and can homogenize intraspecific differences, such as wind-pollination and self-incompatibility. To better understand how historical and contemporary factors shape variation in switchgrass, we use genotyping-by-sequencing to characterize switchgrass from across its range at 98 042 SNPs. Population structuring reflects biogeographic and ploidy differences within and between switchgrass ecotypes and indicates that biogeographic history, ploidy incompatibilities and differential adaptation each have important roles in shaping ecotypic differentiation in switchgrass. At one extreme, we determine that two Panicum taxa are not separate species but are actually conspecific, ecologically divergent types of switchgrass adapted to the extreme conditions of coastal sand dune habitats. Conversely, we identify natural hybrids among lowland and upland ecotypes and visualize their genome-wide patterns of admixture. Furthermore, we determine that genetic differentiation between primarily tetraploid and octoploid lineages is not caused solely by ploidy differences. Rather, genetic diversity in primarily octoploid lineages is consistent with a history of admixture. This suggests that polyploidy in switchgrass is promoted by admixture of diverged lineages, which may be important for maintaining genetic differentiation between switchgrass ecotypes where they are sympatric. These results provide new insights into the mechanisms shaping variation in widespread species and provide a foundation for dissecting the genetic basis of adaptation in switchgrass.

A footprint of past climate change on the diversity and population structure of Miscanthus sinensis

BACKGROUND AND AIMS

Miscanthus is a perennial C4 grass that is a leading potential feedstock crop for the emerging bioenergy industry in North America, Europe and China. However, only a single, sterile genotype of M. × giganteus (M×g), a nothospecies derived from diploid M. sinensis (Msi) and tetraploid M. sacchariflorus (Msa), is currently available to farmers for biomass production. To facilitate breeding of Miscanthus, this study characterized genetic diversity and population structure of Msi in its native range of East Asia.

METHODS

A total of 767 accessions were studied, including 617 Msi from most of its native range in China, Japan and South Korea, and 77 ornamental cultivars and 43 naturalized individuals from the USA. Accessions were evaluated with 21 207 restriction site-associated DNA sequencing single nucleotide polymorphism (SNP) markers, 424 GoldenGate SNPs and ten plastid microsatellite markers.

KEY RESULTS

Six genetic clusters of Msi from geographically distinct regions in Asia were identified. Genetic data indicated that (1) south-eastern China was the origin of Msi populations found in temperate eastern Asia, which is consistent with this area probably having been a refugium during the last glacial maximum (LGM); (2) Msi migrated directly from south-eastern China to Japan before migrating to the same latitudes in China and Korea, which is consistent with the known sequence of warming post-LGM; (3) ornamental Msi cultivars were derived from the southern Japan population, and US naturalized populations were derived from a sub-set of the ornamental cultivars; and (4) many ornamental cultivars previously described as Msi have hybrid ancestry from Msa and Msi, whereas US naturalized populations of Msi do not.

CONCLUSIONS

Population structure of Msi was driven by patterns of warming since the LGM, and secondarily by geographical barriers. This study will facilitate germplasm conservation, association analyses and identification of potential heterotic groups for the improvement of Miscanthus as a bioenergy crop.

Adaptations between ecotypes and along environmental gradients in Panicum virgatum

Determining the patterns and mechanisms of natural selection in the wild is of fundamental importance to understanding the differentiation of populations and the evolution of new species. However, it is often unknown the extent to which adaptive genetic variation is distributed among ecotypes between distinct habitats versus along large-scale geographic environmental gradients, such as those that track latitude. Classic studies of selection in the wild in switchgrass, Panicum virgatum, tested for adaptation at both of these levels of natural variation. Here we review what these field experiments and modern agronomic field trials have taught us about natural variation and selection at both the ecotype and environmental gradient levels in P. virgatum. With recent genome sequencing efforts in P. virgatum, it is poised to become an excellent system for understanding the adaptation of grassland species across the eastern half of North America. The identification of genetic loci involved in different types of adaptations will help to understand the evolutionary mechanisms of diversification within P. virgatum and provide useful information for the breeding of high-yielding cultivars for different ecoregions.

A population genetic transect of Panicum hallii (Poaceae)

PREMISE OF STUDY

Understanding the relationship between climate, adaptation, and population structure is of fundamental importance to botanists because these factors are crucial for the evolution of biodiversity and the response of species to future climate change. Panicum hallii is an emerging model system for perennial grass and bioenergy research, yet very little is known about the relationship between climate and population structure in this system. •

METHODS

We analyzed geographic population differentiation across 39 populations of P. hallii along a longitudinal transect from the savannas of central Texas through the deserts of Arizona and New Mexico. A combination of morphological and genetic (microsatellite) analysis was used to explore patterns of population structure. •

KEY RESULTS

We found strong differentiation between high elevation western desert populations and lower elevation eastern populations of P. hallii, with a pronounced break in structure occurring in western Texas. In addition, we confirmed that there are high levels of morphological and genetic structure between previous recognized varieties (var. hallii and var. filipes) within this species. •

CONCLUSIONS

The results of this study suggest that patterns of population structure within P. hallii may be driven by climatic variation over space. Overall, this study lays the groundwork for future studies on the genetics of local adaptation and reproductive isolation in this system.

Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol

Switchgrass (Panicum virgatum L.) is a perennial grass that has been designated as an herbaceous model biofuel crop for the United States of America. To facilitate accelerated breeding programs of switchgrass, we developed both an association panel and linkage populations for genome-wide association study (GWAS) and genomic selection (GS). All of the 840 individuals were then genotyped using genotyping by sequencing (GBS), generating 350 GB of sequence in total. As a highly heterozygous polyploid (tetraploid and octoploid) species lacking a reference genome, switchgrass is highly intractable with earlier methodologies of single nucleotide polymorphism (SNP) discovery. To access the genetic diversity of species like switchgrass, we developed a SNP discovery pipeline based on a network approach called the Universal Network-Enabled Analysis Kit (UNEAK). Complexities that hinder single nucleotide polymorphism discovery, such as repeats, paralogs, and sequencing errors, are easily resolved with UNEAK. Here, 1.2 million putative SNPs were discovered in a diverse collection of primarily upland, northern-adapted switchgrass populations. Further analysis of this data set revealed the fundamentally diploid nature of tetraploid switchgrass. Taking advantage of the high conservation of genome structure between switchgrass and foxtail millet (Setaria italica (L.) P. Beauv.), two parent-specific, synteny-based, ultra high-density linkage maps containing a total of 88,217 SNPs were constructed. Also, our results showed clear patterns of isolation-by-distance and isolation-by-ploidy in natural populations of switchgrass. Phylogenetic analysis supported a general south-to-north migration path of switchgrass. In addition, this analysis suggested that upland tetraploid arose from upland octoploid. All together, this study provides unparalleled insights into the diversity, genomic complexity, population structure, phylogeny, phylogeography, ploidy, and evolutionary dynamics of switchgrass.

Recent advances in sequencing technologies have enabled large-scale surveys of genetic diversity in model species with a wholly or partly sequenced reference genome. However, thousands of key species, which are essential for food, health, energy, and ecology, do not have reference genomes. To accelerate their breeding cycle via marker assisted selection, high-throughput genotyping is required for these valuable species, in spite of the absence of reference genomes. Based on genotyping by sequencing (GBS) technology, we developed a new single nucleotide polymorphism (SNP) discovery protocol, the Universal Network-Enabled Analysis Kit (UNEAK), which can be widely used in any species, regardless of genome complexity or the availability of a reference genome. Here we test this protocol on switchgrass, currently the prime energy crop species in the United States of America. In addition to the discovery of over a million SNPs and construction of high-density linkage maps, we provide novel insights into the genome complexity, ploidy, phylogeny, and evolution of switchgrass. This is only the beginning: we believe UNEAK offers the key to the exploration and exploitation of the genetic diversity of thousands of non-model species.

Development of a genome-wide multiple duplex-SSR protocol and its applications for the identification of selfed progeny in switchgrass

Background

Switchgrass (Panicum virgatum) is a herbaceous crop for the cellulosic biofuel feedstock development in the USA and Europe. As switchgrass is a naturally outcrossing species, accurate identification of selfed progeny is important to producing inbreds, which can be used in the production of heterotic hybrids. Development of a technically reliable, time-saving and easily used marker system is needed to quantify and characterize breeding origin of progeny plants of targeted parents.

Results

Genome-wide screening of 915 mapped microsatellite (simple sequence repeat, SSR) markers was conducted, and 842 (92.0%) produced clear and scorable bands on a pooled DNA sample of eight switchgrass varieties. A total of 166 primer pairs were selected on the basis of their relatively even distribution in switchgrass genome and PCR amplification quality on 16 tetraploid genotypes. Mean polymorphic information content value for the 166 markers was 0.810 ranging from 0.116 to 0.959. From them, a core set of 48 loci, which had been mapped on 17 linkage groups, was further tested and optimized to develop 24 sets of duplex markers. Most of (up to 87.5%) targeted, but non-allelic amplicons within each duplex were separated by more than 10-bp. Using the established duplex PCR protocol, selfing ratio (i.e., selfed/all progeny x100%) was identified as 0% for a randomly selected open-pollinated ‘Kanlow’ genotype grown in the field, 15.4% for 22 field-grown plants of bagged inflorescences, and 77.3% for a selected plant grown in a growth chamber.

Conclusions

The study developed a duplex SSR-based PCR protocol consisting of 48 markers, providing ample choices of non-tightly-linked loci in switchgrass whole genome, and representing a powerful, time-saving and easily used method for the identification of selfed progeny in switchgrass. The protocol should be a valuable tool in switchgrass breeding efforts.

Karyotype variation is indicative of subgenomic and ecotypic differentiation in switchgrass

BACKGROUND

Karyotypes can provide information about taxonomic relationships, genetic aberrations, and the evolutionary origins of species. However, differentiation of the tiny chromosomes of switchgrass (Panicum virgatum L.) and creation of a standard karyotype for this bioenergy crop has not been accomplished due to lack of distinguishing features and polyploidy.

RESULTS

A cytogenetic study was conducted on a dihaploid individual (2n = 2X = 18) of switchgrass to establish a chromosome karyotype. Size differences, condensation patterns, and arm-length ratios were used as identifying features and fluorescence in-situ hybridization (FISH) assigned 5S and 45S rDNA loci to chromosomes 7 and 2 respectively. Both a maize CentC and a native switchgrass centromeric repeat (PviCentC) that shared 73% sequence identity demonstrated a strong signal on chromosome 3. However, only the PviCentC probe labeled the centromeres of all chromosomes. Unexpected PviCentC and 5S rDNA hybidization patterns were consistent with severe reduction or total deletion of these repeats in one subgenome. These patterns were maintained in tetraploid and octoploid individuals. The 45S rDNA repeat produced the expected number of loci in dihaploid, tetraploid and octoploid individuals. Differences observed at the 5S rDNA loci between the upland and lowland ecotypes of switchgrass provided a basis for distinguishing these subpopulations.

CONCLUSION

Collectively, these results provide a quantitative karyotype of switchgrass chromosomes. FISH analyses indicate genetic divergence between subgenomes and allow for the classification of switchgrass plants belonging to divergent genetic pools. Furthermore, the karyotype structure and cytogenetic analysis of switchgrass provides a framework for future genetic and genomic studies.

Transcriptome analysis and gene expression atlas for Panicum hallii var. filipes, a diploid model for biofuel research

Panicum hallii is an emerging model for genetic studies of agronomic traits in Panicum, presenting a tractable diploid alternative study system to the tetra- or octaploid biofuel crop switchgrass (Panicum virgatum). To characterize the gene complement in P. hallii var. filipes and enable gene expression analysis in this system we sequenced, assembled, and annotated the transcriptome. Over 300 Mb of normalized cDNA prepared from multiple tissues and treatments was sequenced using 454-Titanium, producing an annotated assembly including 15 422 unique gene names. Comparison with other grass genomes identified putative P. hallii homologs for >14 000 previously characterized genes. We also developed an atlas of gene expression across tissues and stages using RNA-Seq (the quantitative analysis of short cDNA reads). SOLiD sequencing and quantitative analysis of more than 40 million cDNA tags identified substantial variation in expression profiles among tissues, consistent with known functional differences. Putative homologs were found for all enzymes in the phenylpropanoid pathway leading to lignin biosynthesis, including genes with known effects on biomass conversion efficiency. The resources developed here will enable studies of the genes underlying variation in cell wall composition, drought tolerance, and biomass production in Panicum.

A high-density simple sequence repeat-based genetic linkage map of switchgrass

Switchgrass (Panicum virgatum) has been identified as a promising cellulosic biofuel crop in the United States. Construction of a genetic linkage map is fundamental for switchgrass molecular breeding and the elucidation of its genetic mechanisms for economically important traits. In this study, a novel population consisting of 139 selfed progeny of a northern lowland genotype, NL 94 LYE 16X13, was used to construct a linkage map. A total of 2493 simple sequence repeat markers were screened for polymorphism. Of 506 polymorphic loci, 80.8% showed a goodness-of-fit of 1:2:1 segregation ratio. Among 469 linked loci on the framework map, 241 coupling vs. 228 repulsion phase linkages were detected that conformed to a 1:1 ratio, confirming disomic inheritance. A total of 499 loci were mapped to 18 linkage groups (LG), of which the cumulative length was 2085.2 cM, with an average marker interval of 4.2 cM. Nine homeologous LG pairs were identified based on multi-allele markers and comparative genomic analysis. Two clusters of segregation-distorted loci were identified on LG 5b and 9b, respectively. Comparative analysis indicated a one-to-one relationship between nine switchgrass homeologous groups and nine foxtail millet (Setaria italica) chromosomes, suggesting strong homology between the two species. The linkage map derived from selfing a heterozygous parent, instead of two separate maps usually constructed for a cross-fertilized species, provides a new genetic framework to facilitate genomics research, quantitative trait locus (QTL) mapping, and marker-assisted breeding.

Genomic diversity in switchgrass (Panicum virgatum): from the continental scale to a dune landscape

Connecting broad-scale patterns of genetic variation and population structure to genetic diversity on a landscape is a key step towards understanding historical processes of migration and adaptation. New genomic approaches can be used to increase the resolution of phylogeographic studies while reducing locus sampling effects and circumventing ascertainment bias. Here, we use a novel approach based on high-throughput sequencing to characterize genetic diversity in complete chloroplast genomes and >10,000 nuclear loci in switchgrass, at continental and landscape scales. Switchgrass is a North American tallgrass species, which is widely used in conservation and perennial biomass production, and shows strong ecotypic adaptation and population structure across the continental range. We sequenced 40.9 billion base pairs from 24 individuals from across the species' range and 20 individuals from the Indiana Dunes. Analysis of plastome sequence revealed 203 variable SNP sites that define eight haplogroups, which are differentiated by 4-127 SNPs and confirmed by patterns of indel variation. These include three deeply divergent haplogroups, which correspond to the previously described lowland-upland ecotypic split and a novel upland haplogroup split that dates to the mid-Pleistocene. Most of the plastome haplogroup diversity present in the northern switchgrass range, including in the Indiana Dunes, originated in the mid- or upper Pleistocene prior to the most recent postglacial recolonization. Furthermore, a recently colonized landscape feature (approximately 150 ya) in the Indiana Dunes contains several deeply divergent upland haplogroups. Nuclear markers also support a deep lowland-upland split, followed by limited gene flow, and show extensive gene flow in the local population of the Indiana Dunes.