Leaf shape evolution has a similar genetic architecture in three edaphic specialists within the Mimulus guttatus species complex

Local adaptation to an altitudinal gradient: the interplay between mean phenotypic trait variation and phenotypic plasticity in Mimulus laciniatus

Organisms can adapt to environmental heterogeneity through two mechanisms: (1) expression of population genetic variation or (2) phenotypic plasticity. In this study we investigated whether patterns of variation in both trait means and phenotypic plasticity along elevational and latitudinal clines in a North American endemic plant, Mimulus laciniatus, were consistent with local adaptation. We grew inbred lines of M. laciniatus from across the species' range in two common gardens varying in day length to measure mean and plastic trait expression in several traits previously shown to be involved in adaptation to M. laciniatus's rocky outcrop microhabitat: flowering time, size-related traits, and leaf shape. We examined correlations between the mean phenotype and phenotypic plasticity, and tested for a relationship between trait variation and population elevation and latitude. We did not find a strong correlation between mean and plastic trait expression at the individual genotype level suggesting that they operate under independent genetic controls. We identified multiple traits that show patterns consistent with local adaptation to elevation: critical photoperiod, flowering time, flower size, mean leaf lobing, and leaf lobing plasticity. These trends occur along multiple geographically independent altitudinal clines indicating that selection is a more likely cause of this pattern than gene flow among nearby populations with similar trait values. We also found that population variation in mean leaf lobing is associated with latitude. Our results indicate that both having more highly lobed leaves and greater leaf shape plasticity may be adaptive at high elevation within M. laciniatus. Our data strongly suggest that traits known to be under divergent selection between M. laciniatus and close relative Mimulus guttatus are also under locally varying selection within M. laciniatus.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Spatially and temporally varying selection influence species boundaries in two sympatric Mimulus

Spatially and temporally varying selection can maintain genetic variation within and between populations, but it is less well known how these forces influence divergence between closely related species. We identify the interaction of temporal and spatial variation in selection and their role in either reinforcing or eroding divergence between two closely related Mimulus species. Using repeated reciprocal transplant experiments with advanced generation hybrids, we compare the strength of selection on quantitative traits involved in adaptation and reproductive isolation in Mimulus guttatus and Mimulus laciniatus between two years with dramatically different water availability. We found strong divergent habitat-mediated selection on traits in the direction of species differences during a drought in 2013, suggesting that spatially varying selection maintains species divergence. However, a relaxation in divergent selection on most traits in an unusually wet year (2019), including flowering time, which is involved in pre-zygotic isolation, suggests that temporal variation in selection may weaken species differences. Therefore, we find evidence that temporally and spatially varying selection may have opposing roles in mediating species boundaries. Given our changing climate, future growing seasons are expected to be more similar to the dry year, suggesting that in this system climate change may actually increase species divergence.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

Promoter variation in a homeobox gene, CpDll, is associated with deeply lobed leaf in Cucurbita pepo L

CpDll, encoding an HD-Zip I transcription factor, positively regulates formation of deeply lobed leaf shape in zucchini, Cucurbita pepo, which is associated with sequence variation in its promoter region. Leaf shape is an important horticultural trait in zucchini (Cucurbita pepo L.). Deeply lobed leaves have potential advantages for high-density planting and hybrid production. However, little is known about the molecular basis of deeply lobed leaf formation in this important vegetable crop. Here, we conducted QTL analysis and fine mapping of the deeply lobed leaf (CpDll) locus using recombinant inbred lines and large F₂ populations developed from crosses between the deeply lobed leaf HM-S2, and entire leaf Jin-GL parental lines. We show that CpDll exhibited incomplete dominance for the deeply lobed leaf shape in HM-S2. Map-based cloning provided evidence that CpCll encodes a type I homeodomain (HD)- and Leu zipper (Zip) element-containing transcription factor. Sequence analysis between HM-S2 and Jin-GL revealed no sequence variations in the coding sequences, whereas a number of variations were identified in the promoter region between them. DUAL-LUC assays revealed significantly stronger promoter activity in HM-S2 than that in Jin-GL. There was also significantly higher expression of CpDll in the leaf base of deeply lobed leaves of HM-S2 compared with entire leaf Jin-GL. Comparative analysis of CpDll gene homologs in nine cucurbit crop species (family Cucurbitaceae) revealed conservation in both structure and function of this gene in regulation of deeply lobed leaf formation. Our work provides new insights into the molecular basis of leaf lobe formation in pumpkin/squash and other cucurbit crops. This work also facilitates marker-assisted selection for leaf shape in zucchini breeding.

tealeaves: an R package for modelling leaf temperature using energy budgets

Plants must regulate leaf temperature to optimize photosynthesis, control water loss and prevent damage caused by overheating or freezing. Physical models of leaf energy budgets calculate the energy fluxes and leaf temperatures for a given set leaf and environmental parameters. These models can provide deep insight into the variation in leaf form and function, but there are few computational tools available to use these models. Here I introduce a new R package called tealeaves to make complex leaf energy budget models accessible to a broader array of plant scientists. This package enables novice users to start modelling leaf energy budgets quickly while allowing experts to customize their parameter settings. The code is open source, freely available and readily integrates with other R tools for scientific computing. This paper describes the current functionality of tealeaves, but new features will be added in future releases. This software tool will advance new research on leaf thermal physiology to advance our understanding of basic and applied plant science.

Parallel genomic architecture underlies repeated sexual signal divergence in Hawaiian Laupala crickets

When the same phenotype evolves repeatedly, we can explore the predictability of genetic changes underlying phenotypic evolution. Theory suggests that genetic parallelism is less likely when phenotypic changes are governed by many small-effect loci compared to few of major effect, because different combinations of genetic changes can result in the same quantitative outcome. However, some genetic trajectories might be favoured over others, making a shared genetic basis to repeated polygenic evolution more likely. To examine this, we studied the genetics of parallel male mating song evolution in the Hawaiian cricket Laupala. We compared quantitative trait loci (QTL) underlying song divergence in three species pairs varying in phenotypic distance. We tested whether replicated song divergence between species involves the same QTL and whether the likelihood of QTL sharing is related to QTL effect size. Contrary to theoretical predictions, we find substantial parallelism in polygenic genetic architectures underlying repeated song divergence. QTL overlapped more frequently than expected based on simulated QTL analyses. Interestingly, QTL effect size did not predict QTL sharing, but did correlate with magnitude of phenotypic divergence. We highlight potential mechanisms driving these constraints on cricket song evolution and discuss a scenario that consolidates empirical quantitative genetic observations with micro-mutational theory.

Monkeyflowers (Mimulus): new model for plant developmental genetics and evo-devo

Contents Summary 694 I. Introduction 694 II. The system 695 III. Regulation of carotenoid pigmentation 695 IV. Formation of periodic pigmentation patterns 696 V. Developmental genetics of corolla tube formation and elaboration 697 VI. Molecular basis of floral trait variation underlying pollinator shift 698 VII. Outlook 699 Acknowledgements 699 References 699 SUMMARY: Monkeyflowers (Mimulus) have long been recognized as a classic ecological and evolutionary model system. However, only recently has it been realized that this system also holds great promise for studying the developmental genetics and evo-devo of important plant traits that are not found in well-established model systems such as Arabidopsis. Here, I review recent progress in four different areas of plant research enabled by this new model, including transcriptional regulation of carotenoid biosynthesis, formation of periodic pigmentation patterns, developmental genetics of corolla tube formation and elaboration, and the molecular basis of floral trait divergence underlying pollinator shift. These examples suggest that Mimulus offers ample opportunities to make exciting discoveries in plant development and evolution.

Estimation of QTL heritability based on pooled sequencing data

Motivation

Bulked segregant analysis combined with next generation sequencing has proven to be a simple and efficient approach for fast mapping of quantitative trait loci (QTLs). However, how to estimate the proportion of phenotypic variance explained by a QTL (or termed QTL heritability) in such pooled QTL mapping is an unsolved problem.

Results

In this paper, we propose a method called PQHE to estimate QTL heritability using pooled sequencing data obtained under different experimental designs. Simulation studies indicated that our method is correct and feasible. Four practical examples from rice and yeast are demonstrated, each representing a different situation.

Availability and implementation

The R scripts of our method are open source under GPLv3 license at http://genetics.fafu.edu.cn/PQHE or https://github.com/biotangweiqi/PQHE. The R scripts require the R package rootSolve.

Contact

wuwr@fafu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic architecture of quantitative flower and leaf traits in a pair of sympatric sister species of Primulina

Flowers and leaves each represent suites of functionally interrelated traits that are often involved in species divergence and local adaptation. However, a major unresolved issue is how the individual component traits that make up a complex trait such as a flower evolve in a coordinated fashion to retain a high degree of functionality. We use a quantitative trait loci (QTL) approach to elucidate the genetic architecture of divergence in flower and leaf traits between the sister species Primulina depressa and Primulina danxiaensis, which grow sympatrically but in contrasting microhabitats. We found that flower traits were controlled by multiple QTL of small effect, while leaf physiological and morphological traits tended to be controlled by QTL of larger effect. The observed floral integration, manifested by a high degree overlap in both individual trait QTL and QTL for principal component scores (PCA QTL), may have been critical for evolutionary divergence of floral morphology in relation to their pollinators. This overlap suggests that direct selection on only one or a few of the component traits could have caused substantial divergence in other floral traits due to genetic correlations, while the low QTL overlap between floral and vegetative traits suggests that these trait suites are genetically unlinked and can evolve independently in response to different selective pressures corresponding to their distinct functions.

The role of cold cues at different life stages on germination and flowering phenology

PREMISE OF THE STUDY

The timing of major phenological transitions is critical to lifetime fitness, and life history theory predicts differences for annual and perennial plants. To correctly time these transitions, many plants rely on environmental cues such as exposure to extended periods of cold, which may occur at different stages throughout their lifetime.

METHODS

We studied the role of cold at different life stages, by jointly exposing seed (stratification) and rosettes (vernalization) to cold. We used 23 populations of Mimulus guttatus, which vary from annuals to perennials, and investigated how cold at one or both stages affected germination, flowering, growth, and biomass.

KEY RESULTS

We found that stratification and vernalization interact to affect life cycle transitions, and that cold at either stage could synchronize flowering phenology. For perennials, either stratification or vernalization is necessary for maximum flowering. We also found that germination timing covaried with later traits. Moreover, plants from environments with dissimilar climates displayed different phenological responses to stratification or vernalization.

CONCLUSIONS

In general, cold is more important for seed germination in annuals and plants from environments with warm temperatures and variable precipitation. In contrast, cold is more important for flowering in perennials: it accelerates flowering in plants from lower precipitation environments, and it increases flowering proportion in plants from cooler, more stable precipitation environments. We discuss our findings in the context of the variable environments plants experience within a population and the variation encountered across the biogeographic native range of the species.

Genetic mapping of the LOBED LEAF 1 (ClLL1) gene to a 127.6-kb region in watermelon (Citrullus lanatus L.)

The lobed leaf character is a unique morphologic trait in crops, featuring many potential advantages for agricultural productivity. Although the majority of watermelon varieties feature lobed leaves, the genetic factors responsible for lobed leaf formation remain elusive. The F2:3 leaf shape segregating population offers the opportunity to study the underlying mechanism of lobed leaf formation in watermelon. Genetic analysis revealed that a single dominant allele (designated ClLL1) controlled the lobed leaf trait. A large-sized F3:4 population derived from F2:3 individuals was used to map ClLL1. A total of 5,966 reliable SNPs and indels were identified genome-wide via a combination of BSA and RNA-seq. Using the validated SNP and indel markers, the location of ClLL1 was narrowed down to a 127.6-kb region between markers W08314 and W07061, containing 23 putative ORFs. Expression analysis via qRT-PCR revealed differential expression patterns (fold-changes above 2-fold or below 0.5-fold) of three ORFs (ORF3, ORF11, and ORF18) between lobed and non-lobed leaf plants. Based on gene annotation and expression analysis, ORF18 (encoding an uncharacterized protein) and ORF22 (encoding a homeobox-leucine zipper-like protein) were considered as most likely candidate genes. Furthermore, sequence analysis revealed no polymorphisms in cDNA sequences of ORF18; however, two notable deletions were identified in ORF22. This study is the first report to map a leaf shape gene in watermelon and will facilitate cloning and functional characterization of ClLL1 in future studies.

A New Advanced Backcross Tomato Population Enables High Resolution Leaf QTL Mapping and Gene Identification

Quantitative Trait Loci (QTL) mapping is a powerful technique for dissecting the genetic basis of traits and species differences. Established tomato mapping populations between domesticated tomato (Solanum lycopersicum) and its more distant interfertile relatives typically follow a near isogenic line (NIL) design, such as the S. pennellii Introgression Line (IL) population, with a single wild introgression per line in an otherwise domesticated genetic background. Here, we report on a new advanced backcross QTL mapping resource for tomato, derived from a cross between the M82 tomato cultivar and S. pennellii This so-called Backcrossed Inbred Line (BIL) population is comprised of a mix of BC2 and BC3 lines, with domesticated tomato as the recurrent parent. The BIL population is complementary to the existing S. pennellii IL population, with which it shares parents. Using the BILs, we mapped traits for leaf complexity, leaflet shape, and flowering time. We demonstrate the utility of the BILs for fine-mapping QTL, particularly QTL initially mapped in the ILs, by fine-mapping several QTL to single or few candidate genes. Moreover, we confirm the value of a backcrossed population with multiple introgressions per line, such as the BILs, for epistatic QTL mapping. Our work was further enabled by the development of our own statistical inference and visualization tools, namely a heterogeneous hidden Markov model for genotyping the lines, and by using state-of-the-art sparse regression techniques for QTL mapping.

The genetic architecture of local adaptation and reproductive isolation in sympatry within the Mimulus guttatus species complex

The genetic architecture of local adaptation has been of central interest to evolutionary biologists since the modern synthesis. In addition to classic theory on the effect size of adaptive mutations by Fisher, Kimura and Orr, recent theory addresses the genetic architecture of local adaptation in the face of ongoing gene flow. This theory predicts that with substantial gene flow between populations local adaptation should proceed primarily through mutations of large effect or tightly linked clusters of smaller effect loci. In this study, we investigate the genetic architecture of divergence in flowering time, mating system-related traits, and leaf shape between Mimulus laciniatus and a sympatric population of its close relative M. guttatus. These three traits are probably involved in M. laciniatus' adaptation to a dry, exposed granite outcrop environment. Flowering time and mating system differences are also reproductive isolating barriers making them 'magic traits'. Phenotypic hybrids in this population provide evidence of recent gene flow. Using next-generation sequencing, we generate dense SNP markers across the genome and map quantitative trait loci (QTLs) involved in flowering time, flower size and leaf shape. We find that interspecific divergence in all three traits is due to few QTL of large effect including a highly pleiotropic QTL on chromosome 8. This QTL region contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in other Mimulus species. Our results are consistent with theory, indicating that local adaptation and reproductive isolation with gene flow should be due to few loci with large and pleiotropic effects.

Breaking the mold: understanding the evolution and development of lateral organs in diverse plant models

The formation of complex three-dimensional shape differs significantly between plants and animals due to the presence of the cell wall in the former, which prevents all cell migration. Instead, in lateral plant organs such as leaves or petals, shape is controlled by a series of developmental phases in which the organ acquires polarity, cells undergo proliferation, and, lastly, cells expand to their final shape and size. Although these processes were first described based on mutagenesis approaches in major model systems like Arabidopsis thaliana, further insight into their complexity is best provided by studies of natural variation in organ shape in alternative model systems that sample a broader range of plant form. Weaving together work from both forward and evolutionary genetics, this review focuses on how modification in polarity establishment, cell proliferation and cell expansion drives modifications in the fundamental lateral organ developmental program to create diversity in shape.

Accessibility, constraint, and repetition in adaptive floral evolution

Adaptive phenotypic evolution is shaped by natural selection on multiple organismal traits as well as by genetic correlations among traits. Genetic correlations can arise through pleiotropy and can bias the production of phenotypic variation to certain combinations of traits. This phenomenon is referred to as developmental bias or constraint. Developmental bias may accelerate or constrain phenotypic evolution, depending on whether selection acts parallel or in opposition to genetic correlations among traits. We discuss examples from floral evolution where genetic correlations among floral traits contribute to rapid, coordinated evolution in multiple floral organ phenotypes and suggest future research directions that will explore the relationship between the genetic basis of adaptation and the pre-existing structure of genetic correlations. On the other hand, natural selection may act perpendicular to a strong genetic correlation, for example when two traits are encoded by a subset of the same genes and natural selection favors change in one trait and stability in the second trait. In such cases, adaptation is constrained by the availability of genetic variation that can influence the focal trait with minimal pleiotropic effects. Examples from plant diversification suggest that the origin of certain adaptations depends on the prior evolution of a gene copy with reduced pleiotropic effects, generated through the process of gene duplication followed by subfunctionalization or neofunctionalization. A history of gene duplication in some developmental pathways appears to have allowed particular flowering plant linages to have repeatedly evolved adaptations that might otherwise have been developmentally constrained.

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep-bodied benthic form adapted to forage near the lake substrate, and a narrow-bodied limnetic form adapted to forage in open water. Here, we use genome-wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.