Genotypic variation of cell wall composition and its conversion efficiency in Miscanthus sinensis,a potential biomass feedstock crop in China

Plant cell walls are composed of cellulose microfibrils embedded in a cross‐linked‐net of matrix polysaccharides and co‐polymerized with lignin. The study presented the genotypic variations of cell wall composition, biohydrogen production, and lignocellulose degradation ratio in a collection of 102 Miscanthus sinensis (M. Sinensis, hereafter) accessions collected from a wide geographical range in China. Significant variations were observed for the determined traits, cellulose content, hemicellulose content, cellulose and hemicellulose degradation efficiency, and biohydrogen yield. The cellulose, hemicellulose, and lignin contents ranged from 30.20–44.25, 28.97–42.65, and 6.96–20.75%, respectively. The degradation ratio of cellulose and hemicellulose varied from 2.08% to 37.87% and from 14.71% to 52.50%, respectively. The feedstock was fermented to produce biohydrogen, and the production varied from 14.59 to 40.66 ml per gram of Miscanthus biomass. The expression profile of three cellulose synthase (MsCesA) genes was initially established to indicate the genotypic difference among the M. sinensis accessions. Pearson's correlations were conducted to reveal the perplexing relationship between the tested traits, biohydrogen yield, cell wall composition and their degradation efficiency. In addition, the relationship pattern, between the test traits and the geographic factors corresponding with the original place, was investigated. The result showed that the significant variation among the M. sinensis genotypes is the result of natural selection in different environments of their original habitats. Improvement in cell wall composition and structure and enhancement of lignocellulose degradation ratio could significantly increase sustainable bioenergy production.     

Influence of salt stress on C4 photosynthesis in Miscanthus sinensis Anderss.

• Miscanthus sinensis Anderss. is a good candidate for C₄ bioenergy crop development for marginal lands. As one of the characteristics of marginal lands, salinization is a major limitation to agricultural production. The present work aimed to investigate the possible factors involved in the tolerance of M. sinensis C₄ photosynthesis to salinity stress.
• Seedlings of two accessions (salt‐tolerant ‘JM0119’ and salt‐sensitive ‘JM0099’) were subjected to 0 mm NaCl (control) or 250 mm NaCl (salt stress treatment) for 2 weeks. The chlorophyll content, parameters of photosynthesis and chlorophyll a fluorescence, activity of C₄ enzymes and expression of C₄ genes were measured.
• The results showed that photosynthesis rate, transpiration rate, chlorophyll content, PSII operating efficiency, coefficient of photochemical quenching, activity of phosphoenolpyruvate carboxylase (PEPC) and pyruvate, orthophosphate dikinase (PPDK) and gene expression of PEPC and PPDK under salinity were higher after long‐term salinity exposure in ‘JM0119’ than in ‘JM0099’, while activity of NADP‐malate dehydrogenase (NADP‐MDH) and NADP‐malic enzyme (NADP‐ME), together with expression of NADP‐MDH and NADP‐ME, were much higher in ‘JM0099’ than in ‘JM0119’.
• In conclusion, the increased photosynthetic capacity under long‐term salt stress in the salt‐tolerant relative to the salt‐sensitive M. sinensis accession was mainly associated with non‐stomatal factors, such as reduced chlorophyll loss, higher PSII operating efficiency, enhanced activity of PEPC and PPDK and relatively lower activity of NADP‐ME.

     

Genetic variation in Miscanthus × giganteus and the importance of estimating genetic distance thresholds for differentiating clones

Miscanthus × giganteus (Mxg) is an important bioenergy feedstock crop, however, genetic diversity among legacy cultivars may be severely constrained. Only one introduction from Japan to Denmark of this sterile, triploid, vegetatively propagated crop was recorded in the 1930s. We sought to determine if the Mxg cultivars in North America were all synonyms, and if they were derived from the European introduction. We used 64 nuclear and five chloroplast simple sequence repeat (SSR) markers to estimate genetic similarity for 27 Mxg accessions from North America, and compared them with six accessions from Europe, including the species’ type‐specimen. A subset of accessions was also evaluated by restriction‐site associated DNA sequencing (RAD‐seq). In addition, we assessed the potential of new crosses to increase Mxg genetic diversity by comparing eight new triploid Mxg progeny grown from seed, along with samples of the parental species M. sacchariflorus and M. sinensis. Estimates of genotyping error rates were essential for distinguishing between experimental error and true genotypic differences among accessions. Given differences in estimated error rates and costs per marker for SSRs and RAD‐seq, the former is currently more cost‐effective for determining if two accessions are genetically identical. We concluded that all of the Mxg legacy cultivars were derived via vegetative propagation from a single genet. In contrast with the Mxg legacy cultivars, genetic similarity to the type‐specimen of eight new triploid Mxg progeny ranged from 0.46 to 0.56. Though genetic diversity among the Mxg legacy cultivars is critically low, new crosses can provide much‐needed variation to growers.     

Genome‐wide association and genomic prediction for biomass yield in a genetically diverse Miscanthus sinensis germplasm panel phenotyped at five locations in Asia and North America

To improve the efficiency of breeding of Miscanthus for biomass yield, there is a need to develop genomics‐assisted selection for this long‐lived perennial crop by relating genotype to phenotype and breeding value across a broad range of environments. We present the first genome‐wide association (GWA) and genomic prediction study of Miscanthus that utilizes multilocation phenotypic data. A panel of 568 Miscanthus sinensis accessions was genotyped with 46,177 single nucleotide polymorphisms (SNPs) and evaluated at one subtropical and five temperate locations over 3 years for biomass yield and 14 yield‐component traits. GWA and genomic prediction were performed separately for different years of data in order to assess reproducibility. The analyses were also performed for individual field trial locations, as well as combined phenotypic data across groups of locations. GWA analyses identified 27 significant SNPs for yield, and a total of 504 associations across 298 unique SNPs across all traits, sites, and years. For yield, the greatest number of significant SNPs was identified by combining phenotypic data across all six locations. For some of the other yield‐component traits, greater numbers of significant SNPs were obtained from single site data, although the number of significant SNPs varied greatly from site to site. Candidate genes were identified. Accounting for population structure, genomic prediction accuracies for biomass yield ranged from 0.31 to 0.35 across five northern sites and from 0.13 to 0.18 for the subtropical location, depending on the estimation method. Genomic prediction accuracies of all traits were similar for single‐location and multilocation data, suggesting that genomic selection will be useful for breeding broadly adapted M. sinensis as well as M. sinensis optimized for specific climates. All of our data, including DNA sequences flanking each SNP, are publicly available. By facilitating genomic selection in M. sinensis and Miscanthus × giganteus, our results will accelerate the breeding of these species for biomass in diverse environments.     

Uncovering the genetic diversity of yams (Dioscorea spp.) in China by combining phenotypic trait and molecular marker analyses

Yam is an important edible tuber and root plant worldwide; China as one of the native places of yams has many diverse local resources. The goal of this study was to clarify the genetic diversity of the commonly cultivated yam landraces and the genetic relationship between the main yam species in China. In this study, 26 phenotypic traits of 112 yam accessions from 21 provinces in China were evaluated, and 24 simple sequence repeat (SSR) and 29 sequence‐related amplified polymorphism (SRAP) markers were used for the genetic diversity analysis. Phenotypic traits revealed that Dioscorea opposita had the highest genetic diversity, followed by D. alata, D. persimilis, D. fordii, and D. esculenta. Among the 26 phenotypic traits, the Shannon diversity indexes of leaf shape, petiole color, and stem color were high, and the range in the variation of tuber‐related traits in the underground part was higher than that in the aboveground part. All accessions were divided into six groups by phenotypic trait clustering, which was also supported by principal component analysis (PCA). Molecular marker analysis showed that SSR and SRAP markers had good amplification effects and could effectively and accurately evaluate the genetic variation of yam. The unweighted pair‐group method with arithmetic means analysis based on SSR‐SRAP marker data showed that the 112 accessions were also divided into six groups, similar to the phenotypic trait results. The results of PCA and population structure analysis based on SSR‐SRAP data also produced similar results. In addition, the analysis of the origin and genetic relationship of yam indicated that the species D. opposita may have originated from China. These results demonstrate the genetic diversity and distinctness among the widely cultivated species of Chinese yam and provide a theoretical reference for the classification, breeding, germplasm innovation, utilization, and variety protection of Chinese yam resources.     

Genetic variation and bidirectional gene flow in the riparian plant Miscanthus lutarioriparius,across its endemic range: implications for adaptive potential

Miscanthus lutarioriparius is an endemic species that grows along the middle and lower reaches of the Yangtze River and is a valuable source of germplasm for the development of second‐generation energy crops. The plant that propagates via seeds, stem nodes, and rhizomes shows high phenotypic variation and strong local adaptation. Here, we examined the magnitude and spatial distribution of genetic variation in M. lutarioriparius across its entire distributional range and tested underlying factors that shaped its genetic variation. Population genetic analyses were conducted on 644 individuals from 25 populations using 16 microsatellite markers. M. lutarioriparius exhibited a high level of genetic variation (H_E = 0.682–0.786; A_r = 4.74–8.06) and a low differentiation (F_ST = 0.063; D_est = 0.153). Of the total genetic variation, 10% was attributed to the differences among populations (df = 24, P < 0.0001), whereas 90% was attributed to the differences among individuals (df = 619, P ≤ 0.0001). Genetic diversity did not differ significantly across longitudes and did not increase in the populations growing downstream of the Yangtze River. However, significant associations were found between genetic differentiation and spatial distance. Six genetic discontinuities were identified, which mostly distributed among downstream populations. We conclude that anthropogenic factors and landscape features both contributed to shaping the pattern of gene flow in M. lutarioriparius, including long‐distance bidirectional dispersal. Our results explain the genetic basis of the high degree of adaptability in M. lutarioriparius and identify potential sources of new germplasm for the domestication of this potential second‐generation energy crop.     

Regional and local patterns of genetic variation and structure in yellow‐necked mice ‐ the roles of geographic distance,population abundance,and winter severity

The goal of this study, conducted in seven large woodlands and three areas with small woodlots in northeastern Poland in 2004–2008, was to infer genetic structure in yellow‐necked mouse Apodemus flavicollis population and to evaluate the roles of environmental and population ecology variables in shaping the spatial pattern of genetic variation using 768 samples genotyped at 13 microsatellite loci. Genetic variation was very high in all studied regions. The primal genetic subdivision was observed between the northern and the southern parts of the study area, which harbored two major clusters and the intermediate area of highly admixed individuals. The probability of assignment of individual mice to the northern cluster increased significantly with lower temperatures of January and July and declined in regions with higher proportion of deciduous and mixed forests. Despite the detected structure, genetic differentiation among regions was very low. Fine‐scale structure was shaped by the population density, whereas higher level structure was mainly shaped by geographic distance. Genetic similarity indices were highly influenced by mouse abundance (which positively correlated with the share of deciduous forests in the studied regions) and exhibited the greatest change between 0 and 1 km in the forests, 0 and 5 km in small woodlots. Isolation by distance pattern, calculated among regions, was highly significant but such relationship between genetic and geographic distance was much weaker, and held the linearity at very fine scale (~1.5 km), when analyses were conducted at individual level.     

Genetic diversity and phylogeography of Daphnia similoides sinensis located in the middle and lower reaches of the Yangtze River

Geographical patterns, climate, and environmental change have important influences on the distribution and spread of aquatic organisms. However, the relationships between the geographical pattern and phylogenetics of Daphnia as well as environmental change are not well known. The genetic diversity and phylogeography of seven D. similoides sinensis populations located in the middle and lower reaches of the Yangtze River were investigated based on the combination of mitochondrial (COI gene) and nuclear (14 microsatellite primers) markers. Based on the mitochondrial gene markers, D. similoides sinensis from the middle and lower reaches of the Yangtze River had one ancestral haplotype and two evolutionary clades. In addition, D. similoides sinensis population deviated from neutral evolution, showing signs of a bottleneck effect followed by population expansion. Based on the microsatellite markers, the seven D. similoides sinensis populations formed three main groups. The dendrogram (NJ/ME) showed that D. similoides sinensis based on the mitochondrial genes marker were obviously clustered two main clades, whereas there were three clades based on the microsatellite markers. Our results suggested that the habitat fragmentation due to the barrier of the dams and sluices promoted the genetic differentiation and phylogeography of D. similoides sinensis populations in the middle and lower reaches of the Yangtze River.     

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337‐species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed‐mating, non‐woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST. Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.     

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high‐amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model‐averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.     

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P. a. axillaris, P. a. parodii and P. a. subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P. axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100 000 years, which is compatible with a divergence time between 35 000 and 107 000 years ago between P. a. axillaris and P. a. parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P. a. axillaris and P. a. parodii in their contact zone.     

No evidence of divergence at neutral genetic markers between the two morphologically different subspecies of the most numerous Arctic seabird

Identifying natural populations that might be considered separate units using morphology, genotype or both is important in understanding the process of speciation and for conservation. We examined the relationships between the only two subspecies of the most numerous Arctic seabird, the Little Auk Alle alle, using both morphological (wing and head‐bill lengths) and genetic data (482 base pairs of the mitochondrial control region and seven nuclear microsatellite loci). We found significant morphological differences between the subspecies, A. a. polaris being significantly larger than the nominate A. a. alle. However, we did not find the subspecies to be differentiated at either mitochondrial DNA or at microsatellite loci. Consequently, one evolutionary significant unit is proposed. The similarity of the two subspecies at neutral genetic markers may be due to contemporary gene flow between populations, as well as large population sizes both in the present and in the past, combined with recent post‐glacial colonization of the Artic.     

Geographical variation of genetic and phenotypic traits in the Mexican sailfin mollies, Poecilia velifera and P. petenensis

Comparing the patterns of population divergence using both neutral genetic and phenotypic traits provides an opportunity to examine the relative importance of evolutionary mechanisms in shaping population differences. We used microsatellite markers to examine population genetic structure in the Mexican sailfin mollies Poecilia velifera and P. petenensis. We compared patterns of genetic structure and divergence to that in two types of phenotypic traits: morphological characters and mating behaviours. Populations within each species were genetically distinct, and conformed to a model of isolation by distance, with populations within different geographical regions being more genetically similar to one another than were populations from different regions. Bayesian clustering and barrier analyses provided additional support for population separation, especially between geographical regions. In contrast, none of the phenotypic traits showed any type of geographical pattern, and population divergence in these traits was uncorrelated with that found in neutral markers. There was also a weaker pattern of regional differences among geographical regions compared to neutral genetic divergence. These results suggest that while divergence in neutral traits is likely a product of population history and genetic drift, phenotypic divergence is governed by different mechanisms, such as natural and sexual selection, and arises at spatial scales independent from those of neutral markers.     

Evidence for adaptation from standing genetic variation on an antimicrobial peptide gene in the mussel Mytilus edulis

Genome scans of population differentiation identify candidate loci for adaptation but provide little information on how selection has influenced the genetic structure of these loci. Following a genome scan, we investigated the nature of the selection responsible for the outlying differentiation observed between populations of the marine mussel Mytilus edulis at a leucine/arginine polymorphism (L31R) in the antimicrobial peptide MGD2. We analysed DNA sequence polymorphisms, allele frequencies and population differentiation of polymorphisms closely linked to L31R, and pairwise and third‐order linkage disequilibria. An outlying level of population differentiation was observed at L31R only, while no departure from panmixia was observed at linked loci surrounding L31R, as in most of the genome. Selection therefore seems to affect L31R directly. Three hypotheses can explain the lack of differentiation in the chromosomal region close to L31R: (i) hitchhiking has occurred but migration and recombination subsequently erased the signal, (ii) selection was weak enough and recombination strong enough to limit the hitchhiking effect to a very small chromosomal region or (iii) selection acted on a pre‐existing polymorphism (i.e. standing variation) at linkage equilibrium with its background. Linkage equilibrium was observed between L31R and linked polymorphisms in every population analysed, as expected under the three hypotheses. However, linkage disequilibrium was observed in some populations between pairs of loci located upstream and downstream to L31R, generating a complex pattern of third‐order linkage disequilibria which is best explained by the hypothesis of selection on a pre‐existing polymorphism. We hypothesise that selection could be either balanced, maintaining alleles at different frequencies depending on the pathogen community encountered locally by mussels, or intermittent, resulting in sporadic fluctuations in allele frequency.     

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black‐spotted pond frog (Pelophylax nigromaculatus), an IUCN‐listed species, on 24 land‐bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen‐binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen‐driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host–pathogen associations in conservation planning and management.     

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species’ evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.     

Balancing selection and genetic drift create unusual patterns of MHCIIβ variation in Galápagos mockingbirds

The extracellular subunit of the major histocompatibility complex MHCIIβ plays an important role in the recognition of pathogens and the initiation of the adaptive immune response of vertebrates. It is widely accepted that pathogen‐mediated selection in combination with neutral micro‐evolutionary forces (e.g. genetic drift) shape the diversity of MHCIIβ, but it has proved difficult to determine the relative effects of these forces. We evaluated the effect of genetic drift and balancing selection on MHCIIβ diversity in 12 small populations of Galápagos mockingbirds belonging to four different species, and one larger population of the Northern mockingbird from the continental USA. After genotyping MHCIIβ loci by high‐throughput sequencing, we applied a correlational approach to explore the relationships between MHCIIβ diversity and population size by proxy of island size. As expected when drift predominates, we found a positive effect of population size on the number of MHCIIβ alleles present in a population. However, the number of MHCIIβ alleles per individual and number of supertypes were not correlated with population size. This discrepancy points to an interesting feature of MHCIIβ diversity dynamics: some levels of diversity might be shaped by genetic drift while others are independent and possibly maintained by balancing selection.