Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

Disentangling the effects of breakdown of self-incompatibility and transition to selfing in North American Arabidopsis lyrata

A breakdown of self‐incompatibility (SI) followed by a shift to selfing is commonly observed in the evolution of flowering plants. Both are expected to reduce the levels of heterozygosity and genetic diversity. However, breakdown of SI should most strongly affect the region of the SI locus (S‐locus) because of the relaxation of balancing selection that operates on a functional S‐locus, and a potential selective sweep. In contrast, a transition to selfing should affect the whole genome. We set out to disentangle the effects of breakdown of SI and transition to selfing on the level and distribution of genetic diversity in North American populations of Arabidopsis lyrata. Specifically, we compared sequence diversity of loci linked and unlinked to the S‐locus for populations ranging from complete selfing to fully outcrossing. Regardless of linkage to the S‐locus, heterozygosity and genetic diversity increased with population outcrossing rate. High heterozygosity of self‐compatible individuals in outcrossing populations suggests that SI is not the only factor preventing the evolution of self‐fertilization in those populations. There was a strong loss of diversity in selfing populations, which was more pronounced at the S‐locus. In addition, selfing populations showed an accumulation of derived mutations at the S‐locus. Our results provide evidence that beyond the genome‐wide consequences of the population bottleneck associated with the shift to selfing, the S‐locus of A. lyrata shows a specific signal either reflecting the relaxation of balancing selection or positive selection.     

Multiple origins of self-fertilization in tristylous Eichhornia paniculata (Pontederiaceae): Inferences from style morph and isozyme variation

Populations of Eichhornia paniculata (Pontederiaceae) exhibit a wide range of mating systems, from predominant outcrossing to high levels of self-fertilization. The origin of self-fertilization in this tristylous species is associated with the loss of style-length morphs from populations and the spread of self-pollinating, floral variants. We examined geographic variation in style morph and allozyme frequencies to determine whether the loss of style morphs and transition to selfing could have multiple origins in E. paniculata. Surveys of floral variation in 167 populations from six states in northeastern Brazil revealed that at least one style morph was absent from 29.3%. Non-trimorphic populations occurred in all states and ranged in frequency from 9% in Ceará to 68% in Alagoas. Selfing variants occurred in 8.5% and 55% of trimorphic and non-trimorphic populations, respectively, and were distributed among five of six states with primary concentrations in Alagoas and Pernambuco. A comparison of electrophoretic variation at 24 isozyme loci in 28 trimorphic, 13 dimorphic and 3 monomorphic populations indicated that non-trimorphic populations contained 84% of the allelic variation present in trimorphic populations and were markedly differentiated from one another. Analyses of genetic distance and the distribution of rare alleles indicated that non-trimorphic populations were often more similar to neighbouring trimorphic populations than to one another. Populations with selfing variants occurred at low frequency in three genetically distinct parts of the range. These results, in combination with genetic and morphological evidence suggest that style morphs are lost repeatedly from populations of E. paniculata and that selfing variants may have originated on at least three separate occasions in northeastern Brazil.     

Genetic analysis of traits distinguishing outcrossing and self-pollinating forms of currant tomato,Lycopersicon pimpinellifolium (Jusl.) Mill

The evolution of inbreeding is common throughout the angiosperms, although little is known about the developmental and genetic processes involved. Lycopersicon pimpinellifolium (currant tomato) is a self-compatible species with variation in outcrossing rate correlated with floral morphology. Mature flowers from inbreeding and outcrossing populations differ greatly in characters affecting mating behavior (petal, anther, and style lengths); other flower parts (sepals, ovaries) show minimal differences. Analysis of genetic behavior, including quantitative trait locus (QTL) mapping, was performed on representative selfing and outcrossing plants derived from two contrasting natural populations. Six morphological traits were analyzed: flowers per inflorescence; petal, anther, and style lengths; and lengths of the fertile and sterile portions of anthers. All traits were smaller in the selfing parent and had continuous patterns of segregation in the F(2). Phenotypic correlations among traits were all positive, but varied in strength. Quantitative trait locus mapping was done using 48 RFLP markers. Five QTL total were found involving four of the six traits: total anther length, anther sterile length, style length, and flowers per inflorescence. Each of these four traits had a QTL of major (>25%) effect on phenotypic variance.     

Genetic linkage analysis of X-ray hypersensitivity in the LEC mutant rat

The LEC rat has been reported to exhibit X-ray hypersensitivity and deficiency in DNA double-strand break (DSB) repair. The present study was performed to map the locus responsible for this phenotype, the xhs (X-ray hypersensitivity), as the first step in identifying the responsible gene. Analysis of the progeny of (BN × LEC)F₁× LEC backcrosses indicated that the X-ray hypersensitive phenotype was controlled by multiple genetic loci in contrast to the results reported previously. Quantitative trait loci (QTL) linkage analysis revealed two responsible loci located on Chromosomes (Chr) 4 and 1. QTL on Chr 4 exhibited very strong linkage to the X-ray hypersensitive phenotype, while QTL on Chr 1 showed weak linkage. The Rad52 locus, mutation of which results in hypersensitivity to ionizing radiation and impairment of DNA DSB repair in yeast, was reported to be located on the synteneic regions of mouse Chr 6 and human Chr 12. However, mapping of the rat Rad52 locus indicated that it was located 23 cM distal to the QTL on Chr 4. Furthermore, none of the radio-sensitivity-related loci mapped previously in the rat chromosome were identical to the QTL on Chrs 4 and 1 in the LEC rat. Thus, it seems that X-ray hypersensitivity in the LEC rat is caused by mutation(s) in as-yet-undefined genes. Received: 14 February 2000 / Accepted: 17 May 2000     

A quantitative trait locus for body fat on chromosome 1q43 in French Canadians: linkage and association studies

Objective: To explore a quantitative trait locus (QTL) on human chromosome 1q affecting BMI, adiposity, and fat‐free mass phenotypes in the Quebec Family Study cohort. Research Methods and Procedures: Non‐parametric sibpair and variance component linkage analyses and family‐based association studies were performed with a dense set of chromosome 1q43 microsatellites and single‐nucleotide polymorphism markers in 885 adult individuals. Results: Linkage was observed between marker D1S184 and BMI (p = 0.0004) and with body fat mass or percentage body fat (p ≤ 0.0003), but no linkage was detected with fat‐free mass. Furthermore, significant linkages (p < 0.0001) were achieved with subsamples of sibpairs at both ends of phenotype distributions. Association studies with quantitative transmission disequilibrium tests refined the linkage to a region overlapping the regulator of G‐protein signaling 7 (RGS7) gene and extending to immediate upstream gene loci. Discussion: The present study indicates that the QTL on chromosome 1q43 specifically affects total adiposity and provides a genetic mapping framework for the dissection of this adiposity locus.     

High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. et Zucc)

High-density genetic map is a valuable tool for fine mapping locus controlling a specific trait especially for perennial woody plants. In this study, we firstly constructed a high-density genetic map of mei (Prunus mume) using SLAF markers, developed by specific locus amplified fragment sequencing (SLAF-seq). The linkage map contains 8,007 markers, with a mean marker distance of 0.195 cM, making it the densest genetic map for the genus Prunus. Though weeping trees are used worldwide as landscape plants, little is known about weeping controlling gene(s) (Pl). To test the utility of the high-density genetic map, we did fine-scale mapping of this important ornamental trait. In total, three statistic methods were performed progressively based on the result of inheritance analysis. Quantitative trait loci (QTL) analysis initially revealed that a locus on linkage group 7 was strongly responsible for weeping trait. Mutmap-like strategy and extreme linkage analysis were then applied to fine map this locus within 1.14 cM. Bioinformatics analysis of the locus identified some candidate genes. The successful localization of weeping trait strongly indicates that the high-density map constructed using SLAF markers is a worthy reference for mapping important traits for woody plants.     

QTL analysis of heterostyly in Primula sieboldii and its application for morph identification in wild populations

Background and Aims

Primula sieboldii is a perennial clonal herb that is distributed around the Sea of Japan and is endangered in Japan. Its breeding system is characterized by heteromorphic self-incompatibility, and the morph ratio within a population is very important for reproductive success. The aims of this study were to construct a linkage map, map the S locus as a qualitative trait and quantitative trait loci (QTLs) for floral morphological traits related to heterostyly, and predict the morph type in wild populations by using molecular markers for devising a conservation strategy.

Methods

A linkage map was constructed with 126 markers. The QTLs for four floral traits and the S locus were mapped. Using the genotypes of loci that were located near both the S locus and the QTLs with large effects, morphs of 59 wild genets were predicted.

Key Results

The linkage map consisted of 14 linkage groups (LGs). The S locus was mapped to LG 7. Major QTLs for stigma and anther heights were detected in the same region as the S locus. These QTLs exhibited high logarithm of the odds scores and explained a high percentage of the phenotypic variance (>85 %). By analysing these two traits within each morph, additional QTLs for each trait were detected. Using the four loci linked to the S locus, the morphs of 43 genets in three wild populations could be predicted.

Conclusions

This is the first report of a linkage map and QTL analysis for floral morphology related to heterostyly in P. sieboldii. Floral morphologies related to heterostyly are controlled by the S locus in LG 7 and by several QTLs in other LGs. Additionally, this study showed that molecular markers are effective tools for investigating morph ratios in a population containing the non-flowering individuals or during the non-flowering seasons.     

Replication and refinement of a quantitative trait locus influencing milk protein percentage on ovine chromosome 3

A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half‐sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (P_c < 0.001; P_exp < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E‐11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.     

Evolutionary pathways to self-fertilization in a tristylous plant species

Evolutionary transitions from outcrossing to selfing occur commonly in heterostylous genera. The morphological polymorphisms that characterize heterostyly provide opportunities for different pathways for selfing to evolve. Here, we investigate the origins and pathways by which selfing has evolved in tristylous Eichhornia paniculata by providing new evidence based on morphology, DNA sequences and genetic analysis. The primary pathway from outcrossing to selfing involves the stochastic loss of the short-styled morph (S-morph) from trimorphic populations, followed by the spread of selfing variants of the mid-styled morph (M-morph). However, the discovery of selfing variants of the long-styled morph (L-morph) in Central America indicates a secondary pathway and distinct origin for selfing. Comparisons of multi-locus nucleotide sequences from 27 populations sampled from throughout the geographical range suggest multiple transitions to selfing. Genetic analysis of selfing variants of the L- and M-morphs demonstrates recessive control of the loss of herkogamy, although the number of factors appears to differ between the forms. Early stages in the establishment of selfing involve developmental instability in the formation of flowers capable of autonomous self-pollination. The relatively simple genetic control of herkogamy reduction and frequent colonizing episodes may often create demographic conditions favouring transitions to selfing in E. paniculata .     

Genetic basis of tristyly in tetraploid Oxalis alpina (Oxalidaceae)

The inheritance of style‐morphs was investigated in tetraploid populations of tristylous Oxalis alpina (Oxalidaceae) to determine if alleles controlling style‐morphs are expressed at duplicated loci. In tetraploid populations, a dominant S allele leads to expression of the short‐styled phenotype at the short/non‐short locus and is epistatic to the M allele at the mid/long locus. The M allele results in expression of the mid‐styled phenotype but only if the S allele is absent. Long‐styled morphs are homozygous recessive at the short and mid loci. Test crosses of many tetraploid short‐styled individuals resulted in segregations of short‐, mid‐ and long‐styled individuals which, because of linkage between the short and mid loci, can only occur with polyploidy and expression of alleles at duplicated loci. Segregation patterns from three crosses suggest the possibility of disomic inheritance via preferential pairing of chromosomes in tetraploid populations of O. alpina. Segregation patterns in the progeny of mid‐styled individuals indicated that only a few individuals had more than one copy of the M allele, despite the potential for accumulation of M alleles via self‐fertilization of partially self‐compatible mid‐styled morphs in some populations. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 308–318.     

Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of Arabidopsis thaliana to drought

Drought stress was imposed on two sets of Arabidopsis thaliana genotypes grown in sand under short‐day conditions and analysed for several shoot and root growth traits. The response to drought was assessed for quantitative trait locus (QTL) mapping in a genetically diverse set of Arabidopsis accessions using genome‐wide association (GWA) mapping, and conventional linkage analysis of a recombinant inbred line (RIL) population. Results showed significant genotype by environment interaction (G×E) for all traits in response to different watering regimes. For the RIL population, the observed G×E was reflected in 17 QTL by environment interactions (Q×E), while 17 additional QTLs were mapped not showing Q×E. GWA mapping identified 58 single nucleotide polymorphism (SNPs) associated with loci displaying Q×E and an additional 16 SNPs associated with loci not showing Q×E. Many candidate genes potentially underlying these loci were suggested. The genes for RPS3C and YLS7 were found to contain conserved amino acid differences when comparing Arabidopsis accessions with strongly contrasting drought response phenotypes, further supporting their candidacy. One of these candidate genes co‐located with a QTL mapped in the RIL population.     

Isolation and characterization of 31 polymorphic microsatellite markers in barfin flounder (Verasper moseri) and the cross-species amplification in spotted halibut (Verasper variegatus)

The construction of a genetic linkage map and evaluation of population genetic diversity both require large numbers of polymorphic molecular markers. In the present study, we report 31 polymorphic microsatellite markers, of which 28 were isolated from two repeat-enriched libraries constructed from genomic DNA, and three were detected in two genes (MCH-R1 and MCH-R2) of barfin flounder (Verasper moseri). A total of 94 alleles were detected with an average of 3.0 alleles per locus. The number of alleles, observed and expected heterozygosity per locus ranged from two to six, from 0.30 to 1.00 and from 0.33 to 0.78, respectively. Three loci significantly deviated from Hardy-Weinberg equilibrium after Bonferroni correction (P < 0.0016) and no significant linkage disequilibrum between pairs of loci was found. Cross-species amplification of these markers was evaluated using the closely related species spotted halibut (Verasper variegatus). This study will potentially be useful for stock management, constructing of a genetic linkage map, mapping economically important quantitative trait loci (QTL), and studying the population genetic diversity of barfin flounder (Verasper moseri).     

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome—that is, the genetic architecture of reproductive isolation—can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male‐transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.     

GENETIC MAPPING OF THE LAMINARIA JAPONICA (LAMINARALES,PHAEOPHYTA) USING AMPLIFIED FRAGMENT LENGTH POLYMORPHISM MARKERS1

To establish a molecular‐marker‐assisted system of breeding and genetic study for Laminaria japonica Aresch., amplified fragment length polymorphism (AFLP) was used to construct a genetic linkage map of L. japonica featuring 230 progeny of F₂ cross population. Eighteen primer combinations produced 370 polymorphic loci and 215 polymorphic loci segregated in a 3:1 Mendelian segregation ratio (P ≤ 0.05). Of the 215 segregated loci, 142 were ordered into 27 linkage groups. The length of the linkage groups ranged from 6.7 to 90.3 centimorgans (cM) with an average length of 49.6 cM, and the total length was 1,085.8 cM, which covered 68.4% of the estimated 1,586.9 cM genome. The number of mapped markers on each linkage group ranged from 2 to 12, averaging 5.3 markers per group. The average density of the markers was 1 per 9.4 cM. Based on the marker density and the resolution of the map, the constructed linkage map can satisfy the need for quantitative trait locus (QTL) location and molecular‐marker‐assisted breeding for Laminaria.     

Identification of quantitative trait loci and candidate genes for an anxiolytic‐like response to ethanol in BXD recombinant inbred strains

Genetic differences in acute behavioral responses to ethanol contribute to the susceptibility to alcohol use disorder and the reduction of anxiety is a commonly reported motive underlying ethanol consumption among alcoholics. Therefore, we studied the genetic variance in anxiolytic‐like responses to ethanol across the BXD recombinant inbred (RI) mouse panel using the light–dark transition model of anxiety. Strain‐mean genetic mapping and a mixed‐model quantitative trait loci (QTL) analysis replicated several previously published QTL for locomotor activity and identified several novel anxiety‐related loci. Significant loci included a chromosome 11 saline anxiety‐like QTL (Salanq1) and a chromosome 12 locus (Etanq1) influencing the anxiolytic‐like response to ethanol. Etanq1 was successfully validated by studies with BXD advanced intercross strains and fine‐mapped to a region comprising less than 3.5 Mb. Through integration of genome‐wide mRNA expression profiles of the mesocorticolimbic reward circuit (prefrontal cortex, nucleus accumbens and ventral midbrain) across the BXD RI panel, we identified high priority candidate genes within Etanq1, the strongest of which was Ninein (Nin), a Gsk3β‐interacting protein that is highly expressed in the brain.     

The identification and mapping of candidate genes and QTL involved in the fatty acid desaturation pathway in <Emphasis Type="Italic">Brassica napus</Emphasis>

We constructed a linkage map for the population QDH, which was derived from a cross between an oilseed rape cultivar and a resynthesised Brassica napus. The linkage map included ten markers linked to loci orthologous to those encoding fatty acid biosynthesis genes in Arabidopsis thaliana. The QDH population contains a high level of allelic variation, particularly in the C genome. We conducted quantitative trait locus (QTL) analyses, using field data obtained over 3 years, for the fatty acid composition of seed oil. The population segregates for the two major loci controlling erucic acid content, on linkage groups A8 and C3, which quantitatively affect the content of other fatty acids and is a problem generally encountered when crossing “wild” germplasm with cultivated “double low” oilseed rape cultivars. We assessed three methods for QTL analysis, interval mapping, multiple QTL mapping and single marker regression analysis of the subset of lines with low erucic acid. We found the third of these methods to be most appropriate for our main purpose, which was the study of the genetic control of the desaturation of 18-carbon fatty acids. This method enabled us to decouple the effect of the segregation of the erucic acid-controlling loci and identify 34 QTL for fatty acid content of seed oil, 14 in the A genome and 20 in the C genome. The QTL indicate the presence of 13 loci with novel alleles inherited from the progenitors of the resynthesised B. napus that might be useful for modulating the content or extent of desaturation of polyunsaturated fatty acids, only one of which coincides with the anticipated position of a candidate gene, an orthologue of FAD2.     

THE DISSOLUTION OF A COMPLEX GENETIC POLYMORPHISM: THE EVOLUTION OF SELF-FERTILIZATION IN TRISTYLOUS EICHHORNIA PANICULATA (PONTEDERIACEAE)

Eichhornia paniculata (Pontederiaceae) displays a wide range of outcrossing levels as a result of the dissolution of the tristylous genetic polymorphism and the evolution of semihomostyly. Population surveys, comparison of fitness components of the style morphs, and computer simulations were used to investigate the breakdown of tristyly and the selective mechanisms responsible for the evolution of self-fertilization. Of 110 populations surveyed in northeast Brazil and Jamaica, 53% were trimorphic, 25% were dimorphic, and 22% were monomorphic for style morph. The short (S) morph was underrepresented in trimorphic populations and absent from nontrimorphic populations. The mid (M) morph predominated in dimorphic populations and was the only morph in monomorphic populations. Stamen modifications promoting selfing, associated with semihomostyle evolution, were largely confined to the M morph. They were rare in trimorphic populations, common in dimorphic populations, and often fixed in monomorphic populations. Stochastic simulations and comparisons of fruit set in natural populations indicate that founder events, population bottlenecks, and lowered fertility of the S morph due to an absence of long-tongued pollinators can each account for loss of the S morph from trimorphic populations. A reduced level of disassortative mating can accentuate the rate at which the S morph is lost by both random and deterministic processes. Nontrimorphic populations occur at the geographical margins of the region surveyed and tend to be smaller and less dense than trimorphic populations. These observations and the higher fruit set of the M morph relative to the L morph in dimorphic populations suggest that reproductive assurance, favoring selfing variants of the M morph under conditions of low pollinator service, has been of primary importance in the origin of most monomorphic populations. Where pollinator service is reliable, however, automatic selection of selfing genes, aided by mating asymmetries between the morphs, can cause the M morph to spread to fixation in dimorphic populations.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.