Software solutions for the livestock genomics SNP array revolution

Since the beginning of the genomic era, the number of available single nucleotide polymorphism (SNP) arrays has grown considerably. In the bovine species alone, 11 SNP chips not completely covered by intellectual property are currently available, and the number is growing. Genomic/genotype data are not standardized, and this hampers its exchange and integration. In addition, software used for the analyses of these data usually requires not standard (i.e. case specific) input files which, considering the large amount of data to be handled, require at least some programming skills in their production. In this work, we describe a software toolkit for SNP array data management, imputation, genome‐wide association studies, population genetics and genomic selection. However, this toolkit does not solve the critical need for standardization of the genotypic data and software input files. It only highlights the chaotic situation each researcher has to face on a daily basis and gives some helpful advice on the currently available tools in order to navigate the SNP array data complexity.     

Genome‐wide association analysis reveals novel loci for hypoxia adaptability in Tibetan chicken

The Tibetan chicken (TBC), an indigenous chicken breed of the Tibetan Plateau, has adapted to its hypoxic, high‐altitude environment over hundreds of years. The objective of this study was to identify the polymorphisms and genes associated with adaptation to hypoxia in this chicken breed. In the present study, samples were collected during days 18–21 of the incubation period from both surviving chicks and dead embryos, all of which were hatched under hypoxic conditions. A genome‐wide association study was conducted using the Illumina iSelect 60K SNP array with a case–control design, in which the case group consisted of the dead chicken embryos (n = 54) and controls were the surviving chicks (n = 82). Four significant SNPs were detected at the genome‐wide level (P < 0.05), and the results indicated that fork head box G1 (FOXG1) was the main candidate gene. The lead SNP NC_006092.4:g.33368893T>C was confirmed with a polymerase chain reaction‐restriction fragment length polymorphism analysis of 122 cases and 212 controls. A chi‐square test showed a significant association between NC_006092.4:g.33368893T>C and hatchability under hypoxic conditions (P < 0.01). Our results revealed novel polymorphisms and a candidate gene associated with hypoxic adaptation, facilitating further study in this field.     

Development and validation of the Axiom®Apple480K SNP genotyping array

Cultivated apple (Malus × domestica Borkh.) is one of the most important fruit crops in temperate regions, and has great economic and cultural value. The apple genome is highly heterozygous and has undergone a recent duplication which, combined with a rapid linkage disequilibrium decay, makes it difficult to perform genome‐wide association (GWA) studies. Single nucleotide polymorphism arrays offer highly multiplexed assays at a relatively low cost per data point and can be a valid tool for the identification of the markers associated with traits of interest. Here, we describe the development and validation of a 487K SNP Affymetrix Axiom^® genotyping array for apple and discuss its potential applications. The array has been built from the high‐depth resequencing of 63 different cultivars covering most of the genetic diversity in cultivated apple. The SNPs were chosen by applying a focal points approach to enrich genic regions, but also to reach a uniform coverage of non‐genic regions. A total of 1324 apple accessions, including the 92 progenies of two mapping populations, have been genotyped with the Axiom^®Apple480K to assess the effectiveness of the array. A large majority of SNPs (359 994 or 74%) fell in the stringent class of poly high resolution polymorphisms. We also devised a filtering procedure to identify a subset of 275K very robust markers that can be safely used for germplasm surveys in apple. The Axiom^®Apple480K has now been commercially released both for public and proprietary use and will likely be a reference tool for GWA studies in apple.     

High‐throughput genomics in sorghum: from whole‐genome resequencing to a SNP screening array

With its small, diploid and completely sequenced genome, sorghum (Sorghum bicolor L. Moench) is highly amenable to genomics‐based breeding approaches. Here, we describe the development and testing of a robust single‐nucleotide polymorphism (SNP) array platform that enables polymorphism screening for genome‐wide and trait‐linked polymorphisms in genetically diverse S. bicolor populations. Whole‐genome sequences with 6× to 12× coverage from five genetically diverse S. bicolor genotypes, including three sweet sorghums and two grain sorghums, were aligned to the sorghum reference genome. From over 1 million high‐quality SNPs, we selected 2124 Infinium Type II SNPs that were informative in all six source genomes, gave an optimal Assay Design Tool (ADT) score, had allele frequencies of 50% in the six genotypes and were evenly spaced throughout the S. bicolor genome. Furthermore, by phenotype‐based pool sequencing, we selected an additional 876 SNPs with a phenotypic association to early‐stage chilling tolerance, a key trait for European sorghum breeding. The 3000 attempted bead types were used to populate half of a dual‐species Illumina iSelect SNP array. The array was tested using 564 Sorghum spp. genotypes, including offspring from four unrelated recombinant inbred line (RIL) and F₂ populations and a genetic diversity collection. A high call rate of over 80% enabled validation of 2620 robust and polymorphic sorghum SNPs, underlining the efficiency of the array development scheme for whole‐genome SNP selection and screening, with diverse applications including genetic mapping, genome‐wide association studies and genomic selection.     

Identifying chromosomal selection‐sweep regions in facial eczema selection‐line animals using an ovine 50K‐SNP array

Facial eczema (FE) is a hepato‐mycotoxicosis found mainly in New Zealand sheep and cattle. When genetics was found to be a factor in FE susceptibility, resistant and susceptible selection lines of Romney sheep were established to enable further investigations of this disease trait. Using the Illumina OvineSNP50 BeadChip, we conducted a selection‐sweep experiment on these FE genetic lines. Two analytical methods were used to detect selection signals, namely the Peddrift test (Dodds & McEwan, 1997) and fixation index F_ST (Weir & Hill, 2002). Of 50 975 single nucleotide polymorphism (SNP) markers tested, there were three that showed highly significant allele frequency differences between the resistant and susceptible animals (Peddrift nominal P < 0.000001). These SNP loci are located on chromosomes OAR1, OAR11 and OAR12 that coincide precisely with the three highest genomic F_ST peaks. In addition, there are nine less significant Peddrift SNPs (nominal P ≤ 0.000009) on OAR6 (n = 2), OAR9 (n = 2), OAR12, OAR19 (n = 2), OAR24 and OAR26. In smoothed F_ST (five‐SNP moving average) plots, the five most prominent peaks are on OAR1, OAR6, OAR7, OAR13 and OAR19. Although these smoothed F_ST peaks do not coincide with the three most significant Peddrift SNP loci, two (on OAR6 and OAR19) overlap with the set of less significant Peddrift SNPs above. Of these 12 Peddrift SNPs and five smoothed F_ST regions, none is close to the FE candidate genes catalase and ABCG2; however, two on OAR1 and one on OAR13 fall within suggestive quantitative trait locus regions identified in a previous genome screen experiment. The present studies indicated that there are at least eight genomic regions that underwent a selection sweep in the FE lines.     

基于高密度SNP标记的肉牛人工选择痕迹筛查

近年来随着遗传改良工作的实施,人工选择大大提高了肉牛的生产性能并使其遗传基础发生巨大改变。文章基于Illumina BovineSNP50(54K)和BovineHD(770K)两款芯片数据,采用FST检验方法分析牛群的遗传分化,并筛查人工选择在牛的基因组留下的印记。通过全基因组范围内的扫描,共发现47 104个"离群"位点和3064个群体特异的人工选择"候选基因",如CLIC5、TG、CACNA2D1、FSHR等。通过基因注释对基因的生物学过程和分子功能进行富集分析。文章构建了我国肉牛的全基因组的选择信号图谱,为深入研究人工选择和理解生物进化提供线索,且研究结果也显示人工选择对基因组的影响在牛品种遗传改良中发挥了重要作用。     

Highly accurate SNP genotyping from historical and low‐quality samples

Historical and other poor‐quality samples are often necessary for population genetics, conservation, and forensics studies. Although there is a long history of using mtDNA from such samples, obtaining and genotyping nuclear loci have been considered difficult and error‐prone at best, and impossible at worst. The primary issues are the amount of nuclear DNA available for genotyping, and the degradation of the DNA into small fragments. Single nucleotide polymorphisms offer potential advantages for assaying nuclear variation in historical and poor‐quality samples, because the amplified fragments can be very small, varying little or not at all in size between alleles, and can be amplified efficiently by polymerase chain reaction (PCR). We present a method for highly multiplexed PCR of SNP loci, followed by dual‐fluorescence genotyping that is very effective for genotyping poor‐quality samples, and can potentially use very little template DNA, regardless of the number of loci to be genotyped. We genotyped 19 SNP loci from DNA extracted from modern and historical bowhead whale tissue, bone and baleen samples. The PCR failure rate was < 1.5%, and the genotyping error rate was 0.1% when DNA samples contained > 10 copies/µL of a 51‐bp nuclear sequence. Among samples with ≤ 10 copies/µL DNA, samples could still be genotyped confidently with appropriate levels of replication from independent multiplex PCRs.     

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 study reveals five nucleotide sequence variants for carcass traits in beef cattle

Genetic associations of nucleotide sequence variants with carcass traits in beef cattle were investigated using a genome-wide single nucleotide polymorphism (SNP) assay. Three hundred and thirteen Korean cattle were genotyped with the Illumina BovineSNP50 BeadChip, and 39,129 SNPs from 311 animals were analysed for each carcass phenotype after filtering by quality assurance. Five sequence markers were associated with one of the meat quantity or quality traits; rs109593638 on chromosome 3 with marbling score, rs109821175 on chromosome 11 and rs110862496 on chromosome 13 with backfat thickness (BFT), and rs110228023 on chromosome 6 and rs110201414 on chromosome 16 with eye muscle area (EMA) (P < 1.27 × 10(-6) , Bonferonni P < 0.05). The ss96319521 SNP, located within a gene with functions of muscle development, dishevelled homolog 1 (DVL1), would be a desirable candidate marker. Individuals with genotype CC at this gene appeared to have increased both EMA and carcass weight. Fine-mapping would be required to refine each of the five association signals shown in the current study for future application in marker-assisted selection for genetic improvement of beef quality and quantity.     

High‐density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool

In wheat, a lack of genetic diversity between breeding lines has been recognized as a significant block to future yield increases. Species belonging to bread wheat's secondary and tertiary gene pools harbour a much greater level of genetic variability, and are an important source of genes to broaden its genetic base. Introgression of novel genes from progenitors and related species has been widely employed to improve the agronomic characteristics of hexaploid wheat, but this approach has been hampered by a lack of markers that can be used to track introduced chromosome segments. Here, we describe the identification of a large number of single nucleotide polymorphisms that can be used to genotype hexaploid wheat and to identify and track introgressions from a variety of sources. We have validated these markers using an ultra‐high‐density Axiom^® genotyping array to characterize a range of diploid, tetraploid and hexaploid wheat accessions and wheat relatives. To facilitate the use of these, both the markers and the associated sequence and genotype information have been made available through an interactive web site.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

机器学习方法在基因交互作用探测中的研究进展

复杂疾病是基因与基因、基因与环境交互作用的结果,高维基因交互作用的探测给计算带来了极大的挑战。在过去20年间,机器学习方法被用于探测基因-基因交互作用,并取得了一定的效果。本文综述了机器学习方法在基因交互作用探测中的研究进展,系统地介绍了神经网络(neural networks, NN)、随机森林(random forest, RF)、支持向量机(support vector machines, SVM)和多因子降维法(multifactor dimensionality reduction, MDR)等机器学习方法在全基因组关联研究(genome wide association study, GWAS)中探测基因交互作用的原理和局限性,并对未来的研究进行了展望。     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High‐throughput high‐density genotyping arrays continue to be a fast, accurate, and cost‐effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high‐density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole‐genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives—honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance

Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non‐additive genetic variation, we present an approach using genome‐wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein–Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome‐wide markers were tested for associations between heterozygosity and bTB status using marker‐based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy‐Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome‐wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy‐number‐variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large‐scale data.     

Color phenotypes are under similar genetic control in two distantly related species of Timema stick insect

Ecology and genetics are both of general interest to evolutionary biologists as they can influence the phenotypic and genetic response to selection. The stick insects Timema podura and Timema cristinae exhibit a green/melanistic body color polymorphism that is subject to different ecologically based selective regimes in the two species. Here, we describe aspects of the genetics of this color polymorphism in T. podura, and compare this to previous results in T. cristinae. We first show that similar color phenotypes of the two species cluster in phenotypic space. We then use genome‐wide association mapping to show that in both species, color is controlled by few loci, dominance relationships between color alleles are the same, and SNPs associated with color phenotypes colocalize to the same linkage group. Regions within this linkage group that harbor genetic variants associated with color exhibit elevated linkage disequilibrium relative to genome wide expectations, but more strongly so in T. cristinae. We use these results to discuss predictions regarding how the genetics of color could influence levels of phenotypic and genetic variation that segregate within and between populations of T. podura and T. cristinae, drawing parallels with other organisms.