Non-Coding Transcriptome Provides Novel Insights into the Escherichia coli F17 Susceptibility of Sheep Lamb

SIMPLE SUMMARY: Diarrhea and vomiting caused by Escherichia coli (E. coli) F17 are considered significant threats to animal farming. In the present study, RNA-Seq was performed to investigate the potential circRNA and miRNA biomarkers for E. coli F17-antagonism (AN) and -sensitive (SE) lambs. The results indicated that circRNA and miRNA expression is closely associated with the susceptibility of E. coli F17 in lambs. Numbers of circRNAs and miRNAs may serve as potential biomarkers for intestinal inflammatory response against E. coli F17 infection. Our study can provide a preliminary understanding of the underlying mechanisms of intestinal immunity. ABSTRACT: It has long been recognized that enterotoxigenic Escherichia coli (ETEC) is the major pathogen responsible for vomiting and diarrhea. E. coli F17, a main subtype of ETEC, is characterized by high morbidity and mortality in young livestock. However, the transcriptomic basis underlying E. coli F17 infection has not been fully understood. In this study, RNA sequencing was performed to explore the expression profiles of circRNAs and miRNAs in the jejunum of E. coli F17-antagonism (AN) and -sensitive (SE) lambs. A total of 16,534 circRNAs and 271 miRNAs (125 novel miRNAs and 146 annotated miRNAs) were screened, and 214 differentially expressed (DE) circRNAs and 53 DE miRNAs were detected between the AN and SE lambs (i.e., novel_circ_0025840, novel_circ_0022779, novel_miR_107, miR-10b). Functional enrichment analyses showed that source genes of DE circRNAs were mainly involved in metabolic-related pathways, while target genes of DE miRNAs were mainly enriched in the immune response pathways. Then, a two-step machine learning approach combining Random Forest (RF) and XGBoost (candidates were first selected by RF and further assessed by XGBoost) was performed, which identified 44 circRNAs and 39 miRNAs as potential biomarkers (i.e., novel_circ_0000180, novel_circ_0000365, novel_miR_192, oar-miR-496-3p) for E. coli infection. Furthermore, circRNA-related and lncRNA-related ceRNA networks were constructed, containing 46 circRNA-miRNA-mRNA competing triplets and 630 lncRNA-miRNA-mRNA competing triplets, respectively. By conducting a serious of bioinformatic analyses, our results revealed important circRNAs and miRNAs that could be potentially developed as candidate biomarkers for intestinal inflammatory response against E. coli F17 infection; our study can provide novel insights into the underlying mechanisms of intestinal immunity

Defining ovine dermal papilla cell markers and identifying key signaling pathways regulating its intrinsic properties

Dermal papilla cell (DPC), one of the key cell types during hair follicle development and regeneration, specifies hair size, shape and cycling. It is also an important in vitro screening model for hair growth. Although some characteristics of DPCs, such as agglutinative growth and marker genes, have been studied in mice and humans, the intrinsic properties of ovine DPCs and the regulatory mechanism of the intrinsic properties during continued culture in vitro remained unknown. In this study, based on our previous single-cell transcriptome sequencing on sheep lambskin, we verified SOX18 and PDGFRA as the novel marker genes of ovine DPCs through immunofluorescence staining on skin sections and cultured DPCs. Using continued cell culture and alkaline phosphatase staining, we found that different from mice and humans, ovine DPCs exhibit particularly robust and stable aggregation with unbated alkaline phosphatase activity till 30 passages during continued culture in vitro. Also, we found that the expression of some marker genes and the activity of Wnt/β-catenin signaling differ between early passaged DPCs and multiple passaged DPCs. Further, using Wnt/β-catenin agonist and antagonist, we demonstrated that Wnt/β-catenin signaling could regulate cell aggregation and alkaline phosphatase activity of ovine DPCs through regulating FGF and IGF signaling. This study provides the basis for isolating ovine DPCs and defines their intrinsic properties, which contribute to improving wool performance and medicine of hair regeneration

MicroRNA-181a Regulates the Proliferation and Differentiation of Hu Sheep Skeletal Muscle Satellite Cells and Targets the YAP1 Gene

MicroRNA (miRNA) is of great importance to muscle growth and development, including the regulation of the proliferation and differentiation of skeletal muscle satellite cells (SMSCs). In our research group’s previous study, we found that miR-181a is differentially expressed in the longissimus dorsi muscle of Hu sheep at different stages. We speculated that miR-181a may participate in the growth and development process of Hu sheep. To understand the mechanism of miR-181a regulating the growth and development of Hu sheep skeletal muscle, we extracted skeletal muscle satellite cells from the longissimus dorsi muscle of 3-month-old Hu sheep fetuses and performed a series of experiments. Our results showed that miR-181a suppressed SMSCs’ proliferation using QRT-PCR, Western blot, CCK-8, EDU, and Flow cytometry cycle tests. In addition, QRT-PCR, Western blot, and immunofluorescence indicated that miR-181a facilitated the differentiation of SMSCs. Then, we used dual-luciferase reporter gene detection, QRT-PCR, and Western blot to find that the Yes1-related transcription regulator (YAP1) is the target gene of miR-181a. Our study supplies a research basis for understanding the regulation mechanism of miR-181a on the growth of Hu sheep skeletal muscle

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat

Background Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear. Methods In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs. Results In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition

Differential Expression of <i>KCNJ12</i> Gene and Association Analysis of Its Missense Mutation with Growth Traits in Chinese Cattle

The potassium inwardly rectifying channel, subfamily J, member 12 (KCNJ12) gene is a promising candidate for economic traits because of its crucial roles in myoblast development. Here, a missense mutation (Cys &gt; Arg) was first detected to be located in exon 3 of KCNJ12 from three Chinese cattle breeds by DNA-pool sequencing. Then, we performed an association analysis of this single-nucleotide polymorphism (SNP) with stature in three Chinese cattle populations (n = 820). A significantly positive correlation was revealed by a reduced animal general linear model and the CC genotype was the most favorable in three breeds. Further, we measured the expression profile of the KCNJ12 gene in various cattle tissues and primary bovine skeletal muscle cells. Ubiquitous expression with high abundance in muscle was observed. Further, in primary bovine skeletal muscle cells, the KCNJ12 mRNA expression was gradually up-regulated in differentiation medium (DM) compared with that in growth medium (GM), suggesting that the KCNJ12 gene is involved in bovine myocyte differentiation. Conclusively, the KCNJ12 gene is a functional candidate gene which can be used as a molecular marker for cattle breeding

Integrated Hair Follicle Profiles of microRNAs and mRNAs to Reveal the Pattern Formation of Hu Sheep Lambskin

Hair follicle development is closely associated with wool curvature. Current studies reveal the crucial role of microRNAs (miRNAs) in hair follicle growth and development. However, few studies are known regarding their role in wool curvature. To reveal the potential roles of miRNAs in Hu sheep lambskin with different patterns, a total of 37 differentially expressed (DE) miRNAs were identified in hair follicles between small waves (SM) and straight wool (ST) groups using RNA-seq. Through functional enrichment and miRNA-mRNA co-expression analysis, some key miRNAs (oar-miR-143, oar-miR-200b, oar-miR-10a, oar-miR-181a, oar-miR-10b, oar-miR-125b, etc.) and miRNA-mRNA pairs (miR-125b target CD34, miR-181a target FGF12, LMO3, miR-200b target ZNF536, etc.) were identified. Though direct or indirect ways affecting hair follicle development, these miRNAs and mRNAs may have possible effects on wool curvature, and this study thus provides valuable insight on potential pattern formation

Whole-genome resequencing of Dorper and Hu sheep to reveal selection signatures associated with important traits

Dorper and Hu sheep exhibit different characteristics in terms of reproduction, growth, and meat quality. Comparison of the genomes of two breeds help to reveal important genomic information. In this study, whole genome resequencing of 30 individuals (Dorper, DB and Hu sheep, HY) identified 15,108,125 single nucleotide polymorphisms (SNPs). Population differentiation (Fst) and cross population extended haplotype homozygosity (XP-EHH) were performed for selective signal analysis. In total, 106 and 515 overlapped genes were present in both the Fst results and XP-EHH results in HY vs DB and in DB vs HY, respectively. In HY vs DB, 106 genes were enriched in 12 GO terms and 83 KEGG pathways, such as ATP binding (GO:0005524) and PI3K-Akt signaling pathway (oas04151). In DB vs HY, 515 genes were enriched in 109 GO terms and 215 KEGG pathways, such as skeletal muscle cell differentiation (GO:0035914) and MAPK signaling pathway (oas04010). According to the annotation results, we identified a series of candidate genes associated with reproduction (UNC5C, BMPR1B, and GLIS1), meat quality (MECOM, MEF2C, and MYF6), and immunity (GMDS, GALK1, and ITGB4). Our investigation has uncovered genomic information for important traits in sheep and provided a basis for subsequent studies of related traits

The Effect of EGR1 on the Proliferation of Dermal Papilla Cells

Early growth response factor 1 (EGR1) is a zinc-finger transcription factor that plays a vital role in the development of hair follicles. According to our previous studies, EGR1 is a transcriptional promoter of the bone morphogenetic protein 7 (BMP7), a candidate gene involved in the proliferation of dermal papilla cells. Since hair follicles are the basis of lambskin pattern formation and dermal papilla cells (DPCs) act on hair follicle growth, in order to elucidate the role of EGR1 and hair follicles, this study aimed to investigate the biological role of EGR1 in DPCs. In our study, the EGR1 coding sequence (CDS) region was firstly cloned by polymerase chain reaction, and bioinformatics analysis was performed. Then, the function of EGR1 was detected by 5-ethynyl-2&rsquo;-deoxyuridine (EDU) and Cell Counting Kit-8 (CCK8), and Western blot (WB) was conducted to analyze the cellular effect of EGR1 on DPCs. The proliferative effect of EGR1 on DPCs was also further confirmed by detecting its expression by qPCR and WB on marker genes of proliferation, including PCNA and CDK2. The sequence of the EGR1 CDS region of a lamb was successfully cloned, and its nucleic acid sequence was analyzed and found to be highly homologous to Rattus norvegicus, Mus musculus, Bos taurus and Homo sapiens. Predictive analysis of the protein encoded by EGR1 revealed that it is an extra-membrane protein, and not a secretory protein, with subcellular localization in the nucleus and cytoplasm. The proliferative effect of DPCs was significantly stronger (p &lt; 0.01) in EGR1 up-regulated DPCs compared to the controls, while the opposite result was observed in EGR1 down-regulated DPCs. Markers of proliferation including PCNA and CDK2 also appeared to be differentially upregulated in EGR1 gene overexpression compared to the controls, with the opposite result in EGR1 gene downregulation. In summary, our study revealed that EGR1 promotes the proliferation of DPCs, and we speculate that EGR1 may be closely associated with hair follicle growth and development

<i>BMP7</i> Functions to Regulate Proliferation of Dermal Papilla Cells in Hu Sheep

Bone morphogenetic proteins (BMPs) are the structurally similar and highly conserved type of functional proteins that play an important role in hair follicle growth and development. BMP7 was a differentially expressed gene in different patterns of Hu sheep lambskin identified using Agilent microarray. Since hair follicle is the basis of pattern formation of lambskin, and its growth and development is governed by dermal papilla cells (DPCs), to clarify the role of BMP7 and hair follicle, our study was designed to investigate the regulation between BMP7 and DPCs. Firstly, the CDS region of BMP7 was cloned by 3’Race and PCR in Hu sheep and performed serious of bioinformatic analysis. Then, the effects of BMP7 on DPCs were analyzed after overexpression and interference of BMP7 in dermal papilla cells by CCK8, EdU, and PI assay. Additionally, qPCR was also conducted to clarify the relationship between BMP7 and the TGF-β/Smad signaling pathway. A total of 1296 bp of the BMP7 CDS region sequence was sucessfully cloned in Hu sheep, encoding a signal peptide of 431 amino acids, molecular weight was 49,316.9 Da and the isoelectric point (Pi) was 7.75. Nucleotide sequencing analysis of BMP7 revealed that Hu sheep had high homology with Bos taurus, Homo sapiens, and Canis lupus familiaris. Structure domain prediction showed that TGF-β superfamily domain exist between 330th–431th amino acid, BMP7 protein is a secreted protein. In BMP7 up-regulated DPCs, DPCs proliferation rate and cell cycle were significantly higher than that of NC group (p Smad3, Smad4, Samd6, and TGF-β1 in TGF-β/Smad signaling pathway were significantly lower than that in NC group (p BMP7 down-regulated DPCs, it presented the opposite result. In conclusion, our study showed that BMP7 had a positive effect on DPCs by accelerating the proliferation and cell cycle of DPCs, and hypothesized that regulate hair follicles growth and development via TGF-β/Smad signaling pathway. These findings may provide a synergistic target for the subsequent research of hair follicle growth and development

MiR-22-3p Inhibits Proliferation and Promotes Differentiation of Skeletal Muscle Cells by Targeting IGFBP3 in Hu Sheep

The growth and development of skeletal muscle require a series of regulatory factors. MiRNA is a non-coding RNA with a length of about 22 nt, which can inhibit the expression of mRNA and plays an important role in the growth and development of muscle cells. The role of miR-22-3p in C2C12 cells and porcine skeletal muscle has been reported, but it has not been verified in Hu sheep skeletal muscle. Through qPCR, CCK-8, EdU and cell cycle studies, we found that overexpression of miR-22-3p inhibited proliferation of skeletal muscle cells (p &lt; 0.01). The results of qPCR and immunofluorescence showed that overexpression of miR-22-3p promoted differentiation of skeletal muscle cells (p &lt; 0.01), while the results of inhibiting the expression of miR-22-3p were the opposite. These results suggested that miR-22-3p functions in growth and development of sheep skeletal muscle cells. Bioinformatic analysis with mirDIP, miRTargets, and RNAhybrid software suggested IGFBP3 was the target of miR-22-3p, which was confirmed by dual-luciferase reporter system assay. IGFBP3 is highly expressed in sheep skeletal muscle cells. Overexpression of IGFBP3 was found to promote proliferation of skeletal muscle cells indicated by qPCR, CCK-8, EdU, and cell cycle studies (p &lt; 0.01). The results of qPCR and immunofluorescence experiments proved that overexpression of IGFBP3 inhibited differentiation of skeletal muscle cells (p &lt; 0.01), while the results of interfering IGFBP3 with siRNA were the opposite. These results indicate that miR-22-3p is involved in proliferation and differentiation of skeletal muscle cells by targeting IGFBP3