Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Evidence of neutral and adaptive genetic divergence between European trout populations sampled along altitudinal gradients

Species with a wide geographical distribution are often composed of distinct subgroups which may be adapted to their local environment. European trout (Salmo trutta species complex) provide an example of such a complex consisting of several genetically and ecologically distinct forms. However, trout populations are strongly influenced by human activities, and it is unclear to what extent neutral and adaptive genetic differences have persisted. We sampled 30 Swiss trout populations from heterogeneous environments along replicated altitudinal gradients in three major European drainages. More than 850 individuals were genotyped at 18 microsatellite loci which included loci diagnostic for evolutionary lineages and candidate markers associated with temperature tolerance, reproductive timing and immune defence. We find that the phylogeographic structure of Swiss trout populations has not been completely erased by stocking. Distinct genetic clusters corresponding to the different drainages could be identified, although nonindigenous alleles were clearly present, especially in the two Mediterranean drainages. We also still detected neutral genetic differentiation within rivers which was often associated with the geographical distance between populations. Five loci showed evidence of divergent selection between populations with several drainage-specific patterns. Lineage-diagnostic markers, a marker linked to a quantitative trait locus for upper temperature tolerance in other salmonids and a marker linked to the major histocompatibility class I gene were implicated in local adaptation and some patterns were associated with altitude. In contrast, tentative evidence suggests a signal of balancing selection at a second immune relevant gene (TAP2). Our results confirm the persistence of both neutral and potentially adaptive genetic differences between trout populations in the face of massive human-mediated dispersal.

Evolution of dominance under frequency-dependent intraspecific competition

A population-genetic analysis is performed of a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under frequency-dependent disruptive selection caused by intraspecific competition for a continuum of resources. The modifier locus determines the degree of dominance at the trait level. We establish the conditions when a modifier allele can invade and when it becomes fixed if sufficiently frequent. In general, these are not equivalent because an unstable internal equilibrium may exist and the condition for successful invasion of the modifier is more restrictive than that for eventual fixation from already high frequency. However, successful invasion implies global fixation, i.e., fixation from any initial condition. Modifiers of large effect can become fixed, and also invade, in a wider parameter range than modifiers of small effect. We also study modifiers with a direct, frequency-independent deleterious fitness effect. We show that they can invade if they induce a sufficiently high level of dominance and if disruptive selection on the ecological trait is strong enough. For deleterious modifiers, successful invasion no longer implies global fixation because they can become stuck at an intermediate frequency due to a stable internal equilibrium. Although the conditions for invasion and for fixation if sufficiently frequent are independent of the linkage relation between the two loci, the rate of spread depends strongly on it. The present study provides further support to the view that evolution of dominance may be an efficient mechanism to remove unfit heterozygotes that are maintained by balancing selection. It also demonstrates that an invasion analysis of mutants of very small effect is insufficient to obtain a full understanding of the evolutionary dynamics under frequency-dependent selection.

The genetics of reproductive organ morphology in two Petunia species with contrasting pollination syndromes

Main conclusion Switches between pollination syndromes have happened frequently during angiosperm evolution. Using QTL mapping and reciprocal introgressions, we show that changes in reproductive organ morphology have a simple genetic basis. In animal-pollinated plants, flowers have evolved to optimize pollination efficiency by different pollinator guilds and hence reproductive success. The two Petunia species, P. axillaris and P. exserta, display pollination syndromes adapted to moth or hummingbird pollination. For the floral traits color and scent, genetic loci of large phenotypic effect have been well documented. However, such large-effect loci may be typical for shifts in simple biochemical traits, whereas the evolution of morphological traits may involve multiple mutations of small phenotypic effect. Here, we performed a quantitative trait locus (QTL) analysis of floral morphology, followed by an in-depth study of pistil and stamen morphology and the introgression of individual QTL into reciprocal parental backgrounds. Two QTLs, on chromosomes II and V, are sufficient to explain the interspecific difference in pistil and stamen length. Since most of the difference in organ length is caused by differences in cell number, genes underlying these QTLs are likely to be involved in cell cycle regulation. Interestingly, conservation of the locus on chromosome II in a different P. axillaris subspecies suggests that the evolution of organ elongation was initiated on chromosome II in adaptation to different pollinators. We recently showed that QTLs for pistil and stamen length on chromosome II are tightly linked to QTLs for petal color and volatile emission. Linkage of multiple traits will enable major phenotypic change within a few generations in hybridizing populations. Thus, the genomic architecture of pollination syndromes in Petunia allows for rapid responses to changing pollinator availability.

Molecular characterization of five porcine candidate genes for drip loss in pork

Drip loss is the loss of fluid from a piece of meat without mechanical force and represents an important meat quality trait. Previous work revealed a quantitative trait locus (QTL) for drip loss in pork in an experimental Duroc x Pietrain (DUPI) F2 family on SSC 5. Based on functional data indicating their possible involvement in water holding capacity and their expression in skeletal muscle, we selected five positional candidates (ACO2, ADSL, CBY1, KCNJ4, PLA2AG6) out of 130 predicted genes in the QTL interval for further analysis. We performed a mutation analysis of all coding exons and discovered 204 polymorphisms. We genotyped 39 single nucleotide polymorphisms (SNPs) in 192 Pietrain pigs with extreme drip loss phenotypes and detected a possible association with drip loss for one non-coding SNP in the ADSL gene (ss107793818, p(raw) = 0.021). Correspondingly, ADSL diplotypes were associated with drip loss and pH1 of M. longissimus dorsi. However, after correction for multiple testing, none of the tested SNPs were significantly associated with drip loss. One possible explanation for these results is that one of the QTL-alleles from the experimental DUPI family may be fixed or nearly fixed in the tested Pietrain population.

Molecular characterization of the equine collagen, type IX, alpha 2 (COL9A2) gene on horse chromosome 2p16-->p15

The mammalian collagen, type IX, alpha 2 gene (COL9A2) encodes the alpha-2 chain of type IX collagen and is located on horse chromosome 2p16-->p14 harbouring a quantitative trait locus for osteochondrosis. We isolated a bacterial artificial chromosome (BAC) clone containing the equine COL9A2 gene and determined the complete genomic sequence of this gene. Cloning and characterization of equine COL9A2 revealed that the equine gene consists of 32 exons spanning approximately 15 kb. The COL9A2 transcript encodes a single protein of 688 amino acids. Thirty two single nucleotide polymorphisms (SNPs) equally distributed in the gene were detected in a mutation scan of eight unrelated Hanoverian warmblood stallions, including one SNP that affects the amino acid sequence of COL9A2. Comparative analyses between horse, human, mouse and rat indicate that the chromosomal location of equine COL9A2 is in agreement with known chromosomal synteny relationships. The comparison of the gene structure and transcript revealed a high degree of conservation towards the other mammalian COL9A2 genes. We chose three informative SNPs for association and linkage disequilibrium tests in three to five paternal half-sib families of Hanoverian warmblood horses consisting of 44 to 75 genotyped animals. The test statistics did not reach the significance threshold of 5% and so we could not show an association of COL9A2 with equine osteochondrosis.

DNA methylation patterns at the IGF2-H19 locus in sperm of Swiss Landrace and Swiss Large White boars

DNA methylation patterns at the IGF2-H19 locus were investigated in sperm DNA from Swiss Landrace (SL) and Swiss Large White (LW) boars. The putative IGF2 differentially methylated regions (DMR) 0, 1 and 2, a quantitative trait nucleotide (QTN) region in the intron 3 and a CpG island in the intron 4 of the IGF2 gene as well as three regions around porcine CTCF binding sites within the H19 differentially methylated domain (DMD) were selected for the DNA methylation analysis. In both breeds putative IGF2 DMR0, 1, 2 and H19 DMD were hypermethylated. Significant differences in DNA methylation content were found between the two breeds in the two DMD regions proximal to the H19 gene. The IGF2 QTN region and the CpG island in the IGF2 intron 4 were hypomethylated in sperm DNA of both breeds. The methylation analysis revealed significantly more methylated CpG sites in the intron 4 of sperm from the LW breed than in that from SL. No difference was found in global DNA methylation between the two breeds. These results indicate differences in DNA methylation patterns between breeds and it remains to be established whether variation in DNA methylation patterns impacts on phenotypic traits.

Tilapia sex determination: Where temperature and genetics meet

This review deals with the complex sex determining system of Nile tilapia, Oreochromis niloticus, governed by the interactions between a genetic determination and the influence of temperature, shown in both domestic and wild populations. Naturally sex reversed individuals are strongly suggested in two wild populations. This can be due to the masculinising temperatures which some fry encounter during their sex differentiation period when they colonise shallow waters, and/or to the influence of minor genetic factors. Differences regarding a) thermal responsiveness of sex ratios between and within Nile tilapia populations, b) maternal and paternal effects on temperature dependent sex ratios and c) nearly identical results in offspring of repeated matings, demonstrate that thermosensitivity is under genetic control. Selection experiments to increase the thermosensitivity revealed high responses in the high and low sensitive lines. The high-line showed ~ 90% males after 2 generations of selection whereas the weakly sensitive line had 54% males. This is the first evidence that a surplus of males in temperature treated groups can be selected as a quantitative trait. Expression profiles of several genes (Cyp19a, Foxl2, Amh, Sox9a,b) from the gonad and brain were analysed to define temperature action on the sex determining/differentiating cascade in tilapia. The coexistence of GSD and TSD is discussed.

EVOLUTION OF DOMINANCE UNDER FREQUENCY-DEPENDENT INTRASPECIFIC COMPETITION IN AN ASSORTATIVELY MATING POPULATION

We study the evolution of higher levels of dominance as a response to negative frequency-dependent selection. In contrast to previous studies, we focus on the effect of assortative mating on the evolution of dominance under frequency-dependent intraspecific competition. We analyze a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under a mixture of frequency-independent stabilizing selection, density-dependent selection, and frequency-dependent selection caused by intraspecific competition for a continuum of resources. The second (modifier) locus determines the degree of dominance at the trait level. Additionally, the population mates assortatively with respect to similarities in the ecological trait. Our analysis shows that the parameter region in which dominance can be established decreases if small levels of assortment are introduced. In addition, the degree of dominance that can be established also decreases. In contrast, if assortment is intermediate, sexual selection for extreme types can be established, which leads to evolution of higher levels of dominance than under random mating. For modifiers with large effects, intermediate levels of assortative mating are most favorable for the evolution of dominance. For large modifiers, the speed of fixation can even be higher for intermediate levels of assortative mating than for random mating.

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G x E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G x E interaction is not necessary for the maintenance of more than two alleles in two demes.

Whole genome scan identifies several chromosomal regions linked to equine sarcoids.

Despite the evidence for a genetic predisposition to develop equine sarcoids (ES), no whole genome scan for ES has been performed to date. The objective of this explorative study was to identify chromosome regions associated with ES. The studied population was comprised of two half-sibling sire families, involving a total of 222 horses. Twenty-six of these horses were affected with ES. All horses had been previously genotyped with 315 microsatellite markers. Quantitative trait locus (QTL) signals were suggested where the F statistic exceeded chromosome-wide significance at P < 0.05. The QTL analyses revealed significant signals reaching P < 0.05 on equine chromosome (ECA) 20, 23 and 25, suggesting a polygenic character for this trait. The candidate regions identified on ECA 20, 23 and 25 include genes regulating virus replication and host immune response. Further investigation of the chromosome regions associated with ES and of genes potentially responsible for the development of ES could form the basis for early identification of susceptible animals, breeding selection or the development of new therapeutic targets.

Natural Variation in Maize Aphid Resistance Is Associated with 2,4-Dihydroxy-7-Methoxy-1,4-Benzoxazin-3-One Glucoside Methyltransferase Activity

Plants differ greatly in their susceptibility to insect herbivory, suggesting both local adaptation and resistance tradeoffs. We used maize (Zea mays) recombinant inbred lines to map a quantitative trait locus (QTL) for the maize leaf aphid (Rhopalosiphum maidis) susceptibility to maize Chromosome 1. Phytochemical analysis revealed that the same locus was also associated with high levels of 2-hydroxy-4,7-dimethoxy-1,4-benzoxazin-3-one glucoside (HDMBOA-Glc) and low levels of 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one glucoside (DIMBOA-Glc). In vitro enzyme assays with candidate genes from the region of the QTL identified three O-methyltransferases (Bx10a-c) that convert DIMBOA-Glc to HDMBOA-Glc. Variation in HDMBOA-Glc production was attributed to a natural CACTA family transposon insertion that inactivates Bx10c in maize lines with low HDMBOA-Glc accumulation. When tested with a population of 26 diverse maize inbred lines, R. maidis produced more progeny on those with high HDMBOA-Glc and low DIMBOA-Glc. Although HDMBOA-Glc was more toxic to R. maidis than DIMBOA-Glc in vitro, BX10c activity and the resulting decline of DIMBOA-Glc upon methylation to HDMBOA-Glc were associated with reduced callose deposition as an aphid defense response in vivo. Thus, a natural transposon insertion appears to mediate an ecologically relevant trade-off between the direct toxicity and defense-inducing properties of maize benzoxazinoids.

No evidence for a genetic association between female mating preference and male secondary sexual trait in a Lake Victoria cichlid fish

Sexual selection by female mating preference for male nuptial coloration has been suggested as a driving force in the rapid speciation of Lake Victoria cichlid fish. This process could have been facilitated or accelerated by genetic associations between female preference loci and male coloration loci. Preferences, as well as coloration, are heritable traits and are probably determined by more than one gene. However, little is known about potential genetic associations between these traits. In turbid water, we found a population that is variable in male nuptial coloration from blue to yellow to red. Males at the extreme ends of the phenotype distribution resemble a reproductively isolated species pair in clear water that has diverged into one species with blue-grey mates and one species with bright red males. Females of the turbid water population vary in mating preference coinciding with the male phenotype distribution. For the current study, these females were mated to blue males. We measured the coloration of the sires and male offspring. Parents-offspring regression showed that the sires did not affect male offspring coloration, which confirms earlier findings that the blue species breeds true. In contrast, male offspring coloration was determined by the identity of the dams, which suggests that there is heritable variation in male color genes between females. However, we found that mating preferences of the dams were not correlated with male offspring coloration. Thus, there is no evidence for strong genetic linkage between mating preference and the preferred trait in this population [Current Zoology 56 (1): 57-64 2010].