Hierarchical classification of switchgrass genotypes using SSR and chloroplast sequences: ecotypes,ploidies, gene pools,and cultivars

Switchgrass (Panicum virgatum L.) is an important crop for bioenergy feedstock development. Switchgrass has two main ecotypes: the lowland ecotype being exclusively tetraploid (2n = 4x = 36) and the upland ecotype being mainly tetraploid and octaploid (2n = 8x = 72). Because there is a significant difference in ploidy, morphology, growth pattern, and zone of adaptation between and within the upland and lowland ecotypes, it is important to discriminate switchgrass plants belonging to different genetic pools. We used 55 simple sequence repeats (SSR) loci and six chloroplast sequences to identify patterns of variation between and within 18 switchgrass cultivars representing seven lowland and 11 upland cultivars from different geographic regions and of varying ploidy levels. We report consistent discrimination of switchgrass cultivars into ecotype membership and demonstrate unambiguous molecular differentiation among switchgrass ploidy levels using genetic markers. Also, SSR and chloroplast markers identified genetic pools related to the geographic origin of the 18 cultivars with respect to ecotype, ploidy, and geographical, and cultivar sources. SSR loci were highly informative for cultivar fingerprinting and to classify plants of unknown origin. This classification system is the first step toward developing switchgrass complementary gene pools that can be expected to provide a significant heterotic increase in biomass yield.     

Effects of ecotypes and morphotypes in feedstock composition of switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) is a C4 grass with high biomass yield potential and is now a model species for the Bioenergy Feedstock Development Program. Two distinct ecotypes (e.g., upland and lowland) and a range of plant morphotypes (e.g., leafy and stemmy) have been observed in switchgrass. The objective of this study was to determine the influence of ecotype and morphotype on biomass feedstock quality. Leaf and stem tissues of leafy and stemmy morphotypes from both lowland and upland ecotypes were analyzed for key feedstock traits. The leaf : stem ratio of leafy morphotype was more than 40% higher than the stemmy morphotype in both upland and lowland ecotypes. Therefore, the stemmy morphotype has significant advantages over leafy morphotype during harvesting, storage, transportation and finally the feedstock quality. Remarkable differences in feedstock quality and mineral composition were observed in switchgrass genotypes with distinct ecotypic origins and variable plant morphotypes. Lignin, hemicelluloses and cellulose concentrations were higher in stems than in the leaves, while ash content was notably high in leaves. A higher concentration of potassium was found in the stems compared to the leaves. In contrast, calcium was higher and magnesium was generally higher in the leaves compared to stems. The upland genotypes demonstrated considerably higher lignin (14.4%) compared with lowland genotypes (12.4%), while hemicellulose was higher in lowland compared with upland. The stemmy type demonstrated slightly higher lignin compared with leafy types (P < 0.1). Differences between the ecotypes and morphotypes for key quality traits demonstrated the potential for improving feedstock composition of switchgrass through selection in breeding programs.     

Tolerance of switchgrass to extreme soil moisture stress: Ecological implications

Switchgrass (Panicum virgatum L.), a native of eastern and central North America, is a leading candidate as a dedicated biofuel feedstock in the US due to its broad adaptability, rapid growth rate, and ability to grow in low production soils. To begin to characterize the important agronomic and ecological traits related to environmental tolerance of switchgrass, we evaluated fitness under stressful growing conditions. We assessed the germination, establishment, performance, and reproductive potential of four common accessions, both upland and lowland ecotypes, at various levels of soil moisture availability (moisture deficit to flooded) in the greenhouse. Seeds emerged and established (55–90% survival) under all soil moisture conditions (−0.3 MPa to flooded). Transplants of lowland ecotypes performed as well in flooded conditions as in field capacity controls, though flooding reduced performance of upland ecotypes. Drought treatments (−4.0 and −11.0 MPa) reduced tiller length and number, leaf area, and biomass production by up to 80%. However, once established, all plants survived at −4.0 MPa and had the same proportion of tillers in flower as at field capacity. The ability of switchgrass to germinate, establish, and flower in low moisture and flooded conditions, particularly lowland ecotypes, may increase the range of environments suitable for biofuel cultivation, and can serve as a baseline for further ecological studies and genetic improvement.     

Genetic diversity in tetraploid switchgrass revealed by AFLP marker polymorphisms

Switchgrass (Panicum virgatum) is a perennial warm-season grass native to North America that has been identified as a dedicated cellulosic biofuel crop. We quantified genetic diversity in tetraploid switchgrass germplasm collected at Oklahoma State University and characterized genetic relatedness among the collections from distinct regions. Fifty-six tetraploid accessions, including seven upland and 49 lowland genotypes from throughout the US, were examined. The amplified fragment length polymorphism (AFLP) procedure was utilized to generate DNA profiling patterns that were scored visually. Sixteen selective AFLP primer combinations were used to amplify 452 polymorphic bands. The accessions' genetic similarity coefficients, UPGMA (unweighted pair-group method with arithmetic averaging) cluster analysis and principle coordinate analysis, were performed. The upland and lowland accessions clustered according to ecotypes, with one exception (TN104). Genetic similarity coefficients among the accessions ranged from 0.73 to 0.95. Analysis of molecular variance (AMOVA) was performed, showing significant differences between the upland and lowland genotypes. The trnL marker confirmed that TN104 was a lowland genotype, but the trnL marker identification of upland and lowland genotypes was not consistent with the AFLP analysis in two germplasms (Miami and AR4).     

Genomic landscape of parallel domestication of upland rice and its implications

Parallel domestication has been widely acknowledged but itsgenetic basis remains largely unclear. As an important rice ecotype, upland rice was assumedly domesticated multiple times in two rice subspecies (Indica and Japonica) and provides a feasible system to explore the genetic basis of parallel domestication. To uncover the genome‐wide pattern of genetic differentiation between upland and lowland rice and explore the parallelism of genetic changes during upland rice domestication, we obtained whole‐genome sequences of 95 rice landraces and yielded genome‐wide expression data for five tissues of representative accessions of upland and lowland rice. Our phylogenetic analyses confirmed multiple domestications of the upland ecotype in two rice subspecies. Genomic scans based on resequencing data identified substantial differentiation between the upland and lowland ecotypes with 11.4% and 14.8% of the genome diverged between the two ecotypes in Indica and Japonica, respectively. Further genome‐wide gene expression analyses found that 30% of effectively expressed genes were significantly differentiated between two ecotypes, indicating the importance of regulation changes in the domestication of upland rice. Importantly, we found that only 1.8% of differentiated genomes and 1.6% of differentially expressed genes were shared by upland Indica and upland Japonica, suggestive of largely unparallel genetic alterations during upland rice domestication. These findings not only provide new insights into the genetic basis of parallel domestication at the genome scale but could also facilitate geneticimprovement and breeding of rice and crops in general.     

Hybridization,polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum)

In plant species, variation in levels of clonality, ploidy and interspecific hybridization can interact to influence geographic patterns of genetic diversity. These factors commonly vary in plants that specialize on saline habitats (halophytes) and may play a role in how they adapt to salinity variation across their range. One such halophyte is the turfgrass and emerging genomic model system seashore paspalum (Paspalum vaginatum Swartz). To investigate how clonal propagation, ploidy variation, and interspecific hybridization vary across ecotypes and local salinity levels in wild P. vaginatum, we employed genotyping‐by‐sequencing, cpDNA sequencing and flow cytometry in 218 accessions representing > 170 wild collections from throughout the coastal southern United States plus USDA germplasm. We found that the two morphologically distinct ecotypes of P. vaginatum differ in their adaptive strategies. The fine‐textured ecotype is diploid and appears to reproduce in the wild both sexually and by clonal propagation; in contrast, the coarse‐textured ecotype consists largely of clonally‐propagating triploid and diploid genotypes. The coarse‐textured ecotype appears to be derived from hybridization between fine‐textured P. vaginatum and an unidentified Paspalum species. These clonally propagating hybrid genotypes are more broadly distributed than clonal fine‐textured genotypes and may represent a transition to a more generalist adaptive strategy. Additionally, the triploid genotypes vary in whether they carry one or two copies of the P. vaginatum subgenome, indicating multiple evolutionary origins. This variation in subgenome composition shows associations with local ocean salinity levels across the sampled populations and may play a role in local adaptation.     

SNP Discovery with EST and NextGen Sequencing in Switchgrass (Panicum virgatum L.)

Although yield trials for switchgrass (Panicum virgatum L.), a potentially high value biofuel feedstock crop, are currently underway throughout North America, the genetic tools for crop improvement in this species are still in the early stages of development. Identification of high-density molecular markers, such as single nucleotide polymorphisms (SNPs), that are amenable to high-throughput genotyping approaches, is the first step in a quantitative genetics study of this model biofuel crop species. We generated and sequenced expressed sequence tag (EST) libraries from thirteen diverse switchgrass cultivars representing both upland and lowland ecotypes, as well as tetraploid and octoploid genomes. We followed this with reduced genomic library preparation and massively parallel sequencing of the same samples using the Illumina Genome Analyzer technology platform. EST libraries were used to generate unigene clusters and establish a gene-space reference sequence, thus providing a framework for assembly of the short sequence reads. SNPs were identified utilizing these scaffolds. We used a custom software program for alignment and SNP detection and identified over 149,000 SNPs across the 13 short-read sequencing libraries (SRSLs). Approximately 25,000 additional SNPs were identified from the entire EST collection available for the species. This sequencing effort generated data that are suitable for marker development and for estimation of population genetic parameters, such as nucleotide diversity and linkage disequilibrium. Based on these data, we assessed the feasibility of genome wide association mapping and genomic selection applications in switchgrass. Overall, the SNP markers discovered in this study will help facilitate quantitative genetics experiments and greatly enhance breeding efforts that target improvement of key biofuel traits and development of new switchgrass cultivars.     

Inter-Annual Precipitation Variability Decreases Switchgrass Productivity from Arid to Mesic Environments

Cellulosic biofuels are an important source of renewable biomass within the alternative energy portfolio. Switchgrass (Panicum virgatum L.), a perennial C₄ grass native to North America, is widely studied as a biofuel feedstock for its consistently high yields and minimal input requirements. The influences of precipitation amount and temporal variability on the fertilizer response of switchgrass productivity are not fully understood. Moreover, global climate models predict changes in rainfall patterns towards lower and increasingly variable soil water availability in several productive areas worldwide, which may impact net primary production of biofuel crops. We conducted a meta-analysis of aboveground net primary production of switchgrass from 48 publications encompassing 82 different locations, 11 soil types, 52 switchgrass cultivars, fertilizer inputs between 0 to 896 kg N ha^?1 year^?1, and 1 to 6 years of annual productivity measures repeated on the same stand. Productivity of the lowland ecotype doubled with N rates >?131 kg N ha^?1 year^?1, but upland ecotype productivity increased only by 50%. Results showed an optimum N rate of 30 to 60 kg N ha^?1 year^?1 for both ecotypes, after which biomass gain per unit of N added decreased. Growing season precipitation (GSPPT) and inter-annual precipitation variability (inter-PPTvar) affected both ecotypes similarly. Long-term mean annual precipitation (MAP) differentially affected lowland and upland productivity, depending on the N level. Productivity responses to MAP and GSPPT were similar for both upland and lowland ecotypes at none or low N rates. When N increased beyond 60 kg N ha^?1 year^?1, lowland cultivars had a greater growth response to MAP than uplands. Productivity increased with increasing GSPPT and MAP and had a positive linear response to MAP ranging from 600 to 1200 mm year^?1. One third of the variability in switchgrass production was accounted for by inter-PPTvar. After accounting for MAP, sites with higher inter-PPTvar had lower switchgrass productivity than sites with lower inter-PPTvar. Increased inter-annual variation in precipitation reduced production of both ecotypes. Predicted changes in the amount and timing of precipitation thus likely will exert greater influence on production of upland than lowland ecotypes of switchgrass.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Copy number variants in Italian Large White pigs detected using high‐density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance‐tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with penncnv software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by penncnv were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.     

Mismatch in the distribution of floral ecotypes and pollinators: insights into the evolution of sexually deceptive orchids

Plants are predicted to show floral adaptation to geographic variation in the most effective pollinator, potentially leading to reproductive isolation and genetic divergence. Many sexually deceptive orchids attract just a single pollinator species, limiting opportunities to experimentally investigate pollinator switching. Here, we investigate Drakaea concolor, which attracts two pollinator species. Using pollinator choice tests, we detected two morphologically similar ecotypes within D. concolor. The common ecotype only attracted Zaspilothynnus gilesi, whereas the rare ecotype also attracted an undescribed species of Pogonothynnus. The rare ecotype occurred at populations nested within the distribution of the common ecotype, with no evidence of ecotypes occurring sympatrically. Surveying for pollinators at over 100 sites revealed that ecotype identity was not correlated with wasp availability, with most orchid populations only attracting the rare Z. gilesi. Using microsatellite markers, genetic differentiation among populations was very low (G_ST = 0.011) regardless of ecotype, suggestive of frequent gene flow. Taken together, these results may indicate that the ability to attract Pogonothynnus has evolved recently, but this ecotype is yet to spread. The nested distribution of ecotypes, rather than the more typical formation of ecotypes in allopatry, illustrates that in sexually deceptive orchids, pollinator switching could occur throughout a species' range, resulting from multiple potentially suitable but unexploited pollinators occurring in sympatry. This unusual case of sympatric pollinators highlights D. concolor as a promising study system for further understanding the process of pollinator switching from ecological, chemical and genetic perspectives.     

Postglacial climate changes and rise of three ecotypes of harbour porpoises,Phocoena phocoena,in western Palearctic waters

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base‐pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23–19 kilo‐years ago, kyrbp ). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp ) than either does with the European continental shelf ecotype (c. 28 kyrbp ), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post‐glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.     

Copy number variation associates with mortality in long‐lived individuals: a genome‐wide assessment

Copy number variants (CNVs) represent a significant source of genetic variation in the human genome and have been implicated in numerous diseases and complex traits. To date, only a few studies have investigated the role of CNVs in human lifespan. To investigate the impact of CNVs on prospective mortality at the extreme end of life, where the genetic component of lifespan appears most profound, we analyzed genomewide CNV data in 603 Danish nonagenarians and centenarians (mean age 96.9 years, range 90.0–102.5 years). Replication was performed in 500 long‐lived individuals from the Leiden Longevity Study (mean age 93.2 years, range 88.9–103.4 years). First, we assessed the association between the CNV burden of each individual (the number of CNVs, the average CNV length, and the total CNV length) and mortality and found a significant increase in mortality per 10 kb increase in the average CNV length, both for all CNVs (hazard ratio (HR) = 1.024, P = 0.002) and for duplications (HR = 1.011, P = 0.005), as well as per 100 kb increase in the total length of deletions (HR = 1.009, P = 0.0005). Next, we assessed the relation between specific deletions and duplications and mortality. Although no genome–wide significant associations were discovered, we identified six deletions and one duplication that showed consistent association with mortality in both or either of the sexes across both study populations. These results indicate that the genome–wide CNV burden, specifically the average CNV length and the total CNV length, associates with higher mortality in long‐lived individuals.     

Local adaptation,genetic divergence,and experimental selection in a foundation grass across the US Great Plains’ climate gradient

Many prior studies have uncovered evidence for local adaptation using reciprocal transplant experiments. However, these studies are rarely conducted for a long enough time to observe succession and competitive dynamics in a community context, limiting inferences for long‐lived species. Furthermore, the genetic basis of local adaptation and genetic associations with climate has rarely been identified. Here, we report on a long‐term (6‐year) experiment conducted under natural conditions focused on Andropogon gerardii, the dominant grass of the North American Great Plains tallgrass ecosystem. We focus on this foundation grass that comprises 80% of tallgrass prairie biomass and is widely used in 20,000 km² of restoration. Specifically, we asked the following questions: (a) Whether ecotypes are locally adapted to regional climate in realistic ecological communities. (b) Does adaptive genetic variation underpin divergent phenotypes across the climate gradient? (c) Is there evidence of local adaptation if the plants are exposed to competition among ecotypes in mixed ecotype plots? Finally, (d) are local adaptation and genetic divergence related to climate? Reciprocal gardens were planted with 3 regional ecotypes (originating from dry, mesic, wet climate sources) of Andropogon gerardii across a precipitation gradient (500–1,200 mm/year) in the US Great Plains. We demonstrate local adaptation and differentiation of ecotypes in wet and dry environments. Surprisingly, the apparent generalist mesic ecotype performed comparably under all rainfall conditions. Ecotype performance was underpinned by differences in neutral diversity and candidate genes corroborating strong differences among ecotypes. Ecotype differentiation was related to climate, primarily rainfall. Without long‐term studies, wrong conclusions would have been reached based on the first two years. Further, restoring prairies with climate‐matched ecotypes is critical to future ecology, conservation, and sustainability under climate change.     

Retention of a chromosomal inversion from an anadromous ancestor provides the genetic basis for alternative freshwater ecotypes in rainbow trout

Migratory behaviour patterns in animals are controlled by a complex genetic architecture. Rainbow trout (Oncorhynchus mykiss) is a salmonid fish that spawns in streams but exhibits three primary life history pathways: stream‐resident (fluvial), lake‐migrant (adfluvial) and ocean‐migrant (anadromous). Previous studies examining fluvial and anadromous O. mykiss have identified several genes associated with life history divergence including the presence of an inversion complex within chromosome 5 (Omy05) that appears to maintain a suite of linked genes controlling migratory behaviour. However, adfluvial trout are migratory without being anadromous, and the genetic basis for this life history has not been investigated from evolutionary perspectives. We sampled wild, native nonanadromous rainbow trout occupying connected stream and lake habitats in a southwest Alaskan watershed to determine whether these fish exhibit genetic divergence between fluvial and adfluvial ecotypes, and whether that divergence parallels that documented in fluvial and anadromous O. mykiss. Data from restriction site‐associated DNA (RAD) sequencing revealed an association between frequencies of both the Omy05 inversion complex and other single nucleotide polymorphisms (SNPs) with habitat type (stream or lake), supporting the genetic divergence of fluvial and adfluvial individuals in sympatry. The presence of a genetic basis for migration into lakes, analogous to that documented for anadromy, indicates that the adfluvial ecotype must be recognized separately from the fluvial form of O. mykiss even though neither is anadromous. These results highlight the genetic architecture underlying migration and the importance of chromosomal inversions in promoting and sustaining intraspecific diversity.     

Chloroplast genome variation in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) exists at multiple ploidies and two phenotypically distinct ecotypes. To facilitate interploidal comparisons and to understand the extent of sequence variation within existing breeding pools, two complete switchgrass chloroplast genomes were sequenced from individuals representative of the upland and lowland ecotypes. The results demonstrated a very high degree of conservation in gene content and order with other sequenced plastid genomes. The lowland ecotype reference sequence (Kanlow Lin1) was 139,677 base pairs while the upland sequence (Summer Lin2) was 139,619 base pairs. Alignments between the lowland reference sequence and short-read sequence data from existing sequence datasets identified as either upland or lowland confirmed known polymorphisms and indicated the presence of other differences. Insertions and deletions principally occurred near stretches of homopolymer simple sequence repeats in intergenic regions while most Single Nucleotide Polymorphisms (SNPs) occurred in intergenic regions and introns within the single copy portions of the genome. The polymorphism rate between upland and lowland switchgrass ecotypes was found to be similar to rates reported between chloroplast genomes of indica and japonica subspecies of rice which were believed to have diverged 0.2-0.4 million years ago.     

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves—an asthmalike inflammatory disease in horses. Analysis of 16 RAO‐susceptible (cases) and six RAO‐resistant (control) horses on a custom‐made whole‐genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO‐affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.