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.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

A genome‐wide association study suggests several novel candidate genes for carcass traits in Chinese Simmental beef cattle

A genome‐wide association study (GWAS) was conducted for two carcass traits in Chinese Simmental beef cattle. The experimental population consisted of 1301 individuals genotyped with the Illumina BovineHD SNP BeadChip (770K). After quality control, 671 990 SNPs and 1217 individuals were retained for the GWAS. The phenotypic traits included carcass weight and bone weight, which were measured after the cattle were slaughtered at 16 to 18 months of age. Three statistical models—a fixed polygene model, a random polygene model and a composite interval mapping polygene model—were used for the GWAS. The genome‐wide significance threshold after Bonferroni correction was 7.44E‐08 (= 0.05/671 990). In this study, we detected eight and seven SNPs significantly associated with carcass weight and bone weight respectively. In total, 11 candidate genes were identified within or close to these significant SNPs. Of these, we found several novel candidate genes, including PBX1, GCNT4, ALDH1A2, LCORL and WDFY3, to be associated with carcass weight and bone weight in Chinese Simmental beef cattle, and their functional roles need to be verified in further studies.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.     

Genome‐wide association study for reproductive traits in a Large White pig population

Using the PorcineSNP80 BeadChip, we performed a genome‐wide association study for seven reproductive traits, including total number born, number born alive, litter birth weight, average birth weight, gestation length, age at first service and age at first farrowing, in a population of 1207 Large White pigs. In total, we detected 12 genome‐wide significant and 41 suggestive significant SNPs associated with six reproductive traits. The proportion of phenotypic variance explained by all significant SNPs for each trait ranged from 4.46% (number born alive) to 11.49% (gestation length). Among them, 29 significant SNPs were located within known QTL regions for swine reproductive traits, such as corpus luteum number, stillborn number and litter size, of which one QTL region associated with litter size contained the ALGA0098819 SNP for total number born. Subsequently, we found that 376 functional genes contained or were near these significant SNPs. Of these, 14 genes—BHLHA15, OCM2, IL1B2, GCK, SMAD2, HABP2, PAQR5, GRB10, PRELID2, DMKN, GPI, GPIHBP1, ADCY2 and ACVR2B—were considered important candidates for swine reproductive traits based on their critical roles in embryonic development, energy metabolism and growth development. Our findings contribute to the understanding of the genetic mechanisms for reproductive traits and could have a positive effect on pig breeding programs.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.     

Genome‐wide association study to identify genomic regions affecting prolificacy in Lori‐Bakhtiari sheep

Several causative mutations in candidate genes affecting prolificacy have been detected in various sheep breeds. A genome‐wide association study was performed on estimated breeding values for litter size in Lori‐Bakhtiari sheep. Prolific ewes with twinning records and others with only singleton records were genotyped using the medium‐density Illumina Ovine SNP50 array. Four single nucleotide polymorphisms (SNPs) associated with litter size were identified on chromosomes 3, 6 and 22. The region on sheep chromosome 3 between 75 739 167 and 75 745 152 bp included two significant SNPs (s52383.1 and OAR3_80038014_X.1) in high linkage disequilibrium with each other. The region that surrounds these SNPs contains a novel putative candidate gene: luteinizing hormone/choriogonadotropin receptor (LHCGR), known to be involved in ovarian steroidogenesis and organism‐specific biosystem pathways in sheep. Known prolificacy genes BMPR1B, BMP15 and GDF9 were not associated with litter size in Lori‐Bakhtiari sheep, suggesting that other biological mechanisms could be responsible for the trait's variation in this breed.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

Genome‐wide screening for DNA variants associated with reading and language traits

Reading and language abilities are heritable traits that are likely to share some genetic influences with each other. To identify pleiotropic genetic variants affecting these traits, we first performed a genome‐wide association scan (GWAS) meta‐analysis using three richly characterized datasets comprising individuals with histories of reading or language problems, and their siblings. GWAS was performed in a total of 1862 participants using the first principal component computed from several quantitative measures of reading‐ and language‐related abilities, both before and after adjustment for performance IQ. We identified novel suggestive associations at the SNPs rs59197085 and rs5995177 (uncorrected P ≈ 10^–7 for each SNP), located respectively at the CCDC136/FLNC and RBFOX2 genes. Each of these SNPs then showed evidence for effects across multiple reading and language traits in univariate association testing against the individual traits. FLNC encodes a structural protein involved in cytoskeleton remodelling, while RBFOX2 is an important regulator of alternative splicing in neurons. The CCDC136/FLNC locus showed association with a comparable reading/language measure in an independent sample of 6434 participants from the general population, although involving distinct alleles of the associated SNP. Our datasets will form an important part of on‐going international efforts to identify genes contributing to reading and language skills.     

Genome‐wide association study for the level of serum electrolytes in Italian Large White pigs

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome‐wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance‐tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome‐wide efficient mixed‐model association. For the level of Ca²⁺, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (P_{nominal value} < 5.00E‐05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.     

Development of EST derived SSRs and SNPs as a genomic resource in Indian catfish, <Emphasis Type="Italic">Clarias batrachus</Emphasis>

Clarias batrachus, an Indian catfish species, is endemic to the Indian subcontinent and potential cultivable species. The genomic resources in C. batrachus in the form of ESTs containing microsatellite repeats (EST-SSR) and single nucleotide polymorphisms (SNPs) that are associated with the expressed genes from spleen were mined. From a total of 1,937 ESTs generated, 1,698 unique sequences were obtained, out of which 221 EST-SSRs were identified and 54% could be functionally annotated by similarity searches. A total of 23 contigs containing 3 or more ESTs were found to contain 31 SNP loci, out of which 8 ESTs showed similarity to genes of known function and 1 for hypothetical protein. Nine ESTs with SSRs and/or SNPs identified in this study were reported to be associated with diseases in human and animals. These identified loci can be developed into markers in C. batrachus, which can be useful in linkage mapping, comparative genomics studies and for its genetic improvement programmes.     

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.     

Identification of genes related to intramuscular fat independent of backfat thickness in Duroc pigs using single‐step genome‐wide association

Intramuscular fat (IMF) is an important meat‐quality trait of pigs, which influences pork’s shearing force, hydraulics, tenderness and juicy flavor. However, to achieve a higher percentage of lean meat, pigs with lower backfat thickness (BF) are intensively selected for, which may lead to a reduction in pork quality. Therefore, the objective of this study was to locate loci that affect IMF without changing BF. A single‐step GWAS was performed on 950 Duroc pigs genotyped by a 50K SNP chip in order to detect genomic variants relevant to IMF and BF. The significant SNPs detected were afterwards divided into a BF subset (seven SNPs), an IMF subset (11 SNPs) and a subset of both traits (12 SNPs), according to their P‐value and LD. After SNP and QTL annotation, our results indicated that SSC1: 167938652, 166363826, 164829874 and 167171587 might be associated with IMF without changing BF. In the subset of both traits, we found that the combined effect of ALGA0006602 (SSC1: 159538854) and 12784636 (SSC1: 160773437) might improve the IMF without changing BF. Our gene annotation result showed that TLE3, ITGA11, SMAD6, PAQR5 and [RNF152^A/G × MC4R^A/A] genes might affect IMF independently of BF. We believe that the SNPs and genes identified in this study will be valuable for the future molecular breeding of IMF in Duroc pigs.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

Genome-Wide Identification of Candidate Genes Associated with <i>β</i>-glucan Traits in a Hulled and Hulless Barley (<i>Hordeum vulgare</i> L.) Population

Barley grain is a valuable source of β-glucan, which is an important component of dietary fiber with significant human health benefits. Although the genetic basis of β-glucan biosynthesis has been widely studied, a genome-wide association study (GWAS) is still required for a scan of the candidate genes related to the complex quantitative trait based on the high-quality barley reference genome. In this study, a GWAS was conducted using a population composed of 87 barley landraces (39 hulled and 48 hulless, β-glucan from 2.07% to 6.56%) with 191,098 nucleotide polymorphisms (SNPs) markers to cover the chromosomes with the highest density. The population was divided into four sub-populations (POP1~POP4), and the β-glucan content in POP2 was significantly higher than that in other groups, in which most of the hulless barley landraces are from Qinghai-Tibet Plateau in China. Among seven SNP markers identified by GWAS, two (SNP2 and SNP3) of them showed positive correlation to β-glucan trait and the remaining five (SNP1, SNP4, SNP5, SNP6 and SNP7) showed the negative relationship. Two candidate genes linked to SNP7, HORVU7Hr1G000320 and HORVU7Hr1G000040, belong to the nucleotide triphosphate hydrolase superfamily which is probable to affect the activities of β-glucan synthase. Another candidate gene associated with SNP1, HORVU1Hr1G000010, is possibly involved in sugar response. In conclusion, our results provide new insights into the genetic basis of β-glucan accumulation in barley grains, and the discovery of new SNP markers distributed in each chromosome and the associated candidate genes will be valuable for the breeding of functional barley varieties with high β-glucan.