Cheng-Ruei Lee

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Summary Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

3D phenotyping and quantitative trait locus mapping identify core regions of the rice genome controlling root architecture

Significance Improving the efficiency of root systems should result in crop varieties with better yields, requiring fewer chemical inputs, and that can grow in harsher environments. Little is known about the genetic factors that condition root growth because of roots’ complex shapes, the opacity of soil, and environmental influences. We designed a 3D root imaging and analysis platform and used it to identify regions of the rice genome that control several different aspects of root system growth. The results of this study should inform future efforts to enhance root architecture for agricultural benefit. Identification of genes that control root system architecture in crop plants requires innovations that enable high-throughput and accurate measurements of root system architecture through time. We demonstrate the ability of a semiautomated 3D in vivo imaging and digital phenotyping pipeline to interrogate the quantitative genetic basis of root system growth in a rice biparental mapping population, Bala × Azucena. We phenotyped >1,400 3D root models and >57,000 2D images for a suite of 25 traits that quantified the distribution, shape, extent of exploration, and the intrinsic size of root networks at days 12, 14, and 16 of growth in a gellan gum medium. From these data we identified 89 quantitative trait loci, some of which correspond to those found previously in soil-grown plants, and provide evidence for genetic tradeoffs in root growth allocations, such as between the extent and thoroughness of exploration. We also developed a multivariate method for generating and mapping central root architecture phenotypes and used it to identify five major quantitative trait loci (r2 = 24–37%), two of which were not identified by our univariate analysis. Our imaging and analytical platform provides a means to identify genes with high potential for improving root traits and agronomic qualities of crops.

Quantifying effects of environmental and geographical factors on patterns of genetic differentiation

Elucidating the factors influencing genetic differentiation is an important task in biology, and the relative contribution from natural selection and genetic drift has long been debated. In this study, we used a regression‐based approach to simultaneously estimate the quantitative contributions of environmental adaptation and isolation by distance on genetic variation in Boechera stricta, a wild relative of Arabidopsis. Patterns of discrete and continuous genetic differentiation coexist within this species. For the discrete differentiation between two major genetic groups, environment has larger contribution than geography, and we also identified a significant environment‐by‐geography interaction effect. Elsewhere in the species range, we found a latitudinal cline of genetic variation reflecting only isolation by distance. To further confirm the effect of environmental selection on genetic divergence, we identified the specific environmental variables predicting local genotypes in allopatric and sympatric regions. Water availability was identified as the possible cause of differential local adaptation in both geographical regions, confirming the role of environmental adaptation in driving and maintaining genetic differentiation between the two major genetic groups. In addition, the environment‐by‐geography interaction is further confirmed by the finding that water availability is represented by different environmental factors in the allopatric and sympatric regions. In conclusion, this study shows that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, which only produced a gradual, clinal pattern of genetic variation. These findings emphasize the importance of environmental selection in shaping patterns of species‐wide genetic variation in the natural environment.

Genetic trade-offs and conditional neutrality contribute to local adaptation

Divergent natural selection promotes local adaptation and can lead to reproductive isolation of populations in contrasting environments; however, the genetic basis of local adaptation remains largely unresolved in natural populations. Local adaptation might result from antagonistic pleiotropy, where alternate alleles are favoured in distinct habitats, and polymorphism is maintained by selection. Alternatively, under conditional neutrality some alleles may be favoured in one environment but neutral at other locations. Antagonistic pleiotropy maintains genetic variation across the landscape; however, there is a systematic bias against discovery of antagonistic pleiotropy because the fitness benefits of local alleles need to be significant in at least two environments. Here, we develop a generally applicable method to investigate polygenic local adaptation and identify loci that are the targets of selection. This approach evaluates allele frequency changes after selection at loci across the genome to distinguish antagonistic pleiotropy from conditional neutrality and deleterious variation. We investigate local adaptation at the qualitative trait loci (QTL) level in field experiments, in which we expose 177 F6 recombinant inbred lines and parental lines of Boechera stricta (Brassicaceae) to their parental environments over two seasons. We demonstrate polygenic selection for native alleles in both environments, with 2.8% of the genome exhibiting antagonistic pleiotropy and 8% displaying conditional neutrality. Our study strongly supports antagonistic pleiotropy at one large‐effect flowering phenology QTL (nFT): native homozygotes had significantly greater probabilities of flowering than foreign homozygotes in both parental environments. Such large‐scale field studies are essential to elucidate the genetic basis of adaptation in natural populations.

A Gain-of-Function Polymorphism Controlling Complex Traits and Fitness in Nature

Natural Selection at Work Catching the evolution of a novel function and determining its selective parameters in nature remains an extremely difficult task. Prasad et al. (p. 1081) have undertaken this quest documenting the molecular basis of a natural allelic polymorphism and its effects on herbivory and survival in the Arabidopsis relative, Boechera stricta, living in the Rocky Mountains. Positive selection for a mutation that enhances resistance to herbivory in the model plant Boechera is described. Identification of the causal genes that control complex trait variation remains challenging, limiting our appreciation of the evolutionary processes that influence polymorphisms in nature. We cloned a quantitative trait locus that controls plant defensive chemistry, damage by insect herbivores, survival, and reproduction in the natural environments where this polymorphism evolved. These ecological effects are driven by duplications in the BCMA (branched-chain methionine allocation) loci controlling this variation and by two selectively favored amino acid changes in the glucosinolate-biosynthetic cytochrome P450 proteins that they encode. These changes cause a gain of novel enzyme function, modulated by allelic differences in catalytic rate and gene copy number. Ecological interactions in diverse environments likely contribute to the widespread polymorphism of this biochemical function.

LIFE-HISTORY QTLS AND NATURAL SELECTION ON FLOWERING TIME IN BOECHERA STRICTA, A PERENNIAL RELATIVE OF ARABIDOPSIS

Plants must precisely time flowering to capitalize on favorable conditions. Although we know a great deal about the genetic basis of flowering phenology in model species under controlled conditions, the genetic architecture of this ecologically important trait is poorly understood in nonmodel organisms. Here, we evaluated the transition from vegetative growth to flowering in Boechera stricta, a perennial relative of Arabidopsis thaliana. We examined flowering time QTLs using 7920 recombinant inbred individuals, across seven laboratory and field environments differing in vernalization, temperature, and photoperiod. Genetic and environmental factors strongly influenced the transition to reproduction. We found directional selection for earlier flowering in the field. In the growth chamber experiment, longer winters accelerated flowering, whereas elevated ambient temperatures delayed flowering. Our analyses identified one large effect QTL (nFT), which influenced flowering time in the laboratory and the probability of flowering in the field. In Montana, homozygotes for the native allele at nFT showed a selective advantage of 6.6%. Nevertheless, we found relatively low correlations between flowering times in the field and the growth chambers. Additionally, we detected flowering‐related QTLs in the field that were absent across the full range of laboratory conditions, thus emphasizing the need to conduct experiments in natural environments.

Boechera, a model system for ecological genomics

The selection and development of a study system for evolutionary and ecological functional genomics (EEFG) depend on a variety of factors. Here, we present the genus Boechera as an exemplary system with which to address ecological and evolutionary questions. Our focus on Boechera is based on several characteristics as follows: (i) native populations in undisturbed habitats where current environments reflect historical conditions over several thousand years; (ii) functional genomics benefitting from its close relationship to Arabidopsis thaliana; (iii) inbreeding tolerance enabling development of recombinant inbred lines, near‐isogenic lines and positional cloning; (iv) interspecific crosses permitting mapping for genetic analysis of speciation; (v) apomixis (asexual reproduction by seeds) in a genetically tractable diploid; and (vi) broad geographic distribution in North America, permitting ecological genetics for a large research community. These characteristics, along with the current sequencing of three Boechera species by the Joint Genome Institute, position Boechera as a rapidly advancing system for EEFG studies.

Genome diversity in Brachypodium distachyon: deep sequencing of highly diverse inbred lines

Brachypodium distachyon is small annual grass that has been adopted as a model for the grasses. Its small genome, high-quality reference genome, large germplasm collection, and selfing nature make it an excellent subject for studies of natural variation. We sequenced six divergent lines to identify a comprehensive set of polymorphisms and analyze their distribution and concordance with gene expression. Multiple methods and controls were utilized to identify polymorphisms and validate their quality. mRNA-Seq experiments under control and simulated drought-stress conditions, identified 300 genes with a genotype-dependent treatment response. We showed that large-scale sequence variants had extremely high concordance with altered expression of hundreds of genes, including many with genotype-dependent treatment responses. We generated a deep mRNA-Seq dataset for the most divergent line and created a de novo transcriptome assembly. This led to the discovery of >2400 previously unannotated transcripts and hundreds of genes not present in the reference genome. We built a public database for visualization and investigation of sequence variants among these widely used inbred lines.

Evolution of flux control in the glucosinolate pathway in Arabidopsis thaliana

Network characteristics of biochemical pathways are believed to influence the rate of evolutionary change in constituent enzymes. One characteristic that may affect rate heterogeneity is control of the amount of product produced by a biochemical pathway or flux control. In particular, theoretical analyses suggest that adaptive substitutions should be concentrated in the enzyme(s) that exert the greatest control over flux. Although a handful of studies have found a correlation between position in a pathway and evolutionary rate, these investigations have not examined the relationship between evolutionary rate and flux control. Given that genes with greater control will experience stronger selection and that the probability of fixation is proportional to the selective advantage, we ask the following: 1) do upstream enzymes have majority flux control, 2) do enzymes with majority flux control accumulate adaptive substitutions, and 3) are upstream enzymes under higher selective constraint? First, by perturbing the enzymes in the aliphatic glucosinolate pathway in Arabidopsis thaliana with gene insertion lines, we show that flux control is focused in the first enzyme in the pathway. Next, by analyzing several sequence signatures of selection, we also show that this enzyme is the only one in the pathway that shows convincing evidence of selection. Our results support the hypothesis that natural selection preferentially acts on enzymes with high flux control.

Molecular Evolution and Functional Diversification of Fatty Acid Desaturases after Recurrent Gene Duplication in Drosophila

Frequent gene duplications in the genome incessantly supply new genetic materials for functional innovation presumably driven by positive Darwinian selection. This mechanism in the desaturase gene family has been proposed to be important in triggering the pheromonal diversification in insects. With the recent completion of a dozen Drosophila genomes, a genome-wide perspective is possible. In this study, we first identified homologs of desaturase genes in 12 Drosophila species and noted that while gene duplication events are relatively frequent, gene losses are not scarce, especially in the desat1–desat2–desatF clade. By reconciling the gene tree with species phylogeny and the chromosomal synteny of the sequenced Drosophila genomes, at least one gene loss in desat2 and a minimum of six gene gains (resulting in seven desatF homologs, α-η), three gene losses and one relocation in desatF were inferred. Upon branching off the ancestral desat1 lineage, both desat2 and desatF gained novel functions through accelerating protein evolution. The amino acid residues under positive selection located near the catalytic sites and the C-terminal region might be responsible for altered substrate selectivity between closely related species. The association between the expression pattern of desatF-α and the chemical composition of cuticular hydrocarbons implies that the ancestral function of desatF-α is the second desaturation at the four carbons after the first double bond in diene synthesis, and the shift from bisexual to female-specific expression in desatF-α occurred in the ancestral lineage of Drosophila melanogaster subgroup. A relationship between the number of expressed desatF homologs and the diene diversification has also been observed. These results suggest that the molecular diversification of fatty acid desaturases after recurrent gene duplication plays an important role in pheromonal diversity in Drosophila.