Genome-wide association study and pathway analysis identify NTRK2 as a novel candidate gene for litter size in sheep

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveals that those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

NTRK2 Promotes Sheep Granulosa Cells Proliferation and Reproductive Hormone Secretion and Activates the PI3K/AKT Pathway

Neurotrophin receptor B (NTRK2), also named TRKB, belongs to the neurotrophic factor family. Previous studies have shown that NTRK2 is associated with high fertility in mammals. However, the molecular mechanism and regulatory pathway of this neurotrophic factor remain unclear. In this study, NTRK2 overexpression and NTRK2-siRNA were constructed to detect the effects of NTRK2 on the proliferation and hormone secretion of the ovarian granulosa cells (GCs) of sheep. We successfully isolated follicular phase granulosa cells in vitro from the ovaries of sheep in simultaneous estrus, and the immunofluorescence results confirmed that NTRK2 was expressed in the collected cells. Subsequently, the effect of NTRK2 on the proliferation of sheep granulosa cells was examined via cell transfection experiments. The results showed that the expression of CDK4 and CyclinD2 was significantly increased after NTRK2 overexpression, while the opposite trend was observed after the inhibition of NTRK2 expression (p < 0.05). The EdU and CCK-8 assays showed that the proliferation rate of sheep GCs was significantly increased after NTRK2 overexpression, while the opposite trend was observed after the inhibition of NTRK2 expression (p < 0.05). Moreover, NTRK2 significantly increased the expression of steroidogenesis-related genes, including steroidogenic acute regulatory protein (STAR) and hydroxy-δ-5-steroid dehydrogenase (HSD3B1), and cytochrome P450 family 19 subfamily A member 1 (CYP19A1). The ELISA results showed that the secretion levels of E2 and P4 significantly increased after NTRK2 overexpression, while the opposite trend was observed after the inhibition of NTRK2 expression (p < 0.05). Previous studies had confirmed that NTRK2 gene belongs to the PI3K-AKT signaling pathway and participates in the signaling of this pathway. This was demonstrated by protein-protein interaction analysis and NTRK2 belongs to the PI3K-AKT pathway. The modification of PI3K and AKT, markers of the PI3K-AKT pathway, via phosphorylation was increased after NTRK2 overexpression in the sheep GCs, while the opposite trend was observed after the inhibition of NTRK2 expression (p < 0.05). Overall, these results suggest that the NTRK2 gene regulates the proliferation of GCs and the secretion of steroid hormones in sheep, and that it influences the phosphorylation level of the PI3K/AKT signaling pathway. These findings provided a theoretical basis and new perspectives for exploring the regulation of NTRK2 gene in the development of ovine follicles.

Characterization and Comparative Analysis of Whole-Transcriptome Sequencing in High- and Low-Fecundity Chongming White Goat Ovaries during the Estrus Phase

Reproductive performance is one of the most important economic traits in the goat industry. Increasing the number of goats is an effective measure to improve production efficiency and reduce production costs. Ovaries are important reproductive organs in female mammals that directly affect the estrous cycle and reproductive abilities. Understanding the complex transcription network of non-coding RNAs (lncRNAs, circRNAs, and miRNAs) and messenger RNA (mRNA) could lead to significant insights into the ovarian regulation of the reproductive processes of animals. However, the whole-transcriptome analysis of the non-coding RNAs and mRNA of the ovaries in Chongming white goats between high-fecundity (HP) and low-fecundity (LP) groups is limited. In this study, a whole-transcriptome sequencing approach was used to identify lncRNA, circRNA, miRNA, and mRNA expression in the ovaries of Chongming white goats during the estrus phase using RNA-Seq technology. More than 20,000 messenger RNAs (mRNAs), 10,000 long non-coding RNAs (lncRNAs), 3500 circular RNAs (circRNAs), and 1000 micro RNAs (miRNAs) were identified. A total of 1024 differential transcripts (724 mRNAs, 112 lncRNAs, 178 circRNAs, and 10 miRNAs) existing between the HP and the LP groups were revealed through a bioinformatics analysis. They were enriched in the prolactin signaling pathway, the Jak-STAT signaling pathway, and the GnRH signaling pathway, as well as various metabolic pathways. Differentially expressed mRNAs (such as LYPD6, VEGFA, NOS3, TNXB, and EPHA2) and miRNAs (such as miR-10a-5p) play key roles in the regulation of goat ovaries during the estrus phase. The enrichment of pathways related to reproduction, such as the Hippo, Hedgehog, PI3K-AKT, and MAPK signaling pathways, suggests that they might be involved in the prolificacy of goat ovaries. Overall, we identified several gene modules associated with goat fecundity and provided a basis for a molecular mechanism in the ovaries of Chongming white goats.

Validation and Analysis of COIL, a Gene Associated with Multiple Lambing Traits in Sheep

In a past study, the team used specific-locus amplified fragment sequencing (SLAF sequencing) to detect single-nucleotide polymorphisms (SNPs) contributing to the differences in lambing numbers in Xinjiang sheep. This study verified the correlation between the COIL gene and lambing number characters in sheep and explored its possible mechanism of action. In this study, three SNPs in the COIL gene, namely COILSNP1 (rs7321466), COILSNP2 (rs7314134), and COILSNP3 (rs7321563), were explored in terms of their possible mechanism of action. A tissue expression profiling analysis revealed that the COIL gene was significantly more expressed in the uterus and ovaries than in other tissues (p < 0.05), whereas an association analysis revealed that the number of lambs born was significantly different among individuals with different genotypes of this COILSNP1 (p < 0.05). The Cell Counting Kit-8(CCK-8) revealed that the overexpression of the COIL gene significantly increased the proliferation of mouse ovarian fibroblasts and sheep fibroblasts (p < 0.05). The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) revealed that the overexpression of the COIL gene significantly increased the activity of sheep fibroblasts (p < 0.01) and mouse ovarian fibroblasts (p < 0.05). The overexpression of the COIL gene affected the biogenesis pathway of spliceosomal U snRNPs by validating protein network connections. This activity affects ovulation, embryonic development, and changes in lambing size in sheep.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

Genome-Wide Association Analysis of Muscle pH in Texel Sheep × Altay Sheep F2 Resource Population

pH was one of the important meat quality traits, which was an important factor affecting the storage/shelf life and quality of meat in meat production. In order to find a way to extend the storage/shelf life, the pH values (pH_45min, pH_24h, pH_48h and pH_72h) of the longissimus dorsi muscles in F₂ individuals of 462 Texel sheep × Altay sheep were determined, genotyping was performed using Illumina Ovine SNP 600 K BeadChip and whole genome resequencing technology, a genome-wide association analysis (GWAS) was used to screen the candidate genes and molecular markers for pH values related to the quality traits of mutton, and the effects of population stratification were detected by Q-Q plots. The results showed that the pH population stratification analysis did not find significant systemic bias, and there was no obvious population stratification effect. The results of the association analysis showed that 28 SNPs significantly associated with pH reached the level of genomic significance. The candidate gene associated with pH_45min was identified as the CCDC92 gene by gene annotation and a search of the literature. Candidate genes related to pH_24h were KDM4C, TGFB2 and GOT2 genes. The candidate genes related to pH_48h were MMP12 and MMP13 genes. The candidate genes related to pH_72h were HILPDA and FAT1 genes. Further bioinformatics analyses showed 24 gene ontology terms and five signaling pathways that were significantly enriched (p ≤ 0.05). Many terms and pathways were related to cellular components, processes of protein modification, the activity of protein dimerization and hydrolase activity. These identified SNPs and genes could provide useful information about meat and the storage/shelf life of meat, thereby extending the storage/shelf life and quality of meat.

Comparative Hypothalamic Transcriptome Analysis Reveals Crucial mRNAs, lncRNAs, and circRNAs Affecting Litter Size in Goats

Litter size is an important indicator to measure the reproductive performance of goats, which is affected by the reproductive function of animals. The hypothalamus, as the regulatory center of the endocrine system, plays an important role in the reproduction of female animals. Here, we performed high-throughput RNA sequencing using hypothalamic tissue from high-fecundity and low-fecundity Leizhou goats to explore critical functional genes associated with litter size. Differentially expressed mRNA, lncRNA, and circRNAs were screened using DESeq and were enriched, and then analyzed by Gene Ontology and Kyoto Encyclopedia of Gene and Genome. Results showed that some of these differentially expressed mRNAs could be enriched in reproductive processes, jak-STAT, prolactin signaling pathway, and other signaling pathways related to reproduction, such as SOCS3. Furthermore, the central proteins POSTN, MFAP5, and DCN from protein-protein interaction may regulate animal reproductive activity by affecting cell proliferation and apoptosis. lncRNA MSTRG.33887.2 as well as circRNAs chicirc_098002, chicirc_072583, and chicirc_053531 may be able to influence animal reproduction by participating in folate metabolism and energy metabolism homeostasis through their respective target genes. Our results expand the molecular mechanism of hypothalamic regulation on animal reproduction.

Role of SNPs located in the exon 9 of ATAPA1 gene on goose egg production

The meat and egg of goose is one of the main components of human food supply. The improvement of goose egg production is particularly important for the increasing human population. However, limited information is available about the effective molecular markers and mechanisms of egg production in goose. In this study, we jointly utilized the data of genome resequencing in different egg production Sichuan white goose and transcriptome at different follicle development stages to identified the molecular markers and mechanisms of egg production. The coefficient of variation of individual egg production in Sichuan white goose population is 0.42 to 0.49. Fifty individuals with the highest (laying 365 days egg number, LEN365 = 79-145) and 50 individuals with the lowest (LEN365 = 8-48) egg production were divided into high and low egg production groups. Based on whole-genome sequencing data of the selected samples, 36 SNPs (annotation novel.12.470, CELF2, ATP1A1, KCNJ6, RAB4A, UST, REV3L, DHX15, CAVN2, SLC5A9, Cldn5, MRPS23, and Tspan2) associated with the LEN365 were identified, involving multiple pathways such as metabolism and endocrinology. Notably, 5 SNPs located in the exon9 of ATP1A1 were identified by GWAS analysis. The association analysis with LEN365 showed the phenotypic variance explained of this haplotype consisting of 5 SNPs is 20.51%. Through transcriptome data analysis, we found the expression of ATP1A1 in the granular layers was increased in the stage of small yellow follicle to large yellow follicle (LYF) and LYF to F5, while decreased in F2 to F1. For the first time, we report the haplotype region formed by 5 SNPS on exon9 of ATP1A1 is associated with egg production in goose and involved in follicle selection and maturation processes.

A Novel c.100C > G Mutation in the FST Gene and Its Relation With the Reproductive Traits of Awassi Ewes

Reproductive traits are affected by many factors, including ovarian function, hormones, and genetics. Genetic polymorphisms of candidate genes are associated with reproductive traits. Several candidate genes are associated with economic traits, including the follistatin (FST) gene. Thus, this study aimed to evaluate whether the genetic variations in the FST gene are associated with the reproductive traits in Awassi ewes. The genomic DNA was extracted from 109 twin ewes and 123 single-progeny ewes. Therefore, 4 sequence fragments from the FST gene were amplified using polymerase chain reaction (PCR) (exon 2/240, exon 3/268, exon 4/254, and exon 5/266 bp, respectively). For a 254 bp amplicon, 3 genotypes were identified: CC, CG, and GG. Sequencing revealed a novel mutation in CG genotypes c.100C > G. The statistical analysis of c.100C > G showed an association with reproductive characteristics. Ewes carrying the c.100C > G had significantly (P ⩽ .01) lower litter sizes, twinning rates, lambing rates, and more days to lambing compared with CG and CC genotypes. Logistic regression analysis confirmed that the c.100C > G single-nucleotide polymorphism (SNP) is responsible for decreasing litter size. According to these results, the variant c.100C > G negatively affects the traits of interest and is associated with lower reproductive traits in Awassi sheep. As a result of this study, ewes carrying the c.100C > G SNP have lower litter size and are less prolific.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Expression of long noncoding RNAs in the ovarian granulosa cells of women with diminished ovarian reserve using high-throughput sequencing

Background

Infertility is a global reproductive-health problem, and diminished ovarian reserve (DOR) is one of the common causes of female infertility. Long noncoding RNAs (lncRNAs) are crucial regulators of numerous physiological and pathological processes in humans. However, whether lncRNAs are involved in the development of DOR remains to be elucidated.

Methods

Ovarian granulosa cells (OGCs) extracted from infertile women with DOR and from women with normal ovarian reserve (NOR) were subjected to high-throughput sequencing. Comprehensive bioinformatics analysis was conducted to identify the differential expression of messenger RNAs (mRNAs) and lncRNAs. Sequencing results were validated by the selection of lncRNAs and mRNAs using real-time reverse transcription-quantitative polymerase chain reaction (RT-qPCR).

Results

Compared with the NOR group, a total of 244 lncRNAs were upregulated (53 known and 191 novel), and 222 lncRNAs were downregulated (36 known and 186 novel) in the DOR group. Similarly, 457 mRNAs had differential expression between the two groups. Of these, 169 were upregulated and 288 were downregulated. Bioinformatics analysis revealed that the differentially expressed genes of mRNA and lncRNAs were considerably enriched in “cell adhesion and apoptosis”, “steroid biosynthesis”, and “immune system”. A co-expression network comprising lncRNAs and their predicted target genes revealed the possible involvement of the “thyroid hormone signaling pathway” and “protein binding, digestion and absorption” in DOR pathogenesis. The expression of SLC16A10 was positively regulated by multiple lncRNAs. After RT-qPCR validation of seven differentially expressed lncRNAs and mRNAs, respectively, the expression of lncRNA NEAT1, GNG12, ZEB2-AS1, and mRNA FN1, HAS3, RGS4, SUOX were in accordance with RNA-sequencing.

Conclusions

We presented the first data showing that the expression profiles of lncRNA and mRNA in OGCs between NOR and DOR patients using RNA sequencing. The lncRNAs and mRNAs that we identified may serve as novel diagnostic biomarkers for patients with DOR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-022-01053-6.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Genome-Wide Identification of mRNAs, lncRNAs, and Proteins, and Their Relationship With Sheep Fecundity

The exploration of multiple birth-related genes has always been a significant focus in sheep breeding. This study aimed to find more genes and proteins related to the litter size in sheep. Ovarian specimens of Small Tail Han sheep (multiple births) and Xinji Fine Wool sheep (singleton) were collected during the natural estrus cycle. Transcriptome and proteome of ovarian specimens were analyzed. The transcriptome results showed that "steroid hormone biosynthesis" and "ovarian steroidogenesis" were significantly enriched, in which HSD17B1 played an important role. The proteome data also confirmed that the differentially expressed proteins (DEPs) were enriched in the ovarian steroidogenesis pathway, and the CYP17A1 was the candidate DEP. Furthermore, lncRNA MSTRG.28645 was highly expressed in Small Tailed Han sheep but lowly expressed in Xinji fine wool sheep. In addition, MSTRG.28645, a hub gene in the co-expression network between mRNAs and lncRNAs, was selected as one of the candidate genes for subsequent verification. Expectedly, the overexpression and interference of HSD17B1 and MSTRG.28645 showed a significant effect on hormone secretion in granulosa cells. Therefore, this study confirmed that HSD17B1 and MSTRG.28645 might be potential genes related to the fecundity of sheep. It was concluded that both HSD17B1 and MSTRG.28645 were critical regulators in the secretion of hormones that affect the fecundity of the sheep.

Meta-analysis of genome-wide association studies for litter size in sheep

Litter size and ovulation rate are important reproduction traits in sheep and have important impacts on the profitability of farm animals. To investigate the genetic architecture of litter size, we report the first meta-analysis of genome-wide association studies (GWAS) using 522 ewes and 564,377 SNPs from six sheep breeds. We identified 29 significant associations for litter size which 27 of which have not been reported in individual GWAS for each population. However, we could confirm the role of BMPR1B in prolificacy. Our gene set analysis discovered biological pathways related to cell signaling, communication, and adhesion. Functional clustering and enrichment using protein databases identified epidermal growth factor-like domain affecting litter size. Through analyzing protein-protein interaction data, we could identify hub genes like CASK, PLCB4, RPTOR, GRIA2, and PLCB1 that were enriched in most of the significant pathways. These genes have a role in cell proliferation, cell adhesion, cell growth and survival, and autophagy. Notably, identified SNPs were scattered on several different chromosomes implying different genetic mechanisms underlying variation of prolificacy in each breed. Given the different layers that make up the follicles and the need for communication and transfer of hormones and nutrients through these layers to the oocyte, the significance of pathways related to cell signaling and communication seems logical. Our results provide genetic insights into the litter size variation in different sheep breeds.

Genes and Pathways Affecting Sheep Productivity Traits: Genetic Parameters, Genome-Wide Association Mapping, and Pathway Enrichment Analysis

Ewe productivity is a composite and maternal trait that is considered the most important economic trait in sheep meat production. The objective of this study was the application of alternative genome-wide association study (GWAS) approaches followed by gene set enrichment analysis (GSEA) on the ewes’ genome to identify genes affecting pregnancy outcomes and lamb growth after parturition in Iranian Baluchi sheep. Three maternal composite traits at birth and weaning were considered. The traits were progeny birth weight, litter mean weight at birth, total litter weight at birth, progeny weaning weight, litter mean weight at weaning, and total litter weight at weaning. GWASs were performed on original phenotypes as well as on estimated breeding values. The significant SNPs associated with composite traits at birth were located within or near genes RDX, FDX1, ARHGAP20, ZC3H12C, THBS1, and EPG5. Identified genes and pathways have functions related to pregnancy, such as autophagy in the placenta, progesterone production by the placenta, placental formation, calcium ion transport, and maternal immune response. For composite traits at weaning, genes (NR2C1, VEZT, HSD17B4, RSU1, CUBN, VIM, PRLR, and FTH1) and pathways affecting feed intake and food conservation, development of mammary glands cytoskeleton structure, and production of milk components like fatty acids, proteins, and vitamin B-12, were identified. The results show that calcium ion transport during pregnancy and feeding lambs by milk after parturition can have the greatest impact on weight gain as compared to other effects of maternal origin.