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

Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep

The types and morphology of sheep horns have been extensively researched, yet the genetic foundation underlying the emergence of diverse horn characteristics during the breeding of polled Tibetan sheep has remained elusive. Genome-wide association analysis (GWAS) was performed on 103 subtypes (normal large horn, scurs, and polled) differentiated from G2 (offspring (G2) of parent (G1) of polled) of the polled core herd. Six single nucleotide polymorphisms (SNPs) located on chromosome 10 of the relaxin family peptide receptor 2 (RXFP2) gene exhibited positive correlations with horn length, horn base circumference, and horn base interval. Furthermore, in genotyping 382 G2 individuals, significant variations were observed for each specific horn type. Three additional mutations were identified near the target SNP upstream of the amplification product. Finally, the RXFP2-specific haplotype associated with the horned trait effectively maintained horn length, horn base circumference, and horn base interval in Tibetan sheep, as confirmed by population validation of nine loci in a sample size of 1125 individuals. The present study offers novel insights into the genetic differentiation of the horned type during improvement breeding and evolution, thereby establishing a robust theoretical foundation for polled Tibetan sheep breeding and providing valuable guidance for practical production.

Construction of a high-density genetic map and QTL localization of body weight and wool production related traits in Alpine Merino sheep based on WGR

BACKGROUND

The Alpine Merino is a new breed of fine-wool sheep adapted to the cold and arid climate of the plateau in the world. It has been popularized in Northwest China due to its superior adaptability as well as excellent production performance. Those traits related to body weight, wool yield, and wool fiber characteristics, which are economically essential traits in Alpine Merino sheep, are controlled by QTL (Quantitative Trait Loci). Therefore, the identification of QTL and genetic markers for these key economic traits is a critical step in establishing a MAS (Marker-Assisted Selection) breeding program.

RESULTS

In this study, we constructed the high-density genetic linkage map of Alpine Merino sheep by sequencing 110 F1 generation individuals using WGR (Whole Genome Resequencing) technology. 14,942 SNPs (Single Nucleotide Polymorphism) were identified and genotyped. The map spanned 2,697.86 cM, with an average genetic marker interval of 1.44 cM. A total of 1,871 high-quality SNP markers were distributed across 27 linkage groups, with an average of 69 markers per LG (Linkage Group). Among them, the smallest genetic distance is 19.62 cM for LG2, while the largest is 237.19 cM for LG19. The average genetic distance between markers in LGs ranged from 0.24 cM (LG2) to 3.57 cM (LG17). The marker density in the LGs ranged from LG14 (39 markers) to LG1 (150 markers).

CONCLUSIONS

The first genetic map of Alpine Merino sheep we constructed included 14,942 SNPs, while 46 QTLs associated with body weight, wool yield and wool fiber traits were identified, laying the foundation for genetic studies and molecular marker-assisted breeding. Notably, there were QTL intervals for overlapping traits on LG4 and LG8, providing potential opportunities for multi-trait co-breeding and further theoretical support for selection and breeding of ultra-fine and meaty Alpine Merino sheep.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

Three Copies of zbed1 Specific in Chromosome W Are Essential for Female-Biased Sexual Size Dimorphism in Cynoglossus semilaevis

The sex chromosome, especially specific in one sex, generally determines sexual size dimorphism (SSD), a phenomenon with dimorphic sexual difference in the body size. For Cynoglossus semilaevis, a flatfish in China, although the importance of chromosome W and its specific gene zbed1 in female-biased SSD have been suggested, its family members and regulation information are still unknown. At present, three zbed1 copies gene were identified on chromosome W, with no gametologs. Phylogenetic analysis for the ZBED family revealed an existence of ZBED9 in the fish. Nine members were uncovered from C. semilaevis, clustering into three kinds, ZBED1, ZBED4 and ZBEDX, which is less than the eleven kinds of ZBED members in mammals. The predominant expression of zbed1 in the female brain and pituitary tissues was further verified by qPCR. Transcription factor c/ebpα could significantly enhance the transcriptional activity of zbed1 promoter, which is opposite to its effect on the male determinant factor-dmrt1. When zbed1 was interfered with, piwil1, esr2 and wnt7b were up-regulated, while cell-cycle-related genes, including cdk4 and ccng1, were down-regulated. Thus, zbed1 is involved in cell proliferation by regulating esr2, piwil1, cell cycle and the Wnt pathway. Further research on their interactions would be helpful to understand fish SSD.

Identification of Critical Genes for Ovine Horn Development Based on Transcriptome during the Embryonic Period

Horns, also known as headgear, are a unique structure of ruminants. As ruminants are globally distributed, the study of horn formation is critical not only for increasing our understanding of natural and sexual selection but also for the breeding of polled sheep breeds to facilitate modern sheep farming. Despite this, a significant number of the underlying genetic pathways in sheep horn remain unclear. In this study, to clarify the gene expression profile of horn buds and investigate the key genes in horn bud formation, RNA-sequencing (RNA-seq) technology was utilized to investigate differential gene expression in the horn buds and adjacent forehead skin of Altay sheep fetuses. There were only 68 differentially expressed genes (DEGs) identified, consisting of 58 up-regulated genes and 10 down-regulated genes. RXFP2 was differentially up-regulated in the horn buds and had the highest significance (p-value = 7.42 × 10^-14). In addition, 32 DEGs were horn-related genes identified in previous studies, such as RXFP2, FOXL2, SFRP4, SFRP2, KRT1, KRT10, WNT7B, and WNT3. Further, Gene Ontology (GO) analysis showed that the DEGs were mainly enriched with regard to growth, development, and cell differentiation. Pathway analysis revealed that the Wnt signaling pathway may be responsible for horn development. Further, through combining the protein-protein interaction networks of the DEGs, it was found that the top five hub genes, namely, ACAN, SFRP2, SFRP4, WNT3, and WNT7B, were also associated with horn development. Our results suggest that only a few key genes, including RXFP2, are involved in bud formation. This study not only validates the expression of candidate genes identified at the transcriptome level in previous studies but also provides new possible marker genes for horn development, which may promote our understanding of the genetic mechanisms of horn formation.

A polygenic basis for birth weight in a wild population of red deer (Cervus elaphus)

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait.

Ewe are what ewe wear: bigger horns, better ewes and the potential consequence of trophy hunting on female fitness in bighorn sheep

In polygynous species, secondary sexual traits such as weapons or elaborate ornaments have evolved through intrasexual competition for mates. In some species, these traits are present in both sexes but are underdeveloped in the sex facing lower intrasexual competition for mates. It is often assumed that these underdeveloped sexually selected traits are a vestige of strong sexual selection on the other sex. Here, we challenge this assumption and investigate whether the expression of secondary sexual traits is associated with fitness in female bighorn sheep. Analyses of 45 years of data revealed that female horn length at 2 years, while accounting for mass and environmental variables, is associated with younger age at primiparity, younger age of first offspring weaned, greater reproductive lifespan and higher lifetime reproductive success. There was no association between horn length and fecundity. These findings highlight a potential conservation issue. In this population, trophy hunting selects against males with fast-growing horns. Intersexual genetic correlations imply that intense selective hunting of large-horned males before they can reproduce can decrease female horn size. Therefore, intense trophy hunting of males based on horn size could reduce female reproductive performance through the associations identified here, and ultimately reduce population growth and viability.

Genomic association with pathogen carriage in bighorn sheep (Ovis canadensis)

Genetic composition can influence host susceptibility to, and transmission of, pathogens, with potential population-level consequences. In bighorn sheep (Ovis canadensis), pneumonia epidemics caused by Mycoplasma ovipneumoniae have been associated with severe population declines and limited recovery across North America. Adult survivors either clear the infection or act as carriers that continually shed M. ovipneumoniae and expose their susceptible offspring, resulting in high rates of lamb mortality for years following the outbreak event. Here, we investigated the influence of genomic composition on persistent carriage of M. ovipneumoniae in a well-studied bighorn sheep herd in the Wallowa Mountains of Oregon, USA. Using 10,605 SNPs generated using RADseq technology for 25 female bighorn sheep, we assessed genomic diversity metrics and employed family-based genome-wide association methodologies to understand variant association and genetic architecture underlying chronic carriage. We observed no differences among genome-wide diversity metrics (heterozygosity and allelic richness) between groups. However, we identified two variant loci of interest and seven associated candidate genes, which may influence carriage status. Further, we found that the SNP panel explained ~55% of the phenotypic variance (SNP-based heritability) for M. ovipneumoniae carriage, though there was considerable uncertainty in these estimates. While small sample sizes limit conclusions drawn here, our study represents one of the first to assess the genomic factors influencing chronic carriage of a pathogen in a wild population and lays a foundation for understanding genomic influence on pathogen persistence in bighorn sheep and other wildlife populations. Future research should incorporate additional individuals as well as distinct herds to further explore the genomic basis of chronic carriage.

Opinion: Is gene mapping in wild populations useful for understanding and predicting adaptation to global change?

Changing environmental conditions will inevitably alter selection pressures. Over the long term, populations have to adapt to these altered conditions by evolutionary change to avoid extinction. Quantifying the 'evolutionary potential' of populations to predict whether they will be able to adapt fast enough to forecasted changes is crucial to fully assess the threat for biodiversity posed by climate change. Technological advances in sequencing and high-throughput genotyping have now made genomic studies possible in a wide range of species. Such studies, in theory, allow an unprecedented understanding of the genomics of ecologically relevant traits and thereby a detailed assessment of the population's evolutionary potential. Aimed at a wider audience than only evolutionary geneticists, this paper gives an overview of how gene-mapping studies have contributed to our understanding and prediction of evolutionary adaptations to climate change, identifies potential reasons why their contribution to understanding adaptation to climate change may remain limited, and highlights approaches to study and predict climate change adaptation that may be more promising, at least in the medium term.

Heritability of Horn Size in Thinhorn Sheep

Understanding the genetic basis of fitness-related trait variation has long been of great interest to evolutionary biologists. Secondary sexual characteristics, such as horns in bovids, are particularly intriguing since they can be potentially affected by both natural and sexual selection. Until recently, however, the study of fitness-related quantitative trait variation in wild species has been hampered by a lack of genomic resources, pedigree, and/or phenotype data. Recent innovations in genomic technologies have enabled wildlife researchers to perform marker-based relatedness estimation and acquire adequate loci density, enabling both the “top-down” approach of quantitative genetics and the “bottom-up” approach of association studies to describe the genetic basis of fitness-related traits. Here we combine a cross species application of the OvineHD BeadChip and horn measurements (horn length, base circumference, and volume) from harvested thinhorn sheep to examine the heritability and to perform a genome-wide single-nucleotide polymorphism association study of horn size in the species. Thinhorn sheep are mountain ungulates that reside in the mountainous regions of northwestern North America. Thinhorn sheep males grow massive horns that determine the social rank and mating success. We found horn length, base circumference, and volume to be moderately heritable and two loci to be suggestively associated with horn length.

Genomic prediction of the polled and horned phenotypes in Merino sheep

BACKGROUND

In horned sheep breeds, breeding for polledness has been of interest for decades. The objective of this study was to improve prediction of the horned and polled phenotypes using horn scores classified as polled, scurs, knobs or horns. Derived phenotypes polled/non-polled (P/NP) and horned/non-horned (H/NH) were used to test four different strategies for prediction in 4001 purebred Merino sheep. These strategies include the use of single 'single nucleotide polymorphism' (SNP) genotypes, multiple-SNP haplotypes, genome-wide and chromosome-wide genomic best linear unbiased prediction and information from imputed sequence variants from the region including the RXFP2 gene. Low-density genotypes of these animals were imputed to the Illumina Ovine high-density (600k) chip and the 1.78-kb insertion polymorphism in RXFP2 was included in the imputation process to whole-genome sequence. We evaluated the mode of inheritance and validated models by a fivefold cross-validation and across- and between-family prediction.

RESULTS

The most significant SNPs for prediction of P/NP and H/NH were OAR10_29546872.1 and OAR10_29458450, respectively, located on chromosome 10 close to the 1.78-kb insertion at 29.5 Mb. The mode of inheritance included an additive effect and a sex-dependent effect for dominance for P/NP and a sex-dependent additive and dominance effect for H/NH. Models with the highest prediction accuracies for H/NH used either single SNPs or 3-SNP haplotypes and included a polygenic effect estimated based on traditional pedigree relationships. Prediction accuracies for H/NH were 0.323 for females and 0.725 for males. For predicting P/NP, the best models were the same as for H/NH but included a genomic relationship matrix with accuracies of 0.713 for females and 0.620 for males.

CONCLUSIONS

Our results show that prediction accuracy is high using a single SNP, but does not reach 1 since the causative mutation is not genotyped. Incomplete penetrance or allelic heterogeneity, which can influence expression of the phenotype, may explain why prediction accuracy did not approach 1 with any of the genetic models tested here. Nevertheless, a breeding program to eradicate horns from Merino sheep can be effective by selecting genotypes GG of SNP OAR10_29458450 or TT of SNP OAR10_29546872.1 since all sheep with these genotypes will be non-horned.

Genomic analysis of morphometric traits in bighorn sheep using the Ovine Infinium® HD SNP BeadChip

Elucidating the genetic basis of fitness-related traits is a major goal of molecular ecology. Traits subject to sexual selection are particularly interesting, as non-random mate choice should deplete genetic variation and thereby their evolutionary benefits. We examined the genetic basis of three sexually selected morphometric traits in bighorn sheep (Ovis canadensis): horn length, horn base circumference, and body mass. These traits are of specific concern in bighorn sheep as artificial selection through trophy hunting opposes sexual selection. Specifically, horn size determines trophy status and, in most North American jurisdictions, if an individual can be legally harvested. Using between 7,994–9,552 phenotypic measures from the long-term individual-based study at Ram Mountain (Alberta, Canada), we first showed that all three traits are heritable (h² = 0.15–0.23). We then conducted a genome-wide association study (GWAS) utilizing a set of 3,777 SNPs typed in 76 individuals using the Ovine Infinium^® HD SNP BeadChip. We found suggestive association for body mass at a single locus (OAR9_91647990). The absence of strong associations with SNPs suggests that the traits are likely polygenic. These results represent a step forward for characterizing the genetic architecture of fitness related traits in sexually dimorphic ungulates.

SNP discovery in candidate adaptive genes using exon capture in a free-ranging alpine ungulate

Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identify SNP markers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein-coding and nearby 5' and 3' untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis) exon capture data and directly from the domestic sheep genome (Ovis aries v. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR-based SNP chip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositan and bayescan), we detected 28 SNP loci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease-regulating functions (e.g. Ovar-DRA, APC, BATF2, MAGEB18), cell regulation signalling pathways (e.g. KRIT1, PI3K, ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene-targeted SNP discovery and subsequent SNP chip genotyping using low-quality samples in a nonmodel species.

Intense selective hunting leads to artificial evolution in horn size

The potential for selective harvests to induce rapid evolutionary change is an important question for conservation and evolutionary biology, with numerous biological, social and economic implications. We analyze 39 years of phenotypic data on horn size in bighorn sheep (Ovis canadensis) subject to intense trophy hunting for 23 years, after which harvests nearly ceased. Our analyses revealed a significant decline in genetic value for horn length of rams, consistent with an evolutionary response to artificial selection on this trait. The probability that the observed change in male horn length was due solely to drift is 9.9%. Female horn length and male horn base, traits genetically correlated to the trait under selection, showed weak declining trends. There was no temporal trend in genetic value for female horn base circumference, a trait not directly targeted by selective hunting and not genetically correlated with male horn length. The decline in genetic value for male horn length stopped, but was not reversed, when hunting pressure was drastically reduced. Our analysis provides support for the contention that selective hunting led to a reduction in horn length through evolutionary change. It also confirms that after artificial selection stops, recovery through natural selection is slow.

Whole-genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations

Dissecting the genetic basis of phenotypic variation in natural populations is a long-standing goal in evolutionary biology. One open question is whether quantitative traits are determined only by large numbers of genes with small effects, or whether variation also exists in large-effect loci. We conducted genomewide association analyses of forehead patch size (a sexually selected trait) on 81 whole-genome-resequenced male collared flycatchers with extreme phenotypes, and on 415 males sampled independent of patch size and genotyped with a 50K SNP chip. No SNPs were genomewide statistically significantly associated with patch size. Simulation-based power analyses suggest that the power to detect large-effect loci responsible for 10% of phenotypic variance was <0.5 in the genome resequencing analysis, and <0.1 in the SNP chip analysis. Reducing the recombination by two-thirds relative to collared flycatchers modestly increased power. Tripling sample size increased power to >0.8 for resequencing of extreme phenotypes (N = 243), but power remained <0.2 for the 50K SNP chip analysis (N = 1245). At least 1 million SNPs were necessary to achieve power >0.8 when analysing 415 randomly sampled phenotypes. However, power of the 50K SNP chip to detect large-effect loci was nearly 0.8 in simulations with a small effective population size of 1500. These results suggest that reliably detecting large-effect trait loci in large natural populations will often require thousands of individuals and near complete sampling of the genome. Encouragingly, far fewer individuals and loci will often be sufficient to reliably detect large-effect loci in small populations with widespread strong linkage disequilibrium.

Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long‐term studies across Europe and consequently are considered an ecological ‘model organism’. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long‐term study populations of great tits—one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population‐specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population‐specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free‐living populations.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis)

Background

Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries).

Results

Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition.

Conclusion

Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1618-x) contains supplementary material, which is available to authorized users.

A microsatellite-based linkage map for song sparrows (Melospiza melodia)

Although linkage maps are important tools in evolutionary biology, their availability for wild populations is limited. The population of song sparrows (Melospiza melodia) on Mandarte Island, Canada, is among the more intensively studied wild animal populations. Its long-term pedigree data, together with extensive genetic sampling, have allowed the study of a range of questions in evolutionary biology and ecology. However, the availability of genetic markers has been limited. We here describe 191 new microsatellite loci, including 160 high-quality polymorphic autosomal, 7 Z-linked and 1 W-linked markers. We used these markers to construct a linkage map for song sparrows with a total sex-averaged map length of 1731 cM and covering 35 linkage groups, and hence, these markers cover most of the 38-40 chromosomes. Female and male map lengths did not differ significantly. We then bioinformatically mapped these loci to the zebra finch (Taeniopygia guttata) genome and found that linkage groups were conserved between song sparrows and zebra finches. Compared to the zebra finch, marker order within small linkage groups was well conserved, whereas the larger linkage groups showed some intrachromosomal rearrangements. Finally, we show that as expected, recombination frequency between linked loci explained the majority of variation in gametic phase disequilibrium. Yet, there was substantial overlap in gametic phase disequilibrium between pairs of linked and unlinked loci. Given that the microsatellites described here lie on 35 of the 38-40 chromosomes, these markers will be useful for studies in this species, as well as for comparative genomics studies with other species.

The role of sex chromosomes in sexual dimorphism: discordance between molecular and phenotypic data

In addition to initial sex determination, genes on the sex chromosomes are theorized to play a particularly important role in phenotypic differences between males and females. Sex chromosomes in many species display molecular signatures consistent with these theoretical predictions, particularly through sex-specific gene expression. However, the phenotypic implications of this molecular signature are unresolved, and the role of the sex chromosomes in quantitative genetic studies of phenotypic sex differences is largely equivocal. In this article, we examine molecular and phenotypic data in the light of theoretical predictions about masculinization and feminization of the sex chromosomes. Additionally, we discuss the role of genetic and regulatory complexities in the genome–phenotype relationship, and ultimately how these affect the overall role of the sex chromosomes in sex differences.

Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter-subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.

Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations

Sexual dimorphism can evolve when males and females differ in phenotypic optima. Genetic constraints can, however, limit the evolution of sexual dimorphism. One possible constraint is derived from alleles expressed in both sexes. Because males and females share most of their genome, shared alleles with different fitness effects between sexes are faced with intralocus sexual conflict. Another potential constraint is derived from genetic correlations between developmental stages. Sexually dimorphic traits are often favoured at adult stages, but selected against as juvenile, so developmental decoupling of traits between ontogenetic stages may be necessary for the evolution of sexual dimorphism in adults. Resolving intralocus conflicts between sexes and ages is therefore a key to the evolution of age-specific expression of sexual dimorphism. We investigated the genetic architecture of divergence in the ontogeny of sexual dimorphism between two populations of the Japanese medaka (Oryzias latipes) that differ in the magnitude of dimorphism in anal and dorsal fin length. Quantitative trait loci (QTL) mapping revealed that few QTL had consistent effects throughout ontogenetic stages and the majority of QTL change the sizes and directions of effects on fin growth rates during ontogeny. We also found that most QTL were sex-specific, suggesting that intralocus sexual conflict is almost resolved. Our results indicate that sex- and age-specific QTL enable the populations to achieve optimal developmental trajectories of sexually dimorphic traits in response to complex natural and sexual selection.

Challenges and prospects in genome-wide quantitative trait loci mapping of standing genetic variation in natural populations

A considerable challenge in evolutionary genetics is to understand the genetic mechanisms that facilitate or impede evolutionary adaptation in natural populations. For this, we must understand the genetic loci contributing to trait variation and the selective forces acting on them. The decreased costs and increased feasibility of obtaining genotypic data on a large number of individuals have greatly facilitated gene mapping in natural populations, particularly because organisms whose genetics have been historically difficult to study are now within reach. Here we review the methods available to evolutionary ecologists interested in dissecting the genetic basis of traits in natural populations. Our focus lies on standing genetic variation in outbred populations. We present an overview of the current state of research in the field, covering studies on both plants and animals. We also draw attention to particular challenges associated with the discovery of quantitative trait loci and discuss parallels to studies on crops, livestock, and humans. Finally, we point to some likely future developments in genetic mapping studies.

A quantitative trait locus analysis of personality in wild bighorn sheep

Personality, the presence of persistent behav105ioral differences among individuals over time or contexts, potentially has important ecological and evolutionary consequences. However, a lack of knowledge about its genetic architecture limits our ability to understand its origin, evolution, and maintenance. Here, we report on a genome-wide quantitative trait locus (QTL) analysis for two personality traits, docility and boldness, in free-living female bighorn sheep from Ram Mountain, Alberta, Canada. Our variance component linkage analysis based on 238 microsatellite loci genotyped in 310 pedigreed individuals identified suggestive docility and boldness QTL on sheep chromosome 2 and 6, respectively. A lack of QTL overlap indicated that genetic covariance between traits was not modulated by pleiotropic effects at a major locus and may instead result from linkage disequilibrium or pleiotropic effects at QTL of small effects. To our knowledge, this study represents the first attempt to dissect the genetic architecture of personality in a free-living wildlife population, an important step toward understanding the link between molecular genetic variation in personality and fitness and the evolutionary processes maintaining this variation.