GENETICS OF INCIPIENT SPECIATION IN DROSOPHILA MOJAVENSIS. III. LIFE‐HISTORY DIVERGENCE IN ALLOPATRY AND REPRODUCTIVE ISOLATION

We carried out a three‐tiered genetic analysis of egg‐to‐adult development time and viability in ancestral and derived populations of cactophilic Drosophila mojavensis to test the hypothesis that evolution of these life‐history characters has shaped premating reproductive isolation in this species. First, a common garden experiment with 11 populations from Baja California and mainland Mexico and Arizona reared on two host species revealed significant host plant X region and population interactions for viability and development time, evidence for host plant adaptation. Second, replicated line crosses with flies reared on both hosts revealed autosomal, X chromosome, cytoplasmic, and autosome X cactus influences on development time. Viability differences were influenced by host plants, autosomal dominance, and X chromosomal effects. Many of the F₁, F₂, and backcross generations showed evidence of heterosis for viability. Third, a QTL analysis of male courtship song and epicuticular hydrocarbon variation based on 1688 Baja × mainland F₂ males also revealed eight QTL influencing development time differences. Mainland alleles at six of these loci were associated with longer development times, consistent with population‐level differences. Eight G × E interactions were also detected caused by longer development times of mainland alleles expressed on a mainland host with smaller differences among Baja genotypes reared on the Baja host plant. Four QTL influenced both development time and epicuticular hydrocarbon differences associated with courtship success, and there was a significant QTL‐based correlation between development time and cuticular hydrocarbon variation. Thus, the regional shifts in life histories that evolved once D. mojavensis invaded mainland Mexico from Baja California by shifting host plants were genetically correlated with variation in cuticular hydrocarbon‐based mate preferences.     

Genetic mapping of biomass yield in three interconnected Miscanthus populations

Improving biomass yield is a major goal of Miscanthus breeding. We conducted a study on one interspecific Miscanthus sinensis × Miscanthus sacchariflorus F₁ population and two intraspecific M. sinensis F₁ populations, each of which shared a common parent. A field trial was established at Urbana, IL during spring 2011, and phenotypic data were collected in 2012 and 2013 for fourteen yield traits. Six high‐density parental genetic maps, as well as a consensus genetic map integrating M. sinensis and M. sacchariflorus, were developed via the pseudotestcross strategy for noninbred parents with ≥1214 single‐nucleotide polymorphism markers generated from restriction site‐associated DNA sequencing. We confirmed for the first time a whole‐genome duplication in M. sacchariflorus relative to Sorghum bicolor, similar to that observed previously for M. sinensis. Four quantitative trait locus (QTL) analysis methods for detecting marker‐trait associations were compared: (1) individual parental map composite interval mapping analysis, (2) individual parental map stepwise analysis, (3) consensus map single‐population stepwise analysis and (4) consensus map joint‐population stepwise analysis. These four methods detected 288, 264, 133 and 109 total QTLs, which resolved into 157, 136, 106 and 86 meta‐QTLs based on QTL congruency, respectively, including a set of 59 meta‐QTLs common to all four analysis methods. Composite interval mapping and stepwise analysis co‐identified 118 meta‐QTLs across six parental maps, suggesting high reliability of stepwise regression in QTL detection. Joint‐population stepwise analysis yielded the highest resolution of QTLs compared to the other three methods across all meta‐QTLs. Strong, frequently advantageous transgressive segregation in the three populations indicated a promising future for breeding new higher‐yielding cultivars of Miscanthus.     

QTL for body weight,morphometric traits and stress response in European sea bass Dicentrarchus labrax

Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two‐step experimental design for detecting QTL in mass‐spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full‐sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half‐sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.     

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.     

Quantitative trait loci responsible for variation in sexually dimorphic traits in Drosophila melanogaster

To understand the mechanisms of morphological evolution and species divergence, it is essential to elucidate the genetic basis of variation in natural populations. Sexually dimorphic characters, which evolve rapidly both within and among species, present attractive models for addressing these questions. In this report, we map quantitative trait loci (QTL) responsible for variation in sexually dimorphic traits (abdominal pigmentation and the number of ventral abdominal bristles and sex comb teeth) in a natural population of Drosophila melanogaster. To capture the pattern of genetic variation present in the wild, a panel of recombinant inbred lines was created from two heterozygous flies taken directly from nature. High-resolution mapping was made possible by cytological markers at the average density of one per 2 cM. We have used a new Bayesian algorithm that allows QTL mapping based on all markers simultaneously. With this approach, we were able to detect small-effect QTL that were not evident in single-marker analyses. Our results show that at least for some sexually dimorphic traits, a small number of QTL account for the majority of genetic variation. The three strongest QTL account for >60% of variation in the number of ventral abdominal bristles. Strikingly, a single QTL accounts for almost 60% of variation in female abdominal pigmentation. This QTL maps to the chromosomal region that Robertson et al. have found to affect female abdominal pigmentation in other populations of D. melanogaster. Using quantitative complementation tests, we demonstrate that this QTL is allelic to the bric a brac gene, whose expression has previously been shown to correlate with interspecific differences in pigmentation. Multiple bab alleles that confer distinct phenotypes appear to segregate in natural populations at appreciable frequencies, suggesting that intraspecific and interspecific variation in abdominal pigmentation may share a similar genetic basis.     

Congruence in QTL for adventitious rooting in Pinus elliottii × Pinus caribaea hybrids resolves between and within-species effects

Targeting between-species effects for improvement in synthetic hybrid populations derived from outcrossing parental tree species may be one way to increase the efficacy and predictability of hybrid breeding. We present a comparative analysis of the quantitative trait loci (QTL) which resolved between from within-species effects for adventitious rooting in two populations of hybrids between Pinus elliottii and P. caribaea, an outbred F₁ (n=287) and an inbred-like F₂ family (n=357). Most small to moderate effect QTL (each explaining 2–5% of phenotypic variation, PV) were congruent (3 out of 4 QTL in each family) and therefore considered within-species effects as they segregated in both families. A single large effect QTL (40% PV) was detected uniquely in the F₂ family and assumed to be due to a between-species effect, resulting from a genetic locus with contrasting alleles in each parental species. Oligogenic as opposed to polygenic architecture was supported in both families (60% and 20% PV explained by 4 QTL in the F₂ and F₁ respectively). The importance of adventitious rooting for adaptation to survive water-logged environments was thought in part to explain oligogenic architecture of what is believed to be a complex trait controlled by many hundreds of genes.     

Mapping 49 quantitative trait loci at high resolution through sequencing-based genotyping of rice recombinant inbred lines

Mapping chromosome regions responsible for quantitative phenotypic variation in recombinant populations provides an effective means to characterize the genetic basis of complex traits. We conducted a quantitative trait loci (QTL) analysis of 150 rice recombinant inbred lines (RILs) derived from a cross between two cultivars, Oryza sativa ssp. indica cv. 93-11 and Oryza sativa ssp. japonica cv. Nipponbare. The RILs were genotyped through next-generation sequencing, which accurately determined the recombination breakpoints and provided a new type of genetic markers, recombination bins, for QTL analysis. We detected 49 QTL with phenotypic effect ranging from 3.2 to 46.0% for 14 agronomics traits. Five QTL of relatively large effect (14.6–46.0%) were located on small genomic regions, where strong candidate genes were found. The analysis using sequencing-based genotyping thus offers a powerful solution to map QTL with high resolution. Moreover, the RILs developed in this study serve as an excellent system for mapping and studying genetic basis of agricultural and biological traits of rice.     

Fine-scale genetic analysis of species-specific female preference in Drosophila simulans

Behavioural differences are thought to be the first components to contribute to species isolation, yet the precise genetic basis of behavioural isolation remains poorly understood. Here, we used a combination of behaviour assays and genetic mapping to provide the first refined map locating candidate genes for interspecific female preference isolating Drosophila simulans from D. melanogaster. First, we tested whether two genes identified as affecting D. melanogaster female intraspecific mate choice also affect interspecific mate choice; neither of these genes was found to contribute to species‐specific female preference. Next, we used deficiency mapping to locate genes on the right arm of the third chromosome for species‐specific female preference and identified five small significant regions that contain candidate genes contributing to behavioural isolation. All five regions were located in areas that would have low interspecific recombination, which mirrors the results of other behavioural isolation studies that used quantitative trait locus (QTL) mapping, but without the potential concern of bias towards regions of low recombination that QTL mapping may have. As this model system may be refined to the individual gene level using the same methodology, this initial map we provide may potentially serve as a ready template for the identification and characterization of the first behavioural isolation genes.     

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.     

THREE SELECTIONS ARE BETTER THAN ONE: CLINAL VARIATION OF THERMAL QTL FROM INDEPENDENT SELECTION EXPERIMENTS IN DROSOPHILA

We report the results of two independent selection experiments that have exposed distinct populations of Drosophila melanogaster to different forms of thermal selection. A recombinant population derived from Arvin California and Zimbabwe isofemale lines was exposed to laboratory natural selection at two temperatures (T_AZ: 18°C and 28°C). Microsatellite mapping identified quantitative trait loci (QTL) on the X‐chromosome between the replicate “Hot” and “Cold” populations. In a separate experiment, disruptive selection was imposed on an outbred California population for the “knockdown” temperature (T_KD) in a thermal column. Microsatellite mapping of the “High” and “Low” populations also uncovered primarily X‐linked QTL. Notably, a marker in the shaggy locus at band 3A was significantly differentiated in both experiments. Finer scale mapping of the 3A region has narrowed the QTL to the shaggy gene region, which contains several candidate genes that function in circadian rhythms. The same allele that was increased in frequency in the High T_KD populations is significantly clinal in North America and is more common at the warm end of the cline (Florida vs. Maine; however, the cline was not apparent in Australia). Together, these studies show that independent selection experiments can uncover the same target of selection and that evolution in the laboratory can recapitulate putatively adaptive clinal variation in nature.     

Mapping QTL for production traits in segregating Piétrain pig populations using genome‐wide association study results of F2 crosses

In this study, genome‐wide association study (GWAS) results of porcine F2 crosses were used to map QTL in outcross Piétrain populations. For this purpose, two F2 crosses (Piétrain × Meishan, n = 304; Piétrain × Wild Boar, n = 291) were genotyped with the PorcineSNP60v2 BeadChip and phenotyped for the dressing yield, carcass length, daily gain and drip loss traits. GWASs were conducted in the pooled F2 cross applying single marker mixed linear models. For the investigated traits, between two and five (in total 15) QTL core regions, spanning 250 segregating SNPs around a significant trait‐associated peak SNP, were identified. The SNPs within the QTL core regions were subsequently tested for trait association in two outcross Piétrain populations consisting of 771 progeny‐tested boars and 210 sows with their own performance records. In the sow (boar) dataset, five (eight) of the 15 mapped QTL were validated. Hence, many QTL mapped in the F2 crosses (with Piétrain as a common founder breed) are still segregating in the current Piétrain breed. This confirms the usefulness of existing F2 crosses for mapping QTL that are still segregating in the recent founder breed generation. The approach utilizes the high power of an F2 cross to map QTL in a breeding population for which it is not guaranteed that they would be found using a GWAS in this population.     

More than one way to evolve a weed: parallel evolution of US weedy rice through independent genetic mechanisms

Many different crop species were selected for a common suite of ‘domestication traits’, which facilitates their use for studies of parallel evolution. Within domesticated rice (Oryza sativa), there has also been independent evolution of weedy strains from different cultivated varieties. This makes it possible to examine the genetic basis of parallel weed evolution and the extent to which this process occurs through shared genetic mechanisms. We performed comparative QTL mapping of weediness traits using two recombinant inbred line populations derived from crosses between an indica crop variety and representatives of each of the two independently evolved weed strains found in US rice fields, strawhull (S) and blackhull awned (B). Genotyping‐by‐sequencing provided dense marker coverage for linkage map construction (average marker interval <0.25 cM), with 6016 and 13 730 SNPs mapped in F₅ lines of the S and B populations, respectively. For some weediness traits (awn length, hull pigmentation and pericarp pigmentation), QTL mapping and sequencing of underlying candidate genes confirmed that trait variation was largely attributable to individual loci. However, for more complex quantitative traits (including heading date, panicle length and seed shattering), we found multiple QTL, with little evidence of shared genetic bases between the S and B populations or across previous studies of weedy rice. Candidate gene sequencing revealed causal genetic bases for 8 of 27 total mapped QTL. Together these findings suggest that despite the genetic bottleneck that occurred during rice domestication, there is ample genetic variation in this crop to allow agricultural weed evolution through multiple genetic mechanisms.     

Insights into the genetic architecture of sexual dimorphism from an interspecific cross between two diverging Silene (Caryophyllaceae) species

The evolution of sexual dimorphism in species with separate sexes is influenced by the resolution of sexual conflicts creating sex differences through genetic linkage or sex‐biased expression. Plants with different degrees of sexual dimorphism are thus ideal to study the genetic basis of sexual dimorphism. In this study we explore the genetic architecture of sexual dimorphism between Silene latifolia and Silene dioica. These species have chromosomal sex determination and differ in the extent of sexual dimorphism. To test whether QTL for sexually dimorphic traits have accumulated on the sex chromosomes and to quantify their contribution to species differences, we create a linkage map and performed QTL analysis of life history, flower and vegetative traits using an unidirectional interspecific F2 hybrid cross. We found support for an accumulation of QTL on the sex chromosomes and that sex differences explained a large proportion of the variance between species, suggesting that both natural and sexual selection contributed to species divergence. Sexually dimorphic traits that also differed between species displayed transgressive segregation. We observed a reversal in sexual dimorphism in the F2 population, where males tended to be larger than females, indicating that sexual dimorphism is constrained within populations but not in recombinant hybrids. This study contributes to the understanding of the genetic basis of sexual dimorphism and its evolution in Silene.     

Occasional recombination of a selfish X‐chromosome may permit its persistence at high frequencies in the wild

The sex‐ratio X‐chromosome (SR) is a selfish chromosome that promotes its own transmission to the next generation by destroying Y‐bearing sperm in the testes of carrier males. In some natural populations of the fly Drosophila neotestacea, up to 30% of the X‐chromosomes are SR chromosomes. To investigate the molecular evolutionary history and consequences of SR, we sequenced SR and standard (ST) males at 11 X‐linked loci that span the ST X‐chromosome and at seven arbitrarily chosen autosomal loci from a sample of D. neotestacea males from throughout the species range. We found that the evolutionary relationship between ST and SR varies among individual markers, but genetic differentiation between SR and ST is chromosome‐wide and likely due to large chromosomal inversions that suppress recombination. However, SR does not consist of a single multilocus haplotype: we find evidence for gene flow between ST and SR at every locus assayed. Furthermore, we do not find long‐distance linkage disequilibrium within SR chromosomes, suggesting that recombination occurs in females homozygous for SR. Finally, polymorphism on SR is reduced compared to that on ST, and loci displaying signatures of selection on ST do not show similar patterns on SR. Thus, even if selection is less effective on SR, our results suggest that gene flow with ST and recombination between SR chromosomes may prevent the accumulation of deleterious mutations and allow its long‐term persistence at relatively high frequencies.     

An interspecific QTL study of Drosophila wing size and shape variation to investigate the genetic basis of morphological differences

The Drosophila wing has been used as a model in studies of morphogenesis and evolution; the use of such models can contribute to our understanding of mechanisms that promote morphological divergence among populations and species. We mapped quantitative trait loci (QTL) affecting wing size and shape traits using highly inbred introgression lines between D. simulans and D. sechellia, two sibling species of the melanogaster subgroup. Eighteen QTL peaks that are associated with 12 wing traits were identified, including two principal components. The wings of D. simulans and D. sechellia significantly diverged in size; two of the QTL peaks could account for part of this interspecific divergence. Both of these putative QTLs were mapped at the same cytological regions as other QTLs for intraspecific wing size variation identified in D. melanogaster studies. In these regions, one or more loci could account for intra- and interspecific variation in the size of Drosophila wings. Three other QTL peaks were related to a pattern of interspecific variation in wing size and shape traits that is summarized by one principal component. In addition, we observed that female wings are significantly larger and longer than male wings and the second, fourth and fifth longitudinal veins are closer together at the distal wing area. This pattern was summarized by another principal component, for which one QTL was mapped.     

Genetic basis for rapidly evolved tolerance in the wild: adaptation to toxic pollutants by an estuarine fish species

Atlantic killifish (Fundulus heteroclitus) residing in some urban and industrialized estuaries of the US eastern seaboard demonstrate recently evolved and extreme tolerance to toxic aryl hydrocarbon pollutants, characterized as dioxin‐like compounds (DLCs). Here, we provide an unusually comprehensive accounting (69%) through quantitative trait locus (QTL) analysis of the genetic basis for DLC tolerance in killifish inhabiting an urban estuary contaminated with PCB congeners, the most toxic of which are DLCs. Consistent with mechanistic knowledge of DLC toxicity in fish and other vertebrates, the aryl hydrocarbon receptor (ahr2) region accounts for 17% of trait variation; however, QTL on independent linkage groups and their interactions have even greater explanatory power (44%). QTL interpreted within the context of recently available Fundulus genomic resources and shared synteny among fish species suggest adaptation via interacting components of a complex stress response network. Some QTL were also enriched in other killifish populations characterized as DLC‐tolerant and residing in distant urban estuaries contaminated with unique mixtures of pollutants. Together, our results suggest that DLC tolerance in killifish represents an emerging example of parallel contemporary evolution that has been driven by intense human‐mediated selection on natural populations.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Joint analysis for heading date QTL in small interconnected barley populations

The purpose of the present work is to validate the effect of the main QTL determining heading date in a set of 281 doubled haploid lines of barley, derived from 17 small interconnected populations, whose parents are cultivars commonly used in the Spanish barley breeding program. We used 72 molecular markers distributed across the seven chromosomes, particularly in regions known to contain flowering time genes or QTL. A combined linkage map over the 17 populations was constructed. The lines were evaluated in four field trials: two autumn sowings and two winter sowings, and in two treatments at a greenhouse trial, under controlled conditions of photoperiod and temperature. We have found that it is possible to carry out QTL detection in a complex germplasm set, representative of the materials used in an active breeding programme. In most cases two alleles per QTL were detected, though polymorphism of flanking markers was notably higher. The results revealed that there is a set of QTL that accounts for an important percentage of the phenotypic variation, suitable for marker assisted selection. Also, the role of the regions carrying the photoperiod response genes Ppd-H1 and Ppd-H2, the vernalization response genes Vrn-H1 and Vrn-H2, and the earliness per se locus Eam6, of which allele-specific or closely linked markers were available, was confirmed. These results support the use of this kind of approach for the validation of QTL found in single cross population studies, or to survey allelic diversity in plant breeding sets of materials. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life‐history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half‐sib families to construct nine linkage groups (LG). In total, four growth‐related QTL (two chromosome‐wide and two genome‐wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome‐wide and one genome‐wide). The proportion of phenotypic variation explained by the body‐weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker‐assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.     

Increasing the accuracy of genomic prediction in pure‐bred Limousin beef cattle by including cross‐bred Limousin data and accounting for an F94L variant in MSTN

Explicitly fitting effects for major genes or QTL that account for a large percentage of variation in a whole genomic prediction model may increase prediction accuracy. This study compared approaches to account for a major effect of an F94L variant in the MSTN gene within the genomic prediction using bovine whole‐genomic SNP markers. Among the beef cattle breeds, Limousin have been known to have an F94L variant that is not present in Angus. The reference population in this study consisted of 3060 beef cattle including pure‐bred Limousin (PL), cross‐bred Limousin with Angus (LF) and pure‐bred Angus, genotyped using a BovineSNP50 BeadChip and directly for the MSTN‐F94L variant. We compared prediction accuracies in PL animals using the three datasets from only the PL population, admixed PL and LF (AL) or multibreed analysis using all of the PL, LF and Angus (MB) population according to four‐fold cross‐validation after K‐means clustering. The MSTN‐F94L variant was the most strongly associated with five traits (birth weight, calving ease direct, milk, weaning weight and yield grade) among the 13 measured traits in PL and AL populations. Fitting the MSTN‐F94L variant as a random effect, the genomic prediction accuracies for birth weight increased by 2.7% in PL, by 2.2% in AL and by 3.2% in MB. Prediction accuracies for five traits increased in the MB analysis. Fitting MSTN‐F94L as a fixed effect in PL, AL and MB analyses resulted in increased prediction accuracy in PL for eight traits. Prediction accuracies can be improved by including a causal variant in genomic evaluation compared with simply using whole‐genome SNP markers. Fitting the causal variant as a fixed effect along with markers fitted as random effects resulted in greater prediction accuracies for most traits. Causal variants should be genotyped along with SNP markers.