Association of a single nucleotide polymorphism in the akirin 2 gene with economically important traits in Korean native cattle

The akirin 2 gene, located on chromosome 9 in cattle, was previously reported to be associated with nuclear factor‐kappa B (NF‐κB), involved in immune reactions and marbling of meat. To determine whether a single nucleotide polymorphism (SNP) in akirin 2 is associated with economically important traits of Korean native cattle, the c.*188G>A SNP DNA marker in the 3′‐UTR region of akirin 2 was analyzed for its association with carcass weight, longissimus muscle area and marbling. The c.*188G>A SNP was genotyped by polymerase chain reaction restriction fragment length polymorphism, and the frequency of the AA, AG, and GG genotypes were 6.82%, 71.29% and 21.88% respectively. This SNP was significantly associated with longissimus muscle area (Bonferroni corrected P < 0.05), and marbling score (Bonferroni corrected P < 0.01). These results suggest that the c.*188G>A SNP of akirin 2 might be useful as a DNA marker for longissimus muscle area and marbling scores in Korean native cattle.     

Two mutations in the 5′‐flanking region of the KITLG gene are associated with litter size of dairy goats

In this study, Xinong Saanen (SN) and Guanzhong (GZ) dairy goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the 5′‐flanking region of the KITLG gene by DNA sequencing and primer‐introduced restriction analysis–polymerase chain reaction. Two novel SNPs (g.13090G>T and g.13664C>A) were identified (GenBank Accession no. KM658964). Furthermore, g.13090G>T and g.13664C>A loci were closely linked in SN and GZ breeds (r² > 0.33). Association analysis results showed that g.13090G>T and g.13664C>A SNPs significantly affected litter size (P < 0.05). The litter size of individuals with the combined genotype GG/CC from both dairy goat breeds was greater than that of individuals with TT/AA in average parity (P < 0.05). Known biochemical and physiological functions, along with our results, indicated that GG/CC could be used in marker‐assisted selection to choose individuals with greater litter size from both breeds. These results extend the spectrum of genetic variation in the caprine KITLG gene and may contribute to genetic resources and breeding of goats.     

Single nucleotide polymorphisms in bovine <Emphasis Type="Italic">PRL</Emphasis> gene and their associations with milk production traits in Chinese Holsteins

Prolactin (PRL) plays a crucial role in the initiation and maintenance of lactation in mammals. In this study, seven PCR fragments representing most important functional domains of PRL gene were screened for single nucleotide polymorphisms (SNPs) in Chinese Holsteins by single-strand conformation polymorphisms and amplicons sequencing, and their genetic effects on milk production traits were evaluated. A total of four SNPs, including two in the promoter (−1043A>G and −402A>G), one in intron 1 (+2723C>T) and one in exon 4 (+8398G>A) were identified in PRL gene. Statistical results showed significant associations between the promoter genotypes and the milk performance traits in Chinese Holsteins. Cows with genotype P1-GG showed higher milk yields (P < 0.01), while cows with genotype P1-AA showed higher fat contents (P < 0.01). Haplotype analysis of two SNPs in promoter region revealed that the Hap(AG) was significantly associated with increased milk yields and Hap(AA) was associated with increased fat contents (P < 0.01). This is the second study reporting SNPs in the 5′-regulatory region of PRL gene, which interfere with milk production traits.     

SNP genotypes of olfactory receptor genes associated with olfactory ability in German Shepherd dogs

To find out the relationship between SNP genotypes of canine olfactory receptor genes and olfactory ability, 28 males and 20 females from German Shepherd dogs in police service were scored by odor detection tests and analyzed using the Beckman GenomeLab SNPstream. The representative 22 SNP loci from the exonic regions of 12 olfactory receptor genes were investigated, and three kinds of odor (human, ice drug and trinitrotoluene) were detected. The results showed that the SNP genotypes at the OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR2K2‐like:c.518G>A, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A loci had a statistically significant effect on the scenting abilities (P < 0.001). The kind of odor influenced the performances of the dogs (P < 0.001). In addition, there were interactions between genotype and the kind of odor at the following loci: OR10H1‐like:c.632C>T, OR10H1‐like:c.770A>T, OR4C11‐like:c.511T>G and OR4C11‐like:c.692G>A (P < 0.001). The dogs with genotype CC at the OR10H1‐like:c.632C>T, genotype AA at the OR10H1‐like:c.770A>T, genotype TT at the OR4C11‐like:c.511T>G and genotype GG at the OR4C11‐like:c.692G>A loci did better at detecting the ice drug. We concluded that there was linkage between certain SNP genotypes and the olfactory ability of dogs and that SNP genotypes might be useful in determining dogs' scenting potential.     

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.     

Genetic variation in the β, β‐carotene‐9′, 10′‐dioxygenase gene and association with fat colour in bovine adipose tissue and milk

β, β‐carotene‐9′, 10′‐dioxygenase (BCO2) plays a role in cleaving β‐carotene eccentrically, and may be involved in the control of adipose and milk colour in cattle. The bovine BCO2 gene was sequenced as a potential candidate gene for a beef fat colour QTL on chromosome (BTA) 15. A single nucleotide base change located in exon 3 causes the substitution of a stop codon (encoded by the A allele) for tryptophan⁸⁰ (encoded by the G allele) (c. 240G>A, p.Trp80stop, referred to herein as SNP W80X). Association analysis showed significant differences in subcutaneous fat colour and beta‐carotene concentration amongst cattle with different BCO2 genotypes. Animals with the BCO2 AA genotype had more yellow beef fat and a higher beta‐carotene concentration in adipose tissues than those with the GA or GG genotype. QTL mapping analysis with the BCO2 SNP W80X fitted as a fixed effect confirmed that this SNP is likely to represent the quantitative trait nucleotide (QTN) for the fat colour‐related traits on BTA 15. Moreover, animals with the AA genotype had yellower milk colour and a higher concentration of beta‐carotene in the milk.     

Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

DEFB1 polymorphisms are involved in susceptibility to human papillomavirus infection in Brazilian gynaecological patients

The human beta defensin 1 (hBD-1) antimicrobial peptide is a member of the innate immune system known to act in the first line of defence against microorganisms, including viruses such as human papillomavirus (HPV). In this study, five functional polymorphisms (namely g-52G>A, g-44C>G and g-20G>A in the 5’UTR and c.*5G>A and c.*87A>G in the 3’UTR) in the DEFB1 gene encoding for hBD-1 were analysed to investigate the possible involvement of these genetic variants in susceptibility to HPV infection and in the development of HPV-associated lesions in a population of Brazilian women. The DEFB1 g-52G>A and c.*5G>A single-nucleotide polymorphisms (SNPs) and the GCAAA haplotype showed associations with HPV-negative status; in particular, the c.*5G>A SNP was significantly associated after multiple test corrections. These findings suggest a possible role for the constitutively expressed beta defensin-1 peptide as a natural defence against HPV in the genital tract mucosa.     

Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

The CDKN2A/p16INK4a 5′UTR sequence and translational regulation: impact of novel variants predisposing to melanoma

Many variants of uncertain functional significance in cancer susceptibility genes lie in regulatory regions, and clarifying their association with disease risk poses significant challenges. We studied 17 germline variants (nine of which were novel) in the CDKN2A 5′UTR with independent approaches, which included mono and bicistronic reporter assays, Western blot of endogenous protein, and allelic representation after polysomal profiling to investigate their impact on CDKN2A mRNA translation regulation. Two of the novel variants (c.‐27del23, c.‐93‐91delAGG) were classified as causal mutations (score ≥3), along with the c.‐21C>T, c.‐34G>T, and c.‐56G>T, which had already been studied by a subset of assays. The novel c.‐42T>A as well as the previously described c.‐67G>C were classified as potential mutations (score 1 or 2). The remaining variants (c.‐14C>T, c.‐20A>G, c.‐25C>T+c.‐180G>A, c.‐30G>A, c.‐40C>T, c.‐45G>A, c.‐59C>G, c.‐87T>A, c.‐252A>T) were classified as neutral (score 0). In conclusion, we found evidence that nearly half of the variants found in this region had a negative impact on CDKN2A mRNA translation, supporting the hypothesis that 5′UTR can act as a cellular Internal Ribosome Entry Site (IRES) to modulate p16^INK4a translation.     

Allelic polymorphisms in the <Emphasis Type="Italic">MTHFR</Emphasis>, <Emphasis Type="Italic">MTR</Emphasis> and <Emphasis Type="Italic">MTRR</Emphasis> genes in patients with cleft lip and/or palate and their mothers

The distribution of polymorphic variants in the MTHFR, MTR and MTRR genes in the contingents of patients with nonsyndromic cleft lip and/or palate (NSCL/P), their mothers, and healthy individuals from Ukraine’s western region was evaluated. It was shown that the presence of the homozygous MTHFR 677TT genotype might lead to a threefold increase in the risk of CL/P; mothers carrying the MTHFR 677TT genotype had a twofold increase in the risk of giving birth to a child with CL/P compared to homozygous carriers of MTHFR 677CC (OR = 3.3 and OR = 1.92, respectively). The presence of the heterozygous genotype MTR 2756AG can lead to a 1.5-fold increase in the risk of CL/P compared to the 2765AA genotype (OR = 1.48). The heterozygous genotype MTRR 66AG is associated with a 5-fold risk of CL/P (OR = 5.56); mothers carrying this genotype had a 2.6-fold increase in the risk of giving birth to a child with CL/P (OR = 2.6). The prevalence of the MTRR 66G allelic variant among inhabitants of Ukraine’s western region was higher compared to the MTRR 66A allelic variant (wild type); and the MTRR 66GG genotype frequency among CL/P patients was significantly lower compared to the control group.     

Plasma total homocysteine level and methylenetetrahydrofolate reductase 677C>T genetic polymorphism in Japanese patients with rheumatoid arthritis

Hyperhomocysteinemia is a known risk factor of cardiovascular disease. Homocysteine has been also linked to inflammation in rheumatoid arthritis (RA). In this study, we investigated the relationship between plasma homocysteine levels and single nucleotide polymorphism (SNP) of the gene coding for methylenetetrahydrofolate reductase (MTHFR), an enzyme involved in the biosynthesis of homocysteine, and the correlation between the plasma homocysteine levels and generally used inflammatory markers (C-reactive protein, erythrocyte sedimentation rate and matrix metalloproteinase-3) in 96 Japanese patients with RA. Plasma homocysteine levels in patients with the MTHFR 677TT genotype were significantly higher than in those with the 677CC genotype (p < 0.05). In addition, plasma homocysteine levels were increased along with the elevation of general inflammatory markers. Therefore, we conclude that homocysteine might affect the inflammatory status of patients, and the MTHFR 677C>T SNP could be a predictive factor of hyperhomocysteinemia in patients with RA.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Identification of four SNPs and association analysis with meat quality traits in the porcine <Emphasis Type="Italic">Pitx2c</Emphasis> gene

The association of the porcine Pitx2c gene with meat quality traits was investigated in the present study. A total of eight single nucleotide polymorphisms (SNPs) were found. Allele frequencies of four SNPs were further detected in four commercial breeds and eight Chinese indigenous breeds. Single SNP and meat quality associations were analyzed in a Yorkshire×Meishan F₂ population. The SNPs c.474C>T (P<0.01) and c.636C>T (P<0.05) showed a significant association with meat color (MCV1). The SNPs c.*37G>A and c.*47G>A were significantly associated with drip loss rate (DLR), water holding capacity (WHC) and meat color value (MCV1) consistently (P<0.05). Linkage disequilibrium (LD) analysis revealed that the adjacent SNPs were in LD. Two major haplotypes were identified, and association analysis between haplotype combinations and meat quality indicated that the presence of two copies of haplotype 1 -CCGG- may improve meat quality.     

Effects of genetic variability of the dairy goat growth hormone releasing hormone receptor (<Emphasis Type="Italic">GHRHR</Emphasis>) gene on growth traits

Growth hormone-releasing hormone receptor (GHRHR) plays a critical role in growth hormone (GH) synthesis, release and regulation of pituitary somatotroph expansion in vertebrates. The objective of this study was to investigate variations in goat GHRHR gene and their associations with growth traits in 668 dairy goats. The results showed four novel single nucleotide polymorphisms (SNPs): NC_007302:g.5203C>T, 7307C>G, 9583G>A and 9668A>C. In detail, the novel SNP C>T in the 5203rd nucleotide identified a missense mutation: CCC (Pro)>TCC (Phe) at position 116aa of the goat GHRHR (423aa). Besides, 9583G>A and 9668A>C polymorphism were in complete linkage disequilibrium. The genetic diversity analysis revealed that the Guanzhong dairy goat possessed intermediate genetic diversity in P3 and P7 loci, and the Xinong Sannen dairy goat belonged to poor genetic diversity in P4 locus. Significant associations between the genotypes of P3 locus and body length, body height and chest circumference was observed in Guanzhong goat (P < 0.05). However, in Xinong saanen population, significant statistical difference was only found in body height and body length (P < 0.05). In P4 and P7 loci, no significant associations were detected between any variant sites and body length, body height and chest circumference, as well as for the milk traits (P > 0.05). These results strongly suggested that the goat GHRHR gene is a candidate gene that influences growth traits in dairy goat.     

Associations of polymorphisms in the promoter I of bovine acetyl‐CoA carboxylase‐α gene with beef fatty acid composition

The objectives of this study were to identify single nucleotide polymorphisms (SNPs) in the promoter I (PI) region of the bovine acetyl‐CoA carboxylase‐α (ACACA) gene and to evaluate the extent to which they were associated with lipid‐related traits. Eight novel SNPs were identified, which were AJ276223:g.2064T>A (SNP1), g.2155C>T (SNP2), g.2203G>T (SNP3), g.2268T>C (SNP4), g.2274G>A (SNP5), g.2340A>G (SNP6), g.2350T>C (SNP7) and g.2370A>G (SNP8). Complete linkage disequilibrium was observed among SNP1, 2, 4, 5, 6 and 8. Phenotypic data were collected from 573 cross‐bred steers with six sire breeds, including Hereford, Angus, Brangus, Beefmaster, Bonsmara and Romosinuano. The genotypes of SNP1/2/4/5/6/8 were significantly associated with adjusted backfat thickness. The genotypes of SNP3 were significantly associated with triacylglycerol (TAG) content and fatty acid composition of longissimus dorsi muscle (LM) in Brangus‐, Romosinuano‐ and Bonsmara‐sired cattle. Cattle with g.2203GG genotype had greater concentrations of TAG, total lipid, total saturated fatty acid and total monounsaturated fatty acid than did cattle with g.2203GT genotype. The genotypes of SNP7 were significantly associated with fatty acid composition of LM. Cattle with genotype g.2350TC had greater amounts of several fatty acids in LM than did cattle with genotype g.2350CC. Our results suggested that the SNPs in the PI region of ACACA gene are associated with variations in the fatty acid contents in LM.     

LIN28A gene polymorphisms modify neuroblastoma susceptibility: A four‐centre case‐control study

Neuroblastoma ranks the most common seen solid tumour in childhood. Overexpression of LIN28A gene has been linked to the development of multiple human malignancies, but the relationship between LIN28A single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility is still under debate. Herein, we evaluated the correlation of four potentially functional LIN28A SNPs (rs3811464 G>A, rs3811463 T>C, rs34787247 G>A, and rs11247957 G>A) and neuroblastoma susceptibility in 505 neuroblastoma patients and 1070 controls from four independent hospitals in China. The correlation strengths were determined by using odds ratios (ORs) and corresponding 95% confidence intervals (CIs). Among these SNPs, rs34787247 G>A exhibited a significant association with increased susceptibility in neuroblastoma (GA vs GG: adjusted OR = 1.30, 95% CI = 1.03‐1.64; AA vs GG: adjusted OR = 2.51, 95% CI = 1.36‐4.64, AA/GA vs GG: adjusted OR = 1.42, 95% CI = 1.12‐1.80, AA vs GG/GA: adjusted OR = 2.39, 95% CI = 1.29‐4.42). Furthermore, the combined analysis of risk genotypes revealed that subjects carrying three risk genotypes (adjusted OR = 1.64, 95% CI = 1.02‐2.63) are more inclined to develop neuroblastoma than those without risk genotype, and so do carriers of 1‐4 risk genotypes (adjusted OR = 1.26, 95% CI = 1.01‐1.56). Stratification analysis further revealed risk effect of rs3811464 G>A, rs34787247 G>A and 1‐4 risk genotypes in some subgroups. Haplotype analysis of these four SNPs yields two haplotypes significantly correlated with increased neuroblastoma susceptibility. Overall, our finding indicated that LIN28A SNPs, especially rs34787247 G>A, may increase neuroblastoma risk.