The contribution of epistasis to species differences in annual sunflowers

The contribution of epistasis to 15 morphological traits differentiating two annual sunflowers (Helianthus annuus and H. debilis ssp. cucumerifolius) and to hybrid pollen sterility was estimated in a first generation backcross (BC1) mapping population. Analysis of digenic interactions among quantitative trait loci (QTLs) with significant main effects revealed significant interaction effects for six of the 15 morphological traits and for pollen sterility. Likewise, a genome-wide scan of all possible two-locus combinations detected additional significant interactions for three of the traits with significant epistasis in the original analysis: stem pigmentation, phyllary pubescence, and pollen viability. However, these were the only traits of the 16 examined in which detected interactions explained more than 5% of phenotypic variance. The implications of these findings for adaptive evolution and for the introgression of advantageous morphological QTLs across a natural hybrid zone between these taxa are discussed.     

Genetics of species differences in the wild annual sunflowers, Helianthus annuus and H. petiolaris

Much of our knowledge of speciation genetics stems from quantitative trait locus (QTL) studies. However, interpretations of the size and distribution of QTL underlying species differences are complicated by differences in the way QTL magnitudes are estimated. Also, many studies fail to exploit information about QTL directions or to compare inter- and intraspecific QTL variation. Here, we comprehensively analyze an extensive QTL data set for an interspecific backcross between two wild annual sunflowers, Helianthus annuus and H. petiolaris, interpret different estimates of QTL magnitudes, identify trait groups that have diverged through selection, and compare inter- and intraspecific QTL magnitudes. Our results indicate that even minor QTL (in terms of backcross variance) may be surprisingly large compared to levels of standing variation in the parental species or phenotypic differences between them. Morphological traits, particularly flower morphology, were more strongly or consistently selected than life history or physiological traits. Also, intraspecific QTL were generally smaller than interspecific ones, consistent with the prediction that larger QTL are more likely to spread to fixation across a subdivided population. Our results inform the genetics of species differences in Helianthus and suggest an approach for the simultaneous mapping of inter- and intraspecific QTL.     

Genetic architecture of quantitative trait loci associated with morphological and agronomic trait differences in a wild by cultivated barley cross.

Hordeum vulgare subsp. spontaneum is the progenitor of cultivated barley (Hordeum vulgare L.). Domestication combined with plant breeding has led to the morphological and agronomic characteristics of modern barley cultivars. The objective of this study was to map the genetic factors that morphologically and agronomically differentiate wild barley from modern barley cultivars. To address this objective, we identified quantitative trait loci (QTLs) associated with plant height, flag leaf width, spike length, spike width, glume length in relation to seed length, awn length, fragility of ear rachis, endosperm width and groove depth, heading date, flag leaf length, number of tillers per plant, and kernel color in a Harrington/OUH602 advanced backcross (BC2F8) population. This population was genotyped with 113 simple sequence repeat markers. Thirty QTLs were identified, of which 16 were newly identified in this study. One to 4 QTLs were identified for each of the traits except glume length, for which no QTL was detected. The portion of phenotypic variation accounted for by individual QTLs ranged from about 9% to 54%. For traits with more than one QTL, the phenotypic variation explained ranged from 25% to 71%. Taken together, our results reveal the genetic architecture of morphological and agronomic traits that differentiate wild from cultivated barley.     

Effective population size is positively correlated with levels of adaptive divergence among annual sunflowers

The role of adaptation in the divergence of lineages has long been a central question in evolutionary biology, and as multilocus sequence data sets have become available for a wide range of taxa, empirical estimates of levels of adaptive molecular evolution are increasingly common. Estimates vary widely among taxa, with high levels of adaptive evolution in Drosophila, bacteria, and viruses but very little evidence of widespread adaptive evolution in hominids. Although estimates in plants are more limited, some recent work has suggested that rates of adaptive evolution in a range of plant taxa are surprisingly low and that there is little association between adaptive evolution and effective population size in contrast to patterns seen in other taxa. Here, we analyze data from 35 loci for six sunflower species that vary dramatically in effective population size. We find that rates of adaptive evolution are positively correlated with effective population size in these species, with a significant fraction of amino acid substitutions driven by positive selection in the species with the largest effective population sizes but little or no evidence of adaptive evolution in species with smaller effective population sizes. Although other factors likely contribute as well, in sunflowers effective population size appears to be an important determinant of rates of adaptive evolution.     

Detecting adaptive trait introgression between Iris fulva and I. brevicaulis in highly selective field conditions

The idea that natural hybridization has served as an important force in evolutionary and adaptive diversification has gained considerable momentum in recent years. By combining genome analyses with a highly selective field experiment, we provide evidence for adaptive trait introgression between two naturally hybridizing Louisiana Iris species, flood-tolerant Iris fulva and dry-adapted I. brevicaulis. We planted reciprocal backcross (BC1) hybrids along with pure-species plants into natural settings that, due to a flooding event, favored I. fulva. As expected, I. fulva plants survived at much higher rates than I. brevicaulis plants. Backcross hybrids toward I. fulva (BCIF) also survived at significantly higher rates than the reciprocal backcross toward I. brevicaulis (BCIB). Survivorship of BCIB hybrids was strongly influenced by the presence of a number of introgressed I. fulva alleles located throughout the genome, while survivorship in the reciprocal BCIF hybrids was heavily influenced by two epistatically acting QTL of opposite effects. These results demonstrate the potential for adaptive trait introgression between these two species and may help to explain patterns of genetic variation observed in naturally occurring hybrid zones.     

Assessing morphological differences in an adaptive trait: a landmark-based morphometric approach

East African cichlid fishes have evolved a stunning array of oral jaw morphologies. To better understand the adaptive evolution of this trait, we performed a morphological analysis of the jaws of two closely related species from Lake Malawi that have very different modes of feeding. Labeotropheus fuelleborni forages along the substrate with a "biting" mode of feeding, while Metriaclima zebra feeds in the water column with a "sucking" mode. We analyzed each of the four skeletal elements that make up the oral jaws: the dentary, articular, premaxilla, and maxilla. In addition, we performed the same analysis on the neurocranium, an element closely associated with the oral jaws. We used the thin-plate spline method to quantify morphological differences, which allowed us to relate our results to the functional biology of the species. We find many aspects of shape change that relate directly to the functional design of the cichlid head. The same series of measurements was made on hybrids between Labeotropheus and Metriaclima. For every character, hybrid progeny are statistically different from both parental species. These results suggest an additive mode of action of the alleles responsible for these phenotypes.     

Genetic relationships among the annual species of Cicer L.

Summary Genetic relationships between 7 annual species of the genus Cicer, including the cultivated chickpea, have been studied. These species were assigned to 3 crossability groups. In each group interspecific hybrids could be obtained but their fertility differed considerably in the various cross combinations. Crosses between members of different groups yielded no viable seeds. The possibility of gene transfer from the wild species to the cultivated chickpea C. arietinum was also assessed. Only two species could be considered for this purpose, C. reticulatum, which is the wild progenitor of the cultivated species, and C. echinospermum, which is in the secondary gene pool of C. arietinum. A unique postzygotic reproductive barrier mechanism was found between the members of Group II, C. judaicum, C. pinnatifidum and C. bijugum. It is based on a disharmony in the growth rate of the stigma and the anthers at the time of anthesis of the F₁ interspecific hybrid so that selfpollination is avoided. It is proposed that this kind of mechanism has been involved only when an effective spatial isolation between the three species had been obtained.     

Optimized breeding strategies for multiple trait integration: I. Minimizing linkage drag in single event introgression

From a breeding standpoint, multiple trait integration (MTI) is a four-step process of converting an elite variety/hybrid for value-added traits (e.g. transgenic events) using backcross breeding, ultimately regaining the performance attributes of the target hybrid along with reliable expression of the value-added traits. In the light of the overarching goal of recovering equivalent performance in the finished conversion, this study focuses on the first step of MTI, single event introgression, exploring the feasibility of marker-aided backcross conversion of a target maize hybrid for 15 transgenic events, incorporating eight events into the female hybrid parent and seven into the male parent. Single event introgression is conducted in parallel streams to convert the recurrent parent (RP) for individual events, with the primary objective of minimizing residual non-recurrent parent (NRP) germplasm, especially in the chromosomal proximity to the event (i.e. linkage drag). In keeping with a defined lower limit of 96.66 % overall RP germplasm recovery (i.e. ≤120 cM NRP germplasm given a genome size of 1,788 cM), a breeding goal for each of the 15 single event conversions was developed: <8 cM of residual NRP germplasm across the genome with ~1 cM in the 20 cM region flanking the event. Using computer simulation, we aimed to identify optimal breeding strategies for single event introgression to achieve this breeding goal, measuring efficiency in terms of number of backcross generations required, marker data points needed, and total population size across generations. Various selection schemes classified as three-stage, modified two-stage, and combined selection conducted from BC1 through BC3, BC4, or BC5 were compared. The breeding goal was achieved with a selection scheme involving five generations of marker-aided backcrossing, with BC1 through BC3 selected for the event of interest and minimal linkage drag at population size of 600, and BC4 and BC5 selected for the event of interest and recovery of the RP germplasm across the genome at population size of 400, with selection intensity of 0.01 for all generations. In addition, strategies for choice of donor parent to facilitate conversion efficiency and quality were evaluated. Two essential criteria for choosing an optimal donor parent for a given RP were established: introgression history showing reduction of linkage drag to ~1 cM in the 20 cM region flanking the event and genetic similarity between the RP and potential donor parents. Computer simulation demonstrated that single event conversions with <8 cM residual NRP germplasm can be accomplished by BC5 with no genetic similarity, by BC4 with 30 % genetic similarity, and by BC3 with 86 % genetic similarity using previously converted RPs as event donors. This study indicates that MTI to produce a ‘quality’ 15-event-stacked hybrid conversion is achievable. Furthermore, it lays the groundwork for a comprehensive approach to MTI by outlining a pathway to produce appropriate starting materials with which to proceed with event pyramiding and trait fixation before version testing.     

Genetic architecture and adaptive significance of the selfing syndrome in Capsella

The transition from outcrossing to predominant self-fertilization is one of the most common evolutionary transitions in flowering plants. This shift is often accompanied by a suite of changes in floral and reproductive characters termed the selfing syndrome. Here, we characterize the genetic architecture and evolutionary forces underlying evolution of the selfing syndrome in Capsella rubella following its recent divergence from the outcrossing ancestor C. grandiflora. We conduct genotyping by multiplexed shotgun sequencing and map floral and reproductive traits in a large (N= 550) F2 population. Our results suggest that in contrast to previous studies of the selfing syndrome, changes at a few loci, some with major effects, have shaped the evolution of the selfing syndrome in Capsella. The directionality of QTL effects, as well as population genetic patterns of polymorphism and divergence at 318 loci, is consistent with a history of directional selection on the selfing syndrome. Our study is an important step toward characterizing the genetic basis and evolutionary forces underlying the evolution of the selfing syndrome in a genetically accessible model system.     

Trait means,trait plasticity and trait differences to other species jointly explain species performances in grasslands of varying diversity

Functional traits may help to explain the great variety of species performances in plant communities, but it is not clear whether the magnitude of trait values of a focal species or trait differences to co‐occurring species are key for trait‐based predictions. In addition, trait expression within species is often plastic, but this variation has been widely neglected in trait‐based analyses. We studied functional traits and plant biomass of 59 species in 66 experimental grassland mixtures of varying species richness (Jena Experiment). We related mean species performances (species biomass and relative yield RY) and their plasticities along the diversity gradient to trait‐based pedictors involving mean species traits (T_mean), trait plasticities along the diversity gradient (T_slope), extents of trait variation across communities (T_CV; coefficient of variation) and hierarchical differences (T_diff) and trait distances (absolute values of trait differences T_dist) between focal and co‐occurring species. T_mean (30–55%) and T_diff (30–33%) explained most variation in mean species performances and their plasticities, but T_slope (20–25%) was also important in explaining mean species performances. The mean species traits and the trait differences between focal species and neighbors with the greatest explanatory power were related to plant size and stature (shoot length, mass:height ratios) and leaf photosynthetic capacity (specific leaf area, stable carbon isotopes and leaf nitrogen concentration). The contribution of trait plasticities in explaining species performances varied in direction (positive or negative) and involved traits related to photosynthetic capacity, nitrogen acquisition (nitrogen concentrations and stable isotopes) as well as structural stability (shoot carbon concentrations). Our results suggest that incorporating plasticity in trait expression as well as trait differences to co‐occurring species is critical for extending trait‐based analyses to understand the assembly of plant communities and the contribution of individual species in structuring plant communities.     

Genetic basis of adaptive shape differences in the cichlid head

East African cichlids exhibit an extraordinary level of morphological diversity. Key to their success has been a dramatic radiation in trophic biology, which has occurred rapidly and repeatedly in different lakes. In this report we take the first step in understanding the genetic basis of differences in cichlid oral jaw design. We estimate the effective number of genetic factors that control differences in the cichlid head through a comprehensive morphological assessment of two Lake Malawi cichlid species and their F(1) and F(2) hybrid progeny. We estimate that between one and 11 factors underlie shape difference of individual bony elements. We show that many of the skeletal differences in the head and oral jaw apparatus are inherited together, suggesting a degree of pleiotropy in the genetic architecture of this character complex. Moreover, we find that cosegregation of shape differences in different elements corresponds to developmental, rather than functional, units.     

A population genomic analysis of species boundaries: neutral processes, adaptive divergence and introgression between two hybridizing plant species

Little is known about the nature of species boundaries between closely related plant species and about the extent of introgression as a consequence of hybridization upon secondary contact. To address these topics we analyzed genome-wide differentiation between two closely related Silene species, Silene latifolia and S. dioica , and assessed the strength of introgression in sympatry. More than 300 AFLP markers were genotyped in three allopatric and three sympatric populations of each species. Outlier analyses were performed separately for sympatric and allopatric populations. Both positive and negative outlier loci were found, indicating that divergent and balancing selection, respectively, have shaped patterns of divergence between the two species. Sympatric populations of the two species were found to be less differentiated genetically than allopatric populations, indicating that hybridization has led to gene introgression. We conclude that differentiation between S. latifolia and S. dioica has been shaped by a combination of introgression and selection. These results challenge the view that species differentiation is a genome-wide phenomenon, and instead support the idea that genomes can be porous and that species differentiation has a genic basis.     

Reconciling extremely strong barriers with high levels of gene exchange in annual sunflowers

In several cases, estimates of gene flow between species appear to be higher than we might predict given the strength of interspecific barriers separating these species pairs. However, as far as we are aware, detailed measurements of reproductive isolation have not previously been compared with a coalescent-based assessment of gene flow. Here, we contrast these two measures in two species of sunflower, Helianthus annuus and H. petiolaris. We quantified the total reproductive barrier strength between these species by compounding the contributions of the following prezygotic and postzygotic barriers: ecogeographic isolation, reproductive asynchrony, niche differentiation, pollen competition, hybrid seed formation, hybrid seed germination, hybrid fertility, and extrinsic postzygotic isolation. From this estimate, we calculated the probability that a reproductively successful hybrid is produced: estimates of P(hyb) range from 10(-4) to 10(-6) depending on the direction of the cross and the degree of independence among reproductive barriers. We then compared this probability with population genetic estimates of the per generation migration rate (m). We showed that the relatively high levels of gene flow estimated between these sunflower species (N(e) m= 0.34-0.76) are mainly due to their large effective population sizes (N(e) > 10(6)). The interspecific migration rate (m) is very small (<10(-7)) and an order of magnitude lower than that expected based on our reproductive barrier strength estimates. Thus, even high levels of reproductive isolation (>0.999) may produce genomic mosaics.     

Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics

This study introduces a new multivariate approach for analyzing the effects of quantitative trait loci (QTL) on shape and demonstrates this method for the mouse mandible. We quantified size and shape with the methods of geometric morphometrics, based on Procrustes superimposition of five morphological landmarks recorded on each mandible. Interval mapping for F(2) mice originating from an intercross of the LG/J and SM/J inbred strains revealed 12 QTL for size, 25 QTL for shape, and 5 QTL for left-right asymmetry. Multivariate ordination of QTL effects by principal component analysis identified two recurrent features of shape variation, which involved the positions of the coronoid and angular processes relative to each other and to the rest of the mandible. These patterns are reminiscent of the knockout phenotypes of a number of genes involved in mandible development, although only a few of these are possible candidates for QTL in our study. The variation of shape effects among the QTL showed no evidence of clustering into distinct groups, as would be expected from theories of morphological integration. Further, for most QTL, additive and dominance effects on shape were markedly different, implying overdominance for specific features of shape. We conclude that geometric morphometrics offers a promising new approach to address problems at the interface of evolutionary and developmental genetics.     

Genetic differences in growth of an invasive tree species

Invasive plants are often more vigorous in their introduced ranges than in their native ranges. This may reflect an innate superiority of plants from some habitats or an escape from their enemies. Another hypothesis proposes that invasive plants evolve increased competitive ability in their introduced range. We present the results of a 14-year common garden experiment with the Chinese Tallow Tree ( Sapium sebiferum ) from its native range (Asia), place of introduction to North America (Georgia) and areas colonized a century later (Louisiana and Texas). Invasive genotypes, especially those from recently colonized areas, were larger than native genotypes and more likely to produce seeds but had lower quality, poorly defended leaves. Our results demonstrate significant post-invasion genetic differences in an invasive plant species. Post-introduction adaptation by introduced plants may contribute to their invasive success and make it difficult to predict problem species.