Differential gene expression in ovaries of Qira black sheep and Hetian sheep using RNA-Seq technique

Genomic Regions Associated with Resistance to Gastrointestinal Parasites in Australian Merino Sheep

The objective of this study was to identify genomic regions and genes associated with resistance to gastrointestinal nematodes in Australian Merino sheep in Uruguay, using the single-step GWAS methodology (ssGWAS), which is based on genomic estimated breeding values (GEBVs) obtained from a combination of pedigree, genomic, and phenotypic data. This methodology converts GEBVs into SNP effects. The analysis included 26,638 animals with fecal egg count (FEC) records obtained in two independent parasitic cycles (FEC1 and FEC2) and 1700 50K SNP genotypes. The comparison of genomic regions was based on genetic variances (gVar(%)) explained by non-overlapping regions of 20 SNPs. For FEC1 and FEC2, 18 and 22 genomic windows exceeded the significance threshold (gVar(%) ≥ 0.22%), respectively. The genomic regions with strong associations with FEC1 were located on chromosomes OAR 2, 6, 11, 21, and 25, and for FEC2 on OAR 5, 6, and 11. The proportion of genetic variance attributed to the top windows was 0.83% and 1.9% for FEC1 and FEC2, respectively. The 33 candidate genes shared between the two traits were subjected to enrichment analysis, revealing a marked enrichment in biological processes related to immune system functions. These results contribute to the understanding of the genetics underlying gastrointestinal parasite resistance and its implications for other productive and welfare traits in animal breeding programs.

Genome-wide comparative analysis reveals selection signatures for reproduction traits in prolific Suffolk sheep

The identification of genome-wide selection signatures can reveal the potential genetic mechanisms involved in the generation of new breeds through natural or artificial selection. In this study, we screened the genome-wide selection signatures of prolific Suffolk sheep, a new strain of multiparous mutton sheep, to identify candidate genes for reproduction traits and unravel the germplasm characteristics and population genetic evolution of this new strain of Suffolk sheep. Whole-genome resequencing was performed at an effective sequencing depth of 20× for genomic diversity and population structure analysis. Additionally, selection signatures were investigated in prolific Suffolk sheep, Suffolk sheep, and Hu sheep using fixation index (F ST) and heterozygosity H) analysis. A total of 5,236.338 Gb of high-quality genomic data and 28,767,952 SNPs were obtained for prolific Suffolk sheep. Moreover, 99 selection signals spanning candidate genes were identified. Twenty-three genes were significantly associated with KEGG pathway and Gene Ontology terms related to reproduction, growth, immunity, and metabolism. Through selective signal analysis, genes such as ARHGEF4, CATIP, and CCDC115 were found to be significantly correlated with reproductive traits in prolific Suffolk sheep and were highly associated with the mTOR signaling pathway, the melanogenic pathway, and the Hippo signaling pathways, among others. These results contribute to the understanding of the evolution of artificial selection in prolific Suffolk sheep and provide candidate reproduction-related genes that may be beneficial for the establishment of new sheep breeds.

Assessment of the Effect of Leonurine Hydrochloride in a Mouse Model of PCOS by Gene Expression Profiling

Polycystic ovary syndrome (PCOS) is an endocrine disease commonly associated with metabolic disorders in females. Leonurine hydrochloride (Leo) plays an important role in regulating immunity, tumours, uterine smooth muscle, and ovarian function. However, the effect of Leo on PCOS has not been reported. Here, we used dehydroepiandrosterone to establish a mouse model of PCOS, and some mice were then treated with Leo by gavage. We found that Leo could improve the irregular oestros cycle of PCOS mice, reverse the significantly greater serum testosterone (T) and luteinising hormone (LH) levels, significantly reduce the follicle-stimulating hormone (FSH) level, and significantly increase the LH/FSH ratio of PCOS mice. Leo could also change the phenomenon of ovaries in PCOS mice presented with cystic follicular multiplication and a lacking corpus luteum. Transcriptome analysis identified 177 differentially expressed genes related to follicular development between the model and Leo groups. Notably, the cAMP signalling pathway, neuroactive ligand-receptor interactions, the calcium signalling pathway, the ovarian steroidogenesis pathway, and the Lhcgr, Star, Cyp11a, Hsd17b7, Camk2b, Calml4, and Phkg1 genes may be most related to improvements in hormone levels and the numbers of ovarian cystic follicles and corpora lutea in PCOS mice treated by Leo, which provides a reference for further study of the mechanism of Leo.

Investigating the expression of fertility-regulating LncRNAs in multiparous and uniparous Shal ewe's ovaries

Sheep is the primary source of animal protein in Iran. Birth type is one of the significant features that determine total meat output. Little is known about how long non-coding RNAs (LncRNAs) affect litter size. The purpose of this research is to investigate the DE-LncRNAs in ovarian tissue between multiparous and uniparous Shal ewes. Through bioinformatics analyses, LncRNAs with variable expression levels between ewes were discovered. Target genes were annotated using the DAVID database, and STRING and Cytoscape software were used to evaluate their interactions. The expression levels of 148 LncRNAs were different in the multiparous and uniparous ewe groups (false discovery rate (FDR) < 0.05). Eight biological process terms, nine cellular component terms, 10 molecular function terms, and 38 KEGG pathways were significant (FDR < 0.05) in the GO analysis. One of the most significant processes impacting fertility is mitogen-activated protein kinase (MAPK) signaling pathway, followed by oocyte meiosis, gonadotropin-releasing hormone signaling pathway, progesterone-mediated oocyte maturation, oxytocin signaling pathway, and cAMP signaling pathway. ENSOARG00000025710, ENSOARG00000025667, ENSOARG00000026034, and ENSOARG00000026632 are LncRNAs that may affect litter size and fertility. The most crucial hub genes include MAPK1, BRD2, GAK, RAP1B, FGF2, RAP1B, and RAP1B. We hope that this study will encourage researchers to further investigate the effect of LncRNAs on fertility.

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.

Construction of a circRNA- lincRNA-lncRNA-miRNA-mRNA ceRNA regulatory network identifies genes and pathways linked to goat fertility

Background: There is growing interest in the genetic improvement of fertility traits in female goats. With high-throughput genotyping, single-cell RNA sequencing (scRNA-seq) is a powerful tool for measuring gene expression profiles. The primary objective was to investigate comparative transcriptome profiling of granulosa cells (GCs) of high- and low-fertility goats, using scRNA-seq. Methods: Thirty samples from Ji'ning Gray goats (n = 15 for high fertility and n = 15 for low fertility) were retrieved from publicly available scRNA-seq data. Functional enrichment analysis and a literature mining approach were applied to explore modules and hub genes related to fertility. Then, interactions between types of RNAs identified were predicted, and the ceRNA regulatory network was constructed by integrating these interactions with other gene regulatory networks (GRNs). Results and discussion: Comparative transcriptomics-related analyses identified 150 differentially expressed genes (DEGs) between high- and low-fertility groups, based on the fold change (≥5 and ≤-5) and false discovery rate (FDR <0.05). Among these genes, 80 were upregulated and 70 were downregulated. In addition, 81 mRNAs, 58 circRNAs, 8 lincRNAs, 19 lncRNAs, and 55 miRNAs were identified by literature mining. Furthermore, we identified 18 hub genes (SMAD1, SMAD2, SMAD3, SMAD4, TIMP1, ERBB2, BMP15, TGFB1, MAPK3, CTNNB1, BMPR2, AMHR2, TGFBR2, BMP4, ESR1, BMPR1B, AR, and TGFB2) involved in goat fertility. Identified biological networks and modules were mainly associated with ovary signature pathways. In addition, KEGG enrichment analysis identified regulating pluripotency of stem cells, cytokine-cytokine receptor interactions, ovarian steroidogenesis, oocyte meiosis, progesterone-mediated oocyte maturation, parathyroid and growth hormone synthesis, cortisol synthesis and secretion, and signaling pathways for prolactin, TGF-beta, Hippo, MAPK, PI3K-Akt, and FoxO. Functional annotation of identified DEGs implicated important biological pathways. These findings provided insights into the genetic basis of fertility in female goats and are an impetus to elucidate molecular ceRNA regulatory networks and functions of DEGs underlying ovarian follicular development.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

Comprehensive Analysis of Differentially Expressed CircRNAs in the Ovaries of Low- and High-Fertility Sheep

CircRNAs are essential in regulating follicle growth and development and the female reproductive system at multiple levels. However, the molecular mechanism by which circRNAs regulate reproduction in sheep is unclear and requires further exploration. In this study, RNA sequencing was performed to reveal the circRNA expression profiles in the ovaries of Cele black sheep and Hetian sheep during estrus. Analysis of the number of circRNAs in their host genes revealed that 5031 genes could produce 20,835 circRNAs. Among the differentially expressed circRNAs (DEcircRNA), 75 were upregulated, and 105 were downregulated. Functional enrichment analysis showed that the host genes of DEcircRNA were involved in several pathways, including the MAPK and Hippo signaling pathway. In addition, we constructed a subnetwork of competitive endogenous RNA (ceRNA) containing 4 mRNAs, 4 microRNAs (miRNAs), and 10 circRNAs, potentially related to follicle development. Functional circRNAs (e.g., novel_circ_0003851, novel_circ_0015526, novel_circ_0008117) were found to act as ceRNAs for follicle growth and development-related mRNAs (CUEDC1, KPNB1, ZFPM2) by sponging functional miRNAs (miR-29a, miR-29b, miR-17-5p). Finally, through an RNA pull-down assay, oar-miR-125b was selected and confirmed as the target miRNA of novel-circ-0041512. We analyzed the overall expression of circRNAs in sheep ovaries. Further, we explored the potential mechanisms underlying the circRNA functions, providing a theoretical basis for the genetic progress of reproductive traits in sheep.

Identification of hub genes associated with follicle development in multiple births sheep by WGCNA

Sheep exhibit a distinct estrous cycle that includes four different phases: proestrus, estrus, late estrus, and luteal phase. As the estrous cycle repeats, follicular development regularly alternates. We thus investigated ovarian transcriptome data from each of the four phases using weighted gene co-expression network analysis (WGCNA) to identify modules, pathways, and genes essential to follicle growth and development. We clustered mRNA and long non-coding RNA (lncRNA) into different modules by WGCNA, and calculated correlation coefficients between genes and Stages of the estrous cycle. Co-expression of the black module (cor = 0.81, P<0.001) and the yellow module (cor = 0.61, P<0.04) was found to be critical for follicle growth and development. A total of 2066 genes comprising the black and yellow modules was used for functional enrichment. The results reveal that these genes are mainly enriched in Cell cycle, PI3K-Akt signaling pathway, Oocyte meiosis, Apoptosis, and other important signaling pathways. We also identified seven hub genes (BUB1B, MAD2L1, ASPM, HSD3B1, WDHD1, CENPA, and MXI1) that may play a role in follicle development. Our study may provide several important new markers allowing in depth exploration of the genetic basis for multiparous reproduction in sheep.

Identification of Signatures of Selection for Litter Size and Pubertal Initiation in Two Sheep Populations

Fecundity is an important economic trait in sheep that directly affects their economic and productive efficiency. Our study aimed to identify SNP loci associated with sheep puberty or litter size which could be used in future breeding programs to improve fertility. Genomic DNA was obtained from Hetian and Cele Black sheep breeds and used for reduced-representation genome sequencing to identify SNP loci associated with pubertal initiation and litter size. Selective signatures analysis was performed based on the fixation index and nucleotide diversity, followed by pathway analysis of the genes contained in the selected regions. The selected SNP loci in the genes associated with pubertal initiation and litter size were validated using both sheep breeds. In total, 384,718 high quality SNPs were obtained and 376 genes were selected. Functional annotation of genes and enrichment analysis identified 12 genes associated with pubertal initiation and 11 genes associated with litter size. SNP locus validation showed that two SNP on PAK1 and four on ADCY1 may be associated with pubertal initiation, and one SNP on GNAQ gene (NC_040253.1: g.62677376G > A) was associated with litter size in Cele Black sheep. Our results provide new theoretical support for sheep breeding.

lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA−mRNA and lncRNA−mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein−protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.

Identification of differentially expressed long noncoding RNAs in the ovarian tissue of ewes Shal and Sangsari using RNA-seq

Background

The ovary has an important role in reproductive function. Animal reproduction is dominated by numerous coding genes and noncoding elements. Although long noncoding RNAs (LncRNAs) are important in biological activity, little is known about their role in the ovary and fertility.

Methods

Three adult Shal ewes and three adult Sangsari ewes were used in this investigation. LncRNAs in ovarian tissue from two breeds were identified using bioinformatics analyses, and then target genes of LncRNAs were discovered. Target genes were annotated using the DAVID database, and their interactions were examined using the STRING database and Cytoscape software. The expression levels of seven LncRNAs with their target genes were assessed by real‐time PCR to confirm the RNA‐seq.

Results

Among all the identified LncRNAs, 124 LncRNAs were detected with different expression levels between the two breeds (FDR < 0.05). According to the DAVID database, target genes were discovered to be engaged in one biological process, one cellular component, and 21 KEGG pathways (FDR < 0.05). The PES1, RPS9, EF‐1, Plectin, SURF6, CYC1, PRKACA MAPK1, ITGB2 and BRD2 genes were some of the most crucial target genes (hub genes) in the ovary.

Conclusion

These results could pave the way for future efforts to address sheep prolificacy barriers.

Changes in Transcriptomic Profiles in Different Reproductive Periods in Yaks

Simple Summary

The molecular regulation mechanism of yak ovarian activity has attracted extensive attention. This study investigated the global gene expression profiles in different reproductive stages (anestrus, estrus, and pregnancy) by RNA-seq technology. Enrichment analysis revealed that DEGs were involved in the process of follicular growth, ovulation, and hormone metabolism. This study explored the regulation mechanism of the yak ovary in the reproductive cycle and laid a theoretical foundation for further understanding the reproductive characteristics of yak.

Abstract

Yak reproductive characteristics have received extensive attention, though the molecular regulation mechanism of its ovarian activity remains to be explored. Therefore, this study initially conducted a comparative analysis of yak ovarian activities in anestrus, estrus, and pregnancy regarding their morphology and histology, followed by implementing RNA sequencing (RNA-seq) technology to detect the overall gene expression and biological mechanism in different reproductive stages. H&E staining showed that there were more growing follicles and mature follicles in ovarian tissue sections during estrus than ovarian tissues during non-estrus. The RNA-seq analysis of yak ovary tissues in three periods showed that DEGs related to follicular development and hormone metabolism were screened in the three comparison groups, such as COL1A2, NR4A1, THBS2, PTGS2, SCARB1, STAR, and WNT2B. Bioinformatics analysis showed that these DEGs are involved in ion binding, cell development, metabolic processes, enriched in ECM–receptor interactions, steroid biosynthesis, together with aldosterone generation/discharge and Wnt/PI3K-Akt signaling pathways. In addition, we speculate alternate splice development events to have important role/s in regulating ovarian functional genomic expression profiles. These results provide essential knowledge aimed at scrutinizing pivotal biomarkers for yak ovarian activity, together with paving the way for enhancing researchers’ focus on improving yak reproductive performance.

Genomic variants from RNA-seq for goats resistant or susceptible to gastrointestinal nematode infection

Gastrointestinal nematodes (GIN) are an important constraint in small ruminant production. Genetic selection for resistant animals is a potential sustainable control strategy. Advances in molecular genetics have led to the identification of several molecular genetic markers associated with genes affecting economic relevant traits. In this study, the variants in the genome of Creole goats resistant or susceptible to GIN were discovered from RNA-sequencing. We identified SNPs, insertions and deletions that distinguish the two genotypes, resistant and susceptible and we characterized these variants through functional analysis. The T cell receptor signalling pathway was one of the top significant pathways that distinguish the resistant from the susceptible genotype with 78% of the genes involved in this pathway showing genomic variants. These genomic variants are expected to provide useful resources especially for molecular breeding for GIN resistance in goats.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Identifications of immune-responsive genes for adaptative traits by comparative transcriptome analysis of spleen tissue from Kazakh and Suffolk sheep

Aridity and heat are significant environmental stressors that affect sheep adaptation and adaptability, thus influencing immunity, growth, reproduction, production performance, and profitability. The aim of this study was to profile mRNA expression levels in the spleen of indigenous Kazakh sheep breed for comparative analysis with the exotic Suffolk breed. Spleen histomorphology was observed in indigenous Kazakh sheep and exotic Suffolk sheep raised in Xinjiang China. Transcriptome sequencing of spleen tissue from the two breeds were performed via Illumina high-throughput sequencing technology and validated by RT-qPCR. Blood cytokine and IgG levels differed between the two breeds and IgG and IL-1β were significantly higher in Kazakh sheep than in Suffolk sheep (p < 0.05), though spleen tissue morphology was the same. A total of 52.04 Gb clean reads were obtained and the clean reads were assembled into 67,271 unigenes using bioinformatics analysis. Profiling analysis of differential gene expression showed that 1158 differentially expressed genes were found when comparing Suffolk with Kazakh sheep, including 246 up-regulated genes and 912 down-regulated genes. Utilizing gene ontology annotation and pathway analysis, 21 immune- responsive genes were identified as spleen-specific genes associated with adaptive traits and were significantly enriched in hematopoietic cell lineage, natural killer cell-mediated cytotoxicity, complement and coagulation cascades, and in the intestinal immune network for IgA production. Four pathways and up-regulated genes associated with immune responses in indigenous sheep played indispensable and promoting roles in arid and hot environments. Overall, this study provides valuable transcriptome data on the immunological mechanisms related to adaptive traits in indigenous and exotic sheep and offers a foundation for research into adaptive evolution.

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

Litter size is one of the most important economic traits in sheep. Identification of gene variants that are associated with the prolificacy rate is an important step in breeding program success and profitability of the farm. So, to identify genetic mechanisms underlying the variation in litter size in Iranian Baluchi sheep, a two-step genome-wide association study (GWAS) was performed. GWAS was conducted using genotype data from 91 Baluchi sheep. Estimated breeding values (EBVs) for litter size calculated for 3848 ewes and then used as the response variable. Besides, a pathway analysis using GO and KEGG databases were applied as a complementary approach. A total of three single nucleotide polymorphisms (SNPs) associated with litter size were identified, one each on OAR2, OAR10, and OAR25. The SNP on OAR2 is located within a novel putative candidate gene, Neurotrophic receptor tyrosine kinase 2. This gene product works as a receptor which is essential for follicular assembly, early follicular growth, and oocyte survival. The SNP on OAR25 is located within RAB4A which is involved in blood vessel formation and proliferation through angiogenesis. The SNP on OAR10 was not associated with any gene in the 1Mb span. Moreover, gene-set analysis using the KEGG database identified several pathways, such as Ovarian steroidogenesis, Steroid hormone biosynthesis, Calcium signaling pathway, and Chemokine signaling. Also, pathway analysis using the GO database revealed several functional terms, such as cellular carbohydrate metabolic, biological adhesion, cell adhesion, cell junction, and cell-cell adherens junction, among others. This is the first study that reports the NTRK2 gene affecting litter size in sheep and our study of this gene functions showed that this gene could be a good candidate for further analysis.

Transcriptomic analysis of gene expression in normal goat ovary and intersex goat gonad

Intersexuality is a congenital reproductive disorder that usually occurs in hornless goats, hindering breeding of goats with hornless traits and the development of the goat industry. In this study, we aimed to identify differentially expressed genes in intersex and normal goat gonads by comparing gene transcription profiles of intersex and normal goat gonads. As intersex goats are genetically based on females, we chose female goats as controls. The goats in the control group and the experimental group were both over one-year old. We evaluated the anatomical characteristics of the reproductive organs of five intersex goats using histopathological methods. The gonads were found to be ovarian and testicular types. RNA-Seq technology was used to identify differentially expressed genes in gonads and normal goat ovary tissues. Transcription analysis results were verified by qPCR. The results showed that 2,748 DEGs were upregulated and 3,327 DEGs were downregulated in intersex ovaries unlike in controls, whereas 2006 DEGs were upregulated and 2032 DEGs were downregulated in the interstitial testes. Many of these genes play important roles in mammalian sex determination and sex differentiation, such as SOX9, WT1, GATA4, DMRT1, DHH, AMH, CYP19A1 and FST. We found that many DEGs are involved in biological developmental regulation by GO and KEGG enrichment analyses, and that most genes associated with the steroid synthesis pathway were downregulated. The DEGs identified in this study may be involved in the regulation of intersex goat sex determination and differentiation, and may increase our understanding of the molecular mechanisms of mammalian sex differentiation.

Ovarian transcriptomic analysis and follicular development of Leizhou black duck

This study investigated the factors that caused the differences in egg production during the development of ovarian follicles in Leizhou black ducks. Leizhou black ducks population was divided into 2 groups as high-yield group (HG) and low-yield group (LG). The number of eggs (NE), age at first egg (AFE), weight at first egg, and egg weight (EW) of both groups were recorded, and differences were analyzed using the t test. The logistic model was used to simulate the egg production curves to analyze the production rules. The ovarian follicles of both duck groups were collected to count the number of different grades sized follicles, weigh the ovaries, and observe follicular sections to analyze the developmental differences. Ovarian transcriptomic sequencing was performed to investigate differentially expressed genes and signal pathways in both duck groups. The results revealed a significant difference (P < 0.01) in the NE laid, AFE, and EW between both groups. Comparatively, HG had significantly more (P < 0.01) large yellow follicles (LYF) than LG. The density of medullary layer cells of the follicle section was greater in HG than LG ducks. Transcriptome sequencing revealed a total of 1,027 differentially expressed genes between the HG and LG ducks of which 495 genes were upregulated, and 532 genes were downregulated. Fifty genes were related to reproduction and reproductive processes. Kyoto Encyclopedia of Genes and Genomes–enriched signaling pathways revealed 274 signal pathways enriched in these differentially expressed genes of which the steroid biosynthesis pathway was significantly enriched. Analysis (Q < 0.05) showed that HSD3β → gonadotropin-releasing hormone (GnRH) and estrogen receptor (ESR) → LHβ/ERK1/2 were enriched in the steroid biosynthesis signal pathway. Follicle-stimulating hormone signal pathway mediated by HSD3β → GnRH and ESR → LHβ/ERK1/2 may be involved in ovarian follicle development to regulate LYF reserve process and affect its ovulation cycle, which in turn influence the egg production of Leizhou black ducks.

Inactivation of the htpsA gene affects capsule development and pathogenicity of Streptococcus suis

Streptococcus suis

serotype 2 (S. suis 2) is an important swine pathogen and also an emerging zoonotic agent. HtpsA has been reported as an immunogenic cell surface protein on the bacterium. In the present study, we constructed an isogenic mutant strain of htpsA, namely ΔhtpsA, to study its role in the development and virulence of S. suis 2. Our results showed that the mutant strain lost its typical encapsulated structure with decreased concentrations of sialic acid. Furthermore, the survival rate in whole blood, the anti-phagocytosis by RAW264.7 murine macrophage, and the adherence ability to HEp-2 cells were all significantly affected in the ΔhtpsA. In addition, the deletion of htpsA sharply attenuated the virulence of S. suis 2 in an infection model of mouse. RNA-seq analysis revealed that 126 genes were differentially expressed between the ΔhtpsA and the wild-type strains, including 28 upregulated and 98 downregulated genes. Among the downregulated genes, many were involved in carbohydrate metabolism and synthesis of virulence-associated factors. Taken together, htpsA was demonstrated to play a role in the morphological development and pathogenesis of the highly virulent S. suis 2 05ZYH33 strain.

Vascular morphology and expression of angiogenic factors in sheep ovaries and their relation with the lambing rate

The study aims to present the vascular anatomy of ovarian arterial system in sheep and to explore its relation with the angiogenic factors vascular endothelial growth factor (VEGF) and angiopoietin-2 (Ang-2), environmental factors and the lambing rate. Small-tailed Han sheep, Hu sheep, Aohan finewool sheep and Tibetan sheep who are from different habitats in China and with different lambing rates were investigated. The vascular corrosion casting technique and scanning electron microscopy were used to understand arterial architecture, and real-time quantitative PCR (qRT-PCR) and Western blotting were used for assessing VEGF and Ang-2 expression. Results showed that high-prolificity breeds had significantly larger vessel diameters in ovaries than low-prolificity ones. VEGF gene transcript and protein expression were related with the lambing rate, while Ang-2 gene transcript was negatively related with the lambing rate. Results indicated that high-prolificity breeds have more abundant arterial blood supply, which is associated with VEGF expression.

Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle

Heat stress affects the physiology and production performance of Chinese Holstein dairy cows. As such, the selection of heat tolerance in cows and elucidating its underlying mechanisms are vital to the dairy industry. This study aimed to investigate the heat tolerance associated genes and molecular mechanisms in Chinese Holstein dairy cows using a high-throughput sequencing approach and bioinformatics analysis. Heat-induced physiological indicators and milk yield changes were assessed to determine heat tolerance levels in Chinese Holstein dairy cows by Principal Component Analysis method following Membership Function Value Analysis. Results indicated that rectal temperature (RT), respiratory rate (RR), and decline in milk production were significantly lower (p < 0.05) in heat tolerant (HT) cows while plasma levels of heat shock protein (HSP: HSP70, HSP90), and cortisol were significantly higher (p < 0.05) when compared to non-heat tolerant (NHT) Chinese Holstein dairy cows. By applying RNA-Seq analysis, we identified 200 (81 down-regulated and 119 up-regulated) significantly (|log2fold change| ≥ 1.4 and p ≤ 0.05) differentially expressed genes (DEGs) in HT versus NHT Chinese Holstein dairy cows. In addition, 14 of which were involved in protein–protein interaction (PPI) network. Importantly, several hub genes (OAS2, MX2, IFIT5 and TGFB2) were significantly enriched in immune effector process. These findings might be helpful to expedite the understanding for the mechanism of heat tolerance in Chinese Holstein dairy cows.

Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep (Ovies aries)

Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.

A Transcriptomic Study of the Tail Fat Deposition in Two Types of Hulun Buir Sheep According to Tail Size and Sex

Simple Summary

Based on tail types, Hulun Buir sheep were divided into two lines including small and big fat-tailed, but these two lines have similar genetic background. In this study, we investigated the morphology and transcription level differences of tail fat between these two lines. The RNA-seq analyses indicated several differentially expressed genes when compared between sexes or two tail sizes. Interestingly, we also found an obvious sex difference in the fat metabolism in Hulun Buir sheep. Two different co-expression networks were only shown either in male or in female sheep. Our findings will provide theoretical background in understanding the genetic mechanism of fat deposition in sheep.

Abstract

Hulun Buir sheep of similar genetic background were divided into two lines based on tail types: Small- and big fat-tailed. To explore the molecular mechanism of fat deposition in sheep tails, we firstly evaluated the morphology and transcription level differences of tail fat between these two lines. RNA-Seq technology was used to identify differentially expressed genes (DEGs) in phenotypic extremes of tail sizes. Five comparisons were performed taking into account two factors, sex and tail type. We screened out 373 DEGs between big-tailed and small-tailed Hulun Buir sheep, and 775 and 578 DEGs between two types of tails in male and female sheep, respectively. The results showed an obvious sex difference in the fat metabolism in sheep based on gene ontology (GO), pathway, and network analyses. Intriguingly, there were two different co-expression networks only respectively shown in male and female sheep, which were insulin-related network acting on upstream pathways and PPARG-related network effect in downstream pathways. Furthermore, these two networks were linked by a classic pathway of regulating adipogenesis. This is the first study to investigate the sex differences of fat metabolism in domestic animals, and it demonstrates a new experimental way to study fat metabolism. Our findings will provide theoretical background in understanding the tail-size phenotype in sheep and can be exploited in breeding small-tailed sheep.

Single Nucleotide Polymorphisms in SLC5A1, CCNA1, and ABCC1 and the Association with Litter Size in Small-Tail Han Sheep

Simple Summary

Litter size is one of the most important reproductive traits in sheep. Four single nucleotide polymorphisms (SNPs), g.70067210 T > C in SLC5A1, g.25350431 C > T and g.25360220 T > C in CCNA1, and g.14413132 C > T in ABCC1, were identified by mass spectrometry and may be associated with litter size in sheep. Four SNPs were genotyped in Small-Tail Han, Hu, Cele Black, Suffolk, Sunite, Prairie Tibetan, and Tan sheep, and the expression patterns of SLC5A1, CCNA1, and ABCC1 were determined in Small-Tail Han sheep with different fecundities. Furthermore, we also studied the FecB mutation’s association with litter size in Small-Tail Han sheep. The results indicated that all genes included in this study were differentially expressed in the ovary and uterus of polytocous and monotocous Small-Tail Han sheep. Furthermore, association analysis indicated that both g.70067210 T > C in SLC5A1 and the FecB mutation in BMPR-IB were significantly associated with litter size in Small-Tail Han sheep. Linear regression analysis of the association of multiple markers (FecB and g.70067210 T > C in SCL5A1) with litter size indicated that homozygous ewes carrying the BB/TT genotype had a larger litter size than any other genotype.

Abstract

SLC5A1, CCNA1, and ABCC1 have been extensively studied as candidate genes because of their great influence on the reproductive traits of animals. However, little is known about the association between polymorphisms of the SLC5A1, CCNA1, and ABCC1 genes and litter size in Small-Tail Han sheep. In this study, the expression levels of SLC5A1, CCNA1, and ABCC1 in HPG (hypothalamic–pituitary–gonadal) axis tissues of polytocous and monotocous Small-Tail Han sheep were analyzed by qPCR. To better understand the effects of four single nucleotide polymorphisms (SNPs) comprising of g.70067210 T > C in SLC5A1, g.25350431 C > T and g.25360220 T > C in CCNA1, and g.14413132 C > T in ABCC1, a population genetic analysis was conducted using data obtained from genotyping in 728 sheep from seven breeds. The results indicated that all genes included in this study were differentially expressed in the pituitary and uterus of polytocous and monotocous Small-Tail Han sheep (p < 0.05). The associations of these four SNPs and the FecB mutation with litter size in 384 Small-Tail Han sheep were analyzed, therefore, and it was found that both g.70067210T > C and the FecB mutation were significantly associated with litter size (p < 0.05). The linear regression analysis of the association of multiple markers (FecB and g.70067210 T > C in SCL5A1) with litter size indicated that homozygous ewes carrying the BB/TT genotype had larger litter size than any ewes with any other genotype. In conclusion, the SLC5A1 SNPs significantly affect litter size in sheep and are useful as genetic marker for litter size.

Identification and characterization of mRNAs and lncRNAs in the uterus of polytocous and monotocous Small Tail Han sheep (Ovis aries)

Background

Long non-coding RNAs (lncRNAs) regulate endometrial secretion and uterine volume. However, there is little research on the role of lncRNAs in the uterus of Small Tail Han sheep (FecB++). Herein, RNA-seq was used to comparatively analyze gene expression profiles of uterine tissue between polytocous and monotocous sheep (FecB++) in follicular and luteal phases.

Methods

To identify lncRNA and mRNA expressed in the uterus, the expression of lncRNA and mRNA in the uterus of Small Tail Han sheep (FecB++) from the polytocous group (n = 6) and the monotocous group (n = 6) using RNA-sequencing and real-time polymerase chain reaction (RT-PCR). Identification of differentially expressed lncRNAs and mRNAs were performed between the two groups and two phases . Gene ontology (GO) and pathway enrichment analyses were performed to analyze the biological functions and pathways for the differentially expressed mRNAs. LncRNA-mRNA co-expression network was constructed to further analyses the function of related genes.

Results

In the follicular phase, 473 lncRNAs and 166 mRNAs were differentially expressed in polytocous and monotocous sheep; in the luteal phase, 967 lncRNAs and 505 mRNAs were differentially expressed in polytocous and monotocous sheep. GO and KEGG enrichment analysis showed that the differentially expressed lncRNAs and their target genes are mainly involved in ovarian steroidogenesis, retinol metabolism, the oxytocin signaling pathway, steroid hormone biosynthesis, and the Foxo signaling pathway. Key lncRNAs may regulate reproduction by regulating genes involved in these signaling pathways and biological processes. Specifically, UGT1A1, LHB, TGFB1, TAB1, and RHOA, which are targeted by MSTRG.134747, MSTRG.82376, MSTRG.134749, MSTRG.134751, and MSTRG.134746, may play key regulatory roles. These results offer insight into molecular mechanisms underlying sheep prolificacy.

RNA-seq Transcriptome Analysis in Ovarian Tissue of Pelibuey Breed to Explore the Regulation of Prolificacy

The Pelibuey sheep (Ovis aries) is an indigenous breed distributed in the tropical regions of Mexico. The prolificacy of this sheep is on average from 1 to 1.5 lambs, being an important breeding characteristic that owners seek to increase with the purpose of economic improvements. New-generation RNA sequencing technology has been used to identify the genes that are expressed in the ovarian tissue of sheep that have two or more lambs per parturition, as well as to elucidate the metabolic pathways that are affected by the expression of these genes, with the purpose of better understanding the prolificacy in the sheep. In the present study, the transcriptional expression of multiparous and uniparous sheep was compared using RNA sequencing. Multiparous (M group) and uniparous (U group) sheep that had a genealogical record for three generations (M, n = 5 and U, n = 5) were selected. RNA was extracted from ovarian tissue and subsequently used to prepare the libraries that were sequenced using the Illumina NextSeq500 platform. A total of 31,575 genes were detected from the transcriptomic analysis of which 4908 were significantly expressed (p-value ≤ 0.001) in the ovary of sheep. Subsequently, a second filter was carried out to evaluate the false discovery rate (FDR) and select those genes with p-values ≤ 0.05 and values of expression ≥ 1 (log2), obtaining 354 differential expressed genes (DEG): 120 genes up-regulated and 234 genes down-regulated in the group M with respect to the group U. Through Gene Ontology (GO) and metabolic analysis, we obtained information on the function of differentially expressed genes, and its importance in the reproduction of multiparous sheep. This result suggest that genes identified in the present study participate in the development of the final stages of follicles.

Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq

Background

A multitude of genes tightly regulate ovarian follicular development and hormone secretion. These complex and coordinated biological processes are altered during pregnancy. In order to further understand the regulatory role of these genes during pregnancy, it is important to screen the differentially expressed genes (DEGs) in the ovaries of pregnant and nonpregnant mammals. To detect the genes associated with the development of pregnancy in goats, DEGs from the ovaries from pregnant and nonpregnant Anhui white goats (pAWGs and nAWGs, respectively) were analyzed using RNA sequencing technology (RNA-Seq).

Results

In this study, 13,676,394 and 13,549,560 clean reads were generated from pAWGs and nAWGs, respectively, and 1724 DEGs were identified between the two libraries. Compared with nAWGs, 1033 genes were upregulated and 691 genes were downregulated in pAWGs, including PGR, PRLR, STAR and CYP19A1, which play important roles in goat reproduction. Gene Ontology analysis showed that the DEGs were enriched for 49 functional GO terms. Kyoto Encyclopedia of Genes and Genomes analysis revealed that 397 DEGs were significantly enriched in 13 pathways, including “cell cycle”, “cytokine–cytokine receptor interaction” and “steroid biosynthesis”, suggesting that the genes may be associated with cell cycle regulation, follicular development and hormone secretion to regulate the reproduction process. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data.

Conclusions

The data obtained in this work enrich the genetic resources of goat and provide a further understanding of the complex molecular regulatory mechanisms occurring during the development of pregnancy and reproduction in goats. The DEGs screened in this study may play an important role in follicular development and hormone secretion and they would provide scientific basis for related research in the future.

Electronic supplementary material

The online version of this article (10.1186/s40709-019-0095-9) contains supplementary material, which is available to authorized users.

Streptococcus suis synthesizes deoxyadenosine and adenosine by 5'-nucleotidase to dampen host immune responses

Streptococcus suis is a major porcine bacterial pathogen and emerging zoonotic agent. S. suis 5ʹ-nucleotidase is able to convert adenosine monophosphate to adenosine, resulting in inhibiting neutrophil functions in vitro and it is an important virulence factor. Here, we show that S. suis 5ʹ-nucleotidase not only enables producing 2ʹ-deoxyadenosine from 2ʹ-deoxyadenosine monophosphate by the enzymatic assay and reversed-phase high performance liquid chromatography (RP-HPLC) analysis in vitro, but also synthesizes both 2ʹ-deoxyadenosine and adenosine in mouse blood in vivo by RP-HPLC and liquid chromatography with tandem mass spectrometry analyses. Cellular cytotoxicity assay and Western blot analysis indicated that the production of 2ʹ-deoxyadenosine by 5ʹ-nucleotidase triggered the death of mouse macrophages RAW 264.7 in a caspase-3-dependent way. The in vivo infection experiment showed that 2ʹ-deoxyadenosine synthesized by 5ʹ-nucleotidase caused monocytopenia in mouse blood. The in vivo transcriptome analysis in mouse blood showed the inhibitory effect of 5ʹ-nucleotidase on neutrophil functions and immune responses probably mediated through the generation of adenosine. Taken together, these findings indicate that S. suis synthesizes 2ʹ-deoxyadenosine and adenosine by 5ʹ-nucleotidase to dampen host immune responses, which represents a new mechanism of S. suis pathogenesis.

Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data

Heat stress has a detrimental effect on the physiological and production performance of buffaloes. Elucidating the underlying mechanisms of heat stress is challenging, therefore identifying candidate genes is urgent and necessary. We evaluated the response of buffaloes (n = 30) to heat stress using the physiological parameters, ELISA indexes, and hematological parameters. We then performed mRNA and microRNA (miRNA) expression profiles analysis between heat tolerant (HT, n = 4) and non-heat tolerant (NHT, n = 4) buffaloes, as well as the specific modules, significant genes, and miRNAs related to the heat tolerance identified using the weighted gene co-expression network analysis (WGCNA). The results indicated that the buffaloes in HT had a significantly lower rectal temperature (RT) and respiratory rate (RR) and displayed a higher plasma heat shock protein (HSP70 and HSP90) and cortisol (COR) levels than those of NHT buffaloes. Differentially expressed analysis revealed a total of 753 differentially expressed genes (DEGs) and 16 differentially expressed miRNAs (DEmiRNAs) were identified between HT and NHT. Using the WGCNA analysis, these DEGs assigned into 5 modules, 4 of which were significantly correlation with the heat stress indexes. Interestingly, 158 DEGs associated with heat tolerance in the turquoise module were identified, 35 of which were found within the protein-protein interaction network. Several hub genes (IL18RAP, IL6R, CCR1, PPBP, IL1B, and IL1R1) were identified that significantly enriched in the Cytokine-cytokine receptor interaction. The findings may help further elucidate the underlying mechanisms of heat tolerance in buffaloes.

Comprehensive Expression Profiling Analysis of Pituitary Indicates that circRNA Participates in the Regulation of Sheep Estrus

The pituitary gland is the most important endocrine organ that mainly regulates animal estrus by controlling the hormones synthesis. There is a significant difference between the estrus state and anestrus state of sheep pituitary system. Here, we studied the circular RNA (circRNA) expression profiles of the anterior pituitary of estrus and anestrus sheep using RNA-seq technology. Through this study, we identified a total of 12,468 circRNAs and 9231 differentially expressed circRNAs in the estrus and anestrus pituitary system of sheep. We analyzed some differentially expressed circRNAs by reverse transcription quantitative-PCR (RT-qPCR), and some circRNAs were demonstrated using RNase-R+ resistance experiments. CircRNAs involving the regulation of estrus-related terms and pathways are enriched by using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. In addition, we also predicted partial microRNA-circRNA interaction network for circRNAs that regulate sheep estrus. Overall, this study explored a potential substantial role played by circRNAs involved in pituitary regulation on sheep estrus and proposed new questions for further study.

Genome-wide association study for longevity in the Holstein cattle population

Longevity is one of the most important traits determining dairy cow profitability. In the last decades dairy cows suffered a lowering in the age at culling. With the aim to identify the genes involved in longevity, dates of birth, yields, dates of calving during lifespan and culling dates were collected for 946 culled cows which had been genotyped with the Bovine High Density panel. Using the GenABEL package in R, genome-wide association analysis was performed on three potential traits of longevity: (1) 'days in production,' (2) 'days in herd,' (3) number of calvings over lifespan.' Five genome-wide significant single nucleotide polymorphisms (SNPs) associated with all three longevity traits were detected. Several consecutive SNPs identified on chromosomes 16 and 30 indicated the presence of two suggestive quantitative trait loci (QTL). The genes comprised in the QTL regions had biological functions related to fertility, reproductive disorders, heat stress and welfare of cows. These findings might contribute to improving breeding strategies to improve longevity.

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein-coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.

Ovarian transcriptomic analysis reveals the alternative splicing events associated with fecundity in different sheep breeds

Alternative splicing (AS) is one of the most common mechanisms that accounts for the greater macromolecular and cellular complexity of higher eukaryotic organisms. This study focused on the splicing events in the ovaries of different sheep breeds, namely the Han and Dorset breeds. Of the groups studied, cassette splicing events accounted for the maximum number of the AS events with significant differences, whereas the splicing events that were mutually exclusive with introns accounted for the smallest proportion of splicing events. Greater than 1000 AS events with significant differences were identified between the Han BB and Dorset sheep. The number of AS events with significant differences between Han ++ and Dorset sheep, however, was fewer than that of the comparison of Han BB and Dorset sheep. Seven randomly selected genes with AS events were detected in this study and were validated by reverse transcription-polymerase chain reaction (RT-PCR). Furthermore, there are many genes which were common to the two genotype groups (Han BB and Dorset sheep, as well as the Han ++ and Dorset sheep). In addition, genes detected in the present study were involved in different pathways, including the pathways related with fertility or fecundity. The present study could provide the detailed understanding on the mechanisms of alternative splicing events associated with fecundity in different sheep breeds.

Comparative transcriptomic analysis of high and low egg-producing duck ovaries

The egg-laying rate is an important indicator of egg production of laying ducks. Egg production directly impacts the economic benefits of the duck industry. In order to obtain better insight into the molecular mechanisms associated with the process of egg production, comparative transcriptomic analysis of the ovaries of Jinding ducks with high and low egg production was performed using the Illumina HiSeq 2500 system. A total of 843 differentially expressed genes (DEGs) was identified, 367 that were down-regulated and 476 that were up-regulated in high egg production (HEP) ovaries, as compared with low egg production (LEP) ovaries. Some genes, such as MC5R, APOD, ORAI1, and DYRK4, were more active in HEP ovaries, indicating that these genes may play important roles in regulation of egg production. Among these 843 DEGs, 685 were assigned to gene ontology (GO) categories. Of these, 25 genes were related to reproduction, and 30 were related to the reproductive process, including some associated with ovarian follicle development, circadian regulation of gene expression, circadian rhythm, and estrogen receptor binding. Furthermore, some important functional pathways were revealed, such as the steroid biosynthesis pathway, the endocrine and other factor-regulated calcium reabsorption pathways, circadian rhythm, the neuroactive ligand-receptor interaction pathway, fatty acid biosynthesis, and the calcium-signaling pathway, which appear to be much more active in the HEP group, as compared to those of the LEP group. The results of this study provide very useful information that may contribute to future functional studies of genes involved in bird reproduction.

Identification of Differentially Expressed Genes in Porcine Ovaries at Proestrus and Estrus Stages Using RNA-Seq Technique

Estrus is an important factor for the fecundity of sows, and it is involved in ovulation and hormone secretion in ovaries. To better understand the molecular mechanisms of porcine estrus, the expression patterns of ovarian mRNA at proestrus and estrus stages were analyzed using RNA sequencing technology. A total of 2,167 differentially expressed genes (DEGs) were identified (P ≤ 0.05, |log2  Ratio| ≥ 1), of which 784 were upregulated and 1,383 were downregulated in the estrus compared with the proestrus group. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the cellular process, single-organism process, cell and cell part, and binding and metabolic process. In addition, a pathway analysis showed that these DEGs were significantly enriched in 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including cell adhesion molecules, ECM-receptor interaction, and cytokine-cytokine receptor interaction. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) confirmed the differential expression of 10 selected DEGs. Many of the novel candidate genes identified in this study will be valuable for understanding the molecular mechanisms of the sow estrous cycle.

Integrated ovarian mRNA and miRNA transcriptome profiling characterizes the genetic basis of prolificacy traits in sheep (Ovis aries)

BACKGROUND

The highly prolific breeds of domestic sheep (Ovis aries) are globally valuable genetic resources for sheep industry. Genetic, nutritional and other environmental factors affect prolificacy traits in sheep. To improve our knowledge of the sheep prolificacy traits, we conducted mRNA-miRNA integrated profiling of ovarian tissues from two pure breeds with large (Finnsheep) vs. small (Texel) litter sizes and their F1 crosses, half of which were fed a flushing diet.

RESULTS

Among the samples, 16,402 genes (60.6% known ovine genes) were expressed, 79 novel miRNAs were found, and a cluster of miRNAs on chromosome 18 was detected. The majority of the differentially expressed genes between breeds were upregulated in the Texel with low prolificacy, owing to the flushing diet effect, whereas a similar pattern was not detected in the Finnsheep. F1 ewes responded similarly to Finnsheep rather than displaying a performance intermediate between the two pure breeds.

CONCLUSIONS

The identification and characterization of differentially expressed genes and miRNAs in the ovaries of sheep provided insights into genetic and environmental factors affecting prolificacy traits. The three genes (CST6, MEPE and HBB) that were differentially expressed between the group of Finnsheep and Texel ewes kept in normal diet appeared to be candidate genes of prolificacy traits and will require further validation.

Effect of dietary chitosan oligosaccharide supplementation on the pig ovary transcriptome

Fecundity improvement is one of the most important economic traits for the swine industry. In this study, we identified 486 differentially expressed genes associated with sow prolificacy from COS administrated sow ovaries by RNA-seq.

A comparison of transcriptomic patterns measured in the skin of Chinese fine and coarse wool sheep breeds

We characterised wool traits, and skin gene expression profiles of fine wool Super Merino (SM) and coarse wool Small Tail Han (STH) sheep. SM sheep had a significantly higher total density of wool follicles, heavier fleeces, finer fibre diameter, and increased crimp frequency, staple length and wool grease (lanolin) production. We found 435 genes were expressed at significantly different levels in the skin of the two breeds (127 genes more highly in SM and 308 genes more highly in STH sheep). Classification of the genes more highly expressed in SM sheep revealed numerous lipid metabolic genes as well as genes encoding keratins, keratin-associated proteins, and wool follicle stem cell markers. In contrast, mammalian epidermal development complex genes and other genes associated with skin cornification and muscle function were more highly expressed in STH sheep. Genes identified in this study may be further evaluated for inclusion in breeding programs, or as targets for therapeutic or genetic interventions, aimed at altering wool quality or yield. Expression of the lipid metabolic genes in the skin of sheep may be used as a novel trait with the potential to alter the content or properties of lanolin or the fleece.

Transcriptome analysis of the pectoral muscles of local chickens and commercial broilers using Ribo-Zero ribonucleic acid sequencing

BACKGROUND

The molecular mechanisms underlying meat quality and muscle growth are not clear. The meat quality and growth rates of local chickens and commercial broilers are very different. The Ribo-Zero RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of this study was to provide a reference for studies of the differences in the meat quality and growth of different breeds of chickens.

RESULTS

Ribo-Zero RNA-Seq technology was used to analyze the pectoral muscle transcriptomes of Gushi chickens and AA broilers. Compared with AA broilers, 1649 genes with annotated information were significantly differentially expressed (736 upregulated and 913 downregulated) in Gushi chickens with Q≤0.05 (Q is the P-value corrected by multiple assumptions test) at a fold change ≥2 or ≤0.5. In addition, 2540 novel significantly differentially expressed (SDE) genes (1405 upregulated and 1135 downregulated) were discovered. The results showed that the main signal transduction pathways that differed between Gushi chickens and AA broilers were related to amino acid metabolism. Amino acids are important for protein synthesis, and they regulate key metabolic pathways to improve the growth, development and reproduction of organisms.

CONCLUSION

This study showed that differentially expressed genes in the pectoral tissues of Gushi chickens and AA broilers were related to fat metabolism, which affects meat. Additionally, a large number of novel genes were found that may be involved in fat metabolism and thus may affect the formation of meat, which requires further study. The results of this study provide a reference for further studies of the molecular mechanisms of meat formation.

Genome-wide differential mRNA expression profiles in follicles of two breeds and at two stages of estrus cycle of gilts

Estrus expression by gilts and sows is hereditable and important for heat detection. To better understand the molecular biological mechanisms of estrus expression in gilts, the mRNA expression profiles of follicular tissue from Large White gilts in diestrus (LD, n = 3) and estrus (LE, n = 3), and Chinese indigenous Mi gilts in diestrus (MD, n = 2) and estrus (ME, n = 3) were investigated using RNA sequencing. We detected 122,804-335,295 SNPs, 6,140-14,947 InDel and 12 types of AS events (39.57% TSS, 34.90% TTS) in 11 samples. A total of 2,838 differentially expressed genes (DEGs) were found in LD vs MD, LE vs ME, LE vs LD, or ME vs MD comparisons. Two DEGs (ACP5 and PIGS) were observed in all comparisons. Two new genes (ENSSSCG00000028235 and ENSSSCG00000021903) were exclusively expressed in Mi and Large White gilts, respectively. Bioinformatics analyses indicate that these DEGs are involved in single-organism process, catalytic activity, cell adhesion and enriched in ECM-receptor interaction, olfactory transduction, ovarian steroidogenesis, steroid biosynthesis and CAMs signaling pathways. These results of RNA-Seq have provided important information for screening the key functional genes or molecular markers of estrus expression in gilts.

Expression Analysis of Genes Associated with Prolificacy in FecB Carrier and Noncarrier Indian Sheep

The effect of FecB mutation on the gene expression in FecB carrier and noncarrier estrous synchronized ewes, has been analyzed. For this study the whole ovarian tissues and Graafian follicles were collected from estrus synchronized FecB carrier Garole, and non-carrier Deccani Indian sheep, showing remarkable differences in the numbers of preovulatory follicles among two groups. Eleven potential candidate genes (BMP15, GDF9, BMP4, BMP7, BMPR1B, BMPR1A, SMAD9, LHCGR, FSHR, IGF1R, and STAT5) were selected for their expression analysis by SybrGreen based real-time PCR, across ovaries and Graafian follicles of different fecundity groups, for having better insights into the effect of FecB genotypes on follicular development. Variable expression was observed for almost all the genes included in the present study among high and low fecundity groups that was most significant for the BMP7, BMP4, LHCGR, and FSHR transcripts in the ovarian follicles of high and low fecundity ewes, indicating their importance in governing the fecundity in FecB carrier, Indian Garole sheep. BMP4 expression among the genes studied was significantly higher in FecB carrier Garole sheep. This study confirms the changes in mRNA expression of the genes implicated in follicular development in FecB carrier and noncarrier Indian sheep breeds.

A Streptococcus suis LysM domain surface protein contributes to bacterial virulence

Streptococcus suis (SS) is a major swine pathogen, as well as a zoonotic agent for humans. Numerous factors contribute to SS virulence, but the pathogenesis of SS infection is poorly understood. Here, we show that a novel SS surface protein containing a LysM at the N-terminus (SS9-LysM) contributes to SS virulence. Homology analysis revealed that the amino acid sequence of SS9-LysM from the SS strain GZ0565 shares 99.8-68.7% identity with homologous proteins from other SS strains and 41.2% identity with Group B Streptococcal protective antigen Sip. Immunization experiments showed that 7 out of 30 mice immunized with recombinant SS9-LysM were protected against challenge with the virulent GZ0565 strain, while all of the control mice died within 48h following bacterial challenge. In mouse infection model, the virulence of the SS9-LysM deletion mutant (ΔSS9-LysM) was reduced compared with the wild-type (WT) strain GZ0565 and SS9-LysM complemented strain. In addition, ΔSS9-LysM was significantly more sensitive to killing by pig blood ex vivo and mouse blood in vivo compared with the WT strain and SS9-LysM complemented strain. In vivo transcriptome analysis in mouse blood showed that the WT strain reduced the expression of host genes related to iron-binding by SS9-LysM. Moreover, the total free iron concentration in blood from infected mice was significantly lower for the ΔSS9-LysM strain compared with the WT strain. Together, our data reveal that SS9-LysM facilitates SS survival within blood by releasing more free iron from the host. This represents a new mechanism of SS pathogenesis.

De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.

Transcriptomic Analysis of Ovaries from Pigs with High And Low Litter Size

Litter size is one of the most important economic traits for pig production as it is directly related to the production efficiency. Litter size is affected by interactions between multiple genes and the environment. While recent studies have identified some genes associated with prolificacy in pigs, transcriptomic studies of specific genes affecting litter size in porcine ovaries are rare. In order to identify candidate genes associated with litter size in swine, we assessed gene expression differences between the ovaries of Yorkshire pigs with extremely high and low litter sizes using the RNA-Seq method. A total of 1 243 differentially expressed genes were identified: 897 genes were upregulated and 346 genes were downregulated in high litter size ovary samples compared with low litter size ovary samples. A large number of these genes related to steroid hormone regulation in animal ovaries, including 59 Gene Ontology terms and 27 Kyoto Encyclopedia of Genes and Genomes pathways involved in steroid biosynthesis and ovarian steroidogenesis. From these differentially expressed genes, we identified a total of 11 genes using a bioinformatics screen that may be associated with high litter size in Yorkshire pigs. These results provide a list of new candidate genes for porcine litter size and prolificacy to be further investigated.