A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.     

Arabidopsis TTG2 Regulates TRY Expression through Enhancement of Activator Complex-Triggered Activation

Trichome patterning in Arabidopsis thaliana is regulated by a regulatory feedback loop of the trichome promoting factors TRANSPARENT TESTA GLABRA1 (TTG1), GLABRA3 (GL3)/ENHANCER OF GL3 (EGL3), and GL1 and a group of homologous R3MYB proteins that act as their inhibitors. Together, they regulate the temporal and spatial expression of GL2 and TTG2, which are considered to control trichome cell differentiation. In this work, we show that TTG2 is a specific activator of TRY (but not CPC or GL2). The WRKY protein TTG2 binds to W-boxes in a minimal promoter fragment of TRY, and these W-boxes are essential for rescue of the try mutant phenotype. We further show that TTG2 alone is not able to activate TRY expression, but rather drastically enhances the activation by TTG1 and GL3. As TTG2 physically interacts with TTG1 and because TTG2 can associate with GL3 through its interaction with TTG1, we propose that TTG2 enhances the activity of TTG1 and GL3 by forming a protein complex.     

High‐resolution computational imaging of leaf hair patterning using polarized light microscopy

The leaf hairs (trichomes) on the aerial surface of many plant species play important roles in phytochemical production and herbivore protection, and have significant applications in the chemical and agricultural industries. Trichome formation in the model plant Arabidopsis thaliana also presents a tractable experimental system to study cell differentiation and pattern formation in plants and animals. Studies of this developmental process suggest that trichome positioning may be the result of a self‐forming pattern, emerging from a lateral inhibition mechanism determined by a network of regulatory factors. Critical to the continued success of these studies is the ability to quantitatively characterize trichome pattern phenotypes in response to mutations in the genes that regulate this process. Advanced protocols for the observation of changes in trichome patterns can be expensive and/or time consuming, and lack user‐friendly analysis tools. In order to address some of these challenges, we describe here a strategy based on polarized light microscopy for the quick and accurate measurement of trichome positions, and provide an online tool designed for the quantitative analyses of trichome number, density and patterning.     

Sporulation Genes Associated with Sporulation Efficiency in Natural Isolates of Yeast

Yeast sporulation efficiency is a quantitative trait and is known to vary among experimental populations and natural isolates. Some studies have uncovered the genetic basis of this variation and have identified the role of sporulation genes (IME1, RME1) and sporulation-associated genes (FKH2, PMS1, RAS2, RSF1, SWS2), as well as non-sporulation pathway genes (MKT1, TAO3) in maintaining this variation. However, these studies have been done mostly in experimental populations. Sporulation is a response to nutrient deprivation. Unlike laboratory strains, natural isolates have likely undergone multiple selections for quick adaptation to varying nutrient conditions. As a result, sporulation efficiency in natural isolates may have different genetic factors contributing to phenotypic variation. Using Saccharomyces cerevisiae strains in the genetically and environmentally diverse SGRP collection, we have identified genetic loci associated with sporulation efficiency variation in a set of sporulation and sporulation-associated genes. Using two independent methods for association mapping and correcting for population structure biases, our analysis identified two linked clusters containing 4 non-synonymous mutations in genes – HOS4, MCK1, SET3, and SPO74. Five regulatory polymorphisms in five genes such as MLS1 and CDC10 were also identified as putative candidates. Our results provide candidate genes contributing to phenotypic variation in the sporulation efficiency of natural isolates of yeast.     

Natural variation in the genes responsible for maturity loci E1, E2, E3 and E4 in soybean

Background and Aims

The timing of flowering has a direct impact on successful seed production in plants. Flowering of soybean (Glycine max) is controlled by several E loci, and previous studies identified the genes responsible for the flowering loci E1, E2, E3 and E4. However, natural variation in these genes has not been fully elucidated. The aims of this study were the identification of new alleles, establishment of allele diagnoses, examination of allelic combinations for adaptability, and analysis of the integrated effect of these loci on flowering.

Methods

The sequences of these genes and their flanking regions were determined for 39 accessions by primer walking. Systematic discrimination among alleles was performed using DNA markers. Genotypes at the E1–E4 loci were determined for 63 accessions covering several ecological types using DNA markers and sequencing, and flowering times of these accessions at three sowing times were recorded.

Key Results

A new allele with an insertion of a long interspersed nuclear element (LINE) at the promoter of the E1 locus (e1-re) was identified. Insertion and deletion of 36 bases in the eighth intron (E2-in and E2-dl) were observed at the E2 locus. Systematic discrimination among the alleles at the E1–E3 loci was achieved using PCR-based markers. Allelic combinations at the E1–E4 loci were found to be associated with ecological types, and about 62–66 % of variation of flowering time could be attributed to these loci.

Conclusions

The study advances understanding of the combined roles of the E1–E4 loci in flowering and geographic adaptation, and suggests the existence of unidentified genes for flowering in soybean,     

Retrospective genomic analysis of sorghum adaptation to temperate-zone grain production

Background

Sorghum is a tropical C₄ cereal that recently adapted to temperate latitudes and mechanized grain harvest through selection for dwarfism and photoperiod-insensitivity. Quantitative trait loci for these traits have been introgressed from a dwarf temperate donor into hundreds of diverse sorghum landraces to yield the Sorghum Conversion lines. Here, we report the first comprehensive genomic analysis of the molecular changes underlying this adaptation.

Results

We apply genotyping-by-sequencing to 1,160 Sorghum Conversion lines and their exotic progenitors, and map donor introgressions in each Sorghum Conversion line. Many Sorghum Conversion lines carry unexpected haplotypes not found in either presumed parent. Genome-wide mapping of introgression frequencies reveals three genomic regions necessary for temperate adaptation across all Sorghum Conversion lines, containing the Dw1, Dw2, and Dw3 loci on chromosomes 9, 6, and 7 respectively. Association mapping of plant height and flowering time in Sorghum Conversion lines detects significant associations in the Dw1 but not the Dw2 or Dw3 regions. Subpopulation-specific introgression mapping suggests that chromosome 6 contains at least four loci required for temperate adaptation in different sorghum genetic backgrounds. The Dw1 region fractionates into separate quantitative trait loci for plant height and flowering time.

Conclusions

Generating Sorghum Conversion lines has been accompanied by substantial unintended gene flow. Sorghum adaptation to temperate-zone grain production involves a small number of genomic regions, each containing multiple linked loci for plant height and flowering time. Further characterization of these loci will accelerate the adaptation of sorghum and related grasses to new production systems for food and fuel.     

A Scale-Corrected Comparison of Linkage Disequilibrium Levels between Genic and Non-Genic Regions

The understanding of non-random association between loci, termed linkage disequilibrium (LD), plays a central role in genomic research. Since causal mutations are generally not included in genomic marker data, LD between those and available markers is essential for capturing the effects of causal loci on localizing genes responsible for traits. Thus, the interpretation of association studies requires a detailed knowledge of LD patterns. It is well known that most LD measures depend on minor allele frequencies (MAF) of the considered loci and the magnitude of LD is influenced by the physical distances between loci. In the present study, a procedure to compare the LD structure between genomic regions comprising several markers each is suggested. The approach accounts for different scaling factors, namely the distribution of MAF, the distribution of pair-wise differences in MAF, and the physical extent of compared regions, reflected by the distribution of pair-wise physical distances. In the first step, genomic regions are matched based on similarity in these scaling factors. In the second step, chromosome- and genome-wide significance tests for differences in medians of LD measures in each pair are performed. The proposed framework was applied to test the hypothesis that the average LD is different in genic and non-genic regions. This was tested with a genome-wide approach with data sets for humans (Homo sapiens), a highly selected chicken line (Gallus gallus domesticus) and the model plant Arabidopsis thaliana. In all three data sets we found a significantly higher level of LD in genic regions compared to non-genic regions. About 31% more LD was detected genome-wide in genic compared to non-genic regions in Arabidopsis thaliana, followed by 13.6% in human and 6% chicken. Chromosome-wide comparison discovered significant differences on all 5 chromosomes in Arabidopsis thaliana and on one third of the human and of the chicken chromosomes.     

Ln Is a Key Regulator of Leaflet Shape and Number of Seeds per Pod in Soybean

Narrow leaflet soybean (Glycine max) varieties tend to have more seeds per pod than broad leaflet varieties. Narrow leaflet in soybean is conferred by a single recessive gene, ln. Here, we show that the transition from broad (Ln) to narrow leaflet (ln) is associated with an amino acid substitution in the EAR motif encoded by a gene (designated Gm-JAGGED1) homologous to Arabidopsis JAGGED (JAG) that regulates lateral organ development and the variant exerts a pleiotropic effect on fruit patterning. The genomic region that regulates both the traits was mapped to a 12.6-kb region containing only one gene, Gm-JAG1. Introducing the Gm-JAG1 allele into a loss-of-function Arabidopsis jagged mutant partially restored the wild-type JAG phenotypes, including leaf shape, flower opening, and fruit shape, but the Gm-jag1 (ln) and EAR-deleted Gm-JAG1 alleles in the jagged mutant did not result in an apparent phenotypic change. These observations indicate that despite some degree of functional change of Gm-JAG1 due to the divergence from Arabidopsis JAG, Gm-JAG1 complemented the functions of JAG in Arabidopsis thaliana. However, the Gm-JAG1 homoeolog, Gm-JAG2, appears to be sub- or neofunctionalized, as revealed by the differential expression of the two genes in multiple plant tissues, a complementation test, and an allelic analysis at both loci.     

A Population Genomics Study of the Arabidopsis Core Cell Cycle Genes Shows the Signature of Natural Selection

Large-scale comparison of sequence polymorphism and divergence at numerous genomic loci within and between closely related species can reveal signatures of natural selection. Here, we present a population genomics study based on direct sequencing of 61 mitotic cell cycle genes from 30 Arabidopsis thaliana accessions and comparison of the resulting data to the close relative Arabidopsis lyrata. We found that the Arabidopsis core cell cycle (CCC) machinery is not highly constrained but is subject to different modes of selection. We found patterns of purifying selection for the cyclin-dependent kinase (CDK), CDK subunit, retinoblastoma, and WEE1 gene families. Other CCC gene families often showed a mix of one or two constrained genes and relaxed purifying selection on the other genes. We found several large effect mutations in CDKB1;2 that segregate in the species. We found a strong signature of adaptive protein evolution in the Kip-related protein KRP6 and departures from equilibrium at CDKD;1 and CYCA3;3 consistent with the operation of selection in these gene regions. Our data suggest that within Arabidopsis, the genetic robustness of cell cycle–related processes is more due to functional redundancy than high selective constraint.     

Mutational Analysis of Natural Alleles at the B Incompatibility Factor of SCHIZOPHYLLUM COMMUNE: alpha2 and beta6

The B incompatibility factor of the fungus Schizophyllum commune having allelic specificity α2–β6 was subjected to mutagenesis by X-irradiation. Five types of mutations were recovered, four of them new types previously unreported. All lead to loss of the B-factor regulatory function; three of the five have retained their allelic specificity. The mutations map in three closely linked sites: Bα, Bβ, and between Bα and Bβ.