Comparative evolutionary genetics of deleterious load in sorghum and maize.

Sorghum and maize share a close evolutionary history that can be explored through comparative genomics1,2. To perform a large-scale comparison of the genomic variation between these two species, we analysed ~13 million variants identified from whole-genome resequencing of 499 sorghum lines together with 25 million variants previously identified in 1,218 maize lines. Deleterious mutations in both species were prevalent in pericentromeric regions, enriched in non-syntenic genes and present at low allele frequencies. A comparison of deleterious burden between sorghum and maize revealed that sorghum, in contrast to maize, departed from the domestication-cost hypothesis that predicts a higher deleterious burden among domesticates compared with wild lines. Additionally, sorghum and maize population genetic summary statistics were used to predict a gene deleterious index with an accuracy greater than 0.5. This research represents a key step towards understanding the evolutionary dynamics of deleterious variants in sorghum and provides a comparative genomics framework to start prioritizing these variants for removal through genome editing and breeding.

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor).

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.

A community resource for exploring and utilizing genetic diversity in the USDA pea single plant plus collection.

Globally, pea (Pisum sativum L.) is an important temperate legume crop for food, feed and fodder, and many breeding programs develop cultivars adapted to these end-uses. In order to assist pea development efforts, we assembled the USDA Pea Single Plant Plus Collection (PSPPC), which contains 431 P. sativum accessions with morphological, geographic and taxonomic diversity. The collection was characterized genetically in order to maximize its value for trait mapping and genomics-assisted breeding. To that end, we used genotyping-by-sequencing-a cost-effective method for de novo single-nucleotide polymorphism (SNP) marker discovery-to generate 66 591 high-quality SNPs. These data facilitated the identification of accessions divergent from mainstream breeding germplasm that could serve as sources of novel, favorable alleles. In particular, a group of accessions from Central Asia appear nearly as diverse as a sister species, P. fulvum, and subspecies, P. sativum subsp. elatius. PSPPC genotypes can be paired with new and existing phenotype data for trait mapping; as proof-of-concept, we localized Mendel's A gene controlling flower color to its known position. We also used SNP data to define a smaller core collection of 108 accessions with similar levels of genetic diversity as the entire PSPPC, resulting in a smaller germplasm set for research screening and evaluation under limited resources. Taken together, the results presented in this study along with the release of a publicly available SNP data set comprise a valuable resource for supporting worldwide pea genetic improvement efforts.

Population Structure and Phylogenetic Relationships in a Diverse Panel of Brassica rapa L.

The crop species Brassica rapa L. has significant economic importance around the world. However, the global distribution and complex evolutionary history of the species has made investigating its genetic population structure difficult. Crop domestication and improvement has resulted in extreme phenotypic diversity and subspecies that are used for oilseed, food for human consumption, and fodder for livestock. These subspecies include the oilseed morphotypes. oleifera (turnip rape), ssp. dichotoma (brown sarson/toria), ssp. trilocularis (yellow sarson); ssp. rapa (turnip); and Asian leafy vegetables ssp. pekinensis (Chinese cabbage), ssp. chinensis (bok choy), ssp. nipposinica (mizuna/mibuna), ssp. rapifera (rapini/broccoli rabe), ssp. narinosa (tatsoi), ssp parachinensis (choy sum), and ssp. perviridis (komatsuna). To date, studies have had insufficient sampling to determine the relationship of all morphotypes, especially oilseed morphotypes, and questions remain over the contribution of morphotype and geographic origin to population structure. We used genotyping-by-sequencing to score 18,272 single nucleotide polymorphism markers in a globally diverse panel of 333 B. rapa National Plant Germplasm System accessions that included 10 recognized subspecies. Our population genetic and phylogenetic analyses were broadly congruent and revealed five subpopulations that were largely reflective of morphotype and geography. These subpopulations were 1. European turnips/oilseed, 2. Asian turnips/oilseed, 3. yellow/brown sarson (ssp. trilocularis and ssp. dichotoma), 4. Chinese cabbage (ssp. pekinensis), and 5. bok choy, choy sum, and tatsoi (ssp. chinensis, ssp. parachinensis, ssp. narinosa). Additionally, we found evidence of polyphyly and/or paraphyly, particularly for oilseed morphotypes (ssp. oleifera and ssp. dichotoma) and turnips. The results of this study have provided improved resolution to the genetic and phylogenetic relationships of subspecies within the species B. rapa. Understanding of these relationships is key to the future genetic study and improvement of this globally important crop species.

Temporal Genetic Dynamics of an Experimental, Biparental Field Population of Phytophthora capsici.

Defining the contributions of dispersal, reproductive mode, and mating system to the population structure of a pathogenic organism is essential to estimating its evolutionary potential. After introduction of the devastating plant pathogen, Phytophthora capsici, into a grower's field, a lack of aerial spore dispersal restricts migration. Once established, coexistence of both mating types results in formation of overwintering recombinant oospores, engendering persistent pathogen populations. To mimic these conditions, in 2008, we inoculated a field with two P. capsici isolates of opposite mating type. We analyzed pathogenic isolates collected in 2009-2013 from this experimental population, using genome-wide single-nucleotide polymorphism markers. By tracking heterozygosity across years, we show that the population underwent a generational shift; transitioning from exclusively F1 in 2009-2010, to multi-generational in 2011, and ultimately all inbred in 2012-2013. Survival of F1 oospores, characterized by heterozygosity excess, coupled with a low rate of selfing, delayed declines in heterozygosity due to inbreeding and attainment of equilibrium genotypic frequencies. Large allele and haplotype frequency changes in specific genomic regions accompanied the generational shift, representing putative signatures of selection. Finally, we identified an approximately 1.6 Mb region associated with mating type determination, constituting the first detailed genomic analysis of a mating type region (MTR) in Phytophthora. Segregation patterns in the MTR exhibited tropes of sex-linkage, where maintenance of allele frequency differences between isolates of opposite mating type was associated with elevated heterozygosity despite inbreeding. Characterizing the trajectory of this experimental system provides key insights into the processes driving persistent, sexual pathogen populations.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L.

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

Field-Based High-Throughput Plant Phenotyping Reveals the Temporal Patterns of Quantitative Trait Loci Associated with Stress-Responsive Traits in Cotton.

The application of high-throughput plant phenotyping (HTPP) to continuously study plant populations under relevant growing conditions creates the possibility to more efficiently dissect the genetic basis of dynamic adaptive traits. Toward this end, we employed a field-based HTPP system that deployed sets of sensors to simultaneously measure canopy temperature, reflectance, and height on a cotton (Gossypium hirsutum L.) recombinant inbred line mapping population. The evaluation trials were conducted under well-watered and water-limited conditions in a replicated field experiment at a hot, arid location in central Arizona, with trait measurements taken at different times on multiple days across 2010-2012. Canopy temperature, normalized difference vegetation index (NDVI), height, and leaf area index (LAI) displayed moderate-to-high broad-sense heritabilities, as well as varied interactions among genotypes with water regime and time of day. Distinct temporal patterns of quantitative trait loci (QTL) expression were mostly observed for canopy temperature and NDVI, and varied across plant developmental stages. In addition, the strength of correlation between HTPP canopy traits and agronomic traits, such as lint yield, displayed a time-dependent relationship. We also found that the genomic position of some QTL controlling HTPP canopy traits were shared with those of QTL identified for agronomic and physiological traits. This work demonstrates the novel use of a field-based HTPP system to study the genetic basis of stress-adaptive traits in cotton, and these results have the potential to facilitate the development of stress-resilient cotton cultivars.

Genetic Diversity of the Two Commercial Tetraploid Cotton Species in the Gossypium Diversity Reference Set.

A diversity reference set has been constructed for the Gossypium accessions in the US National Cotton Germplasm Collection to facilitate more extensive evaluation and utilization of accessions held in the Collection. A set of 105 mapped simple sequence repeat markers was used to study the allelic diversity of 1933 tetraploid Gossypium accessions representative of the range of diversity of the improved and wild accessions of G. hirsutum and G. barbadense. The reference set contained 410 G. barbadense accessions and 1523 G. hirsutum accessions. Observed numbers of polymorphic and private bands indicated a greater diversity in G. hirsutum as compared to G. barbadense as well as in wild-type accessions as compared to improved accessions in both species. The markers clearly differentiated the 2 species. Patterns of diversity within species were observed but not clearly delineated, with much overlap occurring between races and regions of origin for wild accessions and between historical and geographic breeding pools for cultivated accessions. Although the percentage of accessions showing introgression was higher among wild accessions than cultivars in both species, the average level of introgression within individual accessions, as indicated by species-specific bands, was much higher in wild accessions of G. hirsutum than in wild accessions of G. barbadense. The average level of introgression within individual accessions was higher in improved G. barbadense cultivars than in G. hirsutum cultivars. This molecular characterization reveals the levels and distributions of genetic diversity that will allow for better exploration and utilization of cotton genetic resources.

Population genomic analysis of 962 whole genome sequences of humans reveals natural selection in non-coding regions.

Whole genome analysis in large samples from a single population is needed to provide adequate power to assess relative strengths of natural selection across different functional components of the genome. In this study, we analyzed next-generation sequencing data from 962 European Americans, and found that as expected approximately 60% of the top 1% of positive selection signals lie in intergenic regions, 33% in intronic regions, and slightly over 1% in coding regions. Several detailed functional annotation categories in intergenic regions showed statistically significant enrichment in positively selected loci when compared to the null distribution of the genomic span of ENCODE categories. There was a significant enrichment of purifying selection signals detected in enhancers, transcription factor binding sites, microRNAs and target sites, but not on lincRNA or piRNAs, suggesting different evolutionary constraints for these domains. Loci in "repressed or low activity regions" and loci near or overlapping the transcription start site were the most significantly over-represented annotations among the top 1% of signals for positive selection.

Neutral genomic regions refine models of recent rapid human population growth.

Human populations have experienced dramatic growth since the Neolithic revolution. Recent studies that sequenced a very large number of individuals observed an extreme excess of rare variants and provided clear evidence of recent rapid growth in effective population size, although estimates have varied greatly among studies. All these studies were based on protein-coding genes, in which variants are also impacted by natural selection. In this study, we introduce targeted sequencing data for studying recent human history with minimal confounding by natural selection. We sequenced loci far from genes that meet a wide array of additional criteria such that mutations in these loci are putatively neutral. As population structure also skews allele frequencies, we sequenced 500 individuals of relatively homogeneous ancestry by first analyzing the population structure of 9,716 European Americans. We used very high coverage sequencing to reliably call rare variants and fit an extensive array of models of recent European demographic history to the site frequency spectrum. The best-fit model estimates ∼ 3.4% growth per generation during the last ∼ 140 generations, resulting in a population size increase of two orders of magnitude. This model fits the data very well, largely due to our observation that assumptions of more ancient demography can impact estimates of recent growth. This observation and results also shed light on the discrepancy in demographic estimates among recent studies.

Population growth inflates the per-individual number of deleterious mutations and reduces their mean effect.

This study addresses the question of how purifying selection operates during recent rapid population growth such as has been experienced by human populations. This is not a straightforward problem because the human population is not at equilibrium: population genetics predicts that, on the one hand, the efficacy of natural selection increases as population size increases, eliminating ever more weakly deleterious variants; on the other hand, a larger number of deleterious mutations will be introduced into the population and will be more likely to increase in their number of copies as the population grows. To understand how patterns of human genetic variation have been shaped by the interaction of natural selection and population growth, we examined the trajectories of mutations with varying selection coefficients, using computer simulations. We observed that while population growth dramatically increases the number of deleterious segregating sites in the population, it only mildly increases the number carried by each individual. Our simulations also show an increased efficacy of natural selection, reflected in a higher fraction of deleterious mutations eliminated at each generation and a more efficient elimination of the most deleterious ones. As a consequence, while each individual carries a larger number of deleterious alleles than expected in the absence of growth, the average selection coefficient of each segregating allele is less deleterious. Combined, our results suggest that the genetic risk of complex diseases in growing populations might be distributed across a larger number of more weakly deleterious rare variants.

Copy number variation analysis in the great apes reveals species-specific patterns of structural variation.

Copy number variants (CNVs) are increasingly acknowledged as an important source of evolutionary novelties in the human lineage. However, our understanding of their significance is still hindered by the lack of primate CNV data. We performed intraspecific comparative genomic hybridizations to identify loci harboring copy number variants in each of the four great apes: bonobos, chimpanzees, gorillas, and orangutans. For the first time, we could analyze differences in CNV location and frequency in these four species, and compare them with human CNVs and primate segmental duplication (SD) maps. In addition, for bonobo and gorilla, patterns of CNV and nucleotide diversity were studied in the same individuals. We show that CNVs have been subject to different selective pressures in different lineages. Evidence for purifying selection is stronger in gorilla CNVs overlapping genes, while positive selection appears to have driven the fixation of structural variants in the orangutan lineage. In contrast, chimpanzees and bonobos present high levels of common structural polymorphism, which is indicative of relaxed purifying selection together with the higher mutation rates induced by the known burst of segmental duplication in the ancestor of the African apes. Indeed, the impact of the duplication burst is noticeable by the fact that bonobo and chimpanzee share more CNVs with gorilla than expected. Finally, we identified a number of interesting genomic regions that present high-frequency CNVs in all great apes, while containing only very rare or even pathogenic structural variants in humans.

Selection upon genome architecture: conservation of functional neighborhoods with changing genes.

An increasing number of evidences show that genes are not distributed randomly across eukaryotic chromosomes, but rather in functional neighborhoods. Nevertheless, the driving force that originated and maintains such neighborhoods is still a matter of controversy. We present the first detailed multispecies cartography of genome regions enriched in genes with related functions and study the evolutionary implications of such clustering. Our results indicate that the chromosomes of higher eukaryotic genomes contain up to 12% of genes arranged in functional neighborhoods, with a high level of gene co-expression, which are consistently distributed in phylogenies. Unexpectedly, neighborhoods with homologous functions are formed by different (non-orthologous) genes in different species. Actually, instead of being conserved, functional neighborhoods present a higher degree of synteny breaks than the genome average. This scenario is compatible with the existence of selective pressures optimizing the coordinated transcription of blocks of functionally related genes. If these neighborhoods were broken by chromosomal rearrangements, selection would favor further rearrangements reconstructing other neighborhoods of similar function. The picture arising from this study is a dynamic genomic landscape with a high level of functional organization.

Allele-specific gene expression is widespread across the genome and biological processes.

Allelic specific gene expression (ASGE) appears to be an important factor in human phenotypic variability and as a consequence, for the development of complex traits and diseases. In order to study ASGE across the human genome, we have performed a study in which genotyping was coupled with an analysis of ASGE by screening 11,500 SNPs using the Mapping 10 K Array to identify differential allelic expression. We found that from the 5,133 SNPs that were suitable for analysis (heterozygous in our sample and expressed in peripheral blood mononuclear cells), 2,934 (57%) SNPs had differential allelic expression. Such SNPs were equally distributed along human chromosomes and biological processes. We validated the presence or absence of ASGE in 18 out 20 SNPs (90%) randomly selected by real time PCR in 48 human subjects. In addition, we observed that SNPs close to -but not included in- segmental duplications had increased levels of ASGE. Finally, we found that transcripts of unknown function or non-coding RNAs, also display ASGE: from a total of 2,308 intronic SNPs, 1510 (65%) SNPs underwent differential allelic expression. In summary, ASGE is a widespread mechanism in the human genome whose regulation seems to be far more complex than expected.

SNP analysis to results (SNPator): a web-based environment oriented to statistical genomics analyses upon SNP data.

UNLABELLED: Single nucleotide polymorphisms (SNPs) are the most widely used marker in studies to assess associations between genetic variants and complex traits or diseases. They are also becoming increasingly important in the study of the evolution and history of humans and other species. The analysis and processing of SNPs obtained thanks to high-throughput technologies imply the time consuming and costly use of different, complex and usually format-incompatible software. SNPator is a user-friendly web-based SNP data analysis suite that integrates, among many other algorithms, the most common steps of a SNP association study. It frees the user from the need to have large computer facilities and an in depth knowledge of genetic software installation and management. Genotype data is directly read from the output files of the usual genotyping platforms. Phenotypic data on the samples can also be easily uploaded. Many different quality control and analysis procedures can be performed either by using built-in SNPator algorithms or by calling standard genetic software. AVAILABILITY: Access is granted from the SNPator webpage http://www.snpator.org.

On the association between chromosomal rearrangements and genic evolution in humans and chimpanzees.

BACKGROUND: The role that chromosomal rearrangements might have played in the speciation processes that have separated the lineages of humans and chimpanzees has recently come into the spotlight. To date, however, results are contradictory. Here we revisit this issue by making use of the available human and chimpanzee genome sequence to study the relationship between chromosomal rearrangements and rates of DNA sequence evolution. RESULTS: Contrary to previous findings for this pair of species, we show that genes located in the rearranged chromosomes that differentiate the genomes of humans and chimpanzees, especially genes within rearrangements themselves, present lower divergence than genes elsewhere in the genome. Still, there are considerable differences between individual chromosomes. Chromosome 4, in particular, presents higher divergence in genes located within its rearrangement. CONCLUSION: A first conclusion of our analysis is that divergence is lower for genes located in rearranged chromosomes than for those in colinear chromosomes. We also report that non-coding regions within rearranged regions tend to have lower divergence than non-coding regions outside them. These results suggest an association between chromosomal rearrangements and lower non-coding divergence that has not been reported before, even if some chromosomes do not follow this trend and could be potentially associated with a speciation episode. In summary, without excluding it, our results suggest that chromosomal speciation has not been common along the human and chimpanzee lineage.

Patterns and rates of intron divergence between humans and chimpanzees.

BACKGROUND: Introns, which constitute the largest fraction of eukaryotic genes and which had been considered to be neutral sequences, are increasingly acknowledged as having important functions. Several studies have investigated levels of evolutionary constraint along introns and across classes of introns of different length and location within genes. However, thus far these studies have yielded contradictory results. RESULTS: We present the first analysis of human-chimpanzee intron divergence, in which differences in the number of substitutions per intronic site (Ki) can be interpreted as the footprint of different intensities and directions of the pressures of natural selection. Our main findings are as follows: there was a strong positive correlation between intron length and divergence; there was a strong negative correlation between intron length and GC content; and divergence rates vary along introns and depending on their ordinal position within genes (for instance, first introns are more GC rich, longer and more divergent, and divergence is lower at the 3' and 5' ends of all types of introns). CONCLUSION: We show that the higher divergence of first introns is related to their larger size. Also, the lower divergence of short introns suggests that they may harbor a relatively greater proportion of regulatory elements than long introns. Moreover, our results are consistent with the presence of functionally relevant sequences near the 5' and 3' ends of introns. Finally, our findings suggest that other parts of introns may also be under selective constraints.

Chromosomal phylogeny of Robertsonian races of the house mouse on the island of Madeira: testing between alternative mutational processes.

The ancestral karyotype of the house mouse (Mus musculus) consists of 40 acrocentric chromosomes, but numerous races exist within the domesticus subspecies characterized by different metacentric chromosomes formed by the joining at the centromere of two acrocentrics. An exemplary case is present on the island of Madeira where six highly divergent chromosomal races have accumulated different combinations of 20 metacentrics in 500-1000 years. Chromosomal cladistic phylogenies were performed to test the relative performance of Robertsonian (Rb) fusions, Rb fissions and whole-arm reciprocal translocations (WARTs) in resolving relationships between the chromosomal races. The different trees yielded roughly similar topologies, but varied in the number of steps and branch support. The analyses using Rb fusions/fissions as characters resulted in poorly supported trees requiring six to eight homoplasious events. Allowance for WARTs considerably increased nodal support and yielded the most parsimonious trees since homoplasy was reduced to a single event. The WART-based trees required five to nine WARTs and 12 to 16 Rb fusions. These analyses provide support for the role of WARTs in generating the extensive chromosomal diversification observed in house mice. The repeated occurrence of Rb fusions and WARTs highlights the contribution of centromere-related rearrangements to accelerated rates of chromosomal change in the house mouse.

Does radial nuclear organisation influence DNA damage?

It has been suggested that chromatin at the nuclear periphery could act to shield DNA sequences in the nuclear interior from damage. To test this hypothesis, we have examined the nuclear distribution of sites of DNA repair induced by oxidation or UV-C. We do not detect more damage (repair) at the nuclear periphery than in the nuclear interior. In fact, contrary to the body guard hypothesis, there is an excess of damage detected in the nuclear interior. This is further supported by sequence comparison between genes on human chromosomes 18 or 19, and their counterparts in the chimpanzee. The synonymous substitution rate for genes on chromosome 19, which is located towards the centre of the human nucleus, was higher than that for genes on chromosome 18, which is located at the nuclear periphery. We conclude that chromatin at the periphery of the human nucleus is not able to protect more internally located sequences from damage and mutation. We suggest that features of the chromatin structure, or base composition, of sequences in the nuclear centre make them more susceptible to damage.