Experimental estimation of mutation rates in a wheat population with a gene genealogy approach

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

Species diversity and phylogeography of Cornus kousa (Asian dogwood) captured by genomic and genic microsatellites

Cornus kousa (Asian dogwood), an East Asia native tree, is the most economically important species of the dogwood genus, owing to its desirable horticultural traits and ability to hybridize with North America-native dogwoods. To assess the species genetic diversity and to better inform the ongoing and future breeding efforts, we assembled an herbarium and arboretum collection of 131 noncultivated C. kousa specimens. Genotyping and capillary electrophoresis analyses of our C. kousa collection with the newly developed genic and published nuclear genomic microsatellites permitted assessment of genetic diversity and evolutionary history of the species. Regardless of the microsatellite type used, the study yielded generally similar insights into the C. kousa diversity with subtle differences deriving from and underlining the marker used. The accrued evidence pointed to the species distinct genetic pools related to the plant country of origin. This can be helpful in the development of the commercial cultivars for this important ornamental crop with increased pyramided utility traits. Analyses of the C. kousa evolutionary history using the accrued genotyping datasets pointed to an unsampled ancestor population, possibly now extinct, as per the phylogeography of the region. To our knowledge, there are few studies utilizing the same gDNA collection to compare performance of genomic and genic microsatellites. This is the first detailed report on C. kousa species diversity and evolutionary history inference.

Chromosomal structural changes and microsatellite variations in newly synthesized hexaploid wheat mediated by unreduced gametes

Allohexaploid wheat was derived from interspecific hybridization, followed by spontaneous chromosome doubling. Newly synthesized hexaploid wheat by crossing Triticum turgidum and Aegilops tauschii provides a classical model to understand the mechanisms of allohexaploidization in wheat. However, immediate chromosome level variation and microsatellite level variation of newly synthesized hexaploid wheat have been rarely reported. Here, unreduced gametes were applied to develop synthesized hexaploid wheat, NA0928, population by crossing T. turgidum ssp. dicoccum MY3478 and Ae. tauschii SY41, and further S0-S3 generations of NA0928 were assayed by sequential cytological and microsatellite techniques. We demonstrated that plentiful chromosomal structural changes and microsatellite variations emerged in the early generations of newly synthesized hexaploid wheat population NA0928, including aneuploidy with whole-chromosome loss or gain, aneuploidy with telosome formation, chromosome-specific repeated sequence elimination (indicated by fluorescence in situ hybridization) and microsatellite sequence elimination (indicated by sequencing), and many kinds of variations have not been previously reported. Additionally, we reported a new germplasm, T. turgidum accession MY3478 with excellent unreduced gametes trait, and then succeeded to transfer powdery mildew resistance from Ae. tauschii SY41 to synthesized allohexaploid wheat population NA0928, which would be valuable resistance resources for wheat improvement.

Late Glacial Demographic Expansion Motivates a Clock Overhaul for Population Genetics

The molecular clock hypothesis is fundamental in evolutionary biology as by assuming constancy of the molecular rate it provides a timeframe for evolution. However, increasing evidence shows time dependence of inferred molecular rates with inflated values obtained using recent calibrations. As recent demographic calibrations are virtually non-existent in most species, older phylogenetic calibration points (>1 Ma) are commonly used, which overestimate demographic parameters. To obtain more reliable rates of molecular evolution for population studies, I propose the calibration of demographic transition (CDT) method, which uses the timing of climatic changes over the late glacial warming period to calibrate expansions in various species. Simulation approaches and empirical data sets from a diversity of species (from mollusk to humans) confirm that, when compared with other genealogy-based calibration methods, the CDT provides a robust and broadly applicable clock for population genetics. The resulting CDT rates of molecular evolution also confirm rate heterogeneity over time and among taxa. Comparisons of expansion dates with ecological evidence confirm the inaccuracy of phylogenetically derived divergence rates when dating population-level events. The CDT method opens opportunities for addressing issues such as demographic responses to past climate change and the origin of rate heterogeneity related to taxa, genes, time, and genetic information content.

Inference of domestication history and differentiation between early- and late-flowering varieties in pearl millet

Pearl millet (Pennisetum glaucum) is a staple crop in Sahelian Africa. Farmers usually grow varieties with different cycle lengths and complementary functions in Sahelian agrosystems. Both the level of genetic differentiation of these varieties and the domestication history of pearl millet have been poorly studied. We investigated the neutral genetic diversity and population genetic structure of early- and late-flowering domesticated and wild pearl millet populations using 18 microsatellite loci and 8 nucleotide sequences. Strikingly, early- and late-flowering domesticated varieties were not differentiated over their whole distribution area, despite a clear difference in their isolation-by-distance pattern. Conversely, our data brought evidence for two well-differentiated genetic pools in wild pearl millet, allowing us to test scenarios with different numbers and origins of domestication using approximate Bayesian computation (ABC). The ABC analysis showed the likely existence of asymmetric migration between wild and domesticated populations. The model choice procedure indicated that a single domestication from the eastern wild populations was the more likely scenario to explain the polymorphism patterns observed in cultivated pearl millet.

Efficiently tracking selection in a multiparental population: the case of earliness in wheat

Multiparental populations are innovative tools for fine mapping large numbers of loci. Here we explored the application of a wheat Multiparent Advanced Generation Inter-Cross (MAGIC) population for QTL mapping. This population was created by 12 generations of free recombination among 60 founder lines, following modification of the mating system from strict selfing to strict outcrossing using the ms1b nuclear male sterility gene. Available parents and a subset of 380 SSD lines of the resulting MAGIC population were phenotyped for earliness and genotyped with the 9K i-Select SNP array and additional markers in candidate genes controlling heading date. We demonstrated that 12 generations of strict outcrossing rapidly and drastically reduced linkage disequilibrium to very low levels even at short map distances and also greatly reduced the population structure exhibited among the parents. We developed a Bayesian method, based on allelic frequency, to estimate the contribution of each parent in the evolved population. To detect loci under selection and estimate selective pressure, we also developed a new method comparing shifts in allelic frequency between the initial and the evolved populations due to both selection and genetic drift with expectations under drift only. This evolutionary approach allowed us to identify 26 genomic areas under selection. Using association tests between flowering time and polymorphisms, 6 of these genomic areas appeared to carry flowering time QTL, 1 of which corresponds to Ppd-D1, a major gene involved in the photoperiod sensitivity. Frequency shifts at 4 of 6 areas were consistent with earlier flowering of the evolved population relative to the initial population. The use of this new outcrossing wheat population, mixing numerous initial parental lines through multiple generations of panmixia, is discussed in terms of power to detect genes under selection and association mapping. Furthermore we provide new statistical methods for use in future analyses of multiparental populations.

Microsatellite mutation rate during allohexaploidization of newly resynthesized wheat

Simple sequence repeats (SSRs, also known as microsatellites) are known to be mutational hotspots in genomes. DNA rearrangements have also been reported to accompany allopolyploidization. A study of the effect of allopolyploidization on SSR mutation is therefore important for understanding the origin and evolutionary dynamics of SSRs in allopolyploids. Three synthesized double haploid (SynDH) populations were made from 241 interspecific F(1 )haploid hybrids between Triticum turgidum L. and Aegilops tauschii (Coss.) through spontaneous chromosome doubling via unreduced gametes. Mutation events were studied at 160 SSR loci in the S(1) generation (the first generation after chromosome doubling) of the three SynDH populations. Of the 148260 SSR alleles investigated in S(1) generation, only one mutation (changed number of repeats) was confirmed with a mutation rate of 6.74 × 10-6. This mutation most likely occurred in the respective F(1) hybrid. In comparison with previously reported data, our results suggested that allohexaploidization of wheat did not increase SSR mutation rate.

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Effectiveness of in situ and ex situ conservation of crop diversity. What a Phaseolus vulgaris L. landrace case study can tell us

The effectiveness of in situ (on-farm) and ex situ conservation strategies to maintain total genetic diversity was assessed in a threatened Phaseolus vulgaris L. landrace. Farmer seed lots (subpopulations) were sampled initially and then after in situ and ex situ multiplication (two locations). The number of plants used in the ex situ multiplications (120) was much larger than that normally used in germplasm bank procedures and the farmer seed lots were kept separate. In situ, the landrace was multiplied by each farmer with the usual population size. Eighty plants from the initial population, the in situ and the two ex situ multiplications were individually tested using 26 microsatellite markers. Most of the genetic parameters showed a consistent decline in the ex situ populations compared with the in situ population, with a notable loss of less frequent alleles. The differentiation among the farmer subpopulations increased when the multiplication took place outside of the adaptation area. Although 120 plants were multiplied in each ex situ cycle, a bottleneck effect was present. In addition, tests for neutrality detected three loci that are involved in pathogen response and are potentially under selective effects. The diversity conservation and the management practices of autogamous landrace crops are discussed.

Using classical population genetics tools with heterochroneous data: time matters!

BACKGROUND

New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.

METHODOLOGY/PRINCIPAL FINDINGS

We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)).

CONCLUSIONS/SIGNIFICANCE

Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.