Candidate circRNAs related to skeletal muscle development in Dazu black goats.

Circular RNA (CircRNA), as a classical noncoding RNA, has been proven to regulate skeletal muscle development (SMD). However, the molecular genetic basis of circRNA regulation in muscle cells remains unclear. In this study, the expression patterns of circRNAs in the longissimus dorsi muscle at embryonic day 75 and postnatal day 1 in DBGs were investigated to identify the key circRNAs that play an important role in SMD in goats. A total of 140 significantly and differentially expressed circRNAs (DEcircRNAs) were identified among the groups at different developmental stages. Among the 116 host genes (HGs) of DEcircRNAs, 76 were significantly and differentially expressed, which was confirmed by previous RNA_seq data. Furthermore, the expression pattern of 10 DEcircRNAs with RT-qPCR was verified, which showed 80% concordance rate with that of RNA_seq datasets. Moreover, the authenticity of seven randomly selected DEcircRNAs was verified by PCR Sanger sequencing. Based on the functional annotation results, among the 76 significantly and differentially expressed HGs, 74 were enriched in 845 GO terms, whereas 35 were annotated to 85 KEGG pathways. The results of this study could provide a comprehensive understanding of the genetic basis of circRNAs involved in SMD and muscle growth.

Comparative Transcriptomic Profiling in Ovarian Tissues of Lohmann Hens and Chengkou Mountain Chicken.

BACKGROUND: As a crucial economic characteristic and a major indicator of reproductive performance in layers, egg production is controlled by a series of complex regulatory heredity basis. In particular, the interacting regulatory function between noncoding RNAs (ncRNAs) and coding RNA plays important roles in regulating laying performance. METHODS: In this study, the RNA sequencing (RNA-seq) of ovarian tissues from Lohmann hens (n = 3) and Chengkou Mountain chicken (n = 3) under the laying peak period was performed to identify RNA transcriptional differences among different laying-performance populations. RESULTS: Results showed that the expression level of 303 mRNAs, 68 long ncRNAs (lncRNAs), 533 circular RNAs (circRNAs), and 79 microRNAs (miRNAs) was significantly different among the groups. Functional enrichment analysis of these differentially expressed (DE) mRNAs revealed that the laying process was implicated in numerous significantly enriched pathways (p < 0.05), such as the neuroactive ligand-receptor interaction, steroid hormone biosynthesis, and calcium-signaling pathway. Furthermore, the lncRNA/circRNA-miRNA-mRNA regulatory networks related to the regulation of laying performance were constructed. Some randomly selective DE RNAs were verified by Real Time Quantitative (RT-qRCR), indicating that the bioinformatics analysis results of RNA-seq data were credible. CONCLUSIONS: This study could increase our understanding of the heredity basis of transcriptome in the laying performance of chicken.

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data.

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

Genome-Wide Selection Signal Analysis to Investigate Wide Genomic Heredity Divergence between Eurasian Wild Boar and Domestic Pig.

As important livestock species, pigs provide essential meat resources for humans, so understanding the genetic evolution behind their domestic history could help with the genetic improvement of domestic pigs. This study aimed to investigate the evolution of convergence and divergence under selection in European and Asian domestic pigs by using public genome-wide data. A total of 164 and 108 candidate genes (CDGs) were obtained from the Asian group (wild boar vs. domestic pig) and the European group (wild boar vs. domestic pig), respectively, by taking the top 5% of intersected windows of a pairwise fixation index (FST) and a cross population extended haplotype homozygosity test (XPEHH). GO and KEGG annotated results indicated that most CDGs were related to reproduction and immunity in the Asian group. Conversely, rich CDGs were enriched in muscle development and digestion in the European group. Eight CDGs were subjected to parallel selection of Eurasian domestic pigs from local wild boars during domestication. These CDGs were mainly involved in olfactory transduction, metabolic pathways, and progesterone-mediated oocyte maturation. Moreover, 36 and 18 haplotypes of INPP5B and TRAK2 were identified in this study, respectively. In brief, this study did not only improve the understanding of the genetic evolution of domestication in pigs, but also provides valuable CDGs for future breeding and genetic improvement of pigs.

A 1.1 Mb duplication CNV on chromosome 17 contributes to skeletal muscle development in Boer goats.

The Boer goat is one of the top meat breeds in modern animal husbandry and has attracted widespread attention for its unique growth performance. However, the genetic basis of muscle development in the Boer goat remains obscure. In this study, we identified specific structural variants in the Boer goat based on genome-wide selection signals and analyzed the basis of the molecular heredity of related candidate genes in muscle development. A total of 9 959 autosomal copy number variations (CNVs) were identified through selection signal analysis in 127 goat genomes. Specifically, we confirmed that the highest signal CNV (HSV) was a chromosomal arrangement containing an approximately 1.11 Mb (CHIR17: 60062304-61171840 bp) duplicated fragment inserted in reverse orientation and a 5 362 bp deleted region (CHIR17:60145940-60151302 bp) with overlapping genes (e.g., ARHGAP10, NR3C2, EDNRA, PRMT9, and TMEM184C). The homozygous duplicated HSV genotype (+/+) was found in 96% of Boer goats but was not detected in Eurasian goats and was only detected in 4% of indigenous African goats. The expression network of three candidate genes ( ARHGAP10, NR3C2, and EDNRA) regulating dose transcription was constructed by RNA sequencing. Results indicated that these genes were involved in the proliferation and differentiation of skeletal muscle satellite cells (SMSCs) and their overexpression significantly increased the expression of SAA3. The HSV of the Boer goat contributed to superior skeletal muscle growth via the dose effects of overlapping genes.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene.

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The FST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

A deletion variant within the FGF5 gene in goats is associated with gene expression levels and cashmere growth.

The FGF5 gene has been associated with the regulation of fibre length in mammals, including cashmere goats. A deletion variant at ~14 kb downstream of the FGF5 gene showed significant divergence between cashmere and non-cashmere goats in previous studies. In this study, we designed specific primers to genotype the deletion variant. The results of gel electrophoresis and Sanger sequencing revealed that a 507-bp deletion mutation is located at 95 454 685-95 455 191 of chromosome 6 in goats. Genotyping data from a large panel of 288 goats showed that the deletion at the FGF5 gene locus appeared to be associated with cashmere length. The deletion variant was close to fixation (frequency 0.97) in cashmere goats. Furthermore, electrophoretic mobility shift assays for evaluating DNA-protein interaction and mRNA expression levels of FGF5 suggested that the deletion variant may serve as a cis-acting element by specifically binding transcription factors to mediate quantitative changes in FGF5 mRNA expression. Our study illustrates how a structural mutation of the FGF5 gene has contributed to the cashmere growth phenotype in domestic goats. The deletion mutation within the FGF5 gene could potentially serve as a molecular marker of cashmere growth in cashmere goat breeding.

Identification of Body Size Determination Related Candidate Genes in Domestic Pig Using Genome-Wide Selection Signal Analysis.

This study aimed to identify the genes related to the body size of pigs by conducting genome-wide selection analysis (GWSA). We performed a GWSA scan on 50 pigs belonging to four small-bodied pig populations (Diannan small-eared pig, Bama Xiang pig, Wuzhishan pig, and Jeju black pig from South Korea) and 124 large-bodied pigs. We used the genetic parameters of the pairwise fixation index (FST) and π ratio (case/control) to screen candidate genome regions and genes related to body size. The results revealed 47,339,509 high-quality SNPs obtained from 174 individuals, while 280 interacting candidate regions were obtained from the top 1% signal windows of both parameters, along with 187 genes (e.g., ADCK4, AMDHD2, ASPN, ASS1, and ATP6V0C). The results of the candidate gene (CG) annotation showed that a series of CGs (e.g., MSTN, LTBP4, PDPK1, PKMYT1, ASS1, and STAT6) was enriched into the gene ontology terms. Moreover, molecular pathways, such as the PI3K-Akt, HIF-1, and AMPK signaling pathways, were verified to be related to body development. Overall, we identified a series of key genes that may be closely related to the body size of pigs, further elucidating the heredity basis of body shape determination in pigs and providing a theoretical reference for molecular breeding.

Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data.

Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion−deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships.

Genetic basis investigation of wattle phenotype in goat using genome-wide sequence data.

In domestic goats, wattles often appear in even numbers, mostly on the neck and a few under the ear. Goat wattle is composed of ectopic cartilage tissue covered by skin and was reported as a dominant inheritance. Thirty-eight goats from two Southwest Chinese breeds were studied to elucidate the genetic basis of wattle phenotype in goat. Their genomes were sequenced for wide-genome selective sweep analysis (WGSA) and a genome-wide association study (GWAS). The WGSA results revealed 500 candidate genes identified by fixation index and π ratio and 261 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched with 195 genes and 38 significantly enriched KEGG items. In particular, three chondrogenesis-related pathways (Wnt, Hippo and MAPK signaling pathways) were found. Among the 500 genes, 474 were enriched to 2855 Gene Ontology items, and four (BMP2, BMP4, RARA and MSX1) were annotated in the regulation and development of chondrogenesis. Four chondrogenesis-related genes (GREM1, NEDD4, ATG7 and ITGA1) were identified from 519 single-nucleotide polymorphisms (SNPs) with a GWAS above the threshold. Six and 11 SNPs on chromosome 10 are located on GREM1 and NEDD4 respectively, and the highest numbers of SNPs on chromosomes 20 and 22 are located on ITGA1 and ATG7 respectively. All of these genes are related to cartilage development. This study identified a series of genes related to chondroplasia by GWAS and WGSA and presented the possibility that wattle inheritance may be influenced by multiple genes. This work provides a new theoretical understanding of the hereditary basis of wattle phenotype.

Genome-Wide Selective Analysis of Boer Goat to Investigate the Dynamic Heredity Evolution under Different Stages.

Boer goats, as kemp in meat-type goats, are selected and bred from African indigenous goats under a long period of artificial selection. Their advantages in multiple economic traits, particularly their plump growth, have attracted worldwide attention. The current study displayed the genome-wide selection signature analyses of South African indigenous goat (AF), African Boer (BH), and Australian Boer (AS) to investigate the hereditary basis of artificial selection in different stages. Four methods (principal component analysis, nucleotide diversity, linkage disequilibrium decay, and neighbor-joining tree) implied the genomic diversity changes with different artificial selection intensities in Boer goats. In addition, the θπ, FST, and XP-CLR methods were used to search for the candidate signatures of positive selection in Boer goats. Consequently, 339 (BH vs. AF) and 295 (AS vs. BH) candidate genes were obtained from SNP data. Especially, 10 genes (e.g., BMPR1B, DNER, ITGAL, and KIT) under selection in both groups were identified. Functional annotation analysis revealed that these genes are potentially responsible for reproduction, metabolism, growth, and development. This study used genome-wide sequencing data to identify inheritance by artificial selection. The results of the current study are valuable for future molecular-assisted breeding and genetic improvement of goats.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions.

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range.

Genome-Wide Selection Sweep between Wild and Local Pigs from Europe for the Investigation of the Hereditary Characteristics of Domestication in Sus Scrofa.

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

Investigating genetic diversity and population phylogeny of five Chongqing local chicken populations autosomal using microsatellites.

The genetic diversity and population structures of five Chongqing local chicken populations were investigated using by 24 microsatellite markers. Results revealed that the mean number of alleles (NA) ranged from 7.08 (Daninghe chicken, DN) to 8.46 (Nanchuan chicken, NC). The highest observed heterozygosity (HO) and expected heterozygosity (HE) were observed in DN (HO = 0.7252; HE = 0.7409) and the lowest HO and HE were observed in XS (Xiushan native chicken [XS], HO = 0.5910 and HE = 0.6697). The inbreeding coefficient (FIS) within population ranged from 0.022 (DN) to 0.119 (XS). Among the 24 microsatellite markers, four loci (MCW0111, MCW0016, ADL0278, and MCW0104) deviated from the Hardy-Weinberg equilibrium in all the studied populations. The results of population polygenetic analysis based on Nei's genetic distance and STRUCTURE software showed that the clustering of the five populations was incomplete consistent with geographical distribution. Moreover, a large number of gene flows were widespread among different populations, suggesting that genetic material exchanges occurred due to human activities and migration which was also verified by PCoA. In summary, this study preliminarily showed that Chongqing local chicken populations had rich genetic diversity and remarkable genetic divergence, but still high risk in conversion. These findings would be useful to the management of conservation strategies and the utilization of local chicken populations in further.

Genetic Basis of Follicle Development in Dazu Black Goat by Whole-Transcriptome Sequencing.

The follicle development (FD) is an important factor determining litter size in animals. Recent studies have found that noncoding RNAs (ncRNAs) play an important role in FD. In particular, the role of the regulatory mechanism of competing endogenous RNAs (ceRNAs) that drive FD has attracted increasing attention. Therefore, this study explored the genetic basis of goat FD by obtaining the complete follicular transcriptome of Dazu black goats at different developmental stages. Results revealed that 128 messenger RNAs (mRNAs), 4 long noncoding RNAs (lncRNAs), 49 microRNAs (miRNAs), and 290 circular RNAs (circRNAs) were significantly differentially expressed (DE) between large and small follicles. Moreover, DEmRNAs were enriched in many signaling pathways related to FD, as well as GO terms related to molecular binding and enzyme activity. Based on the analysis of the ceRNA network (CRN), 34 nodes (1 DElncRNAs, 10 DEcircRNAs, 14 DEmiRNAs, and 9 DEmRNAs) and 35 interactions (17 DEcircRNAs-DEmRNAs, 2 DElncRNAs-DEmiRNAs, and 16 DEmRNA-DEmiRNAs) implied that the CRN could be involved in the FD of goats. In conclusion, we described gene regulation by DERNAs and lncRNA/circRNA-miRNA-mRNA CRNs in the FD of goats. This study provided insights into the genetic basis of FD in precise transcriptional regulation.

Genome-wide selective detection of Mile red-bone goat using next-generation sequencing technology.

The ecotype population of goats (Capra hircus) was created by long-term artificial selection and natural adaptation. Mile red-bone goat is an indigenous breed with visible red bones, and its special bone structure has received extensive attention. This study aimed to identify genetic variants and candidate genes associated with specific bone phenotypes using next-generation sequencing technology (NGS). The results revealed that 31,828,206 single nucleotide polymorphisms (SNPs) were obtained from 72 goats (20 Mile red-bone goats and 52 common goats) by NGS. A total of 100 candidate genes were identified on the basis top 1% window interaction from nucleotide diversity (π), π ratio (π A/π B), and pairwise fixation index (F ST). Exactly 77 known signaling pathways were enriched. Specifically, three coding genes (NMNAT2, LOC102172983, and PNLIP) were annotated in the vitamin metabolism signaling pathways, and NCF2 was annotated to the osteoclast (OC) differentiation pathway. Furthermore, 5862 reliable copy number variations (CNVs) were obtained, and 14 and 24 genes were annotated with the top 1‰ CNV based on F ST (>0.490) and V ST (>0.527), respectively. Several pathways related to bone development and metabolism of exogenous substances in vivo, including calcium signaling pathway, OC differentiation, and glycerophospholipid metabolism, were annotated. Specifically, six genes from 19 candidate CNVs, which were obtained by interaction of the top 1‰ CNVs with F ST and V ST, were annotated to mucin-type O-glycan biosynthesis and metabolic pathways. Briefly, the results implied that pseudopurpurin and specific genetic variants work together to contribute to the red-bone color and specific bone structure of Mile red-bone goat. This study is helpful to understanding the genetic basis of the unique bone phenotype of Mile red-bone goats.

Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle.

Body weight (BW) is an important longitudinal trait that directly described the growth gain of bovine in production. However, previous genome-wide association study (GWAS) mainly focused on the single-record traits, with less attention paid to longitudinal traits. Compared with traditional GWAS models, the association studies based on the random regression model (GWAS-RRM) have better performance in the control of the false positive rate through considering time-stage effects. In this study, the BW trait data were collected from 808 Chinese Simmental beef cattle aged 0, 6, 12, and 18 months, then we performed a GWAS-RRM to fit the time-varied SNP effect. The results showed a total of 37 significant SNPs were associated with BW. Gene functional annotation and enrichment analysis indicated FGF4, ANGPT4, PLA2G4A, and ITGA5 were promising candidate genes for BW. Moreover, these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. These findings will provide prior molecular information for bovine gene-based selection, as well as facilitate the extensive application of GWAS-RRM in domestic animals.

Genome-wide selection signal analysis of Australian Boer goat reveals artificial selection imprinting on candidate genes related to muscle development.

As one of the best-known commercial goat breeds in the world, Boer goat has undergone long-term artificial selection for nearly 100 years, and its excellent growth rate and meat production performance have attracted considerable worldwide attention. Herein, we used single nucleotide polymorphisms (SNPs) called from the whole-genome sequencing data of 46 Australian Boer goats to detect polymorphisms and identify genomic regions related to muscle development in comparison with those of 81 non-specialized meat goat individuals from Europe, Africa, and Asia. A total of 13 795 202 SNPs were identified, and the whole-genome selective signal screen with a π ratio of nucleotide diversity (πcase /πcontrol ) and pairwise fixation index (FST ) was analyzed. Finally, we identified 1741 candidate selective windows based on the top 5% threshold of both parameters; here, 449 candidate genes were only found in 727 of these regions. A total of 433 genes out of the 449 genes obtained were annotated to 2729 gene ontology terms, of which 51 were directly linked to muscle development (e.g., muscle organ development, muscle cell differentiation) by 30 candidate genes (e.g., JAK2, KCNQ1, PDE5A, PDLIM5, TBX5). In addition, 246 signaling pathways were annotated by 178 genes, and two pathways related to muscle contraction, including vascular smooth muscle contraction (ADCY7, PRKCB, PLA2G4E, ROCK2) and cardiac muscle contraction (CACNA2D3, CASQ2, COX6B1), were identified. The results could improve the current understanding of the genetic effects of artificial selection on the muscle development of goat. More importantly, this study provides valuable candidate genes for future breeding of goats.

Three SNPs within exons of INHA and ACVR2B genes are significantly associated with litter size in Dazu black goats.

The aim of this study was to analyse the association between single-nucleotide polymorphisms within INHA and ACVR2B and litter size in Dazu black goats. In total, twenty-two SNPs were genotyped in 190 individuals by SNaPshot and resequencing. The results showed that three SNPs (SNP_1, SNP_12 and SNP_13 in this study) were detected to have significant additive genetic effect on the recorded goat litter size (p < .05). The SNP_1 (NC_030809.1), a non-synonymous substitution of G for T at chr2-g. 28314990 in the exon 2 of INHA gene (NM_001285606.1), resulted in homozygote 2 (HOM2) contributed 0.25 and heterozygote (HET) contributed 0.12 larger litter than homozygote 1 (HOM1). Meanwhile, SNP_12 (Chr22-g. 11721225 A > T) and SNP_13 (Chr22-g. 11721227 A > C) (NC_030829.1) simultaneously mutated at the first and third position of a triplet AAA (lysine, K) in the exon 4 of ACVR2B gene (XM_018066623.1) had estimated genetic effects of HOM1 (0.00) and HOM2 (0.03) larger than HET (-0.12). In conclusion, one SNPs (chr2-g. 28314990 T > G) within the exon 2 of INHA and two SNPs (Chr22-g. 11721225 A > T and Chr22-g. 11721227 A > C) in the exon 4 of ACVR2B gene were highly recommended as candidate markers of litter size in Dazu black goats. A large-scale association study to assess the impact of these variants on litter size is still necessary.

Genome-Wide Association Analysis of Growth Curve Parameters in Chinese Simmental Beef Cattle.

The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight-age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R2 = 0.954). The parameters' mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.