Identification of Genomic Regions and Candidate Genes Associated with Body Weight and Body Conformation Traits in Karachai Goats

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Transcriptome analysis of adipose tissue and muscle of Laiwu and Duroc pigs

Fat content is an important trait in pig production. Adipose tissue and muscle are important sites for fat deposition and affect production efficiency and quality. To regulate the fat content in these tissues, we need to understand the mechanisms behind fat deposition. Laiwu pigs, a Chinese indigenous breed, have significantly higher fat content in both adipose tissue and muscle than commercial breeds such as Duroc. In this study, we analyzed the transcriptomes in adipose tissue and muscle of 21-d-old Laiwu and Duroc piglets. Results showed that there were 828 and 671 differentially expressed genes (DEG) in subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT), respectively. Functional enrichment analysis showed that these DEG were enriched in metabolic pathways, especially carbohydrate and lipid metabolism. Additionally, in the longissimus muscle (LM) and psoas muscle (PM), 312 and 335 DEG were identified, demonstrating enrichment in the cell cycle and metabolic pathways. The protein-protein interaction (PPI) networks of these DEG were analyzed and potential hub genes were identified, such as FBP1 and SCD in adipose tissues and RRM2 and GADL1 in muscles. Meanwhile, results showed that there were common DEG between adipose tissue and muscle, such as LDHB, THRSP, and DGAT2. These findings showed that there are significant differences in the transcriptomes of the adipose tissue and muscle between Laiwu and Duroc piglets (P < 0.05), especially in metabolic patterns. This insight serves to advance our comprehensive understanding of metabolic regulation in these tissues and provide targets for fat content regulation.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

Genome-Wide Association Study of Body Conformation Traits in Tashi Goats (Capra hircus)

Identifying genetic markers of economically valuable traits has practical benefits for the meat goat industry. To better understand the genomic variations influencing body conformation traits, a genome-wide association study was performed on Tashi goats, an indigenous Chinese goat breed. A total of 155 Tashi goats were phenotyped for eight body conformation traits: body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Then, 100 Tashi goats were randomly selected for whole-genome sequencing and genotyped. We obtained 1676.4 Gb of raw data with an average sequencing depth of 6.2X. Clean reads were aligned to the ARS1.2 reference genome, and 11,257,923 single nucleotide polymorphisms (SNPs) were identified. The structure analysis showed that these Tashi goats were almost not genetically related. The 109, 20, 52, 14, 62, 51, 70, and 7 SNPs were significantly associated with body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Within the ±500 kb region of significant SNPs, 183 genes were annotated. The most significantly enriched KEGG pathway was "olfactory transduction", and the most significantly enriched gene ontology (GO) terms were "cellular process", "cellular anatomical entity", and "molecular transducer activity". Interestingly, we found several SNPs on chromosomes 10 and 11 that have been identified multiple times for all eight body conformation traits located in two fragments (114 kb and 1.03 Mb). In chr.10:25988403-26102739, the six SNPs were tightly linked, the TACTAG genotype was the highest at 91.8%, and the FNTB (Farnesyltransferase, CAAX Box Beta) and CHURC1 (Churchill Domain Containing 1) genes were located. In chr.11:88216493-89250659, ten SNPs were identified with several dependent linkage disequilibrium (LD) blocks, and seven related genes were annotated, but no significant SNP was located in them. Our results provide valuable biological information for improving growth performance with practical applications for genomic selection in goats.

Genome-Wide Association Study Revealed Putative SNPs and Candidate Genes Associated with Growth and Meat Traits in Japanese Quail

The search for SNPs and candidate genes that determine the manifestation of major selected traits is one crucial objective for genomic selection aimed at increasing poultry production efficiency. Here, we report a genome-wide association study (GWAS) for traits characterizing meat performance in the domestic quail. A total of 146 males from an F2 reference population resulting from crossing a fast (Japanese) and a slow (Texas White) growing breed were examined. Using the genotyping-by-sequencing technique, genomic data were obtained for 115,743 SNPs (92,618 SNPs after quality control) that were employed in this GWAS. The results identified significant SNPs associated with the following traits at 8 weeks of age: body weight (nine SNPs), daily body weight gain (eight SNPs), dressed weight (33 SNPs), and weights of breast (18 SNPs), thigh (eight SNPs), and drumstick (three SNPs). Also, 12 SNPs and five candidate genes (GNAL, DNAJC6, LEPR, SPAG9, and SLC27A4) shared associations with three or more traits. These findings are consistent with the understanding of the genetic complexity of body weight-related traits in quail. The identified SNPs and genes can be used in effective quail breeding as molecular genetic markers for growth and meat characteristics for the purpose of genetic improvement.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Screening of InDel variants in PRDM6, myostatin and IGF2BP1 genes and association analysis with body measurement traits in Malabari and Attappadi black goats

An insertion/deletion (InDel) polymorphism study of PR/SET domain family 6 (PRDM6), myostatin (MSTN) and insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) genes was conducted in Malabari and Attappady black goats. An association study of identified InDels and body measurement traits was also performed. Body measurements included body length, chest diameter, chest depth, canon circumference, hip width, and hip height at the hip cross. The body trunk index, the body length index, the canon circumference index, and the chest width index were calculated. The Hardy-Weinberg equilibrium (HWE) was tested using a Chi-square test. The observed heterozygosity (Ho), expected heterozygosity (He), and polymorphism information content (PIC) were calculated. A significant difference in body measurements was found across breeds, ages, and breed x age interactions. The PRDM6 InDel was also associated with body measurement traits, such as body height, canon circumference and canon circumference index. In both Malabari and Attappadi black MSTN and PRDM6 InDels were in a state of HWE, while IGF2BP1 InDels were not. Indel markers found in the present study may be used for marker-assisted selection of growth traits among goats.

Functional Characterization of the Almstn2 Gene and Its Association with Growth Traits in the Yellowfin Seabream Acanthopagrus latus (Hottuyn, 1782)

Myostatin (mstn), also known as GDF8, is a growth and differentiation factor of the transforming growth factor-β (TGF-β) superfamily and plays a key inhibitory effect in the regulation of skeletal muscle development and growth in vertebrates. In the present study, to comprehend the role of the mstn2 gene of the yellowfin seabream Acanthopagrus latus (Almstn2b), the genomic sequence of Almstn2b is 2359 bp, which encodes 360 amino acids and is composed of three exons and two introns, was obtained. Two typical regions, a TGF-β propeptide and TGF-β domain, constitute Almstn2b. The topology indicated that Almstn2 was grouped together with other Perciformes, such as the gilthead seabream Sparus aurata. Moreover, Almstn2b was mainly expressed in the brain, fins, and spleen. Furthermore, five SNPs, one in the exons and four in the introns, were identified in the Almstn2b gene. The allele and genotype frequencies of SNP-Almstn2b +1885 A/G were significantly related to the total weight, interorbital distance, stem length, tail length, caudal length, caudal height, body length, and total length (p < 0.05). The allele and genotype frequencies of SNP-Almstn2b +1888 A/G were significantly related to the weight, interorbital distance, long head behind the eyes, body height, tail length, caudal length, and body length. Additionally, the relationship between the SNP-Almstn2b +1915 A/G locus and weight and long head behind the eyes was significant (p < 0.05). Furthermore, the other two SNPs were not significantly associated with any traits. Thus, the SNPs identified in this study could be utilized as candidate SNPs for breeding and marker-assisted selection in A. latus.

Genomic Selection for Live Weight in the 14th Month in Alpine Merino Sheep Combining GWAS Information

Alpine Merino Sheep is a novel breed reared from Australian Merino Sheep as the father and Gansu Alpine Fine-Wool Sheep as the mother, living all year in cold and arid alpine areas with exceptional wool quality and meat performance. Body weight is an important economic trait of the Alpine Merino Sheep, but there is limited research on identifying the genes associated with live weight in the 14th month for improving the accuracy of the genomic prediction of this trait. Therefore, this study's sample comprised 1310 Alpine Merino Sheep ewes, and the Fine Wool Sheep 50K Panel was used for genome-wide association study (GWAS) analysis to identify candidate genes. Moreover, the trial population (1310 ewes) in this study was randomly divided into two groups. One group was used as the population for GWAS analysis and screened for the most significant top 5%, top 10%, top 15%, and top 20% SNPs to obtain prior marker information. The other group was used to estimate the genetic parameters based on the weight assigned by heritability combined with different prior marker information. The aim of this study was to compare the accuracy of genomic breeding value estimation when combined with prior marker information from GWAS analysis with the optimal linear unbiased prediction method for genome selection (GBLUP) for the breeding value of target traits. Finally, the accuracy was evaluated using the five-fold cross-validation method. This research provides theoretical and technical support to improve the accuracy of sheep genome selection and better guide breeding. The results demonstrated that eight candidate genes were associated with GWAS analysis, and the gene function query and literature search results suggested that FAM184B, NCAPG, MACF1, ANKRD44, DCAF16, FUK, LCORL, and SYN3 were candidate genes affecting live weight in the 14th month (WT), which regulated the growth of muscle and bone in sheep. In genome selection analysis, the heritability of GBLUP to calculate the WT was 0.335-0.374, the accuracy after five-fold cross-verification was 0.154-0.190, and after assigning different weights to the top 5%, top 10%, top 15%, and top 20% of the GWAS results in accordance with previous information to construct the G matrix, the accuracy of the WT in the GBLUP model was improved by 2.59-7.79%.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.