Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Candidate genes associated with testicular development, sperm quality, and hormone levels of inhibin, luteinizing hormone, and insulin-like growth factor 1 in brahman bulls

Bull fertility is an important target for genetic improvement, and early prediction using genetic markers is therefore a goal for livestock breeding. We performed genome-wide association studies to identify genes associated with fertility traits measured in young bulls. Data from 1118 Brahman bulls were collected for six traits: blood hormone levels of inhibin (IN) at 4 mo, luteinizing hormone (LH) following a gonadotropin-releasing hormone challenge at 4 mo, and insulin-like growth factor 1 (IGF1) at 6 mo, scrotal circumference (SC) at 12 mo, ability to produce sperm (Sperm) at 18 mo, and percentage of normal sperm (PNS) at 24 mo. All the bulls were genotyped with the BovineSNP50 chip. Sires and dams of the bull population (n = 304) were genotyped with the high-density chip (～800?000 polymorphisms) to allow for imputation, thereby contributing detail on genome regions of interest. Polymorphism associations were discovered for all traits, except for Sperm. Chromosome 2 harbored polymorphisms associated with IN. For LH, associated polymorphisms were located in five different chromosomes. A region of chromosome 14 contained polymorphisms associated with IGF1 and SC. Regions of the X chromosome showed associations with SC and PNS. Associated polymorphisms yielded candidate genes in chromosomes 2, 14, and X. These findings will contribute to the development of genetic markers to help select cattle with improved fertility and will lead to better annotation of gene function in the context of reproductive biology.     

A genome-wide SNP scan reveals novel loci for egg production and quality traits in white leghorn and brown-egg dwarf layers

Availability of the complete genome sequence as well as high-density SNP genotyping platforms allows genome-wide association studies (GWAS) in chickens. A high-density SNP array containing 57,636 markers was employed herein to identify associated variants underlying egg production and quality traits within two lines of chickens, i.e., White Leghorn and brown-egg dwarf layers. For each individual, age at first egg (AFE), first egg weight (FEW), and number of eggs (EN) from 21 to 56 weeks of age were recorded, and egg quality traits including egg weight (EW), eggshell weight (ESW), yolk weight (YW), eggshell thickness (EST), eggshell strength (ESS), albumen height(AH) and Haugh unit(HU) were measured at 40 and 60 weeks of age. A total of 385 White Leghorn females and 361 brown-egg dwarf dams were selected to be genotyped. The genome-wide scan revealed 8 SNPs showing genome-wise significant (P<1.51E-06, Bonferroni correction) association with egg production and quality traits under the Fisher's combined probability method. Some significant SNPs are located in known genes including GRB14 and GALNT1 that can impact development and function of ovary, but more are located in genes with unclear functions in layers, and need to be studied further. Many chromosome-wise significant SNPs were also detected in this study and some of them are located in previously reported QTL regions. Most of loci detected in this study are novel and the follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNPs.     

Genome Wide Analysis of Fertility and Production Traits in Italian Holstein Cattle

A genome wide scan was performed on a total of 2093 Italian Holstein proven bulls genotyped with 50K single nucleotide polymorphisms (SNPs), with the objective of identifying loci associated with fertility related traits and to test their effects on milk production traits. The analysis was carried out using estimated breeding values for the aggregate fertility index and for each trait contributing to the index: angularity, calving interval, non-return rate at 56 days, days to first service, and 305 day first parity lactation. In addition, two production traits not included in the aggregate fertility index were analysed: fat yield and protein yield. Analyses were carried out using all SNPs treated separately, further the most significant marker on BTA14 associated to milk quality located in the DGAT1 region was treated as fixed effect. Genome wide association analysis identified 61 significant SNPs and 75 significant marker-trait associations. Eight additional SNP associations were detected when SNP located near DGAT1 was included as a fixed effect. As there were no obvious common SNPs between the traits analyzed independently in this study, a network analysis was carried out to identify unforeseen relationships that may link production and fertility traits.     

A genome-wide association study of seed protein and oil content in soybean

Background

Association analysis is an alternative to conventional family-based methods to detect the location of gene(s) or quantitative trait loci (QTL) and provides relatively high resolution in terms of defining the genome position of a gene or QTL. Seed protein and oil concentration are quantitative traits which are determined by the interaction among many genes with small to moderate genetic effects and their interaction with the environment. In this study, a genome-wide association study (GWAS) was performed to identify quantitative trait loci (QTL) controlling seed protein and oil concentration in 298 soybean germplasm accessions exhibiting a wide range of seed protein and oil content.

Results

A total of 55,159 single nucleotide polymorphisms (SNPs) were genotyped using various methods including Illumina Infinium and GoldenGate assays and 31,954 markers with minor allele frequency >0.10 were used to estimate linkage disequilibrium (LD) in heterochromatic and euchromatic regions. In euchromatic regions, the mean LD (r ² ) rapidly declined to 0.2 within 360 Kbp, whereas the mean LD declined to 0.2 at 9,600 Kbp in heterochromatic regions. The GWAS results identified 40 SNPs in 17 different genomic regions significantly associated with seed protein. Of these, the five SNPs with the highest associations and seven adjacent SNPs were located in the 27.6-30.0 Mbp region of Gm20. A major seed protein QTL has been previously mapped to the same location and potential candidate genes have recently been identified in this region. The GWAS results also detected 25 SNPs in 13 different genomic regions associated with seed oil. Of these markers, seven SNPs had a significant association with both protein and oil.

Conclusions

This research indicated that GWAS not only identified most of the previously reported QTL controlling seed protein and oil, but also resulted in narrower genomic regions than the regions reported as containing these QTL. The narrower GWAS-defined genome regions will allow more precise marker-assisted allele selection and will expedite positional cloning of the causal gene(s).     

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole‐genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4–5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome‐wise basis, in combined‐sire (seven QTL) or within‐sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome‐wise basis, in combined‐sire (five QTL) or within‐sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin‐derived and Jersey‐derived alleles on these two chromosomes, respectively.     

Fine mapping QTL for female fertility on BTA04 and BTA13 in dairy cattle using HD SNP and sequence data

Background

Female fertility is important for the maintenance of the production in a dairy cattle herd. Two QTL regions on BTA04 and on BTA13 previously detected in Nordic Holstein (NH) and validated in the Danish Jersey (DJ) and Nordic Red (NR) were investigated further in the present study to further refine the QTL locations. Refined QTL regions were imputed to the full sequence data. The genes in the regions were then studied to ascertain their possible effect on fertility traits.

Results

BTA04 was screened for number of inseminations (AIS), 56-day non-return rate (NRR), days from first to last insemination (IFL), and the interval from calving to first insemination (ICF) in the range of 38,257,758 to 40,890,784 bp, whereas BTA13 was screened for ICF only in the range from 21,236,959 to 46,150,079 with the HD bovine SNP array for NH, DJ and NR. No markers in the DJ and NR breeds reached significance. By analyzing imputed sequence data the QTL position on BTA04 was narrowed down to two regions in the NH. In these two regions a total of 9 genes were identified. BTA13 was analyzed using sequence data for the NH breed. The highest –log₁₀(P-value) was 19.41 at 33,903,159 bp. Two regions were identified: Region 1: 33,900,143-33,908,994 bp and Region 2: 34,051,815-34,056,728 bp. SNPs within and between these two regions were annotated as intergenic.

Conclusion

Screening BTA04 and BTA13 for female fertility traits in NH, NR and DJ suggested that the QTL for female fertility were specific for NH. A missense mutation in CD36 showed the strongest association with fertility traits on BTA04. The annotated SNPs on BTA13 were all intergenic variants. It is possible that BTA13 at this stage is poorly annotated such that the associated polymorphisms are located in as-yet undiscovered genes. Fertility traits are complex traits as many different biological and physiological factors determine whether a cow is fertile. Therefore it is not expected that there is a simple explanation with an obvious candidate gene but it is more likely a network of genes and intragenic variants that explain the variation of these traits.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-790) contains supplementary material, which is available to authorized users.     

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

ABSTRACT: BACKGROUND: Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years. RESULTS: QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups. CONCLUSION: Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).     

Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample population size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein–Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with traditional measures of fertility, as well as a fertility phenotype derived from milk progesterone profiles. Traditional fertility measures investigated were days to first heat, days to first service, pregnancy rate to first service, number of services and calving interval (CI); post-partum interval to the commencement of luteal activity (CLA) was derived using routine milk progesterone assays. Phenotypic and genotypic data on 37 590 single nucleotide polymorphisms (SNPs) were available for up to 1570 primiparous cows. Genetic parameters were estimated using linear animal models, and univariate and bivariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. Heritability estimates of the traditional fertility traits varied from 0.03 to 0.16; the heritability for CLA was 0.13. The posterior quantitative trait locus (QTL) probabilities, across the genome, for the traditional fertility measures were all <0.021. Posterior QTL probabilities of 0.060 and 0.045 were observed for CLA on SNPs each on chromosome 2 and chromosome 21, respectively, in the univariate analyses; these probabilities increased when CLA was included in the bivariate analyses with the traditional fertility traits. For example, in the bivariate analysis with CI, the posterior QTL probability of the two aforementioned SNPs were 0.662 and 0.123. Candidate genes in the vicinity of these SNPs are discussed. The results from this study suggest that the power of genome-wide association studies in cattle may be increased by sharing of data and also possibly by using physiological measures of the trait under investigation.     

Long‐term selection for litter size in swine results in shifts in allelic frequency in regions involved in reproductive processes

High‐density genotype data were analyzed in three lines of swine that express substantial variation in sow fertility to uncover regions of the genome potentially influenced during selection for litter size traits. The experimental lines examined include the Nebraska Index Line (NIL), which has been subjected to long‐term selection for litter size; a control line derived from the same population that founded NIL; and a commercial Duroc × Hampshire (D × H) population, in which no selection for litter size was practiced. Regions of the genome potentially affected by selection for litter size traits in NIL were determined by multiple lines of evidence, including altered allelic frequency compared to the other lines, loss of heterozygosity and relative extended haplotype homozygosity. Additionally, a genome‐wide association study for litter size traits was conducted in a population based on NIL and commercial maternal line genetics. Several genomic regions identified as putative signatures of selection overlapped with QTL for litter size traits. One of these regions, located on SSC2 (13–14 Mb), includes the candidate gene P2X3R, which plays a role in implantation and sustained release of hormones associated with reproductive processes. Sequencing identified synonymous SNPs in P2X3R that are fixed in NIL but polymorphic with nearly equal frequencies in the D × H line, indicating a potential role of P2X3R in sow fertility. These results suggest that data derived from these lines can help to uncover and understand a portion of the genetic variance associated with fertility traits in swine.     

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

A genome-wide association study for feed efficiency-related traits in a crossbred pig population

Feed efficiency (FE) is one of the most important traits in pig production. However, it is difficult and costly to measure it, limiting the collection of large amount of data for an accurate selection for better FE. Therefore, the identification of single-nucleotide polymorphisms (SNPs) associated with FE-related traits to be used in the genetic evaluation is of great interest of pig breeding programs for increasing the prediction accuracy and the genetic progress of these traits. The objective of this study was to identify SNPs significantly associated with FE-related traits: average daily gain (ADG), average daily feed intake (ADFI) and feed conversion ratio (FCR). We also aimed to identify potential candidate genes for these traits. Phenotypic information recorded on a population of 2386 three-way crossbreed pigs that were genotyped for 51 468 SNPs was used. We identified three loci of quantitative trait (QTL) regions associated with ADG and three QTL regions associated with ADFI; however, no significant association was found for FCR. A false discovery rate (FDR) ≤ 0.005 was used as the threshold for declaring an association as significant. The QTL regions associated with ADG on Sus scrofa chromosome (SSC) 1 were located between 177.01 and 185.47 Mb, which overlaps with the QTL regions for ADFI on SSC1 (173.26 and 185.47 Mb). The other QTL region for ADG was located on SSC12 (2.87 and 3.22 Mb). The most significant SNPs in these QTL regions explained up to 3.26% of the phenotypic variance of these traits. The non-identification of genomic regions associated with FCR can be explained by the complexity of this trait, which is a ratio between ADG and ADFI. Finally, the genes CDH19, CDH7, RNF152, MC4R, PMAIP1, FEM1B and GAA were the candidate genes found in the 1 Mb window around the QTL regions identified in this study. Among them, the MC4R gene (SSC1) has a well-known function related to ADG and ADFI. In this study, we identified three QTL regions for ADG (SSC1 and SSC12) and three for ADFI (SSC1). These regions were previously described in purebred pig populations; however, to our knowledge, this is the first study to confirm the relevance of these QTL regions in a crossbred pig population. The potential use of the SNPs and genes identified in this study in prediction models that combine genomic selection and marker-assisted selection should be evaluated for increasing the prediction accuracy of these traits in this population.     

High-throughput single nucleotide polymorphism (SNP) identification and mapping in the sesame (<Emphasis Type="Italic">Sesamum indicum</Emphasis> L.) genome with genotyping by sequencing (GBS) analysis

Sesame (Sesamum indicum L. syn. Sesamum orientale L.) is considered to be the first oil seed crop known to man. Despite its versatile use as an oil seed and a leafy vegetable, sesame is a neglected crop and has not been a subject of molecular genetic research until the last decade. There is thus limited knowledge regarding genome-specific molecular markers that are indispensible for germplasm enhancement, gene identification, and marker-assisted breeding in sesame. In this study, we employed a genotyping by sequencing (GBS) approach to a sesame recombinant inbred line (RIL) population for high-throughput single nucleotide polymorphism (SNP) identification and genotyping. A total of 15,521 SNPs were identified with 14,786 SNPs (95.26 %) located along sesame genome assembly pseudomolecules. By incorporating sesame-specific simple sequence repeat (SSR) markers developed in our previous work, 230.73 megabases (99 %) of sequence from the genome assembly were saturated with markers. This large number of markers will be available for sesame geneticists as a resource for candidate polymorphisms located along the physical chromosomes of sesame. Defining SNP loci in genome assembly sequences provides the flexibility to utilize any genotyping strategy to survey any sesame population. SNPs selected through a high stringency filtering protocol (770 SNPs) for improved map accuracy were used in conjunction with SSR markers (50 SSRs) in linkage analysis, resulting in 13 linkage groups that encompass a total genetic distance of 914 cM with 432 markers (420 SNPs, 12 SSRs). The genetic linkage map constitutes the basis for future work that will involve quantitative trait locus (QTL) analyses of metabolic and agronomic traits in the segregating RIL population.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

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

A candidate gene association study for nine economically important traits in Italian Holstein cattle

We genotyped 58 single nucleotide polymorphisms (SNPs) in 25 candidate genes in about 800 Italian Holstein sires. Fifty‐six (minor allele frequency >0.02) were used to evaluate their association with single traits: milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP), milk protein percentage (PP), milk somatic cell count (MSCC); and complex indexes: longevity, fertility and productivity–functionality type (PFT), using deregressed proofs, after adjustment for familial relatedness. Thirty‐two SNPs were significantly associated (proportion of false positives <0.05) with different traits: 16 with MSCC, 15 with PY, 14 with MY, 12 with PFT, eight with longevity, eight with FY, eight with PP, five with FP and two with fertility. In particular, a SNP in the promoter region of the PRLR gene was associated with eight of nine traits. DGAT1 polymorphisms were highly associated with FP and FY. Casein gene markers were associated with several traits, confirming the role of the casein gene cluster in affecting milk yield, milk quality and health traits. Other SNPs in genes located on chromosome 6 were associated with PY, PP, PFT, MY (PPARGC1A) and MSCC (KIT). This latter association may suggest a biological link between the degree of piebaldism in Holstein and immunological functions affecting somatic cell count and mastitis resistance. Other significant SNPs were in the ACACA, CRH, CXCR1, FASN, GH1, LEP, LGB (also known as PAEP), MFGE8, SRC, TG, THRSP and TPH1 genes. These results provide information that can complement QTL mapping and genome‐wide association studies in Holstein.     

Genome-wide scan for bovine twinning rate QTL using linkage disequilibrium

Twinning is a complex trait with negative impacts on health and reproduction, which cause economic loss in dairy production. Several twinning rate quantitative trait loci (QTL) have been detected in previous studies, but confidence intervals for QTL location are broad and many QTL are unreplicated. To identify genomic regions or genes associated with twinning rate, QTL analysis based on linkage combined with linkage disequilibrium (LLD) and individual marker associations was conducted across the genome using high-throughput single nucleotide polymorphism (SNP) genotypes. A total of 9919 SNP markers were genotyped with 200 sires and sons in 19 half-sib North American Holstein dairy cattle families. After SNPs were genotyped, informative markers were selected for genome-wide association tests and QTL searches. Evidence for twinning rate QTL was found throughout the genome. Thirteen markers significantly associated with twinning rate were detected on chromosomes 2, 5 and 14 ( P  < 2.3 × 10⁻⁵). Twenty-six regions on fourteen chromosomes were identified by LLD analysis at P  < 0.0007. Seven previously reported ovulation or twinning rate QTL were supported by results of single marker association or LLD analyses. Single marker association analysis and LLD mapping were complementary tools for the identification of putative QTL in this genome scan.     

Identification of loci governing eight agronomic traits using a GBS‐GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome‐wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping‐by‐sequencing (GBS) approach was used to provide dense genome‐wide marker coverage (>47 000 SNPs) for a panel of 304 short‐season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Genome-wide association study identifies candidate markers for bull fertility in Holstein dairy cattle

The decline in the reproductive efficiency of dairy cattle has become a challenging problem worldwide. Female fertility is now taken into account in breeding goals while generally less attention is given to male fertility. The objective of this study was to perform a genome-wide association study in Holstein bulls to identify genetic variants significantly related to sire conception rate (SCR), a new phenotypic evaluation of bull fertility. The analysis included 1755 sires with SCR data and 38,650 single nucleotide polymorphisms (SNPs) spanning the entire bovine genome. Associations between SNPs and SCR were analyzed using a mixed linear model that included a random polygenic effect and SNP genotype either as a linear covariate or as a categorical variable. A multiple testing correction approach was used to account for the correlation between SNPs because of linkage disequilibrium. After genome-wide correction, eight SNPs showed significant association with SCR. Some of these SNPs are located close to or in the middle of genes with functions related to male fertility, such as the sperm acrosome reaction, chromatin remodeling during the spermatogenesis, and the meiotic process during male germ cell maturation. Some SNPs showed marked dominance effects, which provide more evidence for the relevance of non-additive effects in traits closely related to fitness such as fertility. The results could contribute to the identification of genes and pathways associated with male fertility in dairy cattle.