A near complete genome assembly of the East Friesian sheep genome.

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Norcantharidin Enhances the Antitumor Effect of 5-Fluorouracil by Inducing Apoptosis of Cervical Cancer Cells: Network Pharmacology, Molecular Docking, and Experimental Validation.

The high recurrence rate of cervical cancer is a leading cause of cancer deaths in women. 5-Fluorouracil (5-FU) is an antitumor drug used to treat many types of cancer, but its diminishing effectiveness and side effects limit its use. Norcantharidin (NCTD), a demethylated derivative of cantharidin, exhibits various biological activities. Here, we investigated whether NCTD could potentiate 5-FU to induce cervical cancer cell death. To assess the cell viability and synergistic effects of the drugs, cell counting kit-8 and colony formation assays were performed using HR-HPV-positive cervical cancer cell lines. Annexin V-FITC/PI staining and TUNEL assays were performed to confirm the induction of apoptosis. The synergistic effect of NCTD on the antitumor activity of 5-FU was analyzed using network pharmacology, molecular docking, and molecular dynamics simulations. Apoptosis-related proteins were examined using immunoblotting. The combination of NCTD and 5-FU was synergistic in cervical cancer cell lines. Network pharmacological analysis identified 10 common targets of NCTD and 5-FU for cervical cancer treatment. Molecular docking showed the strong binding affinity of both compounds with CA12, CASP9, and PTGS1. Molecular dynamics simulations showed that the complex system of both drugs with caspase-9 could be in a stable state. NCTD enhanced 5-FU-mediated cytotoxicity by activating apoptosis-related proteins. NCTD acts synergistically with 5-FU to inhibit cervical cancer cell proliferation. NCTD enhances 5-FU-induced apoptosis in cervical cancer cell lines via the caspase-dependent pathway.

Progress on CRISPR/Cas9 system in the genetic improvement of livestock and poultry.

CRISPR/Cas9 gene editing technology, as a highly efficient genome editing method, has been extensively employed in the realm of animal husbandry for genetic improvement. With its remarkable efficiency and precision, this technology has revolutionized the field of animal husbandry. Currently, CRISPR/Cas9-based gene knockout, gene knock-in and gene modification techniques are widely employed to achieve precise enhancements in crucial production traits of livestock and poultry species. In this review, we summarize the operational principle and development history of CRISPR/Cas9 technology. Additionally, we highlight the research advancements utilizing this technology in muscle growth and development, fiber growth, milk quality composition, disease resistance breeding, and animal welfare within the livestock and poultry sectors. Our aim is to provide a more comprehensive understanding of the application of CRISPR/Cas9 technology in gene editing for livestock and poultry.

Multimodal Drug Target Binding Affinity Prediction using Graph Local Substructure.

Predicting the binding affinity of drug target is essential to reduce drug development costs and cycles. Recently, several deep learning-based methods have been proposed to utilize the structural or sequential information of drugs and targets to predict the drug-target binding affinity (DTA). However, methods that rely solely on sequence features do not consider hydrogen atom data, which may result in information loss. Graph-based methods may contain information that is not directly related to the prediction process. Additionally, the lack of structured division can limit the representation of characteristics. To address these issues, we propose a multimodal DTA prediction model using graph local substructures, called MLSDTA. This model comprehensively integrates the graph and sequence modal information from drugs and targets, achieving multimodal fusion through a cross-attention approach for multimodal features. Additionally, adaptive structure aware pooling is applied to generate graphs containing local substructural information. The model also utilizes the DropNode strategy to enhance the distinctions between different molecules. Experiments on two benchmark datasets have shown that MLSDTA outperforms current state-of-the-art models, demonstrating the feasibility of MLSDTA.

Massively Parallel CRISPR-Cas9 Knockout Screening in Sheep Granulosa Cells for FSH Response Genes.

Follicle-stimulating hormone (FSH) regulates ovarian follicle development through specific gene expression programs. Granulosa cells (GCs) are somatic cells surrounding the oocytes, secreting gonadotropins to regulate ovulation and promote follicular development. By analyzing the effects of different doses of FSH on the proliferation of GCs, we found that adding 10 ng/mL of FSH, as the optimal concentration, could promote the growth of GCs. Furthermore, we have successfully constructed the first CRISPR-Cas9 knockout library targeting the genes on chromosomes 2 and 3 and the X chromosomes of the sheep massively parallel coding gene, as well as an ovarian GCs knockout cell library. For the first time, we have exposed the knockout cell library to a concentration of 10 ng/mL FSH to explore the underlying mechanisms. Through this screening, we have identified 836 positive-negative screening genes that are responsive to FSH, thereby revealing the regulatory mechanisms and screening the functionality of candidate genes. Next, RNA-Seq of control (0 ng/mL), low (10 ng/mL), and high (100 ng/mL) doses of FSH revealed 1708 differentially expressed genes, and combined with 836 genes, we obtained 129 FSH dose-dependent genes with extremely significant differences. This enables us to delve deeper into investigating and identifying the mechanisms by which FSH regulates GCs. More generally, we have discovered new regulatory factors and identified reproductivity-associated major effectors. These findings provide novel research directions for further studies on sheep reproduction.

Metagenomic insights into the response of soil microbial communities to pathogenic Ralstonia solanacearum.

Understanding the response of soil microbial communities to pathogenic Ralstonia solanacearum is crucial for preventing bacterial wilt outbreaks. In this study, we investigated the soil physicochemical and microbial community to assess their impact on the pathogenic R.solanacearum through metagenomics. Our results revealed that certain archaeal taxa were the main contributors influencing the health of plants. Additionally, the presence of the pathogen showed a strong negative correlation with soil phosphorus levels, while soil phosphorus was significantly correlated with bacterial and archaeal communities. We found that the network of microbial interactions in healthy plant rhizosphere soils was more complex compared to diseased soils. The diseased soil network had more linkages, particularly related to the pathogen occurrence. Within the network, the family Comamonadaceae, specifically Ramlibacter_tataouinensis, was enriched in healthy samples and showed a significantly negative correlation with the pathogen. In terms of archaea, Halorubrum, Halorussus_halophilus (family: Halobacteriaceae), and Natronomonas_pharaonis (family: Haloarculaceae) were enriched in healthy plant rhizosphere soils and showed negative correlations with R.solanacearum. These findings suggested that the presence of these archaea may potentially reduce the occurrence of bacterial wilt disease. On the other hand, Halostagnicola_larseniia and Haloterrigena_sp._BND6 (family: Natrialbaceae) had higher relative abundance in diseased plants and exhibited significantly positive correlations with R.solanacearum, indicating their potential contribution to the pathogen's occurrence. Moreover, we explored the possibility of functional gene sharing among the correlating bacterial pairs within the Molecular Ecological Network. Our analysis revealed 468 entries of horizontal gene transfer (HGT) events, emphasizing the significance of HGT in shaping the adaptive traits of plant-associated bacteria, particularly in relation to host colonization and pathogenicity. Overall, this work revealed key factors, patterns and response mechanisms underlying the rhizosphere soil microbial populations. The findings offer valuable guidance for effectively controlling soil-borne bacterial diseases and developing sustainable agriculture practices.

Effects of resveratrol on DLD and NDUFB9 decrease in frozen semen of Mongolian sheep.

Mongolian sheep are a breed of sheep in China known for their excellent cold and drought resistance. Sperm from Mongolian sheep are often cryopreserved to improve breeding outcomes. However, cryopreservation of sperm often results in issues such as reduced vitality and altered morphology. Therefore, the objective of this study was to investigate the impact of the cryoprotectant resveratrol on frozen sperm from Mongolian sheep, specifically examining its effects on key proteins during cryopreservation. In this study, sperm samples were obtained from three adult Mongolian rams and processed through semen centrifugation. The sperm motility parameters of Fresh Sperm Group (FR), Resveratrol added before freezing group (FF-Res), Resveratrol-free frozen sperm group (FT), and Resveratrol added after freeze-thawing group (FA-Res) were determined. The tandem mass tags (TMT) peptide labeling combined with LC-MS/MS was used for proteomic analysis of the total proteins in FR and FT groups. A total of 2651 proteins were identified, among which 41 proteins were upregulated and 48 proteins were downregulated after freezing. In-depth bioinformatics analysis of differentially abundant proteins (DAPs) revealed their close association with the tricarboxylic acid cycle (TCA) and oxidative phosphorylation pathway. The energy-related protein dihydrolipoamide dehydrogenase (DLD) and the reactive oxygen species (ROS)-related protein NADH dehydrogenase 1 beta subcomplex subunit 9 (NDUFB9) exhibited significant decreases, indicating their potential role as key proteins contributing to reduced sperm vitality. The study demonstrated that the addition of resveratrol (RES) to semen could elevate the expression levels of DLD and NDUFB9 proteins. This study represents the pioneering proteomic analysis of Mongolian ram sperm before and after cryopreservation, establishing the significance of DLD and NDUFB9 as key proteins influencing the decline in vitality following cryopreservation of Mongolian ram sperm. These findings clarify that resveratrol can enhance the levels of DLD and NDUFB9 proteins in cryopreserved Mongolian ram sperm, consequently enhancing their vitality.

Bottom Electrode Effects on Piezoelectricity of Pb(Zr0.52,Ti0.48)O3 Thin Film in Flexible Sensor Applications.

Piezoelectric thin films grown on a mechanical, flexible mica substrate have gained significant attention for their ability to convert mechanical deformation into electrical energy though a curved surface. To extract the generated charge from the PZT thin films, bottom electrodes are typically grown on mica substrates. However, this bottom electrode also serves as a buffering layer for the growth of PZT films, and its impact on the piezoelectric properties of PZT thin films remains understudied. In this work, the effect of Pt and LaNiO3 bottom electrodes on the piezoelectric effect of a Pb(Zr0.52,Ti0.48)O3 thin film was investigated. It was observed that the PZT thin films on LNO/Mica substrate possessed weaker stress, stronger (100) preferred orientation, and higher remanent polarization, which is beneficial for a higher piezoelectric response theoretically. However, due to insufficient grain growth resulting in more inactive grain boundaries and lattice imperfections, the piezoelectric coefficient of the PZT thin film on LNO/Mica was smaller than that of the PZT thin film on a Pt/Mica substrate. Therefore, it is concluded that, under the current experimental conditions, PZT films grown with Pt as the bottom electrode are better suited for applications in flexible piezoelectric sensor devices. However, when using LNO as the bottom electrode, it is possible to optimize the grain size of PZT films by adjusting the sample preparation process to achieve piezoelectric performance exceeding that of the PZT/Pt/Mica samples.

LILRB3 Modulates Acute Myeloid Leukemia Progression and Acts as an Effective Target for CAR T-cell Therapy.

Identifying novel cell surface receptors that regulate leukemia cell differentiation and can be targeted to inhibit cellular proliferation is crucial to improve current treatment modalities in acute myeloid leukemia (AML), especially for relapsed or chemotherapy-refractory leukemia. Leukocyte immunoglobulin-like receptor type B (LILRB) is an immunomodulatory receptor originally found to be expressed in myeloid cells. In this study, we found that LILRB receptors can be induced under inflammatory stimuli and chemotherapy treatment conditions. Blockade of LILRB3 inhibited leukemia cell proliferation and leukemia progression. In addition, treatment with LILRB3 blocking antibodies upregulated myeloid lineage differentiation transcription factors, including PU.1, C/EBP family, and IRF, whereas phosphorylation of proliferation regulators, for example, AKT, cyclin D1, and retinoblastoma protein, was decreased. Conversely, transcriptomic analysis showed LILRB3 activation by agonist antibodies may enhance leukemia survival through upregulation of cholesterol metabolism, which has been shown to promote leukemia cell survival. Moreover, LILRB3-targeted CAR T cells exhibited potent antitumor effects both in vitro and in vivo. Taken together, our results suggest that LILRB3 is a potentially potent target for multiple treatment modalities in AML. SIGNIFICANCE: LILRB3 regulates differentiation and proliferation in acute myeloid leukemia and can be targeted with monoclonal antibodies and CAR T cells to suppress leukemia growth.

Research on the Strong Generalization of Coal Gangue Recognition Technology Based on the Image and Convolutional Neural Network under Complex Conditions.

A coal gangue image recognition method based on complex conditions is proposed to address the current issue of image-based coal gangue recognition being greatly affected by complex conditions. First, complex conditions such as different shooting backgrounds, shooting distances, and lighting intensities are set to simulate the underground coal mining environment. Then, based on three convolutional neural network algorithms, the coal gangue recognition model is established, and the influence of different complex conditions on coal gangue image recognition is analyzed. At the same time, a network model with a strong generalization ability is determined. The results show that the accuracy of coal gangue image recognition has no obvious regularity under different shooting background conditions, and complex environments should be the primary factor affecting the accuracy of coal gangue image recognition. The accuracy of coal gangue image recognition is negatively correlated with the increase in shooting distance, and strong light conditions are conducive to improving the accuracy of coal gangue image recognition. The LeNet network model has the strongest generalization ability, which can meet the requirements of recognition accuracy and respond quickly. The accuracy of coal gangue image recognition under different complex conditions can reach more than 0.99, and the average single image recognition time is only 177 ms. This article studies the influence law of different complex conditions on the recognition of coal and gangue images and confirms that the LeNet network has strong generalization ability, achieving accurate and fast recognition of coal gangue images under complex conditions and providing a reference basis for the deployment of underground coal gangue sorting.

Multimodal contrastive representation learning for drug-target binding affinity prediction.

In the biomedical field, the efficacy of most drugs is demonstrated by their interactions with targets, meanwhile, accurate prediction of the strength of drug-target binding is extremely important for drug development efforts. Traditional bioassay-based drug-target binding affinity (DTA) prediction methods cannot meet the needs of drug R&D in the era of big data. Recent years we have witnessed significant success on deep learning-based models for drug-target binding affinity prediction task. However, these models only considered a single modality of drug and target information, and some valuable information was not fully utilized. In fact, the information of different modalities of drug and target can complement each other, and more valuable information can be obtained by fusing the information of different modalities. In this paper, we introduce a multimodal information fusion model for DTA prediction that is called FMDTA, which fully considers drug/target information in both string and graph modalities and balances the feature representations of different modalities by a contrastive learning approach. In addition, we exploited the alignment information of drug atoms and target residues to capture the positional information of string patterns, which can extract more useful feature information in SMILES and target sequences. Experimental results on two benchmark datasets show that FMDTA outperforms the state-of-the-art model, demonstrating the feasibility and excellent feature capture capability of FMDTA. The code of FMDTA and the data are available at: https://github.com/bestdoubleLin/FMDTA.

Integrated Multi-Tissue Transcriptome Profiling Characterizes the Genetic Basis and Biomarkers Affecting Reproduction in Sheep (Ovis aries).

The heritability of litter size in sheep is low and controlled by multiple genes, but the research on its related genes is not sufficient. Here, to explore the expression pattern of multi-tissue genes in Chinese native sheep, we selected 10 tissues of the three adult ewes with the highest estimated breeding value in the early study of the prolific Xinggao sheep population. The global gene expression analysis showed that the ovary, uterus, and hypothalamus expressed the most genes. Using the Uniform Manifold Approximation and Projection (UMAP) cluster analysis, these samples were clustered into eight clusters. The functional enrichment analysis showed that the genes expressed in the spleen, uterus, and ovary were significantly enriched in the Ataxia Telangiectasia Mutated Protein (ATM) signaling pathway, and most genes in the liver, spleen, and ovary were enriched in the immune response pathway. Moreover, we focus on the expression genes of the hypothalamic-pituitary-ovarian axis (HPO) and found that 11,016 genes were co-expressed in the three tissues, and different tissues have different functions, but the oxytocin signaling pathway was widely enriched. To further explore the differences in the expression genes (DEGs) of HPO in different sheep breeds, we downloaded the transcriptome data in the public data, and the analysis of DEGs (Xinggao sheep vs. Sunite sheep in Hypothalamus, Xinggao sheep vs. Sunite sheep in Pituitary, and Xinggao sheep vs. Suffolk sheep in Ovary) revealed the neuroactive ligand-receptor interactions. In addition, the gene subsets of the transcription factors (TFs) of DEGs were identified. The results suggest that 51 TF genes and the homeobox TF may play an important role in transcriptional variation across the HPO. Altogether, our study provided the first fundamental resource to investigate the physiological functions and regulation mechanisms in sheep. This important data contributes to improving our understanding of the reproductive biology of sheep and isolating effecting molecular markers that can be used for genetic selection in sheep.

LSM14B controls oocyte mRNA storage and stability to ensure female fertility.

Controlled mRNA storage and stability is essential for oocyte meiosis and early embryonic development. However, how to regulate mRNA storage and stability in mammalian oogenesis remains elusive. Here we showed that LSM14B, a component of membraneless compartments including P-body-like granules and mitochondria-associated ribonucleoprotein domain (MARDO) in germ cell, is indispensable for female fertility. To reveal loss of LSM14B disrupted primordial follicle assembly and caused mRNA reduction in non-growing oocytes, which was concomitant with the impaired assembly of P-body-like granules. 10× Genomics single-cell RNA-sequencing and immunostaining were performed. Meanwhile, we conducted RNA-seq analysis of GV-stage oocytes and found that Lsm14b deficiency not only impaired the maternal mRNA accumulation but also disrupted the translation in fully grown oocytes, which was closely associated with dissolution of MARDO components. Moreover, Lsm14b-deficient oocytes reassembled a pronucleus containing decondensed chromatin after extrusion of the first polar body, through compromising the activation of maturation promoting factor, while the defects were restored via WEE1/2 inhibitor. Together, our findings reveal that Lsm14b plays a pivotal role in mammalian oogenesis by specifically controlling of oocyte mRNA storage and stability.

Estimates of Genomic Heritability and the Marker-Derived Gene for Re(Production) Traits in Xinggao Sheep.

Xinggao sheep are a breed of Chinese domestic sheep that are adapted to the extremely cold climatic features of the Hinggan League in China. The economically vital reproductive trait of ewes (litter size, LS) and productive traits of lambs (birth weight, BWT; weaning weight, WWT; and average daily gain, ADG) are expressed in females and later in life after most of the selection decisions have been made. This study estimated the genetic parameters for four traits to explore the genetic mechanisms underlying the variation, and we performed genome-wide association study (GWAS) tests on a small sample size to identify novel marker trait associations (MTAs) associated with prolificacy and growth. We detected two suggestive significant single-nucleotide polymorphisms (SNPs) associated with LS and eight significant SNPs for BWT, WWT, and ADG. These candidate loci and genes also provide valuable information for further fine-mapping of QTLs and improvement of reproductive and productive traits in sheep.

Clearance of bacteria from lymph nodes in sheep immunized with Brucella suis S2 vaccine is associated with M1 macrophage activation.

Ovine brucellosis is a global zoonotic disease of sheep caused by Brucella melitensis, which inflicts a significant burden on human and animal health. Brucella suis strain S2 (B. suis S2) is a smooth live attenuated vaccine for the prevention of ovine brucellosis in China. However, no previous studies have assessed the immunogenicity of B. suis S2 vaccine after oral immunization in sheep. Here, we attempted to evaluate the ovine immune response over the course of B. suis S2 immunization and to identify in vivo predictors for vaccine development. Body temperature, serum Brucella antibodies, serum cytokines (IL-12p70 and interferon [IFN]-γ), and bacterial load in the mandibular lymph nodes (LN), superficial cervical LN, superficial inguinal LN, and spleen were investigated to determine the safety and efficacy of the vaccine. The abnormal body temperature of sheep occurred within 8 days post-infection (dpi). Brucella suis S2 persisted for a short time (< 21 dpi) in the mandibular LN. The highest level of IL-12p70 was observed at 9 dpi, whereas serum IFN-γ levels peaked at 12 dpi. Transcriptome analysis and quantitative reverse transcription PCR were performed to determine gene expression profiles in the mandibular LN of sheep. Antigen processing and presentation pathway was the dominant pathway related to the dataset. Our studies suggest that the immune response in ovine LN resembled type 1 immunity with the secretion of IL-12p70 and IFN-γ after B.suis S2 immunization and the vaccine may eliminate Brucella via stimulation of M1 macrophages through the course of Th cells.

An mRNA vaccine elicits STING-dependent antitumor immune responses.

Lipid-formulated RNA vaccines have been widely used for disease prevention and treatment, yet their mechanism of action and individual components contributing to such actions remain to be delineated. Here, we show that a therapeutic cancer vaccine composed of a protamine/mRNA core and a lipid shell is highly potent in promoting cytotoxic CD8+ T cell responses and mediating anti-tumor immunity. Mechanistically, both the mRNA core and lipid shell are needed to fully stimulate the expression of type I interferons and inflammatory cytokines in dendritic cells. Stimulation of interferon-ß expression is exclusively dependent on STING, and antitumor activity from the mRNA vaccine is significantly compromised in mice with a defective Sting gene. Thus, the mRNA vaccine elicits STING-dependent antitumor immunity.

Glutathione ameliorates the meiotic defects of copper exposed ovine oocytes via inhibiting the mitochondrial dysfunctions.

Regardless of the essential role of copper (Cu) in the physiological regulation process of mammalian reproduction, excessive exposure to Cu triggers the meiotic defects of porcine oocytes via compromising the mitochondrial functions. However, the connections between the excessive Cu exposure and meiotic defects of ovine oocytes have not been reported. In this study, the effect of copper sulfate (CuSO4) exposure on the meiotic potentials of ovine oocytes was analyzed. Subsequently, the ameliorative effect of glutathione (GSH) supplementation on the meiotic defects of CuSO4 exposed ovine oocytes was investigated. For these purposes, the in vitro maturation (IVM) of ovine cumulus oocyte complexes (COCs) was conducted in the presence of 5, 10, 20 and 40 µg/mL of CuSO4 supplementation. Subsequently, different concentrations of GSH (2, 4 and 8 mM) were added to the IVM medium containing CuSO4 solution. After IVM, the assay, including nuclear maturation, spindle organization, chromosome alignment, cytoskeleton assembly, cortical granule (CGs) dynamics, mitochondrial function, reactive oxygen species (ROS) generation, apoptosis, epigenetic modification and fertilization capacity of ovine oocytes were performed. The results showed that excessive Cu exposure triggered the meiotic defects of ovine oocytes via promoting the mitochondrial dysfunction related oxidative stress damage. Moreover, the GSH supplementation, not only ameliorated the decreased maturation potential and fertilization defect of CuSO4 exposed oocytes, but inhibited the mitochondrial dysfunction related oxidative stress damage, ROS generation, apoptosis and altered H3K27me3 expression in the CuSO4 exposed oocytes. Combined with the gene expression pattern, the finding in the present study provided fundamental bases for the ameliorative effect of GSH supplementation on the meiotic defects of CuSO4 exposed oocytes via inhibiting the mitochondrial dysfunctions, further benefiting these potential applications of GSH supplementation in the mammalian IVM system and livestock breeding suffering from the excessive Cu exposure.

Effects of dietary oat supplementation on carcass traits, muscle metabolites, amino acid profiles, and its association with meat quality of Small-tail Han sheep.

Oat supplementation of the ruminant diet can improve growth performance and meat quality traits, but the role of muscle metabolites has not been evaluated. This study aimed to establish whether oat grass supplementation (OS) of Small-tail Han sheep improved growth performance and muscle tissue metabolites that are associated with better meat quality and flavor. After 90-day, OS fed sheep had higher live-weight and carcass-weight, and lower carcass fat. Muscle metabolomics analysis showed that OS fed sheep had higher levels of taurine, l-carnitine, inosine-5'-monophospgate, cholic acid, and taurocholic acid, which are primarily involved in taurine and hypotaurine metabolism, purine metabolism, and bile acid biosynthesis and secretion, decreased fat accumulation and they promote functional or flavor metabolites. OS also increased muscle levels of amino acids that are attributed to better quality and flavorsome mutton. These findings provided further evidence for supplementing sheep with oat grass to improve growth performance and meat quality.

Effects of novel variants in BMP15 gene on litter size in Mongolia and Ujimqin sheep breeds.

Reproductive traits, such as ovulation rate and litter size, are important factors influencing the sheep industry. The bone morphogenetic protein 15 (BMP15) is a major gene affecting the reproductive traits in sheep, and multiple mutations in BMP15 gene could affect the ovulation rate and litter size in many sheep breeds, showing high breed specificity. However, identification of novel variations and seeking breed-specific markers associated with litter size in other sheep breeds are still important. In this study, we sequenced the BMP15 gene of Mongolia sheep, and 12 novel variants were detected by direct sequencing and whole-genome resequencing. Among them, the g.50985975 G > A polymorphism in intron and synonymous c.755 T > C (Leu252Pro) in exon 2 of BMP15 were significantly associated with the litter sizes of Mongolia ewes (P < 0.01 and P < 0.01, respectively), as well as the g.50988478C > A and g.50987863G > A in the promoter region of BMP15 were significantly associated with the litter sizes of Ujimqin ewes (P < 0.05 and P < 0.01, respectively). The c.755 T > C mutation is predicted to change the tertiary structure of BMP15. Our findings may provide potentially useful genetic markers for increasing litter size in sheep.