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.

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.

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.

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.

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.

Parallel up-regulation of the profilin gene family following independent domestication of diploid and allopolyploid cotton (Gossypium).

Cotton is remarkable among our major crops in that four species were independently domesticated, two allopolyploids and two diploids. In each case thousands of years of human selection transformed sparsely flowering, perennial shrubs into highly productive crops with seeds bearing the vastly elongated and abundant single-celled hairs that comprise modern cotton fiber. The genetic underpinnings of these transformations are largely unknown, but comparative gene expression profiling experiments have demonstrated up-regulation of profilin accompanying domestication in all three species for which wild forms are known. Profilins are actin monomer binding proteins that are important in cytoskeletal dynamics and in cotton fiber elongation. We show that Gossypium diploids contain six profilin genes (GPRF1-GPRF6), located on four different chromosomes (eight chromosomes in the allopolyploid). All but one profilin (GPRF6) are expressed during cotton fiber development, and both homeologs of GPRF1-GPRF5 are expressed in fibers of the allopolyploids. Remarkably, quantitative RT-PCR and RNAseq data demonstrate that GPRF1-GPRF5 are all up-regulated, in parallel, in the three independently domesticated cottons in comparison with their wild counterparts. This result was additionally supported by iTRAQ proteomic data. In the allopolyploids, there This usage of novel should be fine, since it refers to a novel evolutionary process, not a novel discovery has been novel recruitment of the sixth profilin gene (GPRF6) as a result of domestication. This parallel up-regulation of an entire gene family in multiple species in response to strong directional selection is without precedent and suggests unwitting selection on one or more upstream transcription factors or other proteins that coordinately exercise control over profilin expression.

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.

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.

Evolutionary rate variation, genomic dominance and duplicate gene expression evolution during allotetraploid cotton speciation.

Here, we describe the evolution of gene expression among a diversified cohort of five allopolyploid species in the cotton genus (Gossypium). Using this phylogenetic framework and comparisons with expression changes accompanying F(1) hybridization, we provide a temporal perspective on expression diversification following a shared genome duplication. Global patterns of gene expression were studied by the hybridization of petal RNAs to a custom microarray. This platform measures total expression for c. 42 000 duplicated genes, and genome-specific expression for c. 1400 homoeologs (genes duplicated by polyploidy). We report homoeolog expression bias favoring the allopolyploid D genome over the A genome in all species (among five polyploid species, D biases ranging from c. 54 to 60%), in addition to conservation of biases among genes. Furthermore, we find surprising levels of transgressive up- and down-regulation in the allopolyploids, a diminution of the level of bias in genomic expression dominance but not in its magnitude, and high levels of rate variation among allotetraploid species. We illustrate how phylogenetic and temporal components of expression evolution may be partitioned and revealed following allopolyploidy. Overall patterns of expression evolution are similar among the Gossypium allotetraploids, notwithstanding a high level of interspecific rate variation, but differ strikingly from the direction of genomic expression dominance patterns in the synthetic F(1) hybrid.

Gene duplication and evolutionary novelty in plants.

Duplication is a prominent feature of plant genomic architecture. This has led many researchers to speculate that gene duplication may have played an important role in the evolution of phenotypic novelty within plants. Until recently, however, it was difficult to make this connection. We are now beginning to understand how duplication has contributed to adaptive evolution in plants. In this review we introduce the sources of gene duplication and predictions of the various fates of duplicates. We also highlight several recent and pertinent examples from the literature. These examples demonstrate the importance of the functional characteristics of genes and the source of duplication in influencing evolutionary outcome.

Coordinated and fine-scale control of homoeologous gene expression in allotetraploid cotton.

Within polyploid plant species, it has been demonstrated that homoeologous genes (genes duplicated by polyploidy) often display dynamic expression patterns. To determine if chromosomal location plays a role in establishing these expression patterns, we analyzed the relative levels of homoeolog expression among linked genes from 2 locations in the cotton genome. Genes from the region containing the alcohol dehydrogenase A gene show coordinated expression across several tissues, whereas genes from the region containing cellulose synthase A do not. These results indicate that changes in homoeolog expression may be constrained by linkage in some genomic regions, whereas in other regions, homoeolog expression is largely decoupled from physical proximity. Furthermore, these results suggest that both large- and small-scale regulatory mechanisms may control homoeolog expression patterns.

Spotted cotton oligonucleotide microarrays for gene expression analysis.

BACKGROUND: Microarrays offer a powerful tool for diverse applications plant biology and crop improvement. Recently, two comprehensive assemblies of cotton ESTs were constructed based on three Gossypium species. Using these assemblies as templates, we describe the design and creation and of a publicly available oligonucleotide array for cotton, useful for all four of the cultivated species. RESULTS: Synthetic oligonucleotide probes were generated from exemplar sequences of a global assembly of 211,397 cotton ESTs derived from >50 different cDNA libraries representing many different tissue types and tissue treatments. A total of 22,787 oligonucleotide probes are included on the arrays, optimized to target the diversity of the transcriptome and previously studied cotton genes, transcription factors, and genes with homology to Arabidopsis. A small portion of the oligonucleotides target unidentified protein coding sequences, thereby providing an element of gene discovery. Because many oligonucleotides were based on ESTs from fiber-specific cDNA libraries, the microarray has direct application for analysis of the fiber transcriptome. To illustrate the utility of the microarray, we hybridized labeled bud and leaf cDNAs from G. hirsutum and demonstrate technical consistency of results. CONCLUSION: The cotton oligonucleotide microarray provides a reproducible platform for transcription profiling in cotton, and is made publicly available through http://cottonevolution.info.

PHENETIC COMPARISON OF SEVEN Echinacea SPECIES BASED ON IMMUNOMODULATORY CHARACTERISTICS.

The purpose of the present investigation was to compare similarities and differences in immune response among Echinacea species, which are commonly used to treat upper respiratory infections. The investigation involved two components: acquisition of immunomodulatory data reported here for the first time, and combined phenetic analysis of these data along with previous reports. Experimental data were obtained by stimulating human PBMC in vitro with extracts from Echinacea spp. and assaying production of three cytokines (interleukin-1ß [IL-1ß], interleukin-2 [IL-2], and tumor necrosis factor-α [TNF-α]). Phenetic analyses were employed to compare responses across the entire data set, including UPGMA (Unweighted Pair Group Method with Arithmetic Mean) and neighbor-joining methods. In the immune experiments conducted for this investigation, E. angustifolia, E. paradoxa, E. purpurea, E. simulata, and E. tennesseensis extracts significantly augmented IL-1 ß and TNF-α production, whereas no extracts significantly modulated IL-2. All phenetic methods produced similar dendrograms, revealing two species pairs (E. angustifolia + E. simulata and E. pallida + E.sanguinea) where both species cluster tightly and have similar immune-response profiles. These two species-pairs are maximally dissimilar from each other. The remaining species (E. paradoxa, E. purpurea, and E. tennesseensis) occupy intermediate positions in the dendrogram. Our results suggest that Echinacea spp. act heterogeneously on immune function. The utility of these data for science and industry is discussed.

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.

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.

T-DNA locus structure in a large population of soybean plants transformed using the Agrobacterium-mediated cotyledonary-node method.

Designing transformation experiments for either functional genomics or crop improvement requires knowledge of the transgene locus structure, number, transmission and expression resulting from a specific transformation method. We recently reported an improvement to the soybean [Glycine max (L.) Merrill] cotyledonary-node transformation method that resulted in the efficient production of transgenic plants. To characterize the transgene loci resulting from this method, we analysed 270 independent T0 plants and 95 randomly selected T1 progenies for T-DNA locus complexity using Southern analysis. The lines were transformed with Agrobacterium tumefaciens strains LBA4404 or EHA105 carrying the binary plasmids pGPTV, pTOK233, pCAMBIA1303 or pCAMBIA1309, and regenerated in medium supplemented with or without silver nitrate (AgNO3). Analysis in the T0 generation showed that the number of hpt-hybridizing fragments per plant ranged from 1-15, with 31.5% of the lines having a single hpt-hybridizing fragment. Each primary soybean transformant had, on average, 2.0 unlinked transgene loci and that half of the segregating loci in the T1 progenies were single, simple T-DNA insertions. Of the loci containing multiple T-DNA fragments, a low frequency had tandem and inverted repeat T-DNA structures. Integration of binary plasmid backbone sequences occurred in 37% of primary transformants. A. tumefaciens strain, binary plasmid and thiol treatment had no significant effect on transgene locus structure, numbers or expression. Interestingly, exposure of soybean explants to AgNO3 throughout shoot induction and elongation increased T-DNA locus complexity in the primary transformants and decreased silencing of gusA expression in the T1 generation.

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.

Evolutionary genetics of genome merger and doubling in plants.

Polyploidy is a common mode of evolution in flowering plants. The profound effects of polyploidy on gene expression appear to be caused more by hybridity than by genome doubling. Epigenetic mechanisms underlying genome-wide changes in expression are as yet poorly understood; only methylation has received much study, and its importance varies among polyploids. Genetic diploidization begins with the earliest responses to genome merger and doubling; less is known about chromosomal diploidization. Polyploidy duplicates every gene in the genome, providing the raw material for divergence or partitioning of function in homoeologous copies. Preferential retention or loss of genes occurs in a wide range of taxa, suggesting that there is an underlying set of principles governing the fates of duplicated genes. Further studies are required for general patterns to be elucidated, involving different plant families, kinds of polyploidy, and polyploids of different ages.

Phylogenetic, morphological, and chemotaxonomic incongruence in the North American endemic genus Echinacea.

The study of recently formed species is important because it can help us to better understand organismal divergence and the speciation process. However, these species often present difficult challenges in the field of molecular phylogenetics because the processes that drive molecular divergence can lag behind phenotypic divergence. In the current study we show that species of the recently diverged North American endemic genus of purple coneflower, Echinacea, have low levels of molecular divergence. Data from three nuclear loci and two plastid loci provide neither resolved topologies nor congruent hypotheses about species-level relationships. This lack of phylogenetic resolution is likely due to the combined effects of incomplete lineage sorting, hybridization, and backcrossing following secondary contact. The poor resolution provided by molecular markers contrasts previous studies that found well-resolved and taxonomically supported relationships from metabolic and morphological data. These results suggest that phenotypic canalization, resulting in identifiable morphological species, has occurred rapidly within Echinacea. Conversely, molecular signals have been distorted by gene flow and incomplete lineage sorting. Here we explore the impact of natural history on the genetic organization and phylogenetic relationships of Echinacea.