Single-nucleotide polymorphisms g.151435C>T and g.173057T>C in PRLR gene regulated by bta-miR-302a are associated with litter size in goats

Single-nucleotide polymorphisms (SNPs) located at microRNA-binding sites (miR-SNPs) can affect the expression of genes. This study aimed to identify the miR-SNPs associated with litter size. Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect SNPs in the caprine prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction, and polymerase chain reaction-restriction fragment length polymorphism. Three novel SNPs (g.151435C>T, g.151454A>G, and g.173057T>C) were identified in the caprine PRLR gene. Statistical results indicated that the g.151435C>T and g.173057T>C SNPs were significantly associated with litter size in Guanzhong and Boer goat breeds. Further analysis revealed that combinative genotype C6 (TTAACC) was better than the others for litter size in both goat breeds. Furthermore, the PRLR g.173057T>C polymorphism was predicted to regulate the binding activity of bta-miR-302a. Luciferase reporter gene assay confirmed that 173057C to T substitution disrupted the binding site for bta-miR-302a, resulting in the reduced levels of luciferase. Taken together, these findings suggested that bta-miR-302a can influence the expression of PRLR protein by binding with 3′untranslated region, resulting in that the g.173057T>C SNP had significant effects on litter size.     

Polymorphism identification in goat GNRH1 and GDF9 genes and their association analysis with litter size

This study investigated the polymorphisms of GNRH1 and GDF9 genes in 641 goats of three breeds: Xinong Saanen, Guanzhong and Boer. Two allelic variants were identified in the GNRH1 gene (JN645280:g.3548A>G and JN645281:g.3699G>A) and one allelic variant was found in the GDF9 gene (JN655693:g.4093G>A). Furthermore, g.4093G>A was a missense mutation (p.Val397Ile of GDF9). Results of an association analysis indicated that SNPs g.3548A>G and g.4093G>A had significant effects on litter size (P < 0.05). The combination genotypes of SNPs g.3548A>G, g.3699G>A and g.4093G>A also affected litter size with the C₅ genotype having the highest litter size in the first, third, fourth and average parity. Hence, the biochemical and physiological functions, together with the results obtained in our investigation, suggest that the GNRH1 and GDF9 genes could serve as genetic markers for litter size in goat breeding.     

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.     

Identification of a functional SNP in the 3′‐UTR of caprine MTHFR gene that is associated with milk protein levels

Xinong Saanen (n = 305) and Guanzhong (n = 317) dairy goats were used to detect SNPs in the caprine MTHFR 3′‐UTR by DNA sequencing. One novel SNP (c.*2494G>A) was identified in the said region. Individuals with the AA genotype had greater milk protein levels than did those with the GG genotype at the c.*2494 G>A locus in both dairy goat breeds (P < 0.05). Functional assays indicated that the MTHFR:c.2494G>A substitution could increase the binding activity of bta‐miR‐370 with the MTHFR 3′‐UTR. In addition, we observed a significant increase in the MTHFR protein level of AA carriers relative to that of GG carriers. These altered levels of MTHFR protein may account for the association of the SNP with milk protein level.     

Combined effects of four SNPs within goat PRLR gene on milk production traits

Xinong Saanen (SN, n = 323) and Guanzhong (GZ, n = 197) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions with their intron–exon boundaries of prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction (PIRA-PCR) and PCR-restriction fragment length polymorphism (PCR-RFLP). Four novel SNPs (g.40452T > C, g.40471G > A, g.61677G > A and g.61865G > A) were identified. The g.61677G > A and g.61865G > A SNPs caused amino acid variations p.Ser485Asn and p.Val548Met, respectively. Both g.40452T>C and g.40471G>A loci were closely linked in SN and GZ goat breeds (r² > 0.33). In addition, there was also a close linkage between g.61677G>A and g.61865G>A loci in both goat breeds. Statistical results indicated that the g.40452T > C, g.61677G > A and g.61865G > A SNPs were significantly associated with milk production traits in SN and GZ breeds. Further analysis revealed that combinative genotype C1 (TTAAGGGG) was better than the others for milk yield in SN and GZ goat breeds. These results are consistent with the regulatory function of PRLR in mammary gland development, milk secretion, and expression of milk protein genes, and extend the spectrum of genetic variation of the caprine PRLR gene, which might contribute to goat genetic resources and breeding.     

Sequence analysis of three canine adipokine genes revealed an association between TNF polymorphisms and obesity in Labrador dogs

Obesity is an emerging health problem in purebred dogs. Due to their crucial role in energy homeostasis control, genes encoding adipokines are considered candidate genes, and their variants may be associated with predisposition to obesity. Searching for polymorphism was carried out in three adipokine genes (TNF, RETN and IL6). The study was performed on 260 dogs, including lean (n = 109), overweight (n = 88) and obese (n = 63) dogs. The largest cohort was represented by Labrador Retrievers (n = 136). Altogether, 24 novel polymorphisms were identified: 12 in TNF (including one missense SNP), eight in RETN (including one missense SNP) and four in IL6. Distributions of five common SNPs (two in TNF, two in RETN and one in IL6) were further analyzed with regard to body condition score. Two SNPs in the non‐coding parts of TNF (c.‐40A>C and c.233+14G>A) were associated with obesity in Labrador dogs. The obtained results showed that the studied adipokine genes are highly polymorphic and two polymorphisms in the TNF gene may be considered as markers predisposing Labrador dogs to obesity.     

Genome‐wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.     

Targeted high‐throughput growth hormone 1 gene sequencing reveals high within‐breed genetic diversity in South African goats

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction‐targeted gene amplification together with Illumina MiSeq next‐generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27–58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within‐breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise F_ST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within‐population selection programs, particularly with SA village goats.     

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.     

Polymorphisms of caprine <Emphasis Type="Italic">GDF</Emphasis>9 gene and their association with litter size in Jining Grey goats

The exons 1, 2 and flanking region of growth differentiation factor 9 (GDF9) gene in five randomly selected does of Jining Grey, Boer and Liaoning Cashmere goats were amplified and analyzed. Thirteen nucleotide differences were identified in GDF9 gene between sheep (AF078545) and goats. Four SNPs (G3288A in intron 1, G423A, A959C [Gln320Pro] and G1189A [Val397Ile] in exon 2) were detected in four goat breeds with different prolificacy, in which G3288A was a new SNP in goats. The results showed that loci 3288, 423 and 1189 in Boer goats, loci 3288 and 423 in Guizhou White goats, loci 423 and 1189 in Liaoning Cashmere goats were all in complete linkage disequilibrium (D′ = 1, r ² = 1), respectively. In moderate (Boer goat) and low prolificacy (Liaoning Cashmere goat) breeds, linkage analysis indicated that there were more fervent linkage disequilibrium among loci 3288, 423 and 1189 than high prolificacy (Jining Grey and Guizhou White goats) breeds. For the 959 locus, the genotype distribution showed obvious difference between high prolificacy breeds and moderate or low prolificacy breeds (P < 0.05 or P < 0.01). The Jining Grey goat does with genotype CC or AC had 0.81 (P < 0.01) or 0.63 (P < 0.01) kids more than those with genotype AA, respectively. The present study preliminarily showed an association between allele C at 959 locus of GDF9 gene and high litter size in Jining Grey goats. These results provide further evidence that the GDF9 gene may be significantly correlated with high prolificacy in goats.     

GnRHR gene polymorphisms and their effects on reproductive performance in Chinese goats

In this study, polymorphisms in the goat GnRHR gene exon 1 were detected by PCR-SSCP and DNA sequencing methods in 786 individuals from two different goat breeds. Two haplotypes (A and B), two observed genotypes (AA and AB), and two single nucleotide polymorphisms (SNPs) were detected, which resulted in five amino acid substitutions. The frequencies of haplotypes A and B in the two goat breeds were 0.78–0.83 and 0.17–0.22, respectively. The SNP locus was in Hardy–Weinberg disequilibrium in the two goat breeds (P < 0.05). Polymorphisms of the GnRHR gene were shown to be associated with litter size in the two goat breeds. The SNPs in the goat GnRHR gene had significant effects on litter size (P < 0.05). Therefore, these results suggest that the GnRHR gene is a strong candidate gene that affects litter size in goat.     

Analysis of non‐synonymous SNPs of the porcine SERPINA6 gene as potential causal variants for a QTL affecting plasma cortisol levels on SSC7

Recently, the SERPINA6 gene encoding corticosteroid‐binding globulin (CBG) has been proposed as a candidate gene for a quantitative trait locus (QTL) affecting cortisol level on pig chromosome 7. The QTL was repeatedly detected in different lines, including a Piétrain × (German Landrace × German Large White) cross (PiF1) and purebred German Landrace (LR). In this study, we investigated whether the known non‐synonymous polymorphisms c.44G>T, c.622C>T, c.770C>T, c.793G>A, c.832G>A and c.919G>A of SERPINA6 are sufficient to explain the QTL in these two populations. Our investigations revealed that SNPs c.44G>T, c.622C>T, c.793G>A and c.919G>A are associated with cortisol level in PiF1 (P < 0.01). Haplotype analysis showed that these associations are largely attributable to differences between a major haplotype carrying SNPs c.793G>A and c.919G>A and a haplotype carrying SNPs c.44G>T and c.622C>T. Furthermore, some SNPs, particularly c.44G>T and c.622C>T and the carrier haplotype, showed association with meat quality traits including pH and conductivity (P < 0.05). In LR, the non‐synonymous SNPs segregate at very low frequency (<5%) and/or show only weak association with cortisol level (SNPs c.832G>A and c.919G>A; P < 0.05). These findings suggest that the non‐synonymous SNPs are not sufficient to explain the QTL across different breeds. Therefore, we examined whether the expression of SERPINA6 is affected by cis‐regulatory polymorphisms in liver, the major organ for CBG production. We found allelic expression imbalance of SERPINA6, which suggests that its expression is indeed affected by genetic variation in cis‐acting elements. This represents candidate causal variation for future studies of the molecular background of the QTL.     

PRLH and SOD1 gene variations associated with heat tolerance in Chinese cattle

With the proposed global climate change, heat tolerance is becoming increasingly important to the sustainability of livestock production systems. Results from previous studies showed that variants in the prolactin releasing hormone (PRLH) (AC_000160.1:g.11764610G>A) and superoxide dismutase 1 (SOD1) (AC_000158.1:g.3116044T>A) genes play an important role in heat tolerance in African indicine cattle. However, it is unknown whether or not the mutations are associated with heat tolerance in Chinese cattle. In this study, PCR and DNA sequencing were used to genotype two missense mutations in 725 individuals of 30 cattle breeds. Analysis results demonstrated that two classes of base substitution were detected at two loci: AC_000160.1:g.11764610G>A and AC_000158.1:g.3116044T>A or T>C respectively, with amino acid substitutions arginine to histidine and phenylalanine to isoleucine or leucine. The frequencies of the G and T alleles of the two loci gradually diminished from northern groups to southern groups of native Chinese cattle, whereas the frequencies of A and A or C alleles showed a contrary pattern, displaying a significant geographical difference across native Chinese cattle breeds. Additionally, analysis of these two loci in Chinese indigenous cattle revealed that two SNPs were significantly associated with mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) (P < 0.01), suggesting that cattle with A or C alleles were distributed in regions with higher T, RH and THI. Our results suggest that the two mutations of PRLH and SOD1 genes in Chinese cattle were associated with the heat tolerance.     

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.     

Association between METTL3 gene polymorphisms and neuroblastoma susceptibility: A nine‐centre case‐control study

Neuroblastoma ranks as the most commonly seen and deadly solid tumour in infancy. The aberrant activity of m⁶A‐RNA methyltransferase METTL3 is involved in human cancers. Therefore, functional genetic variants in the METTL3 gene may contribute to neuroblastoma risk. In the current nine‐centre case‐control study, we aimed to analyse the association between the METTL3 gene single nucleotide polymorphisms (SNPs) and neuroblastoma susceptibility. We genotyped four METTL3 gene SNPs (rs1061026 T>G, rs1061027 C>A, rs1139130 A>G, and rs1263801 G>C) in 968 neuroblastoma patients and 1814 controls in China. We found significant associations between these SNPs and neuroblastoma risk in neither single‐locus nor combined analyses. Interestingly, in the stratified analysis, we observed a significant risk association with rs1061027 AA in subgroups of children ≤ 18 months of age (adjusted OR = 1.87, 95% CI = 1.03‐3.41, P = .040) and females (adjusted OR = 1.86, 95% CI = 1.07‐3.24, P = .028). Overall, we identified a significant association between METTL3 gene rs1061027 C>A polymorphism and neuroblastoma risk in children ≤18 months of age and females. Our findings provide novel insights into the genetic determinants of neuroblastoma.     

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short‐tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D₅) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.