Genetic variant of the sheep E2F8 gene and its associations with litter size

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

A repertoire of single nucleotide polymorphisms (SNPs) of major fecundity BMPR1B gene among 75 sheep breeds worldwide

The BMPR1B gene is a major determinant of sheep reproductive capacity. Previous studies revealed that Q249R (FecB) is a profound variant of BMPR1B that influences the ovulation rate and litter size in sheep. However, unlike Q249R locus, the full spectrum of single nucleotide polymorphisms (SNPs) within BMPR1B has not been extensively studied. A systematic screen of SNPs in BMPR1B would facilitate the discovery of novel variants that are associated with litter size. This study aimed to investigate SNPs in the BMPR1B gene via whole genome sequence (WGS) data from 2409 individuals of 75 sheep breeds worldwide. Herein, a total of 9688 variants were screened, among which 15 were coding variants and 8 were novel changes. Specifically, we presented the most comprehensive frequency distribution map of the well-known FecB mutation to date. Besides, among the above-mentioned SNPs, one synonymous mutation (g.30050773C > T) was found to be likely under selection and is potentially associated with fecundity in Duolang sheep. Thus, our study greatly expands the variation repertoire of the ovine BMPR1B gene and provides a valuable resource for exploring causative mutations and genetic markers associated with litter size.

Orally Administered Lactobacillus rhamnosus CY12 Alleviates DSS-Induced Colitis in Mice by Restoring the Intestinal Barrier and Inhibiting the TLR4-MyD88-NF-κB Pathway via Intestinal Microbiota Modulation

Oral ingestion of probiotics is a promising approach to relieving inflammatory disease through regulating the gut microbiota. A newly discovered strain, Lactobacillus rhamnosus CY12 (LCY12), obtained from cattle-yak milk, displayed numerous probiotic properties. These included enhanced viability in low pH and bile environments, adhesion capabilities, and potent antimicrobial effects. The research aimed to explore the beneficial impacts of the novel LCY12 strain on colitis in mice induced by dextran sulfate sodium (DSS) and to elucidate the underlying molecular mechanisms. The results of the study showed that administration of LCY12 effectively helped to reduce the negative effects of DSS-induced body weight loss, disease activity index score, colon length shortening, loss of goblet cells, and overall histopathological scores in the intestines. Simultaneously, LCY12 administration significantly alleviated intestinal inflammation and safeguarded intestinal barrier integrity by enhancing IL-10 levels, while dampening IL-6, IL-1β, and TNF-α production. Additionally, LCY12 boosted the presence of tight junction proteins. Furthermore, LCY12 hindered the TLR4/MyD88/NF-κB signaling pathway by downregulating TLR4 and MyD88 expression, inactivating phosphorylated IκBα, and preventing translocation of NF-κB p65 from the cytoplasm to the nucleus. The LCY12 also increased specific intestinal microbial communities and short-chain fatty acid (SCFA) production. Altogether, LCY12 oral administration alleviated colitis induced with DSS in mice by improving intestinal barrier function and regulating inflammatory cytokines, SCFA production, and intestinal microbiota.

Are two beneficial mutations (p.Q249R and 90-bp Indel) within the ovine BMPRIB gene associated with growth traits?

Background

The problem of achieving economic efficiency in sheep breeding can be largely solved by increasing sheep productivity. Recently, the BMPRIB gene has been revealed by GWAS as a potential candidate gene for sheep body morphometric traits. Therefore, the present study aimed to investigate whether genetic polymorphisms (p.Q249R SNP and 90-bp deletion) in the BMPRIB gene are associated with sheep growth traits.

Methods

PCR-based genotyping was performed on 1,875 sheep, including 1,191 Guiqian semi-fine wool (GQSFW), 560 Luxi Blackhead (LXBH), 55 Lanzhou fat-tailed (LZFT), and 69 Weining (WN) sheep. Genotype-phenotype association was assessed using the independent samples t-test and ANOVA. The significance level was set at αoriginal < 0.05. The threshold p-value for significance was adjusted after correction for multiple comparisons using the Bonferroni correction.

Results

After the Bonferroni correction, it was found that individuals with FecB+/FecB+ genotypes of the p.Q249R had significantly better growth traits in LXBH ewe lambs, including the body length, chest width, paunch girth, cannon circumference, and hip width (P<0.0005). Meanwhile, associations were observed between 90-bp deletion polymorphism and several growth traits (body length, body height, chest depth, and canon circumference) in GQSFW ewe adults after the Bonferroni correction (P < 0.0002), and individuals with the "DD" genotypes had greater growth traits.

Conclusion

Our findings align with the experimental observations from GWAS, which identified the BMPRIB gene as a potential candidate gene for body measurement traits. These findings not only confirm the previous study's results but also expand on them. Therefore, further investigations regarding the impact of BMPRIB polymorphisms on growth traits are necessary in other sheep breeds.

Relationships between the mutations of the goat GATA binding protein 4 gene and growth traits

Osteogenesis is a complex multilevel process regulated by multiple genes. The GATA binding protein 4 (GATA4) gene has been extensively studied for its pivotal role in bone genesis and bone differentiation. However, its relationship with the growth traits of Shaanbei white cashmere (SBWC) and Guizhou black (GB) goats remains unclear. This work aims to investigate the potential influence of genetic mutations in the GATA4 gene on the growth traits goats. Thus, two Insertion/deletion (InDel) polymorphisms (8-bp-InDel and 9-bp-InDel) were screened and detected in a total of 1161 goats (including 980 SBWC goats and 181 GB goats) using PCR and agarose gel electrophoresis. The analyses revealed that there were two genotypes (ID and DD) for these two loci. In SBWC goats, 8-bp-InDel and 9-bp-InDel loci were significantly associated with heart girth (HG) and hip width (HW). Notably, individuals with DD genotype of 8-bp-InDel locus were superior while those with DD genotype of 9-bp-InDel locus were inferior. Correlation analyses of the four combined genotypes revealed significant associations with cannon circumference (CC), body height (BH), HG and HW. This work provides a foundation for the application of molecular marker-assisted selection (MAS) in goat breeding programs. Furthermore, the findings highlight the potential of the GATA4 gene and its genetic variations as valuable indicators for selecting goats with desirable growth traits.

Exploration of the Polymorphism Distribution of Bovine HMGA2 Gene in Worldwide Breeds and Its Associations with Ovarian Traits

The high-mobility group AT-hook 2(HMGA2) gene has been widely studied in the context of cancer and animal growth. However, recently, several studies have uncovered its critical role in cell proliferation. A genome-wide association study (GWAS) further suggests that the HMGA2 gene is a candidate gene in fertility, indicating its connection not only to growth traits but also to reproduction, specifically ovarian traits. Thus, this study aimed to analyze the distribution of the HMGA2 gene in 54 bovine breeds worldwide, identify important short fragment variants (indels), and investigate the relationship between HMGA2 and ovarian development. The dataset included genotypic information from a bovine population of 634 individuals (n = 634). After genotyping and analyzing four selected loci, we found that one out of four loci, rs133750033 (P4-D22-bp), was polymorphic. Our results also reveal that this indel of HMGA2 is significantly associated with certain ovarian traits (p < 0.05). Specifically, it has connection with ovarian length (p = 0.004) and ovarian height (p = 0.026) during diestrus. Additionally, we discovered a higher expression of the HMGA2 gene in Asian cattle breeds. In summary, this study suggests that HMGA2 has the potential to serve as an animal fertility testing marker gene. Moreover, these findings contribute to a more promising outlook for the bovine industry.

Transcriptome Analysis of Compensatory Growth and Meat Quality Alteration after Varied Restricted Feeding Conditions in Beef Cattle

Compensatory growth (CG) is a physiological response that accelerates growth following a period of nutrient limitation, with the potential to improve growth efficiency and meat quality in cattle. However, the underlying molecular mechanisms remain poorly understood. In this study, 60 Huaxi cattle were divided into one ad libitum feeding (ALF) group and two restricted feeding groups (75% restricted, RF75; 50% restricted, RF50) undergoing a short-term restriction period followed by evaluation of CG. Detailed comparisons of growth performance during the experimental period, as well as carcass and meat quality traits, were conducted, complemented by a comprehensive transcriptome analysis of the longissimus dorsi muscle using differential expression analysis, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and weighted correlation network analysis (WGCNA). The results showed that irrespective of the restriction degree, the restricted animals exhibited CG, achieving final body weights comparable to the ALF group. Compensating animals showed differences in meat quality traits, such as pH, cooking loss, and fat content, compared to the ALF group. Transcriptomic analysis revealed 57 genes and 31 pathways differentially regulated during CG, covering immune response, acid-lipid metabolism, and protein synthesis. Notably, complement-coagulation-fibrinolytic system synergy was identified as potentially responsible for meat quality optimization in RF75. This study provides novel and valuable genetic insights into the regulatory mechanisms of CG in beef cattle.

Immunization against inhibin DNA vaccine as an alternative therapeutic for improving follicle development and reproductive performance in beef cattle

The objective of the present study was to investigate the potential role of immunization against INH on follicular development, serum reproductive hormone (FSH, E2, and P4) concentrations, and reproductive performance in beef cattle. A total of 196 non-lactating female beef cattle (4-5 years old) with identical calving records (3 records) were immunized with 0.5, 1.0, 1.5, or 2.0 mg [(T1, n = 58), (T2, n = 46), (T3, n = 42) and (T4, n = 36), respectively] of the pcISI plasmid. The control (C) group (n = 14) was immunized with 1.0 mL 0.9% saline. At 21d after primary immunization, all beef cattle were boosted with half of the primary immunization dose. On day 10 after primary immunization, the beef cattle immunized with INH DNA vaccine evidently induced anti-INH antibody except for the T1 group. The T3 group had the greatest P/N value peak among all the groups. The anti-INH antibody positive rates in T2, T3 and T4 groups were significantly higher than that in C and T1 groups. RIA results indicated that serum FSH concentration in T2 group increased markedly on day 45 after booster immunization; the E2 amount in T3 group was significantly increased on day 10 after primary immunization, and the levels of E2 also improved in T2 and T3 groups after booster immunization; the P4 concentration in T2 group was significantly improved on day 21 after primary immunization. Ultrasonography results revealed that the follicles with different diameter sizes were increased, meanwhile, the diameter and growth speed of ovulatory follicle were significantly increased. Furthermore, the rates of estrous, ovulation, conception, and twinning rate were also significantly enhanced. These findings clearly illustrated that INH DNA vaccine was capable of promoting the follicle development, thereby improving the behavioral of estrous and ovulation, eventually leading to an augment in the conception rates and twinning rate of beef cattle.

Bovine FRAS1: mRNA Expression Profile, Genetic Variations, and Significant Correlations with Ovarian Morphological Traits, Mature Follicle, and Corpus Luteum

The amelioration of bovine fertility caused by a multi-factorial problem has always been a hot topic, among which the detection of available target genes is the most crucial. It was hypothesized that the Fraser extracellular matrix complex subunit 1 (FRAS1) gene detected by GWAS is involved in physiological activities such as ovarian development. Herein, unilateral ovaries from 2111 cows were used to examine the mRNA expression profile and polymorphisms of bovine FRAS1 and their associations with fertility-related characteristics. Firstly, it was confirmed that FRAS1 gene transcripts are expressed in various bovine tissues. Then, among five potential insertion-deletion (indel) loci, the 20 bp (named P3-D20-bp) and 15 bp (P4-D15-bp) deletion mutations were confirmed to be polymorphic with linkage equilibrium. Secondly, the P3-D20-bp polymorphism was significantly associated with ovarian weight and corpus luteum diameter in the metaestrus phase and ovarian length in the dioestrum stage. Additionally, both ovarian length and mature follicle diameter in metaestrus are significantly correlated with different genotypes of P4-D15-bp. Thirdly, the transcriptional expression of the FRAS1 gene in groups with a minimum value of ovarian weight or volume was significantly higher than the expression in groups with a maximum value. Instead of that, the more corpus luteum and mature follicles there are, the higher the transcription expression of the FRAS1 gene is. Furthermore, FRAS1 expression in cows with a heterozygous genotype (ID) of P3-D20-bp was significantly higher than others. Eventually, P3-D20-bp deletion could disturb the binding efficiency of WT1-I and Sox2 to FRAS1 sequence according to binding prediction, indicating that mutation may affect gene expression and traits by influencing the binding of transcription factors. Overall, the polymorphisms of P3-D20-bp and P4-D15-bp of the bovine FRAS1 gene significantly correlated to follicle or ovarian traits that could be applied in optimizing female fertility in cow MAS breeding programs.

Exploring the Sheep MAST4 Gene Variants and Their Associations with Litter Size

The economic efficiency of sheep breeding can be improved by enhancing sheep productivity. A recent genome-wide association study (GWAS) unveiled the potential impact of the MAST4 gene on prolificacy traits in Australian White sheep (AUW)). Herein, whole-genome sequencing (WGS) data from 26 different sheep breeds worldwide (n = 1507), including diverse meat, wool, milk, or dual-purpose sheep breed types from China, Europe, and Africa, were used. Moreover, polymerase chain reaction (PCR) genotyping of the MAST4 gene polymorphisms in (n = 566) Australian white sheep (AUW) was performed. The 3 identified polymorphisms were not homogeneously distributed across the 26 examined sheep breeds. Findings revealed prevalent polymorphisms (P3-ins-29 bp and P6-del-21 bp) with varying frequencies (0.02 to 0.97) across 26 breeds, while P5-del-24 bp was presented in 24 out of 26 breeds. Interestingly, the frequency of the P3-ins-29 bp variant was markedly higher in Chinese meat or dual-purpose sheep breeds, while the other two variants also showed moderate frequencies in meat breeds. Notably, association analysis indicated that all InDels were associated with AUW sheep litter size (p < 0.05). These results suggest that these InDels within the MAST4 gene could be useful in marker-assisted selection in sheep breeding.

circUBE3C modulates myoblast development by binding to miR-191 and upregulating the expression of p27

Noncoding RNAs, including miRNAs (microRNAs) and circRNAs (circular RNA), are crucial regulators of myoblast proliferation and differentiation during muscle development. However, the specific roles and molecular mechanisms of circRNAs in muscle development remain poorly understood. Based on the existing circRNA-miRNA-mRNA network, our study focuses on circUBE3C, exploring its differential expression in fetal and adult muscle tissue of the cattle and investigating its impact on myoblast proliferation, apoptosis, and differentiation. The functional analysis of overexpression plasmids and siRNAs (small interfering RNAs) targeting circUBE3C was comprehensively evaluated by employing an array of advanced assays, encompassing CCK-8 (cell counting kit-8), EdU (5-ethynyl-20-deoxyuridine), flow cytometry, western blot analysis, and RT-qPCR. In vivo investigations indicated that overexpression of circUBE3C impedes the process of skeletal muscle regeneration. Mechanistically, we demonstrated that circUBE3C interacts with miR-191 and alleviates the suppression of p27 through cytoplasmic separation, bioinformatics prediction, dual-luciferase reporter assay, and RIP (RNA immunoprecipitation). Our findings indicate that the novel circRNA circUBE3C competitively binds to miR-191, thereby inhibiting proliferation and promoting apoptosis in bovine primary myoblasts and unveiling a regulatory pathway in bovine skeletal muscle development. These findings expand our understanding of circRNA functions in mammals and provide a basis for further exploration of their role in myogenesis and muscle diseases.

The 12-bp indel in the SMAD family member 2 gene is associated with goat growth traits

SMAD family member 2 (SMAD2) is a member of the TGFβ signaling pathway and functions as an essential regulator in the processes of development, cell proliferation, and bone formation. A previous observation reported that a 12-bp deletion of this gene affected the litter size in goats. However, according to our knowledge, no study has reported an association between this polymorphism and goat body measurement traits. The purpose of this study was to investigate the association of the insertion/deletion (indel) within the SMAD2 gene with the growth traits of goats. The indel polymorphism was found to be significantly associated with chest width and bust (p < 0.05), while cannon circumference was significantly the strongest compared to other traits (p < 0.01) and individuals with the DD genotypes were more dominant genotypes than other genotypes. In summary, we found evidence that the 12-bp indel within the SMAD2 gene could improve goat body measurement traits, paving the way for marker-assisted selection in the field of goat genetics and breeding.

Detection of genetic variation in bovine CRY1 gene and its associations with carcass traits

The biological clock (also known as circadian clock) is closely related to growth and development, metabolism, and diseases in animals. As a part of the circadian clock, the cryptochrome circadian regulator 1 (CRY1) gene is involved in the regulation of biological processes such as osteogenesis, energy metabolism and cell proliferation, however, few studies have been reported on the relationship between this gene and animal carcass traits. Herein, a total of four insertion/deletion (InDel) loci within the CRY1 gene were detected in Shandong Black Cattle Genetic Resource (SDBCGR) population (n = 433). Among them, the P1-6-bp-del locus was polymorphic in population of interest. Moreover, the P1-6-bp-del locus showed two genotypes, with a higher insertion/insertion (II) genotype frequency (0.751) than insertion/deletion (ID) genotype frequency (0.249). Correlation analysis showed that the P1-6-bp-del locus polymorphisms were significantly associated with twenty carcass traits (e.g., slaughter weight, limb weight, and belly meat weight). Individuals with II genotype were significantly better than those with ID genotype for eighteen carcass traits. Therefore, the P1-6-bp-del locus of the CRY1 gene can be used as a molecular marker for beef cattle breeding.

Insertion/deletion (Indel) variant of the goat RORA gene is associated with growth traits

RAR related orphan receptor A (RORA), which encodes the retinoid-acid-related orphan receptor alpha (RORα), is a clock gene found in skeletal muscle. Several studies have shown that RORα plays an important role in bone formation, suggesting that RORA gene may take part in the regulation of growth and development. The purpose of this research is to study the insertion/deletion (indel) variations of the RORA gene and investigate the relationship with the growth traits of Shaanbei white cashmere (SBWC) goats. Herein, the current study identified that the P4-11-bp and P11-28-bp deletion sites are polymorphic among 12 pairs of primers within the RORA gene in the SBWC goats (n = 641). Moreover, the P11-28-bp deletion locus was significantly related to the body height (p = 0.046), height at hip cross (p = 0.012), and body length (p = 0.003). Both of P4-11-bp and P11-28-bp indels showed the moderate genetic diversity (0.25 Collapse

Key Words

• RAR related orphan receptor A (RORA)
• clock gene
• goat
• growth traits
• insertion/deletion

Collapse

MESH Headings

• Pregnancy
• Female
• Animals
• Litter Size/genetics
• Goats/physiology
• INDEL Mutation/genetics
• Phenotype

Collapse

Grants Collapse

Affiliation(s)

Qian Zhou College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Huina Hu College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuta Yang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuxin Kang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianyong Lan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianfeng Wu Institute of Animal Husbandry and Veterinary, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian, China
Zhengang Guo College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China Animal Husbandry and Veterinary Science Institute of Bijie city, Bijie, Guizhou, China
Chuanying Pan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China

Collapse

A 14-bp deletion in bovine EPAS1 gene is associated with carcass traits

EPAS1 (Endothelial PAS Domain Protein 1) gene is well-known for its function in plateau hypoxia adaptability. It encodes HIF-2α, which involved in the induction of genes regulated by oxygen and then affects multiple physiological processes such as angiogenesis and energy metabolism. All of these indicate it may affect the development of animals. In this study, a 14-bp deletion in EPAS1 gene was uncovered in Shandong black cattle population (n = 502). Two genotypes (II and ID) were found and the frequency of the homozygous II genotype is higher than the heterozygous ID genotype. This population is consisted with HWE (p > 0.05). And more importantly, the 14-bp deletion was associated with outside flat (p = 0.003), brisket (p = 0.001), and knuckle (p = 0.032). These findings suggested that the 14-bp deletion is significantly associated with carcass traits, which could be served as a molecular marker applied to cow breeding.

Investigation on mRNA expression and genetic variation within goat SMAD2 gene and its association with litter size

The SMAD family member 2 (SMAD2), a member of the TGF-beta superfamily, executes a significant part in the oogenesis and ovulation process. A genome-wide selective sweep analysis also found SMAD2 was different in the fertility groups of Laoshan dairy goats; whether this gene was linked to litter size was unknown. Therefore, SMAD2 was chosen to study its effects on Shaanbei white cashmere goat reproduction and mRNA expression profile. Herein, the mRNA expression level of SMAD2 was firstly determined in female goat tissues, revealing significant differences in mRNA levels of different tissues (p < 0.05), including ovary tissue, indicating a potential role for SMAD2 in goat prolificacy. Then, using six pairs of primers, only one indel locus (P3-Del-12-bp) was found to be polymorphic in goat SMAD2 (n = 501). ANOVA also revealed that a P3-Del-12-bp deletion was significantly related to first-born litter size (p = 0.037). The Chi-square (χ2) test revealed that the ID genotype was significantly more prevalent in mothers with multiple lambs (p = 0.01), indicating that heterozygous individuals (ID) are more likely to produce multiple lambs. Our findings suggest that the SMAD2 gene's P3-Del-12-bp deletion could be used to improve goat breeds by assisting with litter size selection.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

Copy number variation of bovine S100A7 as a positional candidate affected body measurements

Beef production is closely related to the national economy and the attention has been paid to the improvement of beef cattle by molecular markers associated. Copy number variations (CNVs) recently have been gained many researches and recognized as an important source of genetic variation. Extensive studies have indicated that CNVs have effects on a large range of economic traits by a wide range of gene copy number alteration. S100A7 is a member of S100 family which is a famous family of Ca2+-binding proteins. S100A7 plays a crucial role in many important phenotypes (progress) including inflammatory diseases, psoriasis, obesity, etc. The aim of our study was to explore the phenotypic effects of CNV located in the S100A7 gene of bovine chromosome 3. We detected S100A7 CNV by qPCR in different cattle breeds, including Qinchuan cattle, Yunling cattle, Xianan cattle and a crossbred group Pinan. The copy number was identified as gain, normal and loss type, our results showed that the gain type was the main type in three types of S100A7 CNV of the whole tested breeds. After CNV detection, association analysis between S100A7 CNV and growth traits was carried out in four cattle breeds. We found significant effects of the CNV on cattle growth traits with differently preferred CNV types such as gain type with better chest depth (p = 0.043) in QC, loss type with better body length (p = 0.008) and rump width (p = 0.014) in YL, normal with better chest girth (p = 0.001), gain with better waist width (p = 0.001) and rump width (p = 0.044) in PN. These results suggested that the S100A7 CNV could affect the phenotypic traits and be used as a promising genetic marker for cattle molecular breeding.

Investigation of the association between goat DNMT3B gene polymorphism and growth traits

The DNA methyltransferase 3 beta (DNMT3B) gene is key for DNA methylation and has been well recognized in regulating growth and development. A previous observation indicated that an 11-bp indel of DNMT3B affected the reproductive traits in goats, yet the effect of this polymorphism on body measurement traits in goats has not been reported. This study aims to investigate the associations between DNMT3B gene polymorphism and goat growth traits. We investigated this 11-bp indel in 2184 goats and three genotypes have been found in Shaanbei white cashmere goat (SBWC): insertion/insertion (II), deletion/deletion (DD) and insertion/deletion (ID). Only ID and DD genotypes were detected in Nubian goats and Guizhou heima goat (GZHM). The allele frequencies analyzed revealed that the 'D' allele frequencies were higher in all three goat breeds. Further association analysis demonstrated that this indel is markedly associated with the cannon circumference (CC) and cannon circumference index (CCI) of SBWC and cannon circumference (CC) of Nubian goats (p < .05). The CC and CCI are essential indicators to measure the growth status of goats. In summary, our study sheds some light on the potential impact of the 11-bp indel polymorphism of the DNMT3B gene on improving the growth traits in goats.

CircRNF111 Contributes to Adipocyte Differentiation by Elevating PPARγ Expression via miR-27a-3p

The content and distribution of adipocytes is an important factor that affects meat quality. Previous studies showed that circRNAs are involved in various physiological processes. Nevertheless, more research is needed to investigate the function of circRNAs in adipogenesis. The present study examines the effects of circRNF111 on adipogenesis of bovine preadipocyte and aims to elucidate the underlying molecular mechanisms. In our study, the sequence signature of circRNF111 was identified using bioinformatics, RNA-FISH, and sequencing. Mechanistically, knockdown or exogenous expression of circRNF111 in preadipocytes was done to prove the functional significance of circRNF111. Combined with bioinformatics, a dual fluorescein reporter system, and immunoprecipitation, the interaction between circRNF111, miR-27a-3p, and the target gene PPARγ was verified. The results reveal that circRNF111 is positively correlated with adipocyte differentiation. The newly identified bovine circRNF111 functions as a miR-27a-3p sponge to rescue the inhibitory effect of miR-27a-3p on the PPARγ gene, thereby promoting adipogenesis.

Association analysis of indels in the VISFATIN gene with five cattle breeds and their growth traits

VISFATIN is an adipose cytokine that has been proved to correlate with growth and development traits. In a previous study from our lab, two insertion/deletions (indels; including a 35-bp insertion at its intron 4 and a 6-bp deletion in intron 5) were identified within the VISFATIN gene. To validate these indels and evaluate their association with growth traits in Chinese cattle, a total of 413 samples from four Chinese indigenous breeds and 217 samples from Chinese breeds were detected. Three genotypes (WW, WI and II) at intron 4 were detected based on the 35-bp insertion (allele I) or deletion (allele W) and showed moderate polymorphism in all samples. Two genotypes (WW and WD) at intron 5 were detected based on the 6-bp deletion (allele D) or insertion (allele W) in Xianan (XN) cattle and Jinnan (JN) cattle population but showed poor polymorphisms. Association analysis illustrated that the indel at intron 4 is significantly associated with chest girth, rump length and body weight in Ji'an (JA) cattle and the indel at intron 5 can cause a significant difference in rump length in JN cattle. To our knowledge, it is the first time it has been shown that indels within the VISFATIN gene are associated with growth traits in the two Chinese indigenous cattle breeds. These findings suggest that the VISFATIN gene can be used as a molecular marker for JN and JA cattle breeding.

Insertion/deletions within the bovine FoxO1 gene and their association analysis with growth traits in three Chinese cattle breeds

FOXO1 (FKHR) gene, as a transcription factor, plays a vital role in animal growth and development, participating in many biological processes. The aim of this study was to ascertain Insertion/deletions (Indels) polymorphism within bovine FoxO1 gene in 679 Chinese adult cows and associate them with stature traits. Two Indels (named as Indel-3 and Indel-4, recorded as rs383545622 and rs525318770 in NCBI, respectively) were successfully genotyped by the Once PCR method, which was reliable, rapid and cost effective for simultaneous detection of two or more Indels. Indel-3 and Indel-4 were located at the second intron. All four different haplotypes (H1: D3D4, H2: I3D4, H3: D3I4, H4: I3I4) could be identified, and the D (del-) allele, DD (del-/del-) genotype and D3D4 haplotype retained the highest frequency. However, individuals with DI (D3I3, D4I4 or H1H4/H2H3 genotype) showed significantly better phenotypic traits than those with the other genotypes in Nanyang cattle, showing a hybrid vigor. The results implied that this DI genotype can be applied to early selective breeding to improve the productivity of Nanyang cattle. Our results suggested that these two Indels within the bovine FoxO1 gene might be used as genetic markers for marker-assisted selection (MAS) in cattle breeding and genetics.

Polymorphism within the GATA binding protein 4 gene is significantly associated with goat litter size

GATA binding protein 4 (GATA4) is a typical transcription binding factor, and its main functions include regulating the proliferation, differentiation and apoptosis of ovarian granulosa cells, promoting spermatogenesis and sex differentiation, implying that this gene have possibly roles in animal reproduction. This study aims to detect five potential insertion/deletions (indels) of the GATA4 gene in 606 healthy unrelated Shaanbei white cashmere (SBWC) goats and analyze its association with the litter size. The electrophoresis and DNA sequencing identified two polymorphic indels (e.g., P4-Del-8bp and P5-Ins-9bp indel). Then T-test analysis showed that P4-Del-8bp was significantly correlated with litter size (p = 0.022) because of two different genotypes detected, e.g., insertion-deletion (ID) and deletion-deletion (DD), and the average litter size of individuals with DD genotype goats was higher than that of others. However, there was no correlation between P5-Ins-9bp and lambing of goats. Chi-square (X2) test found that the distribution of and P4-Del-8bp genotypes (X2 = 6.475, p = 0.011) was significantly different between single and multiple-lamb groups, while P5-Ins-9bp (X2 = 0.030, p = 0.862) was not. Therefore, these findings revealed that P4-Del-8bp polymorphism of goat GATA4 gene was a potential molecular marker significantly associated with litter size, which can be used for the marker-assisted selection (MAS) breeding to improve goat industry.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

Distribution of DGAT1 copy number variation in Chinese goats and its associations with milk production traits

Goat is an important sector for meat and dairy products. Diacylglycerol O-acyltransferase 1 (DGAT1), which is a key gene in milk production, has been recently detected to overlap with a novel copy number variation (CNV) in goats. CNVs could be genetic markers providing new insights into the genetic basis of phenotypic variation. Up to now, there are no reports on the DGAT1-related CNV (DGAT1 CNV) in Chinese goats. This study first detected the distribution of the DGAT1 CNV in Chinese seven goat breeds, finding substantial differences among dairy, meat, and fiber goats (P < 0.01). The association analysis between the DGAT1 CNV and milk production traits revealed significant associations: Xinong Sannen (XS) dairy goat with copy number loss type had higher freezing point depression (FPD) (P < 0.01) and milk solids-not-fat (SNF) content (P < 0.05). Overall, our study unraveled the distribution of DGAT1 CNV in Chinese goats for the first time and found the potential role of this CNV in the marker-assisted selection of dairy goat breeding.

Phloretin alleviates palmitic acid-induced oxidative stress in HUVEC cells by suppressing the expression of LncBAG6-AS

Oxidative stress (OS) is an important trigger of vascular endothelial injury (VEI), which then leads to cardiovascular disease (CVDs). Phloretin was previously investigated to alleviate OS in human umbilical vein endothelial cells (HUVECs) by activating the AMPK/Nrf2 pathway; however, whether phloretin exerts cardiovascular health benefits by targeting non-coding RNAs (ncRNAs) remains unclear. Herein, the whole transcriptome sequencing and lncRNA library building were performed on HUVECs, a commonly used cell line for CVDs study, from different groups in control (CK), palmitic acid (PA, 100 μM), and PA + phloretin (50 μM, G50). KEGG analysis demonstrated that DE-lncRNAs regulated the pathway related to OS and metabolism in HUVECs. LncBAG6-AS was highly expressed under OS stimulation, which was reversed by phloretin co-treatment. Moreover, the MMP, activities of SOD, GSH-Px, T-AOC and GR were significantly ameliorated after interference of LncBAG6-AS, which were consistent with phloretin recover group. Furthermore, the expression of DE-genes from previously reported mRNA sequencing, including MAPK10, PIK3R1, ATP2B4, AKT2, and ADCY9, were significantly changed with LncBAG6-AS interference, indicating that LncBAG6-AS may participate in the process of OS attenuation by phloretin through regulating gene expression. So, the transcriptome sequencing of HUVECs with LncBAG6-AS knockdown was subsequently performed and DE-genes for "NC vs. si-ASO-LncBAG6-AS" were significantly enriched with GO terms, such as apoptosis, response to OS, ferroptosis, and others, which were similar to those observed from KEGG analysis. Overall, this study provides new insights into the molecular mechanisms by which bioactive substances alleviate OS and potential targets for the early prevention and treatment of VEI.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Structural variation and introgression from wild populations in East Asian cattle genomes confer adaptation to local environment

BACKGROUND

Structural variations (SVs) in individual genomes are major determinants of complex traits, including adaptability to environmental variables. The Mongolian and Hainan cattle breeds in East Asia are of taurine and indicine origins that have evolved to adapt to cold and hot environments, respectively. However, few studies have investigated SVs in East Asian cattle genomes and their roles in environmental adaptation, and little is known about adaptively introgressed SVs in East Asian cattle.

RESULTS

In this study, we examine the roles of SVs in the climate adaptation of these two cattle lineages by generating highly contiguous chromosome-scale genome assemblies. Comparison of the two assemblies along with 18 Mongolian and Hainan cattle genomes obtained by long-read sequencing data provides a catalog of 123,898 nonredundant SVs. Several SVs detected from long reads are in exons of genes associated with epidermal differentiation, skin barrier, and bovine tuberculosis resistance. Functional investigations show that a 108-bp exonic insertion in SPN may affect the uptake of Mycobacterium tuberculosis by macrophages, which might contribute to the low susceptibility of Hainan cattle to bovine tuberculosis. Genotyping of 373 whole genomes from 39 breeds identifies 2610 SVs that are differentiated along a "north-south" gradient in China and overlap with 862 related genes that are enriched in pathways related to environmental adaptation. We identify 1457 Chinese indicine-stratified SVs that possibly originate from banteng and are frequent in Chinese indicine cattle.

CONCLUSIONS

Our findings highlight the unique contribution of SVs in East Asian cattle to environmental adaptation and disease resistance.

Polymorphism of GTF2A1 gene is associated with litter size in sheep

General transcription factor IIA subunit 1 (GTF2A1) is required for transcriptional activation from most promoters in eukaryotic transcription. Previous whole-genome association analyses articles have predicted the effect of this gene on lambing in sheep. In the study, a total of nine insertion/deletion (indel) variants named L1 to L9 in the gene were selected to detect in 550 Australian White sheep (AuW) adult ewes. The polymorphisms were detected in four loci (L1, L2, L3 and L8locus) and the polymorphism information content (PIC) values were 0.270, 0.375, 0.372 and 0.314. Further, our study revealed that L1, L2 and L3 loci of GTF2A1 gene were significantly related to the first parity litter size and the polymorphism of L8 was significantly correlated with litter size in the second parity. In details, for the first parity, individuals with the II genotype of the L1 locus had higher little size than that with the ID genotype, individuals with the ID and DD genotype of the L2 locus had higher little size than that with the II genotype, and individuals with the DD genotype of the L3 locus had higher little size than that with the II genotype. All four loci do not follow to Hardy-Weinberg equilibrium, and have no linkage between them. In conclusion, the polymorphisms of GTF2A1 were confirmed and analysis results demonstrated that there are some relevance between difference genotypes and litter size, and these findings may provide new insights for accelerating sheep molecular breeding through molecular marker-assisted selection (MAS).

InDel and CNV within the AKAP13 Gene Revealing Strong Associations with Growth Traits in Goat

A-kinase-anchoring protein 13 (AKAP13) is a member of the AKAP protein family that has been found to be associated with bone formation. Thus, we investigated the AKAP13 gene as a potential candidate gene for molecular-marker-assisted selection (MAS) in breeding. Our aim was to explore genetic variations (InDel and CNV) within the AKAP13 gene of Shaanbei white cashmere (SBWC) goats and analyze their relationship with growth traits. Ultimately, we identified three InDel loci (16-bp deletion, 15-bp insertion, and 25-bp deletion) and three CNVs, and the 16-bp and 15-bp loci were significantly associated with goat body length (p < 0.05). Both the 16-bp deletion variant and the 15-bp insertion variant facilitated an increase in body length in goats. In addition to this, there was a certain superposition effect between 16-bp and 15-bp loci, although there was no linkage. Additionally, the CNV1 locus was significantly correlated with body height and body length of goats (p < 0.05), and CNV2 was significantly correlated with chest depth, chest circumference, and cannon circumference of goats (p < 0.05). Individuals with gain type showed excellent growth performance. In conclusion, the InDel and CNV loci that we have identified could possibly serve as effective molecular markers in goat breeding, which is very essential for improving efficiency and success of breeding. Moreover, our findings provide a new avenue for further research into the function of the AKAP13 gene.

Goat Pleomorphic Adenoma Gene 1 (PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits

The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.

Assessing the construction of a Healthy City in China: a conceptual framework and evaluation index system

OBJECTIVES

COVID-19 has brought challenges to the health of all mankind. It is particularly important to promote the construction of a 'Healthy China' and build a 'healthy community'. The aims of this study were to construct a reasonable conceptual framework for the Healthy City concept and to assess Healthy City construction in China.

STUDY DESIGN

This study combined qualitative and quantitative research.

METHODS

This study proposes the concept model of 'nature-human body-Healthy City' and accordingly constructs an evaluation index system for the construction of a Healthy City that integrates five dimensions, namely, the medical level, economic basis, cultural development, social services, and ecological environment to explore the spatial and temporal heterogeneity of Healthy City construction in China. Finally, the influencing factors of Healthy City construction patterns are explored using GeoDetector.

RESULTS

(1) The pace of Healthy City construction is generally on the rise; (2) the construction of Healthy Cities exhibits significant global spatial autocorrelation and gradually increasing agglomeration. The spatial distribution of cold hotspot areas was relatively stable; (3) medical and health progress is an important factor; the level of economic development is the leading support; the endowment of resources and environment is the basic condition; public service support provides important support; and scientific and technological innovation capabilities provide technical support for the construction of a Healthy City.

CONCLUSIONS

The spatial heterogeneity of Healthy City construction in China is evident, and the state of spatial distribution is relatively stable. The spatial pattern of Healthy City construction is shaped by a combination of factors. Our research will provide a scientific basis for promoting the construction of Healthy Cities and helping to implement the Health China Strategy.

New insight into copy number variations of goat SMAD2 gene and their associations with litter size and semen quality

Copy number variations (CNV) contribute significantly to genetic variations. Numerous studies have shown that CNV affects phenotypic traits in livestock. The SMAD family member 2 (SMAD2) is a leading candidate gene in reproduction and has a crucial effect on litter size. Additionally, SMAD2 is also required for male reproduction and influences male germ cell development. However, there are no reports on investigating the effect of CNVs in the SMAD2 gene on reproductive traits in goat. Therefore, the goal of this study was to explore associations between CNV of the SMAD2 gene and litter size and semen quality in Shaanbei white cashmere (SBWC) goats. In this study, two CNVs within the SMAD2 were identified in 352 SBWC goats (50 males and 302 females). The association analysis revealed that only CNV2 was significantly associated with female goat first-born litter size (P = 3.59 × 10^-4), male semen concentration (P < 0.01), ejaculation volume, live sperm count, and sperm deformity rate (P < 0.05). In terms of phenotypic performance, the individuals with loss genotypes outperformed those with other genotypes. CNV1 and CNV2 genotype combinations containing their dominant genotypes were also associated with goat litter size (P = 1.7 × 10^-5), but no differences in semen quality were found. In summary, CNV2 of the SMAD2 gene is useful for molecular marker-assisted selection breeding, as it is associated with essential goat reproductive traits.

Insertion/Deletion (InDel) Variants within the Sheep Fat-Deposition-Related PDGFD Gene Strongly Affect Morphological Traits

Platelet-derived growth factor D (PDGFD) is a member of the PDGF gene family, and it plays an important role in the regulation of adipocyte development in mammals. Furthermore, genome-wide association studies (GWAS) have previously identified it as a candidate gene associated with fleece fiber variation, body size, and the fat-tail phenotype in domestic Chinese sheep. In this study, a total of 1919 indigenous Chinese sheep were genotyped to examine the association between nucleotide sequence variations in PDGFD and body morphology. Our results detected both a 14 bp insertion in intron 2 and a 13 bp deletion in intron 4 of PDGFD. Moreover, these two InDel loci had low to moderate polymorphism. Notably, the 13 bp deletion mutation of PDGFD was found to significantly affect sheep body size. Yearling rams in the Luxi black-headed sheep (LXBH) containing a heterozygous genotype (insertion/deletion, ID) were found to have larger body length, chest depth, and body weight than those with wild genotypes. Furthermore, adult ewes in the Guiqian semi-fine wool sheep (GSFW) containing a homozygous mutation (deletion/deletion, DD) were found to have smaller chest width than their peers. Moreover, yearling ewes in this group with the same homozygous mutation were found to have lower body weight, chest width, and cannon circumference compared to those of other individuals. This study demonstrates that PDGFD InDel polymorphisms have the potential to be effective molecular markers to improve morphological traits in domestic Chinese sheep.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5′ untranslated region (5′ UTR) ofHOXB13is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

Convergent changes in MC1R gene are associated with black-headed coat color in sheep

As one of the most obvious phenotypic traits, the coat color of sheep is an ideal model to study the genetic mechanisms underlying coat color varieties of mammals. One distinguishable coat color is the black-headed type, such as the famous black-headed Dorper sheep from Africa and Bayinbuluke sheep from Asia. In this study, we compared the genome sequences of black-headed and all-white sheep to identify causative genes for the black-headed sheep, including black-headed Dorper versus white-headed Dorper, as well as Bayinbuluke (black-headed) versus Small-tailed Han (all-white). The most differentiating region between black-headed sheep and all-white sheep was found to harbor a haplotype covering melanocortin receptor 1 (MC1R) gene. The share of this haplotype by the black-headed sheep from Africa and Asia suggested that the convergent change in the MC1R region is likely to determine this unique coat color. Two missense mutations (g. 14251947T>A and g. 14252090G>A) within this haplotype of MC1R gene were found. We further analyzed whole genome sequence data of 460 world-wide sheep with diverse coat colors and confirmed the association between the MC1R haplotype with pigmentation variations. Our study provides novel insights into coat color genetics in sheep and expands our knowledge of the link between MC1R gene and varying pigmentation patterns in sheep.

A radiomics model development via the associations with genomics features in predicting axillary lymph node metastasis of breast cancer: a study based on a public database and single-centre verification

AIM

To evaluate the predictive performance of the radiomics model in predicting axillary lymph node (ALN) metastasis through the associations between radiomics features and genomic features in patients with breast cancer.

MATERIALS AND METHODS

Patients with breast cancer were enrolled retrospectively from a public database (111 patients as training group) and one hospital (15 patients as external validation group). The genomics features from transcriptome data and radiomics features from dynamic contrast-enhanced magnetic resonance imaging (MRI) were collected. Firstly, overlapping genes were identified using the Kyoto Encyclopedia of Genes and Genomes and differentially expressed gene analysis, while radiomics features were reduced using a data-driven method. Then, the associations between overlapping genes and retained radiomics features were assessed to obtain key pairs of radiomics-genomics features. Furthermore, the least absolute shrinkage and selection operator (LASSO) algorithm was used to detect the key-pairs features. Finally, radiomics and genomics models were constructed to predict ALN metastasis.

RESULTS

After using the hybrid data- and gene-driven selection method, key pairs of features were detected, which consisted of six radiomic features associated with four genomic features. The radiomics model exhibited comparable performance to the genomics model in predicting ALN metastasis (radiomic model: area under the curve [AUC] = 0.71, sensitivity = 77%, specificity = 56%; genomic model: AUC = 0.72, sensitivity = 85%, specificity = 74%). The four genomic features were enriched in six pathways and related to metabolism and human diseases.

CONCLUSION

The radiomics model established using the gene-driven hybrid selection method could predict ALN metastasis in breast cancer, which showed comparable performance to the genomics model.

Ferroptosis contributes to nickel-induced developmental neurotoxicity in zebrafish

Nickel (Ni) is a widely utilized heavy metal that can cause environmental pollution and health hazards. Its safety has attracted the attention of both the environmental ecology and public health fields. While the central nervous system (CNS) is one of the main targets of Ni, its neurotoxicity and the underlying mechanisms remain unclear. Here, by taking advantage of the zebrafish model for live imaging, genetic analysis and neurobehavioral studies, we reveal that the neurotoxic effects induced by exposure to environmentally relevant levels of Ni are closely related to ferroptosis, a newly-described form of iron-mediated cell death. In vivo two-photon imaging, neurobehavioral analysis and transcriptome sequencing consistently demonstrate that early neurodevelopment, neuroimmune function and vasculogenesis in zebrafish larvae are significantly affected by environmental Ni exposure. Importantly, exposure to various concentrations of Ni activates the ferroptosis pathway, as demonstrated by physiological/biochemical tests, as well as the expression of ferroptosis markers. Furthermore, pharmacological intervention of ferroptosis via deferoxamine (DFO), a classical iron chelating agent, strongly implicates iron dyshomeostasis and ferroptosis in these Ni-induced neurotoxic effects. Thus, this study elucidates the cellular and molecular mechanisms underlying Ni neurotoxicity, with implications for our understanding of the physiologically damaging effects of other environmental heavy metal pollutants.

Whole-genome resequencing reveals the genomic diversity and signatures of selection in Romanov sheep

Romanov sheep are adapted to the extremely cold and harsh environment and display a distinctive grey color. Herein, we analyzed the population structure, genetic diversity, and selection signatures of Romanov sheep based on whole-genome sequencing data of 17 Romanov sheep, 114 individuals from other 10 European breeds. The results of PCA, ADMIXTURE, and NJ-tree showed that the Romanov sheep was closely related to other northern European breeds. A relative high level of genetic diversity, low inbreeding coefficient, and large effective population size was observed in Romanov sheep when compared with other European breeds. We then screened the genomic selection signatures of Romanov sheep using FST, XP-XLP, and XP-EHH methods. The most significant region under selection (CHR14:14.2 to 14.3 Mb) harbored a haplotype that contained MC1R gene. Furthermore, this haplotype was also found in other grey-body breeds including Gotland sheep, Grey Tronder Sheep, and German grey heath sheep, suggesting that it was associated with the unique coat color of these breeds. We also found one region (CHR10:40.8Mb- 41.0Mb) harboring PCDH9 gene which was potentially associated with cold environmental adaptation. In summary, this study identified candidate genes that were associated with the unique grey color and environmental adaptation in Romanov sheep, which provided a basis for understanding the genetic background and utilization of this breed.

The ovine HIAT1 gene: mRNA expression, InDel mutations, and growth trait associations

Background

The hippocampal abundant transcript 1 (HIAT1) gene, also known as major facilitator superfamily domain-containing 14A (MFSD14A), encodes for a transmembrane transporter protein and has been previously shown to be associated with milk production in buffalo and sheep breeds, as well as growth traits in chicken and goats. However, tissue level distribution of the ovine HIAT1 gene, as well as its effect on body morphometric traits in sheep, has yet to be studied.

Methods

The HIAT1 mRNA expression profile of Lanzhou fat-tailed (LFT) sheep was determined by quantitative real-time PCR (qPCR). A total of 1498 sheep of three indigenous Chinese sheep breeds were PCR-genotyped for polymorphisms of HIAT1 gene. Student's t-test was used to observe the association between the genotype and sheep morphometric traits.

Results

HIAT1 was widely expressed in all examined tissues, and was particularly abundant in the testis of male LFT sheep. Additionally, a 9-bp insertion mutation (rs1089950828) located within the 5'-upstream region of HIAT1 was investigated in Luxi black-headed (LXBH) sheep and Guiqian semi-fine wool (GSFW) sheep. The wildtype allele frequency 'D' was found to be more prevalent than that of the mutant allele 'I'. Furthermore, low genetic diversity was confirmed in all sampled sheep populations. Subsequent association analyses indicated an association between the 9-bp InDel mutation of interest and the morphometric traits of LXBH and GSFW sheep. Furthermore, yearling ewes with a heterozygous genotype (ID) demonstrated smaller body sizes, while yearling rams and adult ewes with the heterozygous genotype were found to have overall better growth performance.

Conclusion

These findings imply that functional InDel polymorphism (rs1089950828) has the potential to be utilized for marker-assisted selection (MAS) of growth traits in domestic Chinese sheep populations.

Relationships between the Mini-InDel Variants within the Goat CFAP43 Gene and Body Traits

The cilia- and flagella-associated protein 43 (CFAP43) gene encodes a member of the cilia- and flagellum-associated protein family. Cilia on the cell surface influence intercellular signaling and are involved in biological processes such as osteogenesis and energy metabolism in animals. Previous studies have shown that insertion/deletion (InDel) variants in the CFAP43 gene affect litter size in Shaanbei white cashmere (SBWC) goats, and that litter size and body traits are correlated in this breed. Therefore, we hypothesized that there is a significant relationship between InDel variants within the CFAP43 gene and body traits in SBWC goats. Herein, we first investigated the association between three InDel variant loci (L-13, L-16, and L-19 loci) within CFAP43 and body traits in SBWC goats (n = 1827). Analyses revealed that the L-13, L-16, and L-19 loci were significantly associated with chest depth, four body traits, and three body traits, respectively. The results of this study are in good agreement with those previously reported and could provide useful molecular markers for the selection and breeding of goats for body traits.

A 67-bp variable duplication in the promoter region of the ADIPOQ is associated with milk traits in Xinjiang brown cattle

Adiponectin, also known as ADIPOQ, is a hormone protein secreted by adipocytes. The ADIPOQ gene is expressed primarily in adipose tissue, and the encoded protein circulates in the bloodstream and has the potential to regulate both animal fat metabolism and hormone production. Our previous work uncovered a 67-bp variable duplication in the promoter region of ADIPOQ, which reduced the basal transcriptional activity of ADIPOQ in the 3T3_L1 cell and also inhibits the ADIPOQ mRNA expression in adipose tissue. Accordingly, the present study aimed to identify the relationship between the 67-bp structural variations in ADIPOQ promoter region and the milk traits of Xinjiang brown cattle (XJBC). The results revealed two genotypes, DD and ID, in the XJBC, and minor allelic frequency (MAF) for the 'I' allele was more than 1%. Moreover, the association analysis revealed that the 67-bp duplication in the promoter region of the ADIPOQ gene was significantly correlated with the 305 days of milk production volume, fat yield, and milk fat percentage in the XJBC (p < 0.05). These results obtained in this study suggested that the identified variable duplication could be considered as the potential genetic marker for improving milk traits of XJBC.

Indel mutations of sheep PLAG1 gene and their associations with growth traits

Pleiomorphic adenoma gene 1 (PLAG1) is mainly expressed in embryonic development, and it is reported to take an effect on the growth performance of mice, cattle, pigs, and sheep. To explore how conservative the PLAG1 is in different sheep breeds, the effects of the two indel variants on the growth traits of the Chinese Luxi blackhead (LXBH) sheep were firstly detected. The P2-del 30 bp and P4-del 45 bp indel loci of the sheep PLAG1 gene were significantly related to 15 growth traits (P < 0.05). Genotype ID and genotype II were dominant for the P2-del 30 bp and P4-del 45 bp loci, respectively. The above findings indicated that the two indel mutations in the ovine PLAG1 gene were suggested to become the molecular markers for the selection of economic traits in sheep.

A 7-nt nucleotide sequence variant within the sheep KDM3B gene affects female reproduction traits

Lysine demethylase 3B (KDM3B) gene is a histone demethylase, demonstrating specific demethylation of the histone H3 lysine 9. It was detected as a sheep reproductive candidate gene by genome-wide scans, and related studies also showed its significance in female reproductive process. However, rare study researched its polymorphism. Herein, we hypothesized that the polymorphisms of KDM3B gene were associated with sheep reproduction traits. A 7-nt nucleotide sequence variant (rs1088697156) within KDM3B gene was identified in a total of 888 individuals, including the Australian White (AUW) sheep and Lanzhou Fat-tailed (LFT) sheep. II (insertion/insertion) and ID (insertion/deletion) genotypes of 7-nt variant were detected, which were at Hardy-Weinberg equilibrium (HWE) in detected breeds. Association analysis illustrated the 7-nt variant was significantly associated with the litter size, duration of pregnancy, live lamb number, live lamb rate, stillbirth number, stillbirth rate of average and different parity (P < 0.05) in AUW sheep. Moreover, 'ID' was the dominant genotype with excellent consistency in reproductive traits. It is instrumental to select individuals with ID genotype for improving the sheep reproduction traits. These findings suggest that the 7-nt variant within KDM3B gene can be used as a candidate marker of reproduction traits for sheep breeding improvement by marker-assisted selection.

CircRNA Profiling Reveals an Abundant circBDP1 that Regulates Bovine Fat Development by Sponging miR-181b/miR-204 Targeting Sirt1/TRARG1

The proliferation and differentiation of preadipocytes is an important factor determining bovine fat development, which is closely related to the feed conversion ratio, carcass traits, and beef quality. The purpose of this study was to identify the effects of candidate circRNA and miRNA on the proliferation and differentiation of bovine preadipocytes in order to provide basic materials for molecular breeding in cattle. circRNA sequencing was performed on bovine adipocyte samples at different differentiation time points, and a total of 1830 differentially expressed circRNAs were identified. Among them, circBDP1, derived from the bovine BDP1 gene, has potential binding sites for miR-204 (known as a regulator of bovine fat development) and miR-181b, which gives us a hint that circBDP1 may regulate bovine fat development by adsorbing miR-204 and miR-181b. Here, our results revealed that circBDP1 overexpression promoted the proliferation and differentiation of bovine preadipocytes. The miRNA profile of bovine adipocytes at different differentiation time points was also analyzed using the small RNA sequencing method, and a total of 89 differentially expressed miRNAs were identified, including miR-204 and miR-181b. As expected, dual-luciferase reporter results showed that circBDP1 competitively adsorbed miR-181b and miR-204. Overexpression and interference of miR-181b in bovine preadipocytes and 3T3-L1 showed that miR-181b promoted the proliferation and differentiation of preadipocytes. Further results displayed that miR-181b and miR-204 simultaneously targeted the SIRT1 gene, and miR-204 also targeted the 3' UTR region of the TRARG1 gene. In summary, this study found that miR-181b and miR-204 were involved in fat development by targeting SIRT1 and TRARG1. The results of this study will lay a foundation for the research of fat development and beef cattle industry.

Lactobacillus rhamnosus CY12 Enhances Intestinal Barrier Function by Regulating Tight Junction Protein Expression, Oxidative Stress, and Inflammation Response in Lipopolysaccharide-Induced Caco-2 Cells

The intestinal barrier is vital for preventing inflammatory bowel disease (IBD). The objectives of this study were to assess whether the Lactobacillus rhamnosus CY12 could alleviate oxidative stress, inflammation, and the disruption of tight junction (TJ) barrier functions induced by lipopolysaccharide (LPS), and therefore to explore the potential underlying molecular mechanisms. Our results showed that LPS-induced Cancer coli-2 (Caco-2) cells significantly increased the levels of reactive oxygen species (ROS), lactate dehydrogenase, inflammatory cytokines interleukin-1β, interleukin-6, interleukin-8, and tumor necrosis factor-α (IL-1β, IL-6, IL-8, and TNF-α), and the cell apoptosis rate while decreasing the levels of TJ proteins occludin, zonula occludens-1 (ZO-1), and claudin and antioxidant enzymes, such as catalase, superoxide dismutase, and glutathione peroxidase(CAT, SOD, and GSH-Px) (p < 0.05). However, Lactobacillus rhamnosus CY12 could relieve cytotoxicity, apoptosis, oxidative stress, and pro-inflammatory cytokine expressions, and also inhibit the Toll-like receptor 4/nuclear factor kappa-B(TLR4/NF-κB) signaling pathway. Furthermore, the gene expression of antioxidant enzymes, as well as the mRNA and protein expressions of TJ proteins, was improved. Particularly, the concentration of 108 cfu/mL significantly prevented the inflammatory injury induced by LPS in Caco-2 cells (p < 0.05). These findings support a potential application of Lactobacillus rhamnosus CY12 as a probiotic to prevent LPS-induced intestinal injury and treat intestinal barrier dysfunction.