Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity

Rural populations facilitated early SARS-CoV-2 evolution and transmission in Missouri, USA

In the United States, rural populations comprise 60 million individuals and suffered from high COVID-19 disease burdens. Despite this, surveillance efforts are biased toward urban centers. Consequently, how rurally circulating SARS-CoV-2 viruses contribute toward emerging variants remains poorly understood. In this study, we aim to investigate the role of rural communities in the evolution and transmission of SARS-CoV-2 during the early pandemic. We collected 544 urban and 435 rural COVID-19-positive respiratory specimens from an overall vaccine-naïve population in Southwest Missouri between July and December 2020. Genomic analyses revealed 53 SARS-CoV-2 Pango lineages in our study samples, with 14 of these lineages identified only in rural samples. Phylodynamic analyses showed that frequent bi-directional diffusions occurred between rural and urban communities in Southwest Missouri, and that four out of seven Missouri rural-origin lineages spread globally. Further analyses revealed that the nucleocapsid protein (N):R203K/G204R paired substitutions, which were detected disproportionately across multiple Pango lineages, were more associated with urban than rural sequences. Positive selection was detected at N:204 among rural samples but was not evident in urban samples, suggesting that viruses may encounter distinct selection pressures in rural versus urban communities. This study demonstrates that rural communities may be a crucial source of SARS-CoV-2 evolution and transmission, highlighting the need to expand surveillance and resources to rural populations for COVID-19 mitigation.

UDP-glucosyltransferase HvUGT13248 confers type II resistance to Fusarium graminearum in barley

Fusarium head blight (FHB) of barley (Hordeum vulgare) causes yield losses and accumulation of trichothecene mycotoxins (e.g. deoxynivalenol [DON]) in grains. Glucosylation of DON to the nontoxic DON-3-O-glucoside (D3G) is catalyzed by UDP-glucosyltransferases (UGTs), such as barley UGT13248. We explored the natural diversity of UGT13248 in 496 barley accessions and showed that all carried potential functional alleles of UGT13248, as no genotypes showed strongly increased seedling sensitivity to DON. From a TILLING population, we identified 2 mutant alleles (T368I and H369Y) that, based on protein modeling, likely affect the UDP-glucose binding of UGT13248. In DON feeding experiments, DON-to-D3G conversion was strongly reduced in spikes of these mutants compared to controls, and plants overexpressing UGT13248 showed increased resistance to DON and increased DON-to-D3G conversion. Moreover, field-grown plants carrying the T368I or H369Y mutations inoculated with Fusarium graminearum showed increased FHB disease severity and reduced D3G production. Barley is generally considered to have type II resistance that limits the spread of F. graminearum from the infected spikelet to adjacent spikelets. Point inoculation experiments with F. graminearum showed increased infection spread in T368I and H369Y across the spike compared to wild type, while overexpression plants showed decreased spread of FHB symptoms. Confocal microscopy revealed that F. graminearum spread to distant rachis nodes in T368I and H369Y mutants but was arrested at the rachis node of the inoculated spikelet in wild-type plants. Taken together, our data reveal that UGT13248 confers type II resistance to FHB in barley via conjugation of DON to D3G.

Plant Genome Editing and the Relevance of Off-Target Changes

Site-directed nucleases (SDNs) used for targeted genome editing are powerful new tools to introduce precise genetic changes into plants. Like traditional approaches, such as conventional crossing and induced mutagenesis, genome editing aims to improve crop yield and nutrition. Next-generation sequencing studies demonstrate that across their genomes, populations of crop species typically carry millions of single nucleotide polymorphisms and many copy number and structural variants. Spontaneous mutations occur at rates of ∼10-8 to 10-9 per site per generation, while variation induced by chemical treatment or ionizing radiation results in higher mutation rates. In the context of SDNs, an off-target change or edit is an unintended, nonspecific mutation occurring at a site with sequence similarity to the targeted edit region. SDN-mediated off-target changes can contribute to a small number of additional genetic variants compared to those that occur naturally in breeding populations or are introduced by induced-mutagenesis methods. Recent studies show that using computational algorithms to design genome editing reagents can mitigate off-target edits in plants. Finally, crops are subject to strong selection to eliminate off-type plants through well-established multigenerational breeding, selection, and commercial variety development practices. Within this context, off-target edits in crops present no new safety concerns compared to other breeding practices. The current generation of genome editing technologies is already proving useful to develop new plant varieties with consumer and farmer benefits. Genome editing will likely undergo improved editing specificity along with new developments in SDN delivery and increasing genomic characterization, further improving reagent design and application.

The Fate of Deleterious Variants in a Barley Genomic Prediction Population

Targeted identification and purging of deleterious genetic variants has been proposed as a novel approach to animal and plant breeding. This strategy is motivated, in part, by the observation that demographic events and strong selection associated with cultivated species pose a "cost of domestication." This includes an increase in the proportion of genetic variants that are likely to reduce fitness. Recent advances in DNA resequencing and sequence constraint-based approaches to predict the functional impact of a mutation permit the identification of putatively deleterious SNPs (dSNPs) on a genome-wide scale. Using exome capture resequencing of 21 barley lines, we identified 3855 dSNPs among 497,754 total SNPs. We generated whole-genome resequencing data of Hordeum murinum ssp. glaucum as a phylogenetic outgroup to polarize SNPs as ancestral vs. derived. We also observed a higher proportion of dSNPs per synonymous SNPs (sSNPs) in low-recombination regions of the genome. Using 5215 progeny from a genomic prediction experiment, we examined the fate of dSNPs over three breeding cycles. Adjusting for initial frequency, derived alleles at dSNPs reduced in frequency or were lost more often than other classes of SNPs. The highest-yielding lines in the experiment, as chosen by standard genomic prediction approaches, carried fewer homozygous dSNPs than randomly sampled lines from the same progeny cycle. In the final cycle of the experiment, progeny selected by genomic prediction had a mean of 5.6% fewer homozygous dSNPs relative to randomly chosen progeny from the same cycle.

The Role of Deleterious Substitutions in Crop Genomes

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.

Genome-wide patterns of recombination, linkage disequilibrium and nucleotide diversity from pooled resequencing and single nucleotide polymorphism genotyping unlock the evolutionary history of Eucalyptus grandis

We used high-density single nucleotide polymorphism (SNP) data and whole-genome pooled resequencing to examine the landscape of population recombination (ρ) and nucleotide diversity (ϴw ), assess the extent of linkage disequilibrium (r(2) ) and build the highest density linkage maps for Eucalyptus. At the genome-wide level, linkage disequilibrium (LD) decayed within c. 4-6 kb, slower than previously reported from candidate gene studies, but showing considerable variation from absence to complete LD up to 50 kb. A sharp decrease in the estimate of ρ was seen when going from short to genome-wide inter-SNP distances, highlighting the dependence of this parameter on the scale of observation adopted. Recombination was correlated with nucleotide diversity, gene density and distance from the centromere, with hotspots of recombination enriched for genes involved in chemical reactions and pathways of the normal metabolic processes. The high nucleotide diversity (ϴw = 0.022) of E. grandis revealed that mutation is more important than recombination in shaping its genomic diversity (ρ/ϴw = 0.645). Chromosome-wide ancestral recombination graphs allowed us to date the split of E. grandis (1.7-4.8 million yr ago) and identify a scenario for the recent demographic history of the species. Our results have considerable practical importance to Genome Wide Association Studies (GWAS), while indicating bright prospects for genomic prediction of complex phenotypes in eucalypt breeding.

A Male-Specific Genetic Map of the Microcrustacean Daphnia pulex Based on Single-Sperm Whole-Genome Sequencing

Genetic linkage maps are critical for assembling draft genomes to a meaningful chromosome level and for deciphering the genomic underpinnings of biological traits. The estimates of recombination rates derived from genetic maps also play an important role in understanding multiple aspects of genomic evolution such as nucleotide substitution patterns and accumulation of deleterious mutations. In this study, we developed a high-throughput experimental approach that combines fluorescence-activated cell sorting, whole-genome amplification, and short-read sequencing to construct a genetic map using single-sperm cells. Furthermore, a computational algorithm was developed to analyze single-sperm whole-genome sequencing data for map construction. These methods allowed us to rapidly build a male-specific genetic map for the freshwater microcrustacean Daphnia pulex, which shows significant improvements compared to a previous map. With a total of mapped 1672 haplotype blocks and an average intermarker distance of 0.87 cM, this map spans a total genetic distance of 1451 Kosambi cM and comprises 90% of the resolved regions in the current Daphnia reference assembly. The map also reveals the mistaken mapping of seven scaffolds in the reference assembly onto chromosome II by a previous microsatellite map based on F2 crosses. Our approach can be easily applied to many other organisms and holds great promise for unveiling the intragenomic and intraspecific variation in the recombination rates.

Genome-wide linkage-disequilibrium profiles from single individuals

Although the analysis of linkage disequilibrium (LD) plays a central role in many areas of population genetics, the sampling variance of LD is known to be very large with high sensitivity to numbers of nucleotide sites and individuals sampled. Here we show that a genome-wide analysis of the distribution of heterozygous sites within a single diploid genome can yield highly informative patterns of LD as a function of physical distance. The proposed statistic, the correlation of zygosity, is closely related to the conventional population-level measure of LD, but is agnostic with respect to allele frequencies and hence likely less prone to outlier artifacts. Application of the method to several vertebrate species leads to the conclusion that >80% of recombination events are typically resolved by gene-conversion-like processes unaccompanied by crossovers, with the average lengths of conversion patches being on the order of one to several kilobases in length. Thus, contrary to common assumptions, the recombination rate between sites does not scale linearly with distance, often even up to distances of 100 kb. In addition, the amount of LD between sites separated by <200 bp is uniformly much greater than can be explained by the conventional neutral model, possibly because of the nonindependent origin of mutations within this spatial scale. These results raise questions about the application of conventional population-genetic interpretations to LD on short spatial scales and also about the use of spatial patterns of LD to infer demographic histories.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Genetic diversity and ecological niche modelling of wild barley: refugia, large-scale post-LGM range expansion and limited mid-future climate threats?

Describing genetic diversity in wild barley (Hordeum vulgare ssp. spontaneum) in geographic and environmental space in the context of current, past and potential future climates is important for conservation and for breeding the domesticated crop (Hordeum vulgare ssp. vulgare). Spatial genetic diversity in wild barley was revealed by both nuclear- (2,505 SNP, 24 nSSR) and chloroplast-derived (5 cpSSR) markers in 256 widely-sampled geo-referenced accessions. Results were compared with MaxEnt-modelled geographic distributions under current, past (Last Glacial Maximum, LGM) and mid-term future (anthropogenic scenario A2, the 2080s) climates. Comparisons suggest large-scale post-LGM range expansion in Central Asia and relatively small, but statistically significant, reductions in range-wide genetic diversity under future climate. Our analyses support the utility of ecological niche modelling for locating genetic diversity hotspots and determine priority geographic areas for wild barley conservation under anthropogenic climate change. Similar research on other cereal crop progenitors could play an important role in tailoring conservation and crop improvement strategies to support future human food security.

Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins

The levels of diversity and extent of linkage disequilibrium in cultivated species are largely determined by diversity in their wild progenitors. We report a comparison of nucleotide sequence diversity in wild and cultivated barley (Hordeum vulgare ssp. spontaneum and ssp. vulgare) at 7 nuclear loci totaling 9296bp, using sequence from Hordeum bulbosum to infer the ancestral state of mutations. The sample includes 36 accessions of cultivated barley, including 23 landraces (cultivated forms not subject to modern breeding) and 13 cultivated lines and genetic stocks compared to either 25 or 45 accessions of wild barley for the same loci. Estimates of nucleotide sequence diversity indicate that landraces retain >80% of the diversity in wild barley. The primary population structure in wild barley, which divides the species into eastern and western populations, is reflected in significant differentiation at all loci in wild accessions and at 3 of 7 loci in landraces. "Oriental" landraces have slightly higher diversity than "Occidental" landraces. Genetic assignment suggests more admixture from Occidental landraces into Oriental landraces than the converse, which may explain this difference. Based on θπ for silent sites, modern western cultivars have ~73% of the diversity found in landraces and ~71% of the diversity in wild barley.

Tracing the geographic origins of weedy Ipomoea purpurea in the southeastern United States

Ipomoea purpurea (common morning glory) is an annual vine native to Mexico that is well known for its large, showy flowers. Humans have spread morning glories worldwide, owing to the horticultural appeal of morning glory flowers. Ipomoea purpurea is an opportunistic colonizer of disturbed habitats including roadside and agricultural settings, and it is now regarded as a noxious weed in the Southeastern US. Naturalized populations in the Southeastern United States are highly polymorphic for a number of flower color morphs, unlike native Mexican populations that are typically monomorphic for the purple color morph. Although I. purpurea was introduced into the United States from Mexico, little is known about the specific geographic origins of US populations relative to the Mexican source. We use resequencing data from 11 loci and 30 I. purpurea accessions collected from the native range of the species in Central and Southern Mexico and 8 accessions from the Southeastern United States to infer likely geographic origins in Mexico. Based on genetic assignment analysis, haplotype composition, and the degree of shared polymorphism, I. purpurea samples from the Southeastern United States are genetically most similar to samples from the Valley of Mexico and Veracruz State. This supports earlier speculation that I. purpurea in the Southeastern United States was likely to have been introduced by European colonists from sources in Central Mexico.

Molecular evolution of human adenoviruses

The recent emergence of highly virulent human adenoviruses (HAdVs) with new tissue tropisms underscores the need to determine their ontogeny. Here we report complete high quality genome sequences and analyses for all the previously unsequenced HAdV serotypes (n = 20) within HAdV species D. Analysis of nucleotide sequence variability for these in conjunction with another 40 HAdV prototypes, comprising all seven HAdV species, confirmed the uniquely hypervariable regions within species. The mutation rate among HAdV-Ds was low when compared to other HAdV species. Homologous recombination was identified in at least two of five examined hypervariable regions for every virus, suggesting the evolution of HAdV-Ds has been highly dependent on homologous recombination. Patterns of alternating GC and AT rich motifs correlated well with hypervariable region recombination sites across the HAdV-D genomes, suggesting foci of DNA instability lead to formulaic patterns of homologous recombination and confer agility to adenovirus evolution.

Nucleotide sequence diversity of floral pigment genes in Mexican populations of Ipomoea purpurea (morning glory) accord with a neutral model of evolution

The common morning glory (Ipomoea purpurea) is an annual vine native to Central and Southern Mexico. The genetics of flower color polymorphisms and interactions with the biotic environment have been extensively studied in I. purpurea and in its sister species I. nil. In this study, we examine nucleotide sequence polymorphism in 11 loci, 9 of which are known to participate in a pathway that produces floral pigments. A sample of 30 I. purpurea accessions from the native range of Central and Southern Mexico comprise the data, along with one accession from each of the two sister species I. alba and I. nil. We observe moderate levels of nucleotide sequence polymorphism of ~1%. The ratio of recombination to mutation parameter estimates (ρ/θ) of ~2.5 appears consistent with a mixed-mating system. Ipomoea resequencing data from these genic regions are noteworthy in providing a good fit to the standard neutral model of molecular evolution. The derived silent site frequency spectrum is very close to that predicted by coalescent simulations of a drift-mutation process, and Tajima's D values are not significantly different from expectations under neutrality.

Fine-scale population recombination rates, hotspots, and correlates of recombination in the Medicago truncatula genome

Recombination rates vary across the genome and in many species show significant relationships with several genomic features, including distance to the centromere, gene density, and GC content. Studies of fine-scale recombination rates have also revealed that in several species, there are recombination hotspots, that is, short regions with recombination rates 10–100 greater than those in surrounding regions. In this study, we analyzed whole-genome resequence data from 26 accessions of the model legume Medicago truncatula to gain insight into the genomic features that are related to high- and low-recombination rates and recombination hotspots at 1 kb scales. We found that high-recombination regions (1-kb windows among those in the highest 5% of the distribution) on all three chromosomes were significantly closer to the centromere, had higher gene density, and lower GC content than low-recombination windows. High-recombination windows are also significantly overrepresented among some gene functional categories—most strongly NB–ARC and LRR genes, both of which are important in plant defense against pathogens. Similar to high-recombination windows, recombination hotspots (1-kb windows with significantly higher recombination than the surrounding region) are significantly nearer to the centromere than nonhotspot windows. By contrast, we detected no difference in gene density or GC content between hotspot and nonhotspot windows. Using linear model wavelet analysis to examine the relationship between recombination and genomic features across multiple spatial scales, we find a significant negative correlation with distance to the centromere across scales up to 512 kb, whereas gene density and GC content show significantly positive and negative correlations, respectively, only up to 64 kb. Correlations between recombination and genomic features, particularly gene density and polymorphism, suggest that they are scale dependent and need to be assessed at scales relevant to the evolution of those features.

Retention of agronomically important variation in germplasm core collections: implications for allele mining

The primary targets of allele mining efforts are loci of agronomic importance. Agronomic loci typically exhibit patterns of allelic diversity that are consistent with a history of natural or artificial selection. Natural or artificial selection causes the distribution of genetic diversity at such loci to deviate substantially from the pattern found at neutral loci. The germplasm utilized for allele mining should contain maximum allelic variation at loci of interest, in the smallest possible number of samples. We show that the popular core collection assembly procedure "M" (marker allele richness), which leverages variation at neutral loci, performs worse than random assembly for retaining variation at a locus of agronomic importance in sugar beet (Beta vulgaris L. subsp. vulgaris) that is under selection. We present a corrected procedure ("M+") that outperforms M. An extensive coalescent simulation was performed to demonstrate more generally the retention of neutral versus selected allelic variation in core subsets assembled with M+. A negative correlation in level of allelic diversity between neutral and selected loci was observed in 42% of simulated data sets. When core collection assembly is guided by neutral marker loci, as is the current common practice, enhanced allelic variation at agronomically important loci should not necessarily be expected.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

The contribution of recombination to heterozygosity differs among plant evolutionary lineages and life-forms

BACKGROUND

Despite its role as a generator of haplotypic variation, little is known about how the rates of recombination evolve across taxa. Recombination is a very labile force, susceptible to evolutionary and life trait related processes, which have also been correlated with general levels of genetic diversity. For example, in plants, it has been shown that long-lived outcrossing taxa, such as trees, have higher heterozygosity (He) at SSRs and allozymes than selfing or annual species. However, some of these tree taxa have surprisingly low levels of nucleotide diversity at the DNA sequence level, which points to recombination as a potential generator of genetic diversity in these organisms. In this study, we examine how genome-wide and within-gene rates of recombination evolve across plant taxa, determine whether such rates are influenced by the life-form adopted by species, and evaluate if higher genome-wide rates of recombination translate into higher He values, especially in trees.

RESULTS

Estimates of genome-wide (cM/Mb) recombination rates from 81 higher plants showed a significant phylogenetic signal. The use of different comparative phylogenetic models demonstrated that there is a positive correlation between recombination rate and He (0.83 +/- 0.29), and that trees have higher rates of genome-wide recombination than short-lived herbs and shrubs. A significant taxonomic component was further made evident by our models, as conifers exhibited lower recombination rates than angiosperms. This trend was also found at the within-gene level.

CONCLUSIONS

Altogether, our results illustrate how both common ancestry and life-history traits have to be taken into account for understanding the evolution of genetic diversity and genomic rates of recombination across plant species, and highlight the relevance of species life forms to explain general levels of diversity and recombination.

Strong correlation of wild barley (Hordeum spontaneum) population structure with temperature and precipitation variation

In this study, we present the genetic analysis of a new collection of wild barley (Hordeum spontaneum) using 42 simple sequence repeat (SSR) markers that represent the seven chromosomes. The Barley1K (B1K) infrastructure consists of 1020 accessions collected in a hierarchical sampling mode (HSM) from 51 sites across Israel and represents the wide adaptive niche of the modern barley's ancestor. According to the genetic structure analysis, the sampled sites can be divided into seven groups, and sampled microsites located on opposing slopes or in different soil types did not show significant genetic differentiation. Although the genetic analysis indicates a simple isolation-by-distance model among the populations, examination of the genetic populations' structure with abiotic parameters in an ordination analysis revealed that the combination of elevation, mid-day temperature and rainfall explains a high proportion of the variance in the principal components analysis. Our findings demonstrate that the current populations have therefore been shaped and distinguished by non-selective forces such as migration; however, we suggest that aridity and temperature gradients played major roles as selective forces in the adaptation of wild barley in this part of the Fertile Crescent. This unique collection is a prelude for the investigation of the molecular basis underlying plant adaptation and responsiveness to harsh environments.

Multilocus patterns of nucleotide diversity, population structure and linkage disequilibrium in Boechera stricta, a wild relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in approximately 10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have approximately 45% higher recombination per kilobase than outcrossing A. lyrata.

Complex signatures of locus-specific selective pressures and gene conversion on Human Growth Hormone/Chorionic Somatomammotropin genes

Reduced birth weight and slow neonatal growth are risks correlated with the development of common diseases in adulthood. The Human Growth Hormone/Chorionic Somatomammotropin (hGH/CSH) gene cluster (48 kb) at 17q22-24, consisting of one pituitary-expressed postnatal (GH1) and four placental genes (GH2, CSH1, CSH2, and CSHL1) may contribute to common variation in intrauterine and infant growth, and also to the regulation of feto-maternal and adult glucose metabolism. In contrast to GH1, there are limited genetic data on the hGH/CSH genes expressed in utero. We report the first survey of sequence variation encompassing all five hGH/CSH genes. Resequencing identified 113 SNPs/indels (ss86217675-ss86217787 in dbSNP) including 66 novel variants, and revealed remarkable differences in diversity patterns among the homologous duplicated genes as well as between the study populations of European (Estonians), Asian (Han Chinese), and African (Mandenkalu) ancestries. A dominant feature of the hGH/CSH region is hyperactive gene conversion, with the rate exceeding tens to hundreds of times the rate of reciprocal crossing-over and resulting in near absence of linkage disequilibrium. The initiation of gene conversion seems to be uniformly distributed because the data do not predict any recombination hotspots. Signatures of different selective constraints acting on each gene indicate functional specification of the hGH/CSH genes. Most strikingly, the GH2 coding for placental growth hormone shows strong intercontinental diversification (F(ST)=0.41-0.91; p<10(-6)) indicative of balancing selection, whereas the flanking CSH1 exhibits low population differentiation (F(ST)=0.03-0.09), low diversity (non-Africans, pi=8-9 x 10(-5); Africans, pi=8.2 x 10(-4)), and one dominant haplotype worldwide, consistent with purifying selection. The results imply that the success of an association study targeted to duplicated genes may be enhanced by prior resequencing of the study population in order to determine polymorphism distribution and relevant tag-SNPs.

Population genomics of the wild yeast Saccharomyces paradoxus: Quantifying the life cycle

Most microbes have complex life cycles with multiple modes of reproduction that differ in their effects on DNA sequence variation. Population genomic analyses can therefore be used to estimate the relative frequencies of these different modes in nature. The life cycle of the wild yeast Saccharomyces paradoxus is complex, including clonal reproduction, outcrossing, and two different modes of inbreeding. To quantify these different aspects we analyzed DNA sequence variation in the third chromosome among 20 isolates from two populations. Measures of mutational and recombinational diversity were used to make two independent estimates of the population size. In an obligately sexual population these values should be approximately equal. Instead there is a discrepancy of about three orders of magnitude between our two estimates of population size, indicating that S. paradoxus goes through a sexual cycle approximately once in every 1,000 asexual generations. Chromosome III also contains the mating type locus (MAT), which is the most outbred part in the entire genome, and by comparing recombinational diversity as a function of distance from MAT we estimate the frequency of matings to be approximately 94% from within the same tetrad, 5% with a clonemate after switching the mating type, and 1% outcrossed. Our study illustrates the utility of population genomic data in quantifying life cycles.

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.

The new biology: beyond the Modern Synthesis

BACKGROUND

The last third of the 20th Century featured an accumulation of research findings that severely challenged the assumptions of the "Modern Synthesis" which provided the foundations for most biological research during that century. The foundations of that "Modernist" biology had thus largely crumbled by the start of the 21st Century. This in turn raises the question of foundations for biology in the 21st Century.

CONCLUSION

Like the physical sciences in the first half of the 20th Century, biology at the start of the 21st Century is achieving a substantive maturity of theory, experimental tools, and fundamental findings thanks to relatively secure foundations in genomics. Genomics has also forced biologists to connect evolutionary and molecular biology, because these formerly Balkanized disciplines have been brought together as actors on the genomic stage. Biologists are now addressing the evolution of genetic systems using more than the concepts of population biology alone, and the problems of cell biology using more than the tools of biochemistry and molecular biology alone. It is becoming increasingly clear that solutions to such basic problems as aging, sex, development, and genome size potentially involve elements of biological science at every level of organization, from molecule to population. The new biology knits together genomics, bioinformatics, evolutionary genetics, and other such general-purpose tools to supply novel explanations for the paradoxes that undermined Modernist biology.

The extent of linkage disequilibrium in rice (Oryza sativa L.)

Despite its status as one of the world's major crops, linkage disequilibrium (LD) patterns have not been systematically characterized across the genome of Asian rice (Oryza sativa). Such information is critical to fully exploit the genome sequence for mapping complex traits using association techniques. Here we characterize LD in five 500-kb regions of the rice genome in three major cultivated rice varieties (indica, tropical japonica, and temperate japonica) and in the wild ancestor of Asian rice, Oryza rufipogon. Using unlinked SNPs to determine the amount of background linkage disequilibrium in each population, we find that the extent of LD is greatest in temperate japonica (probably >500 kb), followed by tropical japonica (approximately 150 kb) and indica (approximately 75 kb). LD extends over a shorter distance in O. rufipogon (<<40 kb) than in any of the O. sativa groups assayed here. The differences in the extent of LD among these groups are consistent with differences in outcrossing and recombination rate estimates. As well as heterogeneity between groups, our results suggest variation in LD patterns among genomic regions. We demonstrate the feasibility of genomewide association mapping in cultivated Asian rice using a modest number of SNPs.

Error detection in SNP data by considering the likelihood of recombinational history implied by three-site combinations

MOTIVATION

Errors in nucleotide sequence and SNP genotyping data are problematic when inferring haplotypes. Previously published methods for error detection in haplotype data make use of pedigree information; however, for many samples, individuals are not related by pedigree. This article describes a method for detecting errors in haplotypes by considering the recombinational history implied by the patterns of variation, three SNPs at a time.

RESULTS

Coalescent simulations provide evidence that the method is robust to high levels of recombination as well as homologous gene conversion, indicating that patterns produced by both proximate and distant SNPs may be useful for detecting unlikely three-site haplotypes.

AVAILABILITY

The perl script implementing the described method is called EDUT (Error Detection Using Triplets) and is available on request from the authors.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Plant domestication, a unique opportunity to identify the genetic basis of adaptation

Despite the fundamental role of plant domestication in human history and the critical importance of a relatively small number of crop plants to modern societies, we still know little about adaptation under domestication. Here we focus on efforts to identify the genes responsible for adaptation to domestication. We start from a historical perspective, arguing that Darwin's conceptualization of domestication and unconscious selection provides valuable insight into the evolutionary history of crops and also provides a framework to evaluate modern methods used to decipher the genetic mechanisms underlying phenotypic change. We then review these methods, framing the discussion in terms of the phenotype-genotype hierarchy. Top-down approaches, such as quantitative trait locus and linkage disequilibrium mapping, start with a phenotype of interest and use genetic analysis to identify candidate genes. Bottom-up approaches, alternatively, use population genetic analyses to identify potentially adaptive genes and then rely on standard bioinformatics and reverse genetic tools to connect selected genes to a phenotype. We discuss the successes, advantages, and challenges of each, but we conclude that bottom-up approaches to understanding domestication as an adaptive process hold greater promise both for the study of adaptation and as a means to identify genes that contribute to agronomically important traits.

Genetic evidence for a second domestication of barley (Hordeum vulgare) east of the Fertile Crescent

Cereal agriculture originated with the domestication of barley and early forms of wheat in the Fertile Crescent. There has long been speculation that barley was domesticated more than once. We use differences in haplotype frequency among geographic regions at multiple loci to infer at least two domestications of barley; one within the Fertile Crescent and a second 1,500-3,000 km farther east. The Fertile Crescent domestication contributed the majority of diversity in European and American cultivars, whereas the second domestication contributed most of the diversity in barley from Central Asia to the Far East.