What Has Natural Variation Taught Us about Plant Development,Physiology, and Adaptation?

Nearly 100 genes and functional polymorphisms underlying natural variation in plant development and physiology have been identified. In crop plants, these include genes involved in domestication traits, such as those related to plant architecture, fruit and seed structure and morphology, as well as yield and quality traits improved by subsequent crop breeding. In wild plants, comparable traits have been dissected mainly in Arabidopsis thaliana. In this review, we discuss the major contributions of the analysis of natural variation to our understanding of plant development and physiology, focusing in particular on the timing of germination and flowering, plant growth and morphology, primary metabolism, and mineral accumulation. Overall, functional polymorphisms appear in all types of genes and gene regions, and they may have multiple mutational causes. However, understanding this diversity in relation to adaptation and environmental variation is a challenge for which tools are now available.     

Beyond Arabidopsis. Translational biology meets evolutionary developmental biology

Developmental processes shape plant morphologies, which constitute important adaptive traits selected for during evolution. Identifying the genes that act in developmental pathways and determining how they are modified during evolution is the focus of the field of evolutionary developmental biology, or evo-devo. Knowledge of genetic pathways in the plant model Arabidopsis serves as the starting point for investigating how the toolkit of developmental pathways has been used and reused to form different plant body plans. One productive approach is to identify genes in other species that are orthologous to genes known to control developmental pathways in Arabidopsis and then determine what changes have occurred in the protein coding sequence or in the gene's expression to produce an altered morphology. A second approach relies on natural variation among wild populations or crop plants. Natural variation can be exploited to identify quantitative trait loci that underlie important developmental traits and, thus, define those genes that are responsible for adaptive changes. The possibility of applying comparative genomics approaches to Arabidopsis and related species promises profound new insights into the interplay of evolution and development.     

Nucleotide diversity of vernalization and flowering‐time‐related genes in a germplasm collection of meadow fescue (Festuca pratensis Huds. syn. Lolium pratense (Huds.) Darbysh.)

In plant species, control of flowering time is an important factor for adaptation to local natural environments. The Vrn1, CO, FT1 and CK2α genes are key components in the flowering‐specific signaling pathway of grass species. Meadow fescue is an agronomically important forage grass species, which is naturally distributed across Europe and Western Asia. In this study, meadow fescue flowering‐time‐related genes were resequenced to assess nucleotide diversity in European and Western Asian subpopulations. Identified sequence polymorphisms were then converted into PCR‐based molecular genetic markers, and a meadow fescue germplasm collection was genotyped to investigate global allelic variation. Lower nucleotide diversities were observed for the Vrn1 and CO orthologs, while relatively higher values were observed for the FT1 and casein kinase II α‐subunit (CK2α) orthologs. The nucleotide diversity for FT1 orthologs in the Western Asian subpopulation was significantly higher than those of the European subpopulation. Similarly, significant differences in nucleotide diversity for the remaining genes were observed between several combinations of subpopulation. The global allele distribution pattern was consistent with observed level of nucleotide diversity. These results suggested that the degree of purifying selection acting on the genes differs according to geographical location. As previously shown for model plant species, functional specificities of flowering‐time‐related genes may also vary according to environmental conditions.     

How to be early flowering: an evolutionary perspective

In wild and cultivated annual plant species, flowering time is an important life-history trait that coordinates the life cycle with local environmental conditions. Extensive studies on the genetic basis of flowering time in the model species Arabidopsis thaliana have revealed a complex genetic network that can detect environmental and internal signals. Based on this knowledge and on known pleiotropic effects associated with flowering time genes, we suggest that a natural shift towards an early-flowering life cycle might involve only particular functional regions in a limited number of genes. Our predictions are supported by genetic theories of adaptation and by recent data about genes associated with natural variation. We analyse the extent to which these predictions can also apply to crop species.     

Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.     

Lack of genetic diversity across diverse immune genes in an endangered mammal,the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll‐like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome‐level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole‐genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29–220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long‐term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad‐scale immunogenetic diversity analysis in threatened species.     

Molecular population genetics of herbivore-induced protease inhibitor genes in European aspen (Populus tremula L., Salicaceae)

Plants defend themselves against the attack of natural enemies by using an array of both constitutively expressed and induced defenses. Long-lived woody perennials are overrepresented among plant species that show strong induced defense responses, whereas annual plants and crop species are underrepresented. However, most studies of plant defense genes have been performed on annual or short-lived perennial weeds or crop species. Here I use molecular population genetic methods to survey six wound-inducible protease inhibitors (PIs) in a long-lived woody, perennial plant species, the European aspen (Populus tremula), to evaluate the likelihood of either recurrent selective sweeps or balancing selection maintaining amino acid polymorphisms in these genes. The results show that none of the six PI genes have reduced diversities at synonymous sites, as would be expected in the presence of recurrent selective sweeps. However, several genes show some evidence of nonneutral evolution such as enhanced linkage disequilibrium and a large number of high-frequency-derived mutations. A group of at least four Kunitz trypsin inhibitor genes appear to have experienced elevated levels of nonsynonymous substitutions, indicating allelic turnover on an evolutionary timescale. One gene, TI1, has enhanced levels of intraspecific polymorphism at nonsynonymous sites and also has an unusual haplotype structure characterized by two divergent haplotypes occurring at roughly equal frequencies in the sample. One haplotype has very low levels of intraallelic nucleotide diversity, whereas the other haplotype has levels of diversity comparable to other genes in P. tremula. Patterns of sequence diversity at TI1 do not fit a simple model of either balancing selection or recurrent selective sweeps. This suggests that selection at TI1 is more complex, possibly involving allelic cycling.     

The complex genetic architecture of the metabolome

Discovering links between the genotype of an organism and its metabolite levels can increase our understanding of metabolism, its controls, and the indirect effects of metabolism on other quantitative traits. Recent technological advances in both DNA sequencing and metabolite profiling allow the use of broad-spectrum, untargeted metabolite profiling to generate phenotypic data for genome-wide association studies that investigate quantitative genetic control of metabolism within species. We conducted a genome-wide association study of natural variation in plant metabolism using the results of untargeted metabolite analyses performed on a collection of wild Arabidopsis thaliana accessions. Testing 327 metabolites against >200,000 single nucleotide polymorphisms identified numerous genotype-metabolite associations distributed non-randomly within the genome. These clusters of genotype-metabolite associations (hotspots) included regions of the A. thaliana genome previously identified as subject to recent strong positive selection (selective sweeps) and regions showing trans-linkage to these putative sweeps, suggesting that these selective forces have impacted genome-wide control of A. thaliana metabolism. Comparing the metabolic variation detected within this collection of wild accessions to a laboratory-derived population of recombinant inbred lines (derived from two of the accessions used in this study) showed that the higher level of genetic variation present within the wild accessions did not correspond to higher variance in metabolic phenotypes, suggesting that evolutionary constraints limit metabolic variation. While a major goal of genome-wide association studies is to develop catalogues of intraspecific variation, the results of multiple independent experiments performed for this study showed that the genotype-metabolite associations identified are sensitive to environmental fluctuations. Thus, studies of intraspecific variation conducted via genome-wide association will require analyses of genotype by environment interaction. Interestingly, the network structure of metabolite linkages was also sensitive to environmental differences, suggesting that key aspects of network architecture are malleable.     

Levels and patterns of nucleotide variation in domestication QTL regions on rice chromosome 3 suggest lineage-specific selection

Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa--indica and tropical japonica--as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups.     

Genetic diversity, recombination and cryptic clades in Pseudomonas viridiflava infecting natural populations of Arabidopsis thaliana

Species-level genetic diversity and recombination in bacterial pathogens of wild plant populations have been nearly unexplored. Pseudomonas viridiflava is a common natural bacterial pathogen of Arabidopsis thaliana, for which pathogen defense genes and mechanisms are becoming increasing well known. The genetic variation contained within a worldwide sample of P. viridiflava collected from wild populations of A. thaliana was investigated using five genomic sequence fragments totaling 2.3 kb. Two distinct and deeply diverged clades were found within the P. viridiflava sample and in close proximity in multiple populations, each genetically diverse with synonymous variation as high as 9.3% in one of these clades. Within clades, there is evidence of frequent recombination within and between each sequenced locus and little geographic differentiation. Isolates from both clades were also found in a small sample of other herbaceous species in Midwest populations, indicating a possibly broad host range for P. viridiflava. The high levels of genetic variation and recombination together with a lack of geographic differentiation in this pathogen distinguish it from other bacterial plant pathogens for which intraspecific variation has been examined.     

Sequence diversity and haplotype associations with phenotypic responses to crowding: GIGANTEA affects fruit set in Arabidopsis thaliana

Identifying the molecular genetic basis of intraspecific variation in quantitative traits promises to provide novel insight into their evolutionary history as well as genetic mechanisms of adaptation. In an attempt to identify genes responsible for natural variation in competitive responses in Arabidopsis thaliana, we examined DNA sequence diversity at seven loci previously identified as members of the phytochrome B signalling network. For one gene, GIGANTEA (GI), we detected significant haplotype structure. To test for GI haplogroup-phenotype associations, we genotyped 161 A. thaliana accessions at GI and censused the same accessions for total fruit set and the expression of three phenotypic traits (days to flowering, petiole length, and inflorescence height) in a greenhouse experiment where plants were grown in crowded and uncrowded environments. We detected a significant association between GI and total fruit set that resulted in a 14% difference in average fruit set among GI haplogroups. Given that fruit set is an important component of fitness in this species and given the magnitude of the effect, the question arises as to how variation at this locus is maintained. Our observation of frequent and significant epistasis between GI and background single nucleotide polymorphisms (SNP), where the fitness ranking of the GI allele either reverses or does not differ depending on the allele at the interacting SNP, suggests that epistatic selection may actively maintain or at least slow the loss of variation at GI. This result is particularly noteworthy in the light of the ongoing debate regarding the genetic underpinnings of phenotypic evolution and recent observations that epistasis for phenotypic traits and components of fitness is common in A. thaliana.     

Polymorphism of DNA sequences of cryptochrome genes is not associated with the photoperiodic flowering of wild soybean along a latitudinal cline

Both cultivated soybean and its wild relative Glycine soja exhibit strong photoperiodic sensitivity at different latitudes. Recent studies have demonstrated that the blue light-absorbing cryptochrome gene, CRY1a, is involved in the photoperiodic flowering of soybeans. However, no sequence variation was found in the cDNA among cultivars at different latitudinal clines. In the present study, we examined whether positive selection due to polymorphisms in the cryptochrome genes of G. soja occurs. Partial DNA sequences, mainly exons, of cryptochrome genes CRY1a-1d and CRY2a-2c were analyzed for 18 accessions in the Japanese archipelago. The neutral evolutionary pattern of the polymorphisms for all cryptochrome genes except for CRY1a was summarized using Tajima's D test and low nucleotide diversity was shown for all genes. Although CRY1a did not show neutral evolution, balancing selection was recognized in the intron while not in the exon. No geographical pattern of polymorphisms was observed in the cryptochrome genes. These results reject the possibility of cryptochrome genes being involved in the photoperiodic flowering of wild soybeans along a latitudinal cline.     

珍稀濒危飘带兜兰叶绿体全基因组种内变异研究(附录)

飘带兜兰(Paphiopedilum parishii)分布范围狭窄,仅在中国、缅甸、泰国以及老挝有少量分布。近年来,因生境破坏和人为滥采而导致飘带兜兰野生种群极度缩减。为开发种内多态性的分子标记用于保护生物学研究,该研究对飘带兜兰4个野生个体经测序、组装、注释获得的叶绿体基因组序列,与已公布的飘带兜兰2个个体的叶绿体全基因组序列进行比对,分析飘带兜兰叶绿体基因组的种内差异。结果表明:(1)飘带兜兰叶绿体基因组具有典型被子植物叶绿体基因组环状四分体结构,基因组长度为154 403～154 809 bp,共编码129个基因,包括78个蛋白质编码基因、39个tRNA基因、8个rRNA基因,以及4个假基因。(2)在飘带兜兰6个个体叶绿体基因组中检测到103～107个SSRs(simple sequence repeats)位点,其中21个SSR位点具有多态性。此外,在6个个体叶绿体基因组中还检测到60个长序列重复,包括17～21个正向重复、18～29个反向重复、9～16个回文重复、4～9个互补重复。(3)通过比较6个个体叶绿体基因组序列的核苷酸多样性,共发现70处变异,包括10个SNPs(single nucleotide polymorphism)、60个插入缺失(InDels)。其中,有3个SNP位点发生了非同义替换,导致编码功能基因的氨基酸发生改变; 19个插入缺失多态性较高,具有开发为分子标记的潜力。(4)通过计算核苷酸多样性值(P_i)共发现8个有变异的区域,P_i值为0～0.006 32,其中变异度较大的是rps3-rpl22、trnL-UAC-rpl32、rpoB-trnC-GCA以及ycf4,这些高变区可开发为分子标记用于评估飘带兜兰遗传多样性。(5)系统发生分析结果表明,飘带兜兰6个个体叶绿体基因组序列聚在一起,与长瓣兜兰互为姐妹群。综上表明,飘带兜兰叶绿体基因组的SSRs、长序列重复、SNPs、InDels以及核苷酸序列呈现了足够的种内多样性,可开发成分子标记用于该种的系统演化及保护生物学研究。