Genetic contributions to precocity traits in racing Thoroughbreds

Adaptation to early training and racing (i.e. precocity), which is highly variable in racing Thoroughbreds, has implications for the selection and training of horses. We hypothesised that precocity in Thoroughbred racehorses is heritable. Age at first sprint training session (work day), age at first race and age at best race were used as phenotypes to quantify precocity. Using high‐density SNP array data, additive SNP heritability () was estimated to be 0.17, 0.14 and 0.17 for the three traits respectively. In genome‐wide association studies (GWAS) for age at first race and age at best race, a 1.98‐Mb region on equine chromosome 18 (ECA18) was identified. The most significant association was with the myostatin (MSTN) g.66493737C>T SNP (P = 5.46 × 10^?12 and P = 1.89 × 10^?14 respectively). In addition, two SNPs on ECA1 (g.37770220G>A and g.37770305T>C) within the first intron of the serotonin receptor gene HTR7 were significantly associated with age at first race and age at best race. Although no significant associations were identified for age at first work day, the MSTN:g.66493737C>T SNP was among the top 20 SNPs in the GWAS (P = 3.98 × 10^?5). Here we have identified variants with potential roles in early adaptation to training. Although there was an overlap in genes associated with precocity and distance aptitude (i.e. MSTN), the HTR7 variants were more strongly associated with precocity than with distance. Because HTR7 is closely related to the HTR1A gene, previously implicated in tractability in young Thoroughbreds, this suggests that behavioural traits may influence precocity.     

Genome‐wide association QTL mapping for teat number in a purebred population of Duroc pigs

Because of increasing litter size in Western pig breeds, additional teats are desirable to increase the capacity for nursing offspring. We applied genome‐wide SNP markers to detect QTL regions that affect teat number in a Duroc population. We phenotyped 1024 animals for total teat number. A total of 36 588 SNPs on autosomes were used in the analysis. The estimated heritability for teat number was 0.34 ± 0.05 on the basis of a genomic relationship matrix constructed from all SNP markers. Using a BayesC method, we identified a total of 18 QTL regions that affected teat number in Duroc pigs; 9 of the 18 regions were newly detected.     

Genome‐wide association mapping and pathway analysis of leukosis incidence in a US Holstein cattle population

Bovine leukosis virus is an oncogenic virus that infects B cells, causing bovine leukosis disease. This disease is known to have a negative impact on dairy cattle production and, because no treatment or vaccine is available, finding a possible genetic solution is important. Our objective was to perform a comprehensive genetic analysis of leukosis incidence in dairy cattle. Data on leukosis occurrence, pedigree and molecular information were combined into multitrait GBLUP models with milk yield (MY) and somatic cell score (SCS) to estimate genetic parameters and to perform whole‐genome scans and pathway analysis. Leukosis data were available for 11 554 Holsteins daughters of 3002 sires from 112 herds in 16 US states. Genotypes from a 60K SNP panel were available for 961 of those bulls as well as for 2039 additional bulls. Heritability for leukosis incidence was estimated at about 8%, and the genetic correlations of leukosis disease incidence with MY and SCS were moderate at 0.18 and 0.20 respectively. The genome‐wide scan indicated that leukosis is a complex trait, possibly modulated by many genes. The gene set analysis identified many functional terms that showed significant enrichment of genes associated with leukosis. Many of these terms, such as G‐Protein Coupled Receptor Signaling Pathway, Regulation of Nucleotide Metabolic Process and different calcium‐related processes, are known to be related to retrovirus infection. Overall, our findings contribute to a better understanding of the genetic architecture of this complex disease. The functional categories associated with leukosis may be useful in future studies on fine mapping of genes and development of dairy cattle breeding strategies.     

猪的全基因组水平SNP研究现状

SNP（single nucleotide polymorphism,单核苷酸多态）在猪基因组中的分布极其广泛,平均分布间隔为300～400 bp,相关数据库收录已达55万条。猪基因组测序已取得实质性进展,大规模搜索发现基因组及EST（expressed sequence tag）序列中的SNP已展开,应用于猪全基因组水平的SNP芯片已建立。在此基础上,基于猪SNP标记的遗传图谱绘制、QTL（quantitative trait loci）定位、遗传多样性检测及全基因组关联分析等也都相继出现。     

Genome‐wide association study reveals candidate genes associated with body measurement traits in Chinese Wagyu beef cattle

We performed a genome‐wide association study to identify candidate genes for body measurement traits in 463 Wagyu beef cattle typed with the Illumina Bovine HD 770K SNP array. At the genome‐wide level, we detected 18, five and one SNPs associated with hip height, body height and body length respectively. In total, these SNPs are within or near 11 genes, six of which (PENK, XKR4, IMPAD1, PLAG1, CCND2 and SNTG1) have been reported previously and five of which (CSMD3, LAP3, SYN3, FAM19A5 and TIMP3) are novel candidate genes that we found to be associated with body measurement traits. Further exploration of these candidate genes will facilitate genetic improvement in Chinese Wagyu beef cattle.     

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.     

Imputation of 3 million SNPs in the Arabidopsis regional mapping population

Natural variation has become a prime resource to identify genetic variants that contribute to phenotypic variation. The regional mapping (RegMap) population is one of the most important populations for studying natural variation in Arabidopsis thaliana, and has been used in a large number of association studies and in studies on climatic adaptation. However, only 413 RegMap accessions have been completely sequenced, as part of the 1001 Genomes (1001G) Project, while the remaining 894 accessions have only been genotyped with the Affymetrix 250k chip. As a consequence, most association studies involving the RegMap are either restricted to the sequenced accessions, reducing power, or rely on a limited set of SNPs. Here we impute millions of SNPs to the 894 accessions that are exclusive to the RegMap, using the 1135 accessions of the 1001G Project as the reference panel. We assess imputation accuracy using a novel cross‐validation scheme, which we show provides a more reliable measure of accuracy than existing methods. After filtering out low accuracy SNPs, we obtain high‐quality genotypic information for 2029 accessions and 3 million markers. To illustrate the benefits of these imputed data, we reconducted genome‐wide association studies on five stress‐related traits and could identify novel candidate genes.     

Genetic factors controlling wool shedding in a composite Easycare sheep flock

Historically, sheep have been selectively bred for desirable traits including wool characteristics. However, recent moves towards extensive farming and reduced farm labour have seen a renewed interest in Easycare breeds. The aim of this study was to quantify the underlying genetic architecture of wool shedding in an Easycare flock. Wool shedding scores were collected from 565 pedigreed commercial Easycare sheep from 2002 to 2010. The wool scoring system was based on a 10‐point (0–9) scale, with score 0 for animals retaining full fleece and 9 for those completely shedding. DNA was sampled from 200 animals of which 48 with extreme phenotypes were genotyped using a 50‐k SNP chip. Three genetic analyses were performed: heritability analysis, complex segregation analysis to test for a major gene hypothesis and a genome‐wide association study to map regions in the genome affecting the trait. Phenotypes were treated as a continuous or binary variable and categories. High estimates of heritability (0.80 when treated as a continuous, 0.65–0.75 as binary and 0.75 as categories) for shedding were obtained from linear mixed model analyses. Complex segregation analysis gave similar estimates (0.80 ± 0.06) to those above with additional evidence for a major gene with dominance effects. Mixed model association analyses identified four significant (P < 0.05) SNPs. Further analyses of these four SNPs in all 200 animals revealed that one of the SNPs displayed dominance effects similar to those obtained from the complex segregation analyses. In summary, we found strong genetic control for wool shedding, demonstrated the possibility of a single putative dominant gene controlling this trait and identified four SNPs that may be in partial linkage disequilibrium with gene(s) controlling shedding.     

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ～36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.     

A complex structural variant at the KIT locus in cattle with the Pinzgauer spotting pattern

A specific white spotting phenotype, termed finching or line‐backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux‐Zillertaler cattle, two Austrian breeds. The so‐called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome‐wide association study using 27 white spotted and 16 solid‐coloured Tux‐Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole‐genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4‐kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310‐kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KIT^PINZ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux‐Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KIT^PINZ variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT^PINZ variant confirms admixture and the reported historical relationship of these short‐headed breeds with Austrian Tux‐Zillertaler and suggests a mutation event, occurring before breed formation.     

Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle

A putative functional mutation (rs109231213) near PLAG1 (BTA14) associated with stature was studied in beef cattle. Data from 8199 Bos taurus, Bos indicus and Tropical Composite cattle were used to test the associations between rs109231213 and various phenotypes. Further, 23 496 SNPs located on BTA14 were tested for association with these phenotypes, both independently and fitted together with rs109231213. The C allele of rs109231213 significantly increased hip height, weight, net food intake, age at puberty in males and females and decreased IGF‐I concentration in blood and fat depth. When rs109231213 was fitted as a fixed effect in the model, there was an overall reduction in associations between other SNPs and these traits but some SNPs remained associated (P < 10^?4). Frequency of the mutant C allele of rs109231213 differed among B. indicus (0.52), B. taurus (0.96) and Tropical Composite (0.68). Most chromosomes carrying the C allele had the same surrounding 10 SNP haplotype, probably because the C allele was introgressed into Brahman from B. taurus cattle. A region of reduced heterozygosity surrounds the C allele; this is small in B. taurus but 20 Mb long in Brahmans, indicating recent and strong selection for the mutant allele. Thus, the C allele appears to mark a mutation that has been selected almost to fixation in the B. taurus breeds studied here and introduced into Brahman cattle during grading up and selected to a frequency of 0.52 despite its negative effects on fertility.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Genome‐wide association mapping for phenotypic plasticity in rice

Phenotypic plasticity of plants in response to environmental changes is important for adapting to changing climate. Less attention has been paid to exploring the advantages of phenotypic plasticity in resource‐rich environments to enhance the productivity of agricultural crops. Here, we examined genetic variation for phenotypic plasticity in indica rice (Oryza sativa L.) across two diverse panels: (1) a Phenomics of Rice Adaptation and Yield (PRAY) population comprising 301 accessions; and (2) a Multi‐parent Advanced Generation Inter‐Cross (MAGIC) indica population comprising 151 accessions. Altered planting density was used as a proxy for elevated atmospheric CO₂ response. Low planting density significantly increased panicle weight per plant compared with normal density, and the magnitude of the increase ranged from 1.10 to 2.78 times among accessions for the PRAY population and from 1.05 to 2.45 times for the MAGIC population. Genome‐wide‐association studies validate three E nvironmental R esponsiveness (ER) candidate alleles (qER1–3) that were associated with relative response of panicle weight to low density. Two of these alleles were tested in 13 genotypes to clarify their biomass responses during vegetative growth under elevated CO₂ in Japan. Our study provides evidence for polymorphisms that control rice phenotypic plasticity in environments that are rich in resources such as light and CO₂.     

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non‐invasive phenotyping

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase‐specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non‐invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome‐wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker‐trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non‐parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage‐specific growth affecting genes to elucidate important processes operating at different developmental phases.     

A whole genome association analysis identifies loci associated with Mycobacterium avium subsp. paratuberculosis infection status in US holstein cattle

The purpose of this study was to identify loci associated with Mycobacterium avium subspecies paratuberculosis ( Map ) infection status in US Holsteins using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. Two hundred forty-five cows from dairies in New York, Pennsylvania and Vermont enrolled in longitudinal herd studies between January 1999 and November 2007 were assessed for the presence of Map in both faecal and tissue samples. An animal was considered tissue infected if any sample contained at least one colony forming unit of Map per gram of tissue (CFU/g) and the same definition was employed for faecal samples. Each animal was genotyped with the Illumina BovineSNP50 BeadChip and after quality assurance filtering, 218 animals and 45 683 SNPs remained. We sought to identify loci associated with four different case/control classifications: presence of Map in the tissue, presence of Map in faeces, presence of Map in both tissue and faeces and presence of Map in tissue but not faeces. A case–control genome wide association study was conducted to test the four different classifications of Map infection status (cases) when compared with a Map -negative control group (control). Regions on chromosomes 1, 5, 7, 8, 16, 21 and 23 were identified with moderate significance ( P  < 5 × 10⁻⁵). Two regions, one on chromosome 3 (near EDN2 ) and another on chromosome 9 (no positional gene candidates), were identified with a high level of association to the presence of Map in tissue and both tissue and faeces respectively ( P  <   5 × 10⁻⁷, genome-wide Bonferonni P  <   0.05).     

Downy mildew resistance is genetically mediated by prophylactic production of phenylpropanoids in hop

Downy mildew in hop (Humulus lupulus L.) is caused by Pseudoperonospora humuli and generates significant losses in quality and yield. To identify the biochemical processes that confer natural downy mildew resistance (DMR), a metabolome‐ and genome‐wide association study was performed. Inoculation of a high density genotyped F1 hop population (n = 192) with the obligate biotrophic oomycete P. humuli led to variation in both the levels of thousands of specialized metabolites and DMR. We observed that metabolites of almost all major phytochemical classes were induced 48 hr after inoculation. But only a small number of metabolites were found to be correlated with DMR and these were enriched with phenylpropanoids. These metabolites were also correlated with DMR when measured from the non‐infected control set. A genome‐wide association study revealed co‐localization of the major DMR loci and the phenylpropanoid pathway markers indicating that the major contribution to resistance is mediated by these metabolites in a heritable manner. The application of three putative prophylactic phenylpropanoids led to a reduced degree of leaf infection in susceptible genotypes, confirming their protective activity either directly or as precursors of active compounds.     

Heterogeneity of genetic architecture of body size traits in a free‐living population

Knowledge of the underlying genetic architecture of quantitative traits could aid in understanding how they evolve. In wild populations, it is still largely unknown whether complex traits are polygenic or influenced by few loci with major effect, due to often small sample sizes and low resolution of marker panels. Here, we examine the genetic architecture of five adult body size traits in a free‐living population of Soay sheep on St Kilda using 37 037 polymorphic SNPs. Two traits (jaw and weight) show classical signs of a polygenic trait: the proportion of variance explained by a chromosome was proportional to its length, multiple chromosomes and genomic regions explained significant amounts of phenotypic variance, but no SNPs were associated with trait variance when using GWAS. In comparison, genetic variance for leg length traits (foreleg, hindleg and metacarpal) was disproportionately explained by two SNPs on chromosomes 16 (s23172.1) and 19 (s74894.1), which each explained >10% of the additive genetic variance. After controlling for environmental differences, females heterozygous for s74894.1 produced more lambs and recruits during their lifetime than females homozygous for the common allele conferring long legs. We also demonstrate that alleles conferring shorter legs have likely entered the population through a historic admixture event with the Dunface sheep. In summary, we show that different proxies for body size can have very different genetic architecture and that dense SNP helps in understanding both the mode of selection and the evolutionary history at loci underlying quantitative traits in natural populations.     

Non‐replication of genome‐wide‐based associations of efficient food conversion in dairy cows

Animal growth relative to food energy input is of key importance to agricultural production. Several recent studies highlighted genetic markers associated with food conversion efficiency in beef cattle, and there is now a requirement to validate these associations in additional populations and to assess their potential utility for selecting animals with enhanced food‐use efficiency. The current analysis tested a population of dairy cattle using 138 DNA markers previously associated with food intake and growth in a whole‐genome association analysis of beef animals. Although seven markers showed point‐wise significance at P < 0.05, none of the single‐nucleotide polymorphisms tested were significantly associated with food conversion efficiency after correction for multiple testing. These data do not support the involvement of this subset of previously implicated markers in the food conversion efficiency of the physiologically distinct New Zealand Holstein‐Friesian dairy breed.