IntroUNET: Identifying introgressed alleles via semantic segmentation.

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

IntroUNET: identifying introgressed alleles via semantic segmentation.

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge.

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

Variance component estimates, phenotypic characterization, and genetic evaluation of bovine congestive heart failure in commercial feeder cattle.

The increasing incidence of bovine congestive heart failure (BCHF) in feedlot cattle poses a significant challenge to the beef industry from economic loss, reduced performance, and reduced animal welfare attributed to cardiac insufficiency. Changes to cardiac morphology as well as abnormal pulmonary arterial pressure (PAP) in cattle of mostly Angus ancestry have been recently characterized. However, congestive heart failure affecting cattle late in the feeding period has been an increasing problem and tools are needed for the industry to address the rate of mortality in the feedlot for multiple breeds. At harvest, a population of 32,763 commercial fed cattle were phenotyped for cardiac morphology with associated production data collected from feedlot processing to harvest at a single feedlot and packing plant in the Pacific Northwest. A sub-population of 5,001 individuals were selected for low-pass genotyping to estimate variance components and genetic correlations between heart score and the production traits observed during the feeding period. At harvest, the incidence of a heart score of 4 or 5 in this population was approximately 4.14%, indicating a significant proportion of feeder cattle are at risk of cardiac mortality before harvest. Heart scores were also significantly and positively correlated with the percentage Angus ancestry observed by genomic breed percentage analysis. The heritability of heart score measured as a binary (scores 1 and 2 = 0, scores 4 and 5 = 1) trait was 0.356 in this population, which indicates development of a selection tool to reduce the risk of congestive heart failure as an EPD (expected progeny difference) is feasible. Genetic correlations of heart score with growth traits and feed intake were moderate and positive (0.289-0.460). Genetic correlations between heart score and backfat and marbling score were -0.120 and -0.108, respectively. Significant genetic correlation to traits of high economic importance in existing selection indexes explain the increased rate of congestive heart failure observed over time. These results indicate potential to implement heart score observed at harvest as a phenotype under selection in genetic evaluation in order to reduce feedlot mortality due to cardiac insufficiency and improve overall cardiopulmonary health in feeder cattle.

The utility of genomic prediction models in evolutionary genetics.

Variation in complex traits is the result of contributions from many loci of small effect. Based on this principle, genomic prediction methods are used to make predictions of breeding value for an individual using genome-wide molecular markers. In breeding, genomic prediction models have been used in plant and animal breeding for almost two decades to increase rates of genetic improvement and reduce the length of artificial selection experiments. However, evolutionary genomics studies have been slow to incorporate this technique to select individuals for breeding in a conservation context or to learn more about the genetic architecture of traits, the genetic value of missing individuals or microevolution of breeding values. Here, we outline the utility of genomic prediction and provide an overview of the methodology. We highlight opportunities to apply genomic prediction in evolutionary genetics of wild populations and the best practices when using these methods on field-collected phenotypes.

Estimation of Cry3Bb1 resistance allele frequency in field populations of western corn rootworm using a genetic marker.

Following the discovery of western corn rootworm (WCR; Diabrotica virgifera virgifera) populations resistant to the Bacillus thuringiensis (Bt) protein Cry3Bb1, resistance was genetically mapped to a single locus on WCR chromosome 8 and linked SNP markers were shown to correlate with the frequency of resistance among field-collected populations from the US Corn Belt. The purpose of this paper is to further investigate the relationship between one of these resistance-linked markers and the causal resistance locus. Using data from laboratory bioassays and field experiments, we show that one allele of the resistance-linked marker increased in frequency in response to selection, but was not perfectly linked to the causal resistance allele. By coupling the response to selection data with a genetic model of the linkage between the marker and the causal allele, we developed a model that allowed marker allele frequencies to be mapped to causal allele frequencies. We then used this model to estimate the resistance allele frequency distribution in the US Corn Belt based on collections from 40 populations. These estimates suggest that chromosome 8 Cry3Bb1 resistance allele frequency was generally low (<10%) for 65% of the landscape, though an estimated 13% of landscape has relatively high (>25%) resistance allele frequency.

Molecular characterization of Cry1F resistance in fall armyworm, Spodoptera frugiperda from Brazil.

Fall armyworm, Spodoptera frugiperda (J.E. Smith) is a major lepidopteran pest of maize in Brazil and its control particularly relies on the use of genetically engineered crops expressing Bacillus thuringiensis (Bt) toxins such as Cry1F. However, control failures compromising the efficacy of this technology have been reported in many regions in Brazil, but the mechanism of Cry1F resistance in Brazilian fall armyworm populations remained elusive. Here we investigated the molecular mechanism of Cry1F resistance in two field-collected strains of S. frugiperda from Brazil exhibiting high levels of Cry1F resistance. We first rigorously evaluated several candidate reference genes for normalization of gene expression data across strains, larval instars and gut tissues, and identified ribosomal proteins L10, L17 and RPS3A to be most suitable. We then investigated the expression pattern of ten potential Bt toxin receptors/enzymes in both neonates and 2nd instar gut tissue of Cry1F resistant fall armyworm strains compared to a susceptible strain. Next we sequenced the ATP-dependent Binding Cassette subfamily C2 gene (ABCC2) and identified three mutated sites present in ABCC2 of both Cry1F resistant strains: two of them, a GY deletion (positions 788-789) and a P799 K/R amino acid substitution, located in a conserved region of ABCC2 extracellular loop 4 (EC4) and another amino acid substitution, G1088D, but in a less conserved region. We further characterized the role of the novel mutations present in EC4 by functionally expressing both wild type and mutated ABCC2 transporters in insect cell lines, and confirmed a critical role of both sites for Cry1F binding by cell viability assays. Finally, we assessed the frequency of the mutant alleles by pooled population sequencing and pyrosequencing in 40 fall armyworm populations collected from maize fields in different regions in Brazil. We found that the GY deletion being present at high frequency. However we also observed many rare alleles which disrupt residues between sites 783-799, and their diversity and abundance in field collected populations lends further support to the importance of the EC4 domain for Cry1F toxicity.

Evaluating maize phenotypic variance, heritability, and yield relationships at multiple biological scales across agronomically relevant environments.

A challenge to improve an integrative phenotype, like yield, is the interaction between the broad range of possible molecular and physiological traits that contribute to yield and the multitude of potential environmental conditions in which they are expressed. This study collected data on 31 phenotypic traits, 83 annotated metabolites, and nearly 22,000 transcripts from a set of 57 diverse, commercially relevant maize hybrids across three years in central U.S. Corn Belt environments. Although variability in characteristics created a complex picture of how traits interact produce yield, phenotypic traits and gene expression were more consistent across environments, while metabolite levels showed low repeatability. Phenology traits, such as green leaf number and grain moisture and whole plant nitrogen content showed the most consistent correlation with yield. A machine learning predictive analysis of phenotypic traits revealed that ear traits, phenology, and root traits were most important to predicting yield. Analysis suggested little correlation between biomass traits and yield, suggesting there is more of a sink limitation to yield under the conditions studied here. This work suggests that continued improvement of maize yields requires a strong understanding of baseline variation of plant characteristics across commercially-relevant germplasm to drive strategies for consistently improving yield.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity.

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

The Unreasonable Effectiveness of Convolutional Neural Networks in Population Genetic Inference.

Population-scale genomic data sets have given researchers incredible amounts of information from which to infer evolutionary histories. Concomitant with this flood of data, theoretical and methodological advances have sought to extract information from genomic sequences to infer demographic events such as population size changes and gene flow among closely related populations/species, construct recombination maps, and uncover loci underlying recent adaptation. To date, most methods make use of only one or a few summaries of the input sequences and therefore ignore potentially useful information encoded in the data. The most sophisticated of these approaches involve likelihood calculations, which require theoretical advances for each new problem, and often focus on a single aspect of the data (e.g., only allele frequency information) in the interest of mathematical and computational tractability. Directly interrogating the entirety of the input sequence data in a likelihood-free manner would thus offer a fruitful alternative. Here, we accomplish this by representing DNA sequence alignments as images and using a class of deep learning methods called convolutional neural networks (CNNs) to make population genetic inferences from these images. We apply CNNs to a number of evolutionary questions and find that they frequently match or exceed the accuracy of current methods. Importantly, we show that CNNs perform accurate evolutionary model selection and parameter estimation, even on problems that have not received detailed theoretical treatments. Thus, when applied to population genetic alignments, CNNs are capable of outperforming expert-derived statistical methods and offer a new path forward in cases where no likelihood approach exists.

Development and characterization of the first dsRNA-resistant insect population from western corn rootworm, Diabrotica virgifera virgifera LeConte.

The use of dsRNA to control insect pests via the RNA interference (RNAi) pathway is being explored by researchers globally. However, with every new class of insect control compounds, the evolution of insect resistance needs to be considered, and understanding resistance mechanisms is essential in designing durable technologies and effective resistance management strategies. To gain insight into insect resistance to dsRNA, a field screen with subsequent laboratory selection was used to establish a population of DvSnf7 dsRNA-resistant western corn rootworm, Diabrotica virgifera virgifera, a major maize insect pest. WCR resistant to ingested DvSnf7 dsRNA had impaired luminal uptake and resistance was not DvSnf7 dsRNA-specific, as indicated by cross resistance to all other dsRNAs tested. No resistance to the Bacillus thuringiensis Cry3Bb1 protein was observed. DvSnf7 dsRNA resistance was inherited recessively, located on a single locus, and autosomal. Together these findings will provide insights for dsRNA deployment for insect pest control.

Mutational disruption of the ABCC2 gene in fall armyworm, Spodoptera frugiperda, confers resistance to the Cry1Fa and Cry1A.105 insecticidal proteins.

The use of Bt proteins in crops has revolutionized insect pest management by offering effective season-long control. However, field-evolved resistance to Bt proteins threatens their utility and durability. A recent example is field-evolved resistance to Cry1Fa and Cry1A.105 in fall armyworm (Spodoptera frugiperda). This resistance has been detected in Puerto Rico, mainland USA, and Brazil. A S. frugiperda population with suspected resistance to Cry1Fa was sampled from a maize field in Puerto Rico and used to develop a resistant lab colony. The colony demonstrated resistance to Cry1Fa and partial cross-resistance to Cry1A.105 in diet bioassays. Using genetic crosses and proteomics, we show that this resistance is due to loss-of-function mutations in the ABCC2 gene. We characterize two novel mutant alleles from Puerto Rico. We also find that these alleles are absent in a broad screen of partially resistant Brazilian populations. These findings confirm that ABCC2 is a receptor for Cry1Fa and Cry1A.105 in S. frugiperda, and lay the groundwork for genetically enabled resistance management in this species, with the caution that there may be several distinct ABCC2 resistances alleles in nature.

Co-expression network analysis of duplicate genes in maize (Zea mays L.) reveals no subgenome bias.

BACKGROUND: Gene duplication is prevalent in many species and can result in coding and regulatory divergence. Gene duplications can be classified as whole genome duplication (WGD), tandem and inserted (non-syntenic). In maize, WGD resulted in the subgenomes maize1 and maize2, of which maize1 is considered the dominant subgenome. However, the landscape of co-expression network divergence of duplicate genes in maize is still largely uncharacterized. RESULTS: To address the consequence of gene duplication on co-expression network divergence, we developed a gene co-expression network from RNA-seq data derived from 64 different tissues/stages of the maize reference inbred-B73. WGD, tandem and inserted gene duplications exhibited distinct regulatory divergence. Inserted duplicate genes were more likely to be singletons in the co-expression networks, while WGD duplicate genes were likely to be co-expressed with other genes. Tandem duplicate genes were enriched in the co-expression pattern where co-expressed genes were nearly identical for the duplicates in the network. Older gene duplications exhibit more extensive co-expression variation than younger duplications. Overall, non-syntenic genes primarily from inserted duplications show more co-expression divergence. Also, such enlarged co-expression divergence is significantly related to duplication age. Moreover, subgenome dominance was not observed in the co-expression networks - maize1 and maize2 exhibit similar levels of intra subgenome correlations. Intriguingly, the level of inter subgenome co-expression was similar to the level of intra subgenome correlations, and genes from specific subgenomes were not likely to be the enriched in co-expression network modules and the hub genes were not predominantly from any specific subgenomes in maize. CONCLUSIONS: Our work provides a comprehensive analysis of maize co-expression network divergence for three different types of gene duplications and identifies potential relationships between duplication types, duplication ages and co-expression consequences.

Rapid molecular evolution across amniotes of the IIS/TOR network.

The insulin/insulin-like signaling and target of rapamycin (IIS/TOR) network regulates lifespan and reproduction, as well as metabolic diseases, cancer, and aging. Despite its vital role in health, comparative analyses of IIS/TOR have been limited to invertebrates and mammals. We conducted an extensive evolutionary analysis of the IIS/TOR network across 66 amniotes with 18 newly generated transcriptomes from nonavian reptiles and additional available genomes/transcriptomes. We uncovered rapid and extensive molecular evolution between reptiles (including birds) and mammals: (i) the IIS/TOR network, including the critical nodes insulin receptor substrate (IRS) and phosphatidylinositol 3-kinase (PI3K), exhibit divergent evolutionary rates between reptiles and mammals; (ii) compared with a proxy for the rest of the genome, genes of the IIS/TOR extracellular network exhibit exceptionally fast evolutionary rates; and (iii) signatures of positive selection and coevolution of the extracellular network suggest reptile- and mammal-specific interactions between members of the network. In reptiles, positively selected sites cluster on the binding surfaces of insulin-like growth factor 1 (IGF1), IGF1 receptor (IGF1R), and insulin receptor (INSR); whereas in mammals, positively selected sites clustered on the IGF2 binding surface, suggesting that these hormone-receptor binding affinities are targets of positive selection. Further, contrary to reports that IGF2R binds IGF2 only in marsupial and placental mammals, we found positively selected sites clustered on the hormone binding surface of reptile IGF2R that suggest that IGF2R binds to IGF hormones in diverse taxa and may have evolved in reptiles. These data suggest that key IIS/TOR paralogs have sub- or neofunctionalized between mammals and reptiles and that this network may underlie fundamental life history and physiological differences between these amniote sister clades.

Genetic markers for western corn rootworm resistance to Bt toxin.

Western corn rootworm (WCR) is a major maize (Zea mays L.) pest leading to annual economic losses of more than 1 billion dollars in the United States. Transgenic maize expressing insecticidal toxins derived from the bacterium Bacillus thuringiensis (Bt) are widely used for the management of WCR. However, cultivation of Bt-expressing maize places intense selection pressure on pest populations to evolve resistance. Instances of resistance to Bt toxins have been reported in WCR. Developing genetic markers for resistance will help in characterizing the extent of existing issues, predicting where future field failures may occur, improving insect resistance management strategies, and in designing and sustainably implementing forthcoming WCR control products. Here, we discover and validate genetic markers in WCR that are associated with resistance to the Cry3Bb1 Bt toxin. A field-derived WCR population known to be resistant to the Cry3Bb1 Bt toxin was used to generate a genetic map and to identify a genomic region associated with Cry3Bb1 resistance. Our results indicate that resistance is inherited in a nearly recessive manner and associated with a single autosomal linkage group. Markers tightly linked with resistance were validated using WCR populations collected from Cry3Bb1 maize fields showing significant WCR damage from across the US Corn Belt. Two markers were found to be correlated with both diet (R2 = 0.14) and plant (R2 = 0.23) bioassays for resistance. These results will assist in assessing resistance risk for different WCR populations, and can be used to improve insect resistance management strategies.

Evidence of natural selection acting on a polymorphic hybrid incompatibility locus in Mimulus.

As a common cause of reproductive isolation in diverse taxa, hybrid incompatibilities are fundamentally important to speciation. A key question is which evolutionary forces drive the initial substitutions within species that lead to hybrid dysfunction. Previously, we discovered a simple genetic incompatibility that causes nearly complete male sterility and partial female sterility in hybrids between the two closely related yellow monkeyflower species Mimulus guttatus and M. nasutus. In this report, we fine map the two major incompatibility loci-hybrid male sterility 1 (hms1) and hybrid male sterility 2 (hms2)-to small nuclear genomic regions (each <70 kb) that include strong candidate genes. With this improved genetic resolution, we also investigate the evolutionary dynamics of hms1 in a natural population of M. guttatus known to be polymorphic at this locus. Using classical genetic crosses and population genomics, we show that a 320-kb region containing the hms1 incompatibility allele has risen to intermediate frequency in this population by strong natural selection. This finding provides direct evidence that natural selection within plant species can lead to hybrid dysfunction between species.

Speciation and introgression between Mimulus nasutus and Mimulus guttatus.

Mimulus guttatus and M. nasutus are an evolutionary and ecological model sister species pair differentiated by ecology, mating system, and partial reproductive isolation. Despite extensive research on this system, the history of divergence and differentiation in this sister pair is unclear. We present and analyze a population genomic data set which shows that M. nasutus budded from a central Californian M. guttatus population within the last 200 to 500 thousand years. In this time, the M. nasutus genome has accrued genomic signatures of the transition to predominant selfing, including an elevated proportion of nonsynonymous variants, an accumulation of premature stop codons, and extended levels of linkage disequilibrium. Despite clear biological differentiation, we document genomic signatures of ongoing, bidirectional introgression. We observe a negative relationship between the recombination rate and divergence between M. nasutus and sympatric M. guttatus samples, suggesting that selection acts against M. nasutus ancestry in M. guttatus.

Western corn rootworm (Diabrotica virgifera virgifera) transcriptome assembly and genomic analysis of population structure.

BACKGROUND: Western corn rootworm (WCR) is one of the most significant insect pests of maize in North America. WCR has dramatically increased its range in the last century, invading key maize production areas in the US and abroad. In addition, this species has a history of evolving traits that allow it to escape various control options. Improved genetic and genomic resources are crucial tools for understanding population history and the genetic basis of trait evolution. Here we produce and analyze a transcriptome assembly for WCR. We also perform whole genome population resequencing, and combine these resources to better understand the evolutionary history of WCR. RESULTS: The WCR transcriptome assembly presented here contains approximately 16,000 unigenes, many of which have high similarity to other insect species. Among these unigenes we found several gene families that have been implicated in insecticide resistance in other species. We also identified over 500,000 unigene based SNPs among 26 WCR populations. We used these SNPs to scan for outliers among the candidate genes, and to understand how population processes have shaped genetic variation in this species. CONCLUSIONS: This study highlights the utility of transcriptomic and genomic resources as foundational tools for dealing with highly adaptive pest species. Using these tools we identified candidate gene families for insecticide resistance and reveal aspects of WCR population history in light of the species' recent range expansion.

The standing pool of genomic structural variation in a natural population of Mimulus guttatus.

Major unresolved questions in evolutionary genetics include determining the contributions of different mutational sources to the total pool of genetic variation in a species, and understanding how these different forms of genetic variation interact with natural selection. Recent work has shown that structural variants (SVs) (insertions, deletions, inversions, and transpositions) are a major source of genetic variation, often outnumbering single nucleotide variants in terms of total bases affected. Despite the near ubiquity of SVs, major questions about their interaction with natural selection remain. For example, how does the allele frequency spectrum of SVs differ when compared with single nucleotide variants? How often do SVs affect genes, and what are the consequences? To begin to address these questions, we have systematically identified and characterized a large set of submicroscopic insertion and deletion (indel) variants (between 1 and 200 kb in length) among ten inbred lines from a single natural population of the plant species Mimulus guttatus. After extensive computational filtering, we focused on a set of 4,142 high-confidence indels that showed an experimental validation rate of 73%. All but one of these indels were less than 200 kb. Although the largest were generally at lower frequencies in the population, a surprising number of large indels are at intermediate frequencies. Although indels overlapping with genes were much rarer than expected by chance, approximately 600 genes were affected by an indel. Nucleotide-binding site leucine-rich repeat (NBS-LRR) defense response genes were the most enriched among the gene families affected. Most indels associated with genes were rare and appeared to be under purifying selection, though we do find four high-frequency derived insertion alleles that show signatures of recent positive selection.

Duplicate gene evolution, homoeologous recombination, and transcriptome characterization in allopolyploid cotton.

BACKGROUND: Modern allotetraploid cotton contains an "A" and "D" genome from an ancestral polyploidy event that occurred approximately 1-2 million years ago. Diploid A- and D-genome species can be compared to the A- and D-genomes found within these allotetraploids to make evolutionary inferences about polyploidy. In this paper we present a comprehensive EST assembly derived from diploid and model allotetraploid cottons and demonstrate several evolutionary inferences regarding genic evolution that can be drawn from these data. RESULTS: We generated a set of cotton expressed sequence tags (ESTs), comprising approximately 4.4 million Sanger and next-generation (454) transcripts supplemented by approximately 152 million Illumina reads from diploid and allotetraploid cottons. From the EST alignments we inferred 259,192 genome-specific single nucleotide polymorphisms (SNPs). Molecular evolutionary analyses of protein-coding regions demonstrate that the rate of nucleotide substitution has increased among both allotetraploid genomes relative to the diploids, and that the ratio of nonsynonymous to synonymous substitutions has increased in one of the two polyploid lineages we sampled. We also use these SNPs to show that a surprisingly high percentage of duplicate genes (~7 %) show a signature of non-independent evolution in the allotetraploid nucleus, having experienced one or more episodes of nonreciprocal homoeologous recombination (NRHR). CONCLUSIONS: In this study we characterize the functional and mutational properties of the cotton transcriptome, produce a large genome-specific SNP database, and detect illegitimate genetic exchanges between duplicate genomes sharing a common allotetraploid nucleus. Our findings have important implications for our understanding of the consequences of polyploidy and duplicate gene evolution. We demonstrate that cotton genes have experienced an increased rate of molecular evolution following duplication by polyploidy, and that polyploidy has enabled considerable levels of nonreciprocal exchange between homoeologous genes.