Genome-Wide Association Study of Body Weight Traits in Inner Mongolia Cashmere Goats

Objective: Body weight is an important economic trait for a goat, which greatly affects animal growth and survival. The purpose of this study was to identify genes associated with birth weight (BW), weaning weight (WW), and yearling weight (YW).

Materials and Methods: In this study, a genome-wide association study (GWAS) of BW, WW, and YW was determined using the GGP_Goat_70K single-nucleotide polymorphism (SNP) chip in 1,920 Inner Mongolia cashmere goats.

Results: We discovered that 21 SNPs were significantly associated with BW on the genome-wide levels. These SNPs were located in 10 genes, e.g., Mitogen-Activated Protein Kinase 3 (MAPK3), LIM domain binding 2 (LDB2), and low-density lipoprotein receptor-related protein 1B (LRP1B), which may be related to muscle growth and development in Inner Mongolia Cashmere goats. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these genes were significantly enriched in the regulation of actin cytoskeleton and phospholipase D signaling pathway etc.

Conclusion: In summary, this study will improve the marker-assisted breeding of Inner Mongolia cashmere goats and the molecular mechanisms of important economic traits.

Introduction

Body weight is one of the most important economic traits for livestock that can be measured during the entire animal lifetime from birth to slaughter. The early growth rate of a goat has a strong implication on both reproductive and production performances. Early growth performance traits, such as birth weight (BW) and weaning weight (WW), are the basis for selection in genetic improvement programs for meat production (1). To investigate the relationships between body weight and meat quality traits, PCR-restriction fragment length polymorphism and PCR-single strand conformation polymorphism were applied to suggest a low negative relationship between BW, meat quality, and genetic markers (IGF-II and CAPN1) in chickens breeding for meat quality (2). However, conventional breeding methods cannot make significant progress in a short period of time. To a large extent, it hinders the rapid industry development of animal production.

Genome-wide association study (GWAS) is a tool of genome-wide linkage disequilibrium to determine genetic marker information (e.g., single-nucleotide polymorphism, SNP, loci, and copy number) and related genes that genetically affect complex phenotypic traits or important quantitative economic traits throughout the genome (3–6). In 2005, the first GWAS report of macular degeneration of the retina was published in Science (7). With the completion of multi-species genome sequencing, the continuous development of high-throughput SNP chips, and genome typing technologies, GWAS has become a popular solution to the genetic positioning of complex human diseases and animal quantitative traits. For example, great progress has been made in GWAS research on obesity (8), coronary heart disease (9), diabetes (10), and other diseases, successfully identifying complex diseases significant SNPs. GWASs apply SNP arrays to accurately screen and identify the major genes of important economic traits, which becomes a prerequisite for rapid improvement for the molecular breeding of domestic animals, e.g., cattle, pigs, chickens, and sheep. For instance, A GWAS of five meat quality traits in 231 Yorkshire pigs identified 344 significant SNPs associated with five meat quality traits using PorcineSNP60 BeadChip, in which 323 SNPs were located in the reported QTL regions and 21 were novel (11). In a GWAS of the 600 K Affymetrix Axiom Chicken Genotyping Array, KCTD4, LIM domain binding 2 (LDB2), HEP21, and PCASP2 were closely correlated with the spleen weight of Layer chicken (12). Ghasemi et al. carried out a GWAS on the BW of Lori-Bakhtiari sheep and identified three genes (RAB6B, Tf serotransferrin, and GIGYF2) as candidate genes for this trait (13). GWAS of body weight in two cattle populations from the Russian Federation (Siberian region) was determined using the GGP HD150 and identified five statistically significant SNPs and the CCND2 genes as the candidate gene for body weight trait (14).

The Inner Mongolia cashmere goat is a local breed in China with dual-purpose, producing both pleasant meat and excellent cashmere. It is famous for its “white cashmere and meat” in Inner Mongolia. Besides, Inner Mongolia cashmere goats also provide sufficient meat resources and have become one of the most important livestock species for the herdsmen. However, there are no GWAS studies focusing on the body weight traits of Inner Mongolia cashmere goats. Therefore, in this study, we examined BW, WW, and yearling weight (YW) by GGP_Goat_70K SNP chip and GWAS methodology for Inner Mongolia cashmere goat population of 1,920 individuals to identify significant SNPs and the major candidate genes associated with body weight traits in cashmere goats. This study may facilitate the potential use of major genes involving in growth and production traits for genetic improvement of productivity in cashmere goats.

Materials and Methods

Ethics Statement

In this experiment, the breeding environment was in compliance with the standards relevant to an ordinary animal laboratory facility in China National Standard “Laboratory Animal Environment and Facilities” (GB14925-2010). The feeding and the experimental operations on animals were in accordance with the animal welfare requirements.

Phenotypic Measurements and Sample Collection

Goats used in the present study (n = 1,920, 2 years old from six herds; 1 year old from three herds) were selected randomly from the Arbas Stock Farm in Inner Mongolia, China. The BW, WW (3.5 months), and YW (12 months) were measured using an electronic scale for all goats. BW was recorded within 0.5 h after birth, while WW and YW were measured after 12 h of fasting. Samples of ear tissue were collected by ear deficiency forceps and quickly placed in a prepared cryopreservation tube containing 75% alcohol for storage at −80°C until DNA extraction.

Statistical Analysis

Generally, the fixed effect of each trait in the model was identified by General Linear Model (GLM) procedure using the Statistical Analysis System (SAS) program. The influencing factors included the year of production (three levels, 2018-2020), herd (six herds, 1–6), birth status (simple or twin), maternal ages (six levels, 3–8), and sex (male and female, 1–2). Therefore, repeatability and multivariate animal model was used with restricted maximum likelihood method (AIREML) in the WOMBAT software to estimate the variance components of each trait. Then heritability for each trait was obtained under this model. The default convergence criterion was 10⁻⁸. The formula for this model was as follow:

$\begin{array}{l}{y}_i=X_i{b}_i+Z_i{a}_i+Y_i{m}_i+W_i{p}_i+e_i\end{array}$

Where y_i is the vector of observation of animals for trait i; b_i is the vector of fixed effects for trait i; a_i is the vector of direct additive genetic effects of animal for trait i; m_i is the vector of random maternal genetic effects of animal for trait i; and p_i is the vector of individual performant environmental effects for trait i. Y_i is maternal genetic effects for trait i, X_i, Z_i, and W_i are design matrices of the corresponding effects, respectively. e_i is the vector of random residuals effects for trait i. The heritability (h²), was calculated using the formula $h^2=\frac{{\sigma }_a^2}{{\sigma }_p^2}$.

Genotyping and Quality Control

Deoxyribonucleic acid was extracted from ear tissue with the standard phenol-chloroform method according to the protocol of the manufacturer. DNA integrity and purity were tested using 2% agarose gel electrophoresis and a NanoDrop 2000 ultraviolet spectrophotometer (Thermo, Waltham, MA, USA).

The samples were genotyped using the Illumina Goat SNP 70K BeadChip panel that included 67,088 SNPs (Inner Mongolia Agricultural University, China). The samples with call rates <90% were removed from the analysis. The SNPs with GenCall (GC) scores <0.6, genotype call rates <90%, minor allele frequencies (MAF) < 0.01, and significant Hardy-Weinberg (HWE) disequilibrium at 10⁻⁵ were removed from the analysis. Plink 1.90 beta (15) and R software were used for quality control.

GWAS

The correlations between the SNPs and the traits were tested using mixed linear models in GEMMA software version 0.98 (16–18). Therefore, herd and age were combined as fixed effects, as they all significantly related to body weight traits. The statistical analysis model used in this study was y = Xα + Zβ + Wμ + e, where y is the phenotypic trait, X is a matrix of fixed effects, α is the estimation parameter of the fixed effects, Z is a matrix of SNPs, β is the effect of the SNPs, W is a matrix of random effects, μ is the predicted random individuals, and e is the random error, with the distribution e~ N (0, ${\delta }_{\text{e}}^2$). The significance threshold for the GWAS was defined using the Bonferroni correction method. The suggestive genome-wide association significance threshold was P < 9.92 × 10⁻⁷ (0.05/50383), and the chromosome-wide significance level threshold was P < 2.88 × 10⁻⁵ (0.05/50,383/29). Chromosome-wide significance level SNPs were defined to call chromosome-wide significance associations, with a suggestive association corresponds cutoff P < 10⁻⁴ (19). The quantile–quantile (Q–Q) plots were visualized by plotting the distribution of obtained vs. expected log₁₀ (P-value) with inflation factors (λ). The association map and the significant SNPs were visualized in the Manhattan plot with a threshold line. The Manhattan and Q–Q graphics were generated with R v. 3.5.2.

Bioinformatics Analysis

Genes associated with significantly correlated SNP loci were annotated with the goat reference genome (ARS1, GCF_001704415.1). The genes were used as input for Cytoscape that analyzes and visualizes Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (20, 21). GO terms and KEGG pathways with adjusting P < 0.05 were statistically significant.

Results

Descriptive Statistics and Heritability for Weight Traits of Cashmere Goats

Observed phenotypes were analyzed by first descriptive statistics (Table 1). The BW, WW, and YW ranged from 1.7 to 4.2, 16.5 to 29.5, and 24.0 to 49.5 kg, respectively. The coefficients of variation (CV) for BW, WW, and YW were 17.45, 57.93, and 37.98%, respectively. The results indicate that substantial phenotypic variation of these three traits exists in the population of Inner Mongolia cashmere goats. The heritability estimates for BW, WW, and YW were 0.10, 0.27, and 0.10, respectively.

TABLE 1

Table 1. Descriptive statistics and variance components of body weight traits in Inner Mongolia cashmere goats.

Quality Control of Genotyping

The body weight traits of the 1,920 Inner Mongolia cashmere goats measured in this study included BW, WW, and YW. A total of 67,088 SNPs were genotyped. Subsequently, 2,655 SNPs were excluded from our dataset as they did not pass the HWE, while 10,446 SNPs were excluded from our dataset as they did not pass the HWE MAF tests (–MAF 0.01–HWE 1e⁻⁵). After the quality control (QC) was performed on the raw genotypes, a total of 1,909 animals and 50,383 SNPs were obtained and were distributed over the 29 goat chromosomes (Supplementary Figure 1). The descriptive statistics of the studied traits are listed in Table 1. All traits are normally distributed. The PCA result demonstrated that there were no genetic differences between the samples (Supplementary Figure 2).

GWAS

Based on the 50,383 SNPs, GWAS was then performed with a mixed linear model for the BW, WW, and YW. A total of 21 SNPs reached the genome-wide significance levels for BW traits and 50 SNPs reached the chromosome-wide significance levels for the three traits (Figure 1). For BW, 21 genome-wide significance levels correlated SNPs were detected on the chromosomes Chr1, Chr2, Chr3, Chr4, Chr6, Chr7, Chr8, Chr9, Chr10, Chr12, Chr14, Chr18, Chr21, Chr25, and Chr29, these SNPs were electively annotated to 10 genes, e.g., Mitogen-Activated Protein Kinase 3 (MAPK3), LDB2,low-density lipoprotein receptor-related protein 1B (LRP1B; Figure 1A, Table 2).

FIGURE 1

Figure 1. Manhattan plots of body weight traits for Inner Mongolia cashmere goats. (A) Birth weight (BW); (B) weaning weight (WW); and (C) yearling weight (YW).

TABLE 2

Table 2. Significant single-nucleotide polymorphisms (SNPs), associated traits, and candidate genes identified in genome-wide association studies (GWASs).

In addition, 46 correlated SNPs were identified at chromosome-wide significance levels. For WW, two chromosome-wide significance levels correlated SNPs identified on Chr18 were related to two genes, CRISPLD2 and TUB-like protein 2 (TULP2, Figure 1B, Table 2). For YW, two chromosome-wide significance levels correlated SNPs were detected on Chr4 and Chr11, corresponding to one gene, progestogen-associated endometrial protein (PAEP; Figure 1C, Table 2). The Q–Q plot showed that the genomic inflation factors (λ) of BW, WW, and YW were 0.983, 0.990, 1.017, respectively (Figure 2), indicating no genome expansion. Based on the observed and expected P-values in the Q–Q plots of the three body weight traits, there were no population stratification phenomena. However, this phenomenon was significant for the SNPs, which were highly correlated with the three body weight traits.

FIGURE 2

Figure 2. Quantile–quantile plots of body weight traits using a mixed linear model approach. Blue dots represent the –log₁₀ (P-value) of the entire study, and the red line represents the expected values under the null hypothesis of no association. (A) Birth weight (BW); (B) weaning weight (WW); and (C) yearling weight (YW).

Functional Enrichment Analysis of GO and KEGG

In an attempt to better understand the biological functions and signaling pathways of the trait-associated genes, we performed GO and KEGG enrichment analysis for 13 genes near the SNPs of the three body weight traits, and uncovered enriched 23 GO terms, i.e., 11 biological processes, 8 cellular components, and 4 molecular functions (Figures 3, 4), which are closely correlated to growth and organ formation.

FIGURE 3

Figure 3. Gene ontology (GO) enrichment analysis for the regional candidate genes with chromosome-wide significant association.

FIGURE 4

Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for the regional candidate genes with the chromosome-wide significant association.

Discussion

Body weight is an important economic trait in goat production. Therefore, it is important to understand the underlying molecular mechanisms behind goat body weight and identify important functional genes that influence cashmere goat production and breeding. However, compared with reported GWAS results for sheep body weight (6, 13, 22), the goat candidate genes of body weight are not well-studied. Therefore, in this study, genome-wide associations of BW, WW, and YW were determined using the GGP_Goat_70K SNP chip in 1920 Inner Mongolia cashmere goats.

In this study, the heritability of WW is moderately high (0.27), which is similar to the values estimated by Ma et al. (23). The heritabilities of BW and YW (0.11, and 0.10, respectively) are lower than the values estimated by Hossein and Mohammad et al. (24) (0.22 and 0.25, respectively), and Ma et al. (23) (0.647 and 0.571, respectively). The differences between studies are probably caused by the source of goat breeds, model design, and data structure.

Maternal age only has significant effects on BW and WW in early traits, but not YW, indicating great influences of maternal effects on lamb BW and WW. Studies discovered that with the increase of the age of mother, BW and WW also increase, suggesting lower BW and WW of lambs born to primiparous ewes than that of adult ewes. The reason is that the physical maturity and physical maturity of primiparous ewes are not well-developed. In short, lactating lambs have a strong dependence on the ewe, strengthen the feeding and management, improve the nutritional status, and promote lactation, which is conducive to the normal development of lamb.

To correct the effects of population structure and individual kinship relationships, a population genetic structure is considered as a fixed effect while an individual kinship is thought of as a random effect. Therefore, the GWAS of BW, WW, and YW of 1,920 Inner Mongolia cashmere goats was performed using a mixed linear model in this study. The Inner Mongolia cashmere goats that we chose are mainly grazing animals and the nutritional level of their pastures varies among the seasons, which will further affect gene expressions (25, 26).

Body weight growth is closely related to the growth of obesity, fat, and muscle. The candidate genes identified in this study were also closely correlated with growth and tissue development. For instance, the TULP2 gene is located on chromosome 18 and is correlated with WW. TULP2 is a member of the tubby gene family, and gene SNPs in the 19q13.33–13.43 chromosomal region are significantly related to severe obesity in French Caucasians (27). Jackson et al. demonstrate that CRISPLD2 is a circulating adipokine that may regulate adipocyte remodeling during weight loss (28).

The LDB2 gene is located on chromosome 6 and correlated with BW. LDB2, also known as CLIM-1, was identified as a LIM domain-associated cofactor and functions as a transcriptional regulatory factor (29–31). The LIM domain-binding factor 2 (LDB2) gene in the chicken chromosome 4 (GGA4) region ~8.6 Mb in length (71.6–80.2 Mb) had the strongest association with body weight for weeks 7–12 and with average daily gain for weeks 6–12 (32). Liu et al. identify LDB2 genes on chicken (Gallus gallus) chromosome 4 associated with Carcass Weight (CW) and Eviscerated Weight (EW) traits, GWASs were conducted using the Illumina 60 K SNP Bead chip to genotype 724 Beijing-You chickens (33). In Wang's et al. they performed a new GWAS using specific locus amplified fragment sequencing (SLAF-seq) technology to discover that the LDB2 gene in this region had a very strong association with body weight (34). Wei et al. indicated that a 31-bp indel located in the second intron region of the LDB2 gene was significantly correlated with some growth traits and carcass traits of chickens (35).

The LRP1B gene is located on chromosome 2 and correlated with BW. Houde et al. found epigenetic variations at LRP1B, a gene associated with the development of obesity or cardiometabolic complications, are involved in fetal metabolic programming (36). Wang et al. performed a GWAS in 82 sows with an extreme SD of BWs within the first parity to identify LRP1B may contribute to the BW variability trait (37). LRP1B gene encodes LRP1B and mediates cellular cholesterol uptake (38). Dietrich et al. reported that knockout of LRP1B in mice results in early embryonic lethality (39). Association analysis identified LRP1B as a determinant of rat cholesterol concentrations in low-density lipoproteins (LDL) and a significant association with child body mass index (BMI) in humans (40). CNR1 is involved in the growth cone (GO: 0030426). CNR1 can affect the regulation of leptin signaling to regulate BW, food intake, and metabolism (41). Besides, the genes associated with SNP are involved in the biological process. MAPK3 is involved in fundamental cellular processes, such as membrane trafficking, actin cytoskeleton remodeling, cell proliferation and cell survival, animal organ formation, and glucose metabolism. In this study, Prominin 1 (PROM1), FBXL3, and LOC102176015 genes are new candidate genes related to BW. PROM1 is a protein-coding gene, which plays a role in cell differentiation, proliferation, and apoptosis. FBXL3 encodes a member of the F-box protein family, which plays a key role in the maintenance of both the speed and the robustness of the circadian clock oscillation.

Conclusion

In summary, the first GWAS of Inner Mongolia cashmere goats was performed to identify gene-associated to three body traits. We identified 21 SNPs that reached the genome-wide significance levels and 50 SNPs that reached the chromosome-wide significance levels for the three traits in a population of 1,920 Inner Mongolia cashmere goats. MAPK3, LDB2, and LRP1B genes are key candidate genes for BW. PROM1, FBXL3, and LOC102176015 genes are new genes related to BW and can be given to follow-up priority research. The candidate genes in this study were closely correlated to BW development. These findings will make a significant contribution to the understanding of the mechanisms underlying BW traits and facilitate genetics improvement of productivity in Inner Mongolia cashmere goats by providing essential genes related to goat BW traits.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: figshare.com/s/7734276ce6c1239e9474.

Ethics Statement

All samples were collected in accordance with the International Guiding Principles for Biomedical Research involving animals and approved by the Special Committee on Scientific Research and Academic Ethics of Inner Mongolia Agricultural University, responsible for the approval of Biomedical Research Ethics of Inner Mongolia Agricultural University [Approval No. (2020) 056]. No specific permissions were required for these activities, and no endangered or protected species were involved.

Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

This work was financially supported by National Natural Science Foundation of China (31860637); Science and Technology Major Project of Inner Mongolia Autonomous Region 2021SZD0010); Science and Technology Project of Inner Mongolia Autonomous Region (2019GG243); China Agriculture Research System of MOF and MARA (CARS-39); Scientific Project of Inner Mongolia Agricultural University on High-level Introduced Talented Personnel (NDYB2018-1); Central Government Guides Local Science and Technology Development Fund Projects (2020ZY0007); Youth Foundation in Inner Mongolia Agricultural University (No. QN202003). Special Funding Project for the Iconic Achievements of the College of Animal Science, Inner Mongolia Agricultural University (No. BZCG202111).

Conflict of Interest

LZ and LH were employed by company Inner Mongolia Jinlai Livestock Technology.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors are grateful to the staff of Inner Mongolia Jinlai Livestock Technology Co for providing assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.752746/full#supplementary-material

References

1. Buzanskas ME, Grossi DA, Ventura RV, Schenkel FS, Sargolzaei M, Meirelles SLC, et al. Genome-wide association for growth traits in Canchim beef cattle. PLoS ONE. (2014) 9:e94802. doi: 10.1371/journal.pone.0094802

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Molee A, Kuadsantia P, Kaewnakian P. Gene effects on body weight, carcass yield, and meat quality of thai indigenous chicken. J Poult Sci. (2017) 55:94–102. doi: 10.2141/jpsa.0160159

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhang ZW, Xu LY, Yang J, Gao HJ, Zhang LP, Gao X, et al. Weighted single-step genome-wide association study for growth traits in chinese simmental beef cattle. Genes. (2020) 11:189. doi: 10.3390/genes11020189

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Tam V, Patel N, Turcotte M, Bossé Y, Paré G, Meyre D. Benefits and limitations of genome-wide association studies. Nat Rev Genet. (2019) 20:467–84. doi: 10.1038/s41576-019-0127-1

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zhang L, Liu JS, Zhao FP, Ren HX, Xu LY, Lu J, et al. Genome-Wide association studies for growth and meat production traits in sheep. PLoS ONE. (2013) 8:e66569. doi: 10.1371/journal.pone.0066569

PubMed Abstract | CrossRef Full Text | Google Scholar

6. AI-Mamun HA, Kwan P, Clark SA, Ferdosi MH, Tellam R, Gondro C. Genome-wide association study of body weight in Australian merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight. Genet Select Evolut. (2015) 14:66–6. doi: 10.1186/s12711-015-0142-4

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, Haynes C, et al. Complement factor H polymorphism in age-related macular degeneration. Science. (2005) 308:385–9. doi: 10.1016/j.ajo.2005.06.004

CrossRef Full Text | Google Scholar

8. Scuteri A, Sanna S, Chen WM, Uda M, Albai G, Strait J, et al. Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet. (2007) 3:e115. doi: 10.1371/journal.pgen.0030115

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, et al. Genomewide association analysis of coronary artery disease. N Engl J Med. (2007) 357:443–53. doi: 10.1056/NEJMoa072366

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Sladek R, Rocheleau G, Rung J, Dina C, Shen LS, Serre D, et al. A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. (2007) 445:881–5. doi: 10.1038/nature05616

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Dong Q, Liu HY, Li XY, Wei W, Zhao SH, Cao JH. A genome-wide association study of five meat quality traits in Yorkshire pigs. Front Agric Sci Eng. (2014) 1:137–43. doi: 10.15302/J-FASE-2014014

CrossRef Full Text | Google Scholar

12. Shen MM, Liang Q, Dou TC, Meng M, Guo J, Lu J, et al. Genome-Wide association of spleen weight in layer chicken. Sci Agric Sin. (2018) 51:1213–22. doi: 10.3864/j.issn.0578-1752.2018.06.019

CrossRef Full Text | Google Scholar

13. Ghasemi M, Zamani P, Vatankhah M, Abdoli R. Genome-wide association study of birth weight in sheep. Animal. (2019) 13:1797–1803. doi: 10.1017/S1751731118003610

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Igoshin AV, Yudin NS, Belonogova NM, Larkin DM. Genome-wide association study for body weight in cattle populations from Siberia. Anim Genet. (2019) 50:250–3. doi: 10.1111/age.12786

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chang CC, Chow CC, Tellier LCM, Vattikuti S, Purcell SM, Lee JL. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. (2015) 4:7. doi: 10.1186/s13742-015-0047-8

PubMed Abstract | CrossRef Full Text

16. Yu JM, Pressoir G, Briggs WH, Bi IV, Yamasaki M, Doebley JF, et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. (2006) 38:203–8. doi: 10.1038/ng1702

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhou X, Stephens M. Genome-wide efficient mixed model analysis for association studies. Nat Genet. (2012) 44:821–4. doi: 10.1038/ng.2310

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lu ZK, Yue YJ, Yuan C, Liu JB, Chen ZQ, Sun XP, et al. Genome-Wide association study of body weight traits in Chinese fine-wool sheep. Animals. (2020) 10:170. doi: 10.3390/ani10010170

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sahana G, Guldbrandtsen B, Bendixen C, Lund MS. Genome-wide association mapping for female fertility traits in Danish and Swedish Holstein cattle. Animal Genetics. (2010) 41:579–88. doi: 10.1111/j.1365-2052.2010.02064.x

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kramarz B, Huntley RP, Rodríguez-López M, Roncaglia P, Saverimuttu SCC, Parkinson H, et al. Gene ontology curation of neuroinflammation biology improves the interpretation of Alzheimer's disease gene expression data. J Alzheimers Dis. (2020) 75:1417–35. doi: 10.3233/JAD-200207

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Corfitsen HT, Krantz B, Larsen A, Drago A. Molecular pathway analysis associates alterations in obesity-related genes and antipsychotic-induced weight gain. Acta Neuropsychiatr. (2020) 32:72–83. doi: 10.1017/neu.2019.41

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gholizadeh M, Rahimi-Mianji G, Nejati-Javaremi AB. Genomewide association study of body weight traits in baluchi sheep. J Genet. (2015) 94:143–6. doi: 10.1007/s12041-015-0469-1

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ma N, Li YJ, Song YQ, Lou YJ, Song XC. Estimates of genetics parameters of main economic tratits liaoning cachmere goat. J Jilin Agric Univer. (2005) 27:323–7. doi: 10.13327/j.jjlau.2005.04.022

CrossRef Full Text

24. Hossein M, Mohammad MS, Hossein MS. Genetic parameter estimates for growth traits and prolificacy in Raeini Cashmere goats. Trop Anim Health Product. (2012) 44:1213–20. doi: 10.1007/s11250-011-0059-z

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Dervishi E, Joy M, Alvarez-Rodriguez J, Serrano M, Czlvo JH. The forage type (grazing versus hay pasture) fed to ewes and the lamb sex affect fatty acid profile and lipogenic gene expression in the longissimus muscle of suckling lambs. J Anim Sci. (2012) 90:54–66. doi: 10.2527/jas.2011-4057

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sabitzer S, Metzler-Zebeli BU, Klein D, Zebeli Q, Hollmann M, Podstatzky-Lichtenstein L. Epithelial response to high-grain diets involves alteration in nutrient transporters and Na+/K+-ATPase mRNA expression in rumen and colon of goats. J Anim Sci. (2013) 91:4256–66. doi: 10.2527/jas.2012-5570

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bell CG, Benzinou M, Siddiq A, Lecoeur C, Dina C, Lemainque A, et al. Genome-wide linkage analysis for severe obesity in french caucasians finds significant susceptibility locus on chromosome 19q. Diabetes. (2004) 53:1857–65. doi: 10.2337/diabetes.53.7.1857

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Jackson RM, Griesel BA, Short KR, Sparling D, Freeman WM, Olson AL. Weight loss results in increased expression of anti-inflammatory protein CRISPLD2 in mouse adipose tissue. Obesity. (2019) 27:2025–36. doi: 10.1002/oby.22652

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Agulnick AD, Taira M, Breen JJ, Tanaka T, Dawid IB, Westphal H. Interactions of the LIM-domain-binding factor Ldbl with LIM homeodomain proteins. Nature. (1996) 384:270–2. doi: 10.1038/384270a0

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Bach I, Carriere C, Ostendorff HP, Andersen B, Rosenfeld MG. A family of LIM domain-associated cofactors confer transcriptional synergism between LIM and Otx homeodomain proteins. Genes Dev. (1997) 11:1370–80. doi: 10.1101/gad.11.11.1370

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Visvader JE, Mao X, Fujiwara Y, Hahm K, Orkin SH. The LIM-domain binding protein Ldb1 and its partner LMO2 act as negative regulators of erythroid differentiation. Proc Natl Acad Sci USA. (1997) 94:13707–12. doi: 10.1073/pnas.94.25.13707

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gu XR, Feng CG, Ma L, Song C, Wang YQ, Da Y, et al. Genome-Wide association study of body weight in chicken F2 resource population. PLoS ONE. (2011) 6:e21872. doi: 10.1371/journal.pone.0021872

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Liu R, Sun Y, Zhao G, Wang F, Wu D, Zheng M, et al. Genome-Wide association study identifies loci and candidate genes for body composition and meat quality traits in Beijing-you chickens. PLoS ONE. (2013) 8:e61172. doi: 10.1371/journal.pone.0061172

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Wang WH, Wang JY, Zhang T, Wang Y, Zhang Y, Han K. Genome-wide association study of growth traits in Jinghai Yellow chicken hens using SLAF-seq technology. Anim Genet. (2019) 50:175–6. doi: 10.1111/age.12346

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wei CJ, Hou D, Feng Y, Li T, Jing ZZ, Li WY, et al. Molecular characterization and a duplicated 31-bp indel within the LDB2 gene and its associations with production performance in chickens. Gene. (2020) 761:145046. doi: 10.1016/j.gene.2020.145046

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Houde AA, Ruchat SM, Allard C, Baillargeon JP, St-Pierre J, Perron P, et al. LRP1B, BRD2 and CACNA1D: new candidate genes in fetal metabolic programming of newborns exposed to maternal hyperglycemia. Epigenomics. (2015) 15:72. doi: 10.2217/epi.15.72

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wang XM, Liu XL, Deng DD, Yu M, Li XP. Genetic determinants of pig birth weight variability. BMC Genetics. (2016) 17:15. doi: 10.1186/s12863-015-0309-6

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Haas J, Beer AG, Widschwendter P, Oberdanner J, Salzmann K, Sarg B, et al. LRP1b shows restricted expression in human tissues and binds to several extracellular ligands, including fibrinogen and apoE-carrying lipoproteins. Atherosclerosis. (2011) 216:342–7. doi: 10.1016/j.atherosclerosis.2011.02.030

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Dietrich MF, van der Weyden L, Prosser HM, Bradley A, Herz J, Adams DJ. Ectodomains of the LDL receptor-related proteins LRP1b and LRP4 have anchorage independent functions in vivo. PLoS ONE. (2010) 5:e9960. doi: 10.1371/journal.pone.0009960

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Burgdorf KS, Gjesing AP, Grarup N, Justesen JM, Sandholt CH, Witte DR, et al. Association studies of novel obesity-related gene variants with quantitative metabolic phenotypes in a population-based sample of 6,039 Danish individuals. Diabetologia. (2012) 55:105–13. doi: 10.1007/s00125-011-2320-4

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Ehsan PA, Mathilde B, Emmanouil T, Monika S. Transcriptional study of appetite regulating genes in the brain of zebrafish (Danio rerio) with impaired leptin signalling. Sci Rep. (2019) 9:20166. doi: 10.1038/s41598-019-56779-z

PubMed Abstract | CrossRef Full Text | Google Scholar