Admixture mapping of quantitative traits in Populus hybrid zones: power and limitations

Uncovering the genetic architecture of species differences is of central importance for understanding the origin and maintenance of biological diversity. Admixture mapping can be used to identify the number and effect sizes of genes that contribute to the divergence of ecologically important traits, even in taxa that are not amenable to laboratory crosses because of their long generation time or other limitations. Here, we apply admixture mapping to naturally occurring hybrids between two ecologically divergent Populus species. We map quantitative trait loci for eight leaf morphological traits using 77 mapped microsatellite markers from all 19 chromosomes of Populus. We apply multivariate linear regression analysis allowing the modeling of additive and non-additive gene action and identify several candidate genomic regions associated with leaf morphology using an information-theoretic approach. We perform simulation studies to assess the power and limitations of admixture mapping of quantitative traits in natural hybrid populations for a variety of genetic architectures and modes of gene action. Our results indicate that (1) admixture mapping has considerable power to identify the genetic architecture of species differences if sample sizes and marker densities are sufficiently high, (2) modeling of non-additive gene action can help to elucidate the discrepancy between genotype and phenotype sometimes seen in interspecific hybrids, and (3) the genetic architecture of leaf morphological traits in the studied Populus species involves complementary and overdominant gene action, providing the basis for rapid adaptation of these ecologically important forest trees.     

The genetic architecture of quantitative traits: lessons from Drosophila

Understanding the genetic architecture of quantitative traits begins with identifying the genes regulating these traits, mapping the subset of genetically varying quantitative trait loci (QTLs) in natural populations, and pinpointing the molecular polymorphisms defining QTL alleles. Studies in Drosophila have revealed large numbers of pleiotropic genes that interact epistatically to regulate quantitative traits, and large numbers of QTLs with sex-, environment- and genotype-specific effects. Multiple molecular polymorphisms in regulatory regions of candidate genes are often associated with variation for complex traits. These observations offer valuable lessons for understanding the genetic basis of variation for complex traits in other organisms, including humans.     

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.     

QTL mapping for sexually dimorphic fitness-related traits in wild bighorn sheep

Dissecting the genetic architecture of fitness-related traits in wild populations is key to understanding evolution and the mechanisms maintaining adaptive genetic variation. We took advantage of a recently developed genetic linkage map and phenotypic information from wild pedigreed individuals from Ram Mountain, Alberta, Canada, to study the genetic architecture of ecologically important traits (horn volume, length, base circumference and body mass) in bighorn sheep. In addition to estimating sex-specific and cross-sex quantitative genetic parameters, we tested for the presence of quantitative trait loci (QTLs), colocalization of QTLs between bighorn sheep and domestic sheep, and sex × QTL interactions. All traits showed significant additive genetic variance and genetic correlations tended to be positive. Linkage analysis based on 241 microsatellite loci typed in 310 pedigreed animals resulted in no significant and five suggestive QTLs (four for horn dimension on chromosomes 1, 18 and 23, and one for body mass on chromosome 26) using genome-wide significance thresholds (Logarithm of odds (LOD) >3.31 and >1.88, respectively). We also confirmed the presence of a horn dimension QTL in bighorn sheep at the only position known to contain a similar QTL in domestic sheep (on chromosome 10 near the horns locus; nominal P<0.01) and highlighted a number of regions potentially containing weight-related QTLs in both species. As expected for sexually dimorphic traits involved in male-male combat, loci with sex-specific effects were detected. This study lays the foundation for future work on adaptive genetic variation and the evolutionary dynamics of sexually dimorphic traits in bighorn sheep.     

Failure to replicate two mate preference QTLs across multiple strains of Drosophila pseudoobscura

Behavioral genetic mapping studies in model organisms predominantly use crosses originating from a single pair of inbred lines to determine the location of alleles that confer genetic variation in the trait of interest, and they often make sweeping generalizations about the genetic architecture of the trait based on these results. A previous study fine mapped mate preference variation between one pair of Drosophila pseudoobscura lines and identified 2 strong-effect behavioral quantitative trait loci (QTLs). Here, we replicated the previous study's mapping design to examine the extent of variation at these behavioral QTLs across 6 pairs of lines, but we were unable to detect effects of either QTL region in the pairs of lines studied. We suggest that the low-discrimination alleles at these 2 QTLs may occur at low frequency within D. pseudoobscura, although other explanations for the inconsistency are possible. These results underscore the need to examine multiple strains across a species when describing the genetic variation underlying behavioral traits.     

Natural variation in the genetic architecture of a host-parasite interaction in the bumblebee Bombus terrestris

The genetic architecture of fitness-relevant traits in natural populations is a topic that has remained almost untouched by quantitative genetics. Given the importance of parasitism for the host's fitness, we used QTL mapping to study the genetic architecture of traits relevant for host-parasite interactions in the trypanosome parasite, Crithidia bombi and its host, Bombus terrestris. The three traits analysed were the parasite's infection intensity, the strength of the general immune response (measured as the encapsulation of a novel antigen) and body size. The genetic architecture of these traits was examined in three natural, unmanipulated mapping populations of B. terrestris. Our results indicate that the intracolonial phenotypic variation of all three traits is based on a network of QTLs and epistatic interactions. While these networks are similar between mapping populations in complexity and number of QTLs, as well as in their epistatic interactions, the variability in the position of QTL and the interacting loci was high. Only one QTL for body size was plausibly found in at least two populations. QTLs for encapsulation and Crithidia infection intensity were located on the same linkage groups.     

Consequences of near and far between-population crosses for offspring fitness in a rare herb

Crosses between plants from different populations may result in heterosis or outbreeding depression. However, despite its importance for conservation, little is known about the spatial scale over which these effects may arise. To investigate the consequences of between-population crosses at two distinct spatial scales, we conducted reciprocal crosses between four populations from two regions in the rare perennial herb Aster amellus . We assessed seed set and offspring fitness in a common garden experiment. Overall, between-population crosses within regions (10 km) resulted in 8% lower seed set than within-population crosses, while between-region crosses (70 km) resulted in 17% higher seed set than within-population crosses. Moreover, offspring from between-population crosses produced 18% more flower heads than offspring from within-population crosses. We conclude that hybridisation between A. amellus plants from different populations did not lead to immediate outbreeding depression and, thus, could represent a valid conservation option to increase genetic diversity. Moreover, our results suggest that the distance between populations affects the outputs of between-population crosses and therefore needs to be taken into account when promoting gene flow between populations.     

Assessing the complex architecture of polygenic traits in diverged yeast populations

Phenotypic variation arising from populations adapting to different niches has a complex underlying genetic architecture. A major challenge in modern biology is to identify the causative variants driving phenotypic variation. Recently, the baker's yeast, Saccharomyces cerevisiae has emerged as a powerful model for dissecting complex traits. However, past studies using a laboratory strain were unable to reveal the complete architecture of polygenic traits. Here, we present a linkage study using 576 recombinant strains obtained from crosses of isolates representative of the major lineages. The meiotic recombinational landscape appears largely conserved between populations; however, strain-specific hotspots were also detected. Quantitative measurements of growth in 23 distinct ecologically relevant environments show that our recombinant population recapitulates most of the standing phenotypic variation described in the species. Linkage analysis detected an average of 6.3 distinct QTLs for each condition tested in all crosses, explaining on average 39% of the phenotypic variation. The QTLs detected are not constrained to a small number of loci, and the majority are specific to a single cross-combination and to a specific environment. Moreover, crosses between strains of similar phenotypes generate greater variation in the offspring, suggesting the presence of many antagonistic alleles and epistatic interactions. We found that subtelomeric regions play a key role in defining individual quantitative variation, emphasizing the importance of the adaptive nature of these regions in natural populations. This set of recombinant strains is a powerful tool for investigating the complex architecture of polygenic traits.     

Within- and between-population variation for resistance of Daphnia magna to the bacterial endoparasite Pasteuria ramosa

Genetic variation among hosts for resistance to parasites is an important assumption underlying evolutionary theory of host and parasite evolution. Using the castrating bacterial parasite Pasteuria ramosa and its cladoceran host Daphnia magna, we examined both within- and between-population genetic variation for resistance. First, we tested hosts from four populations for genetic variation for resistance to three parasite isolates. Allozyme analysis revealed significant host population divergence and that genetic distance corresponds to geographic distance. Host and parasite fitness components showed strong genetic differences between parasite isolates for host population by parasite interactions and for clones within populations, whereas host population effects were significant for only a few traits. In a second experiment we tested explicitly for within-population differences in variation for resistance by challenging nine host clones from a single population with four different parasite spore doses. Strong clone and dose effects were evident. More susceptible clones also suffered higher costs once infected. The results indicate that within-population variation for resistance is high relative to between-population variation. We speculate that P. ramosa adapts to individual host clones rather than to its host population.     

THE DIVERGENCE OF NEUTRAL QUANTITATIVE CHARACTERS AMONG PARTIALLY ISOLATED POPULATIONS

The neutral model of phenotypic evolution has yielded several simple predictions about the long-term rates of between-population divergence of polygenic traits and about the equilibrium level of within-population variance when mutation and random genetic drift are the sole evolutionary forces. These conclusions must be modified if populations are only partially isolated. A quantitative model is presented for the development of within- and between-population variance for neutral quantitative characters in pairs of populations with arbitrary effective sizes and migration rates. Both the variance in the base population and subsequent variance generated by mutation are considered, and several dynamical and equilibrium properties are shown to be adequately described by simple approximations. The resultant formulations provide some insight into the sensitivity of measures of morphological distance to gene flow, the necessity of isolation for the accumulation of variation between incipient species, and the consequences of gene flow into captive populations of endangered species.     

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.     

The relationship between parental genetic or phenotypic divergence and progeny variation in the maize nested association mapping population

Appropriate selection of parents for the development of mapping populations is pivotal to maximizing the power of quantitative trait loci detection. Trait genotypic variation within a family is indicative of the family's informativeness for genetic studies. Accurate prediction of the most useful parental combinations within a species would help guide quantitative genetics studies. We tested the reliability of genotypic and phenotypic distance estimators between pairs of maize inbred lines to predict genotypic variation for quantitative traits within families derived from biparental crosses. We developed 25 families composed of ~200 random recombinant inbred lines each from crosses between a common reference parent inbred, B73, and 25 diverse maize inbreds. Parents and families were evaluated for 19 quantitative traits across up to 11 environments. Genetic distances (GDs) among parents were estimated with 44 simple sequence repeat and 2303 single-nucleotide polymorphism markers. GDs among parents had no predictive value for progeny variation, which is most likely due to the choice of neutral markers. In contrast, we observed for about half of the traits measured a positive correlation between phenotypic parental distances and within-family genetic variance estimates. Consequently, the choice of promising segregating populations can be based on selecting phenotypically diverse parents. These results are congruent with models of genetic architecture that posit numerous genes affecting quantitative traits, each segregating for allelic series, with dispersal of allelic effects across diverse genetic material. This architecture, common to many quantitative traits in maize, limits the predictive value of parental genotypic or phenotypic values on progeny variance.     

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Abstract Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with Q _ST generally exceeding G _ST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of Q _ST and G _ST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria , regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex .     

Quantitative trait loci controlling leaf venation in Arabidopsis

Leaf veins provide the mechanical support and are responsible for the transport of nutrients and water to the plant. High vein density is a prerequisite for plants to have C4 photosynthesis. We investigated the genetic variation and genetic architecture of leaf venation traits within the species Arabidopsis thaliana using natural variation. Leaf venation traits, including leaf vein density (LVD) were analysed in 66 worldwide accessions and 399 lines of the multi‐parent advanced generation intercross population. It was shown that there is no correlation between LVD and photosynthesis parameters within A. thaliana. Association mapping was performed for LVD and identified 16 and 17 putative quantitative trait loci (QTLs) in the multi‐parent advanced generation intercross and worldwide sets, respectively. There was no overlap between the identified QTLs suggesting that many genes can affect the traits. In addition, linkage mapping was performed using two biparental recombinant inbred line populations. Combining linkage and association mapping revealed seven candidate genes. For one of the candidate genes, RCI2c, we demonstrated its function in leaf venation patterning.     

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

BackgroundIncreasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.ResultsThe populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.ConclusionsThis study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.     

Exploring the genetic characteristics of 93-11 and Nipponbare recombination inbred lines based on the GoldenGate SNP assay

Understanding genetic characteristics in rice populations will facilitate exploring evolutionary mechanisms and gene cloning. Numerous molecular markers have been utilized in linkage map construction and quantitative trait locus (QTL) mappings. However, segregation-distorted markers were rarely considered, which prevented understanding genetic characteristics in many populations. In this study, we designed a 384-marker GoldenGate SNP array to genotype 283 recombination inbred lines (RILs) derived from 93-11 and Nipponbare Oryza sativa crosses. Using 294 markers that were highly polymorphic between parents, a linkage map with a total genetic distance of 1,583.2 cM was constructed, including 231 segregation-distorted markers. This linkage map was consistent with maps generated by other methods in previous studies. In total, 85 significant quantitative trait loci (QTLs) with phenotypic variation explained (PVE) values≥5% were identified. Among them, 34 QTLs were overlapped with reported genes/QTLs relevant to corresponding traits, and 17 QTLs were overlapped with reported sterility-related genes/QTLs. Our study provides evidence that segregation-distorted markers can be used in linkage map construction and QTL mapping. Moreover, genetic information resulting from this study will help us to understand recombination events and segregation distortion. Furthermore, this study will facilitate gene cloning and understanding mechanism of inter-subspecies hybrid sterility and correlations with important agronomic traits in rice.     

Bayesian analysis for genetic architecture of dynamic traits

The dissection of the genetic architecture of quantitative traits, including the number and locations of quantitative trait loci (QTL) and their main and epistatic effects, has been an important topic in current QTL mapping. We extend the Bayesian model selection framework for mapping multiple epistatic QTL affecting continuous traits to dynamic traits in experimental crosses. The extension inherits the efficiency of Bayesian model selection and the flexibility of the Legendre polynomial model fitting to the change in genetic and environmental effects with time. We illustrate the proposed method by simultaneously detecting the main and epistatic QTLs for the growth of leaf age in a doubled-haploid population of rice. The behavior and performance of the method are also shown by computer simulation experiments. The results show that our method can more quickly identify interacting QTLs for dynamic traits in the models with many numbers of genetic effects, enhancing our understanding of genetic architecture for dynamic traits. Our proposed method can be treated as a general form of mapping QTL for continuous quantitative traits, being easier to extend to multiple traits and to a single trait with repeat records.     

Two reproductive traits show contrasting genetic architectures in Plantago lanceolata

In many species, temperature‐sensitive phenotypic plasticity (i.e., an individual's phenotypic response to temperature) displays a positive correlation with latitude, a pattern presumed to reflect local adaptation. This geographical pattern raises two general questions: (a) Do a few large‐effect genes contribute to latitudinal variation in a trait? (b) Is the thermal plasticity of different traits regulated pleiotropically? To address the questions, we crossed individuals of Plantago lanceolata derived from northern and southern European populations. Individuals naturally exhibited high and low thermal plasticity in floral reflectance and flowering time. We grew parents and offspring in controlled cool‐ and warm‐temperature environments, mimicking what plants would encounter in nature. We obtained genetic markers via genotype‐by‐sequencing, produced the first recombination map for this ecologically important nonmodel species, and performed quantitative trait locus (QTL) mapping of thermal plasticity and single‐environment values for both traits. We identified a large‐effect QTL that largely explained the reflectance plasticity differences between northern and southern populations. We identified multiple smaller‐effect QTLs affecting aspects of flowering time, one of which affected flowering time plasticity. The results indicate that the genetic architecture of thermal plasticity in flowering is more complex than for reflectance. One flowering time QTL showed strong cytonuclear interactions under cool temperatures. Reflectance and flowering plasticity QTLs did not colocalize, suggesting little pleiotropic genetic control and freedom for independent trait evolution. Such genetic information about the architecture of plasticity is environmentally important because it informs us about the potential for plasticity to offset negative effects of climate change.     

Towards the genetic architecture of seed lipid biosynthesis and accumulation in Arabidopsis thaliana

We report the quantitative genetic analysis of seed oil quality and quantity in six Arabidopsis thaliana recombinant inbred populations, in which the parent accessions were from diverse geographical origins, and were selected on the basis of variation for seed oil content and lipid composition. Although most of the biochemical steps involved in lipid biosynthesis are known and the key genes have been identified, the regulation of the processes that results in the final oil composition and total amount is not understood. By using physically anchored markers it was possible to compare results across populations. A total of 219 quantitative trait loci (QTLs) were identified, of which 81 were significant at P<0.001. Some of these colocalise with QTLs identified previously, but many novel QTLs were also identified. The results highlight the importance of studying traits in multiple populations, which will lead to a better understanding of the contribution that natural variation makes to the genetic architecture of a phenotype.