Characterization and comparative analyses of muscle transcriptomes in Dorper and small-tailed Han sheep using RNA-Seq technique

Integrative Meta-Analysis: Unveiling Genetic Factors in Meat Sheep Growth and Muscular Development through QTL and Transcriptome Studies

OBJECTIVE

The study aimed to investigate the effects of castration on performance, carcass characteristics, and meat quality in sheep, as well as explore the expression of key genes related to metabolic pathways and muscle growth following castration.

METHODS

A meta-analysis approach was utilized to analyze data from multiple studies to compare the performance, carcass characteristics, and meat quality of castrated sheep (wethers) with intact rams. Additionally, protein-protein interaction (PPI) networks, differential gene expression (DEG) interactions, Gene Ontology (GO) terms, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were examined to identify molecular mechanisms associated with fat metabolism and muscle development in sheep tails.

RESULTS

The analysis revealed that castrated sheep (wethers) exhibited improved average daily gain, increased tenderness, lower backfat thickness, and a tendency for greater loin muscle area compared to intact rams. This suggests that castration promotes faster growth and results in leaner carcasses with potentially higher muscle content. Furthermore, the identification of downregulated DEGs like ACLY, SLC27A2, and COL1A1 and upregulated DEGs such as HOXA9, PGM2L1, and ABAT provides insights into the molecular mechanisms underlying fat deposition and muscle development in sheep.

CONCLUSIONS

The findings support the practice of castration in sheep production as it enhances growth performance, leads to leaner carcasses with higher muscle content, and improves meat tenderness. The identified changes in gene expression offer valuable insights for further research into understanding the impact of castration on muscle development and fat metabolism in sheep. This meta-analysis contributes to the knowledge of molecular mechanisms involved in fat deposition in sheep, opening avenues for future investigations in livestock fat metabolism research.

DIA-Based Proteomic Analysis Reveals MYOZ2 as a Key Protein Affecting Muscle Growth and Development in Hybrid Sheep

Hybridization of livestock can be used to improve varieties, and different hybrid combinations produce unique breeding effects. In this study, male Southdown and Suffolk sheep were selected to hybridize with female Hu sheep to explore the effects of male parentage on muscle growth and the development of offspring. Using data-independent acquisition technology, we identified 119, 187, and 26 differentially abundant proteins (DAPs) between Hu × Hu (HH) versus Southdown × Hu (NH), HH versus Suffolk × Hu (SH), and NH versus SH crosses. Two DAPs, MYOZ2 and MYOM3, were common to the three hybrid groups and were mainly enriched in muscle growth and development-related pathways. At the myoblast proliferation stage, MYOZ2 expression decreased cell viability and inhibited proliferation. At the myoblast differentiation stage, MYOZ2 expression promoted myoblast fusion and enhanced the level of cell fusion. These findings provide new insights into the key proteins and metabolic pathways involved in the effect of male parentage on muscle growth and the development of hybrid offspring in sheep.

Gene Structure, Expression and Function Analysis of MEF2 in the Pacific White Shrimp Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is the most economically important crustacean in the world. The growth and development of shrimp muscle has always been the focus of attention. Myocyte Enhancer Factor 2 (MEF2), a member of MADS transcription factor, has an essential influence on various growth and development programs, including myogenesis. In this study, based on the genome and transcriptome data of L. vannamei, the gene structure and expression profiles of MEF2 were characterized. We found that the LvMEF2 was widely expressed in various tissues, mainly in the Oka organ, brain, intestine, heart, and muscle. Moreover, LvMEF2 has a large number of splice variants, and the main forms are the mutually exclusive exon and alternative 5′ splice site. The expression profiles of the LvMEF2 splice variants varied under different conditions. Interestingly, some splice variants have tissue or developmental expression specificity. After RNA interference into LvMEF2, the increment in the body length and weight decreased significantly and even caused death, suggesting that LvMEF2 can affect the growth and survival of L. vannamei. Transcriptome analysis showed that after LvMEF2 was knocked down, the protein synthesis and immune-related pathways were affected, and the associated muscle protein synthesis decreased, indicating that LvMEF2 affected muscle formation and the immune system. The results provide an important basis for future studies of the MEF2 gene and the mechanism of muscle growth and development in shrimp.

De Novo Transcriptome Dataset Generation of the Swamp Buffalo Brain and Non-Brain Tissues

The sequenced data availability opened new horizons related to buffalo genetic control of economic traits and genomic diversity. The visceral organs (brain, liver, etc.) significantly involved in energy metabolism, docility, or social interactions. We performed swamp buffalo transcriptomic profiling of 24 different tissues (brain and non-brain) to identify novel transcripts and analyzed the differentially expressed genes (DEGs) of brain vs. non-brain tissues with their functional annotation. We obtained 178.57 Gb clean transcriptomic data with GC contents 52.77%, reference genome alignment 95.36%, exonic coverage 88.49%. Totally, 26363 mRNAs transcripts including 5574 novel genes were obtained. Further, 7194 transcripts were detected as DEGs by comparing brain vs. non-brain tissues group, of which 3,999 were upregulated and 3,195 downregulated. These DEGs were functionally associated with cellular metabolic activities, signal transduction, cytoprotection, and structural and binding activities. The related functional pathways included cancer pathway, PI3k-Akt signaling, axon guidance, JAK-STAT signaling, basic cellular metabolism, thermogenesis, and oxidative phosphorylation. Our study provides an in-depth understanding of swamp buffalo transcriptomic data including DEGs potentially involved in basic cellular activities and development that helped to maintain their working capacity and social interaction with humans, and also, helpful to disclose the genetic architecture of different phenotypic traits and their gene expression regulation.

Co-expression of candidate genes regulating growth performance and carcass traits of Barki lambs in Egypt

Sheep are considered one of the main sources of animal protein in Egypt and the producers of sheep mutton eagers to find biological criteria for selecting fast-growing lambs that reach market weight early. Therefore, the present study aimed to find a link between the expression profile of selected candidate genes with growth performance and carcass traits of Barki lambs. Thirty-eight Barki lambs were kept and fed individually after weaning till 12 months of age and were divided into 3 groups according to growth performance (fast, intermediate, and slow-growing). Three samples were taken from different body tissues (eye muscle, liver, and fat tail) of each group, directly during slaughtering and stored at − 80 °C until RNA isolation. Real-time PCR was used to profile selected candidate genes (RPL7, CTP1, FABP4, ADIPOQ, and CAPN3) and GAPDH was used as a housekeeping gene. The results indicated that the final body weight was significantly (P ≤ 0.05) greater in the fast (49.9 kg) and intermediate (40.7 kg) compared to slow-growing animals (30.8 kg). The hot carcass weight was heavier (P ≤ 0.05) in the fast and intermediate-growing (24.57 and 19.07 kg) than slow-growing lambs (15.10 kg). The blood profiles of T3 and T4 hormones in addition to other parameters such as total protein, total lipids, and calcium level showed no clear variations among different experimental groups. At the molecular level, our data demonstrated upregulation of genes involved in protein biosynthesis (RPL7), fatty acid oxidation (CPT1), and lipolysis (FABP4) in the fast and intermediate-growing lambs in all studied tissues which facilitate protein accretion, energy expenditure, and fatty acid partitioning required for muscle building up. Moreover, the expression profile of the gene involved in muscle development (CAPN3) was increased in fast and intermediate-growing compared to slow-growing lambs in order to support muscle proper development. On the other hand, a candidate gene involved in lipogenesis (ADIPOQ) was expressed similarly in fat and liver tissues; however, its expression was increased in muscles of fast and intermediate-growing lambs compared to slow-growing animals. In conclusion, the current study indicated that the expression profile of genes involved in metabolic activities of liver, muscle, and adipose tissue is linked with the growth performance of lambs although no variations were detected in blood parameters. This provides an evidence for the importance of co-expression of these genes in body tissues to determine the final body weight and carcass characteristics of Barki sheep.

Comparative Transcriptome Analyses of Gayal (Bos frontalis), Yak (Bos grunniens), and Cattle (Bos taurus) Reveal the High-Altitude Adaptation

Gayal and yak are well adapted to their local high-altitude environments, yet the transcriptional regulation difference of the plateau environment among them remains obscure. Herein, cross-tissue and cross-species comparative transcriptome analyses were performed for the six hypoxia-sensitive tissues from gayal, yak, and cattle. Gene expression profiles for all single-copy orthologous genes showed tissue-specific expression patterns. By differential expression analysis, we identified 3,020 and 1,995 differentially expressed genes (DEGs) in at least one tissue of gayal vs. cattle and yak vs. cattle, respectively. Notably, we found that the adaptability of the gayal to the alpine canyon environment is highly similar to the yak living in the Qinghai-Tibet Plateau, such as promoting red blood cell development, angiogenesis, reducing blood coagulation, immune system activation, and energy metabolism shifts from fatty acid β-oxidation to glycolysis. By further analyzing the common and unique DEGs in the six tissues, we also found that numerous expressed regulatory genes related to these functions are unique in the gayal and yak, which may play important roles in adapting to the corresponding high-altitude environment. Combined with WGCNA analysis, we found that UQCRC1 and COX5A are the shared differentially expressed hub genes related to the energy supply of myocardial contraction in the heart-related modules of gayal and yak, and CAPS is a shared differential hub gene among the hub genes of the lung-related module, which is related to pulmonary artery smooth muscle contraction. Additionally, EDN3 is the unique differentially expressed hub gene related to the tracheal epithelium and pulmonary vasoconstriction in the lung of gayal. CHRM2 is a unique differentially expressed hub gene that was identified in the heart of yak, which has an important role in the autonomous regulation of the heart. These results provide a basis for further understanding the complex transcriptome expression pattern and the regulatory mechanism of high-altitude domestication of gayal and yak.

Exploring the effects of different male parent crossings on sheep muscles and related regulatory genes using mRNA-Seq

Objective

With improvements in living standards and increase in population globally, the demand for meat products has been increasing; improved meat production from livestock could effectively meet this demand. In this study, we examined the differences in the muscle traits of different male crossbred sheep and attempted to identify key genes that regulate these traits.

Methods

Dubo sheep × Small-tailed Han sheep (DP×STH) and Suffolk × Small-tailed Han sheep (SFK×STH) were selected to determine meat quality and production performance by Masson staining. Transcriptome sequencing and bioinformatic analysis were performed to identify differentially expressed genes (DEGs) related to meat quality. The presence of DEGs was confirmed by real-time polymerase chain reaction (qPCR).

Results

The production performance of SFK×STH sheep was better than that of DP×STH sheep, but the meat quality of DP×STH sheep was better than that of SFK×STH sheep. The muscle fiber diameter of DP×STH sheep was smaller than that of SFK×STH sheep. Twenty-two DEGs were identified. Among them, four Gene Ontology terms were related to muscle traits, and three DEGs were related to muscle or muscle fibers. There were no significant differences in the number of single nucleotide mutations and mutation sites in the different male parent cross combinations.

Conclusion

This study provides genetic resources for future sheep muscle development and cross-breeding research.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Characterization and Comparative Transcriptomic Analysis of Skeletal Muscle in Pekin Duck at Different Growth Stages Using RNA-Seq

Simple Summary

Skeletal muscle is an important tissue and its development is strictly regulated by genes. In this study, in order to understand the muscle-related gene expression in Pekin duck, RNA-seq was performed to analyze and compare skeletal muscle at different growth stages. Alternative splicing, single nucleotide polymorphisms and insertion–deletions were detected, and 299 novel genes were discovered. MYL4, IGF2BP1, CSRP3, SPP1, KLHL31, LAMB2, LAMA2, ITGB1 and OPN played crucial roles in skeletal muscle development. Oxidative phosphorylation, ECM-receptor interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids participated in the regulation of skeletal muscle development in Pekin duck. This study provides an important reference for revealing the developmental mechanisms of pectoral and leg muscles in duck.

Abstract

Skeletal muscle, accounting for approximately 50% of body weight, is the largest and most important tissue. In this study, the gene expression profiles and pathways in skeletal muscle of Pekin duck were investigated and compared at embryonic day 17, 21, and 27 and postnatally at 6 months of age. An average of 49,555,936 reads in each sample was obtained from the transcriptome libraries. Over 70.0% of alternative splicing (AS) in each sample was mainly alternative 5′ first exon (transcription start site)—the first exon splicing (TSS) and alternative 3′ last exon (transcription terminal site)—the last exon splicing (TTS), indicating that TSS and TTS were the most common AS event in Pekin ducks, and these AS events were closely related to the regulation of muscle development at different growth stages. The results provided a valuable genomic resource for selective breeding and functional studies of genes. A total of 299 novel genes with ≥2 exons were obtained. There were 294 to 2806 differentially expressed genes (DEGs) in each pairwise comparison of Pekin duck. Notably, 90 DEGs in breast muscle and 9 DEGs in leg muscle were co-expressed at all developmental points. DEGs were validated by qPCR analysis, which confirmed the tendency of the expression. DEGs related to muscle development were involved in biological processes such as “endodermal cell differentiation”, “muscle cell cellular homeostasis”, “skeletal muscle tissue growth” and “skeletal muscle cell differentiation”, and were involved in pathways such as oxidative phosphorylation, ECM-receptor (extracellular matrix receptor) interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids. Some DEGs, including MYL4, IGF2BP1, CSRP3, SPP1 and KLHL31, as well as LAMB2, LAMA2, ITGB1 and OPN, played crucial roles in muscle growth and development. This study provides valuable information about the expression profile of mRNAs and pathways from duck skeletal muscle at different growth stages, and further functional study of these mRNAs and pathways could provide new ideas for studying the molecular networks of growth and development in duck skeletal muscle.

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.

Transcriptome signature of two lactation stages in Ghezel sheep identifies using RNA-Sequencing

The expression of genes and their regulation during lactation in Ghezel sheep breed remains less understood. To explore the underlying molecular mechanism of the lactation process in the mammary gland, transcriptome profiles of Iranian fat-tailed Ghezel sheep breed milk at two stages, before (BF) and after peak (AF) stages of lactation were investigated. Functional impacts of differentially expressed genes (DEGs) between BF and AF stages were surveyed using Gene Ontology (GO) and Protein-Protein Interaction (PPI) network analysis. Totally, 75 DEGs were identified between BF and AF stages of lactation. The RNA-Seq results were validated by Q-RT-PCR. Gene ontology of DEGs mainly enriched in metabolic process and oxidative phosphorylation. PPI network analysis also highlighted the contribution of peroxisome proliferator-activated receptors (PPAR) signaling, oxidative phosphorylation and metabolic pathways in the lactation process. Intriguingly, the genes involved in fat metabolism dominantly down-regulated at AF stage. Our results provide new insight into transcriptional changes and add to growing body of knowledge on the lactation process in fat-tailed sheep breeds.

Identification of Differentially Expressed Genes in Different Types of Broiler Skeletal Muscle Fibers Using the RNA-seq Technique

The difference in muscle fiber types is very important to the muscle development and meat quality of broilers. At present, the molecular regulation mechanisms of skeletal muscle fiber-type transformation in broilers are still unclear. In this study, differentially expressed genes between breast and leg muscles in broilers were analyzed using RNA-seq. A total of 767 DEGs were identified. Compared with leg muscle, there were 429 upregulated genes and 338 downregulated genes in breast muscle. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in cellular processes, single organism processes, cells, and cellular components, as well as binding and catalytic activity. KEGG analysis shows that a total of 230 DEGs were mapped to 126 KEGG pathways and significantly enriched in the four pathways of glycolysis/gluconeogenesis, starch and sucrose metabolism, insulin signalling pathways, and the biosynthesis of amino acids. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) was used to verify the differential expression of 7 selected DEGs, and the results were consistent with RNA-seq data. In addition, the expression profile of MyHC isoforms in chicken skeletal muscle cells showed that with the extension of differentiation time, the expression of fast fiber subunits (types IIA and IIB) gradually increased, while slow muscle fiber subunits (type I) showed a downward trend after 4 days of differentiation. The differential genes screened in this study will provide some new ideas for further understanding the molecular mechanism of skeletal muscle fiber transformation in broilers.

Magnaporthe oryzae Auxiliary Activity Protein MoAa91 Functions as Chitin-Binding Protein To Induce Appressorium Formation on Artificial Inductive Surfaces and Suppress Plant Immunity

The appressoria that are generated by the rice blast fungus Magnaporthe oryzae in response to surface cues are important for successful colonization. Previous work showed that regulators of G-protein signaling (RGS) and RGS-like proteins play critical roles in appressorium formation. However, the mechanisms by which these proteins orchestrate surface recognition for appressorium induction remain unclear. Here, we performed comparative transcriptomic studies of ΔMorgs mutant and wild-type strains and found that M. oryzae Aa91 (MoAa91), a homolog of the auxiliary activity family 9 protein (Aa9), was required for surface recognition of M. oryzae We found that MoAA91 was regulated by the MoMsn2 transcription factor and that its disruption resulted in defects in both appressorium formation on the artificial inductive surface and full virulence of the pathogen. We further showed that MoAa91 was secreted into the apoplast space and was capable of competing with the immune receptor chitin elicitor-binding protein precursor (CEBiP) for chitin binding, thereby suppressing chitin-induced plant immune responses. In summary, we have found that MoAa91 is a novel signaling molecule regulated by RGS and RGS-like proteins and that MoAa91 not only governs appressorium development and virulence but also functions as an effector to suppress host immunity.IMPORTANCE The rice blast fungus Magnaporthe oryzae generates infection structure appressoria in response to surface cues largely due to functions of signaling molecules, including G-proteins, regulators of G-protein signaling (RGS), mitogen-activated protein (MAP) kinase pathways, cAMP signaling, and TOR signaling pathways. M. oryzae encodes eight RGS and RGS-like proteins (MoRgs1 to MoRgs8), and MoRgs1, MoRgs3, MoRgs4, and MoRgs7 were found to be particularly important in appressorium development. To explore the mechanisms by which these proteins regulate appressorium development, we have performed a comparative in planta transcriptomic study and identified an auxiliary activity family 9 protein (Aa9) homolog that we named MoAa91. We showed that MoAa91 was secreted from appressoria and that the recombinant MoAa91 could compete with a chitin elicitor-binding protein precursor (CEBiP) for chitin binding, thereby suppressing chitin-induced plant immunity. By identifying MoAa91 as a novel signaling molecule functioning in appressorium development and an effector in suppressing host immunity, our studies revealed a novel mechanism by which RGS and RGS-like proteins regulate pathogen-host interactions.

Comparative Transcriptome Analysis Identifying the Different Molecular Genetic Markers Related to Production Performance and Meat Quality in Longissimus Dorsi Tissues of MG × STH and STH Sheep

Crossbred sheep have many prominent traits, such as excellent production performance and high-quality meat, when compared to local sheep breeds. However, the genetic molecular markers related to these characteristics remain unclear. The crossbred MG × STH (small-tailed Han sheep (STH) × Mongolian sheep (MG)) breed and the STH breed were selected to measure production performance and meat quality. We used 14 indexes of production performance and meat quality, which in the MG × STH population showed significant differences compared to the STH breed. Subsequently, the longissimusdorsi from the two sheep were subjected to comparative transcriptomic analyses to identify differentially expressed genes (DEGs) related to production performance and meat quality. A total of 874 DEGs were identified between the two sheep groups. A total of 110 unique DEGs related to sheep production performance and meat quality were selected as the candidate DEGs. We found 6 production-performance-related and 30 meat-quality-related DEGs through a correlation analysis, including SPARC, ACVRL1, FNDC5 and FREM1. The expression levels of 11 DEGs were validated by real-time PCR, and the results were in accordance with the results of the comparative transcriptomic and correlation analyses. These results will assist in understanding sheep heterosis and molecular marker-assisted selection.

miRNA transcriptome comparison between muscle and adipose tissues indicates potential miRNAs associated with intramuscular fat in Chinese swamp buffalo

The amount of intramuscular fat (IMF) affects the tenderness and juiciness of beef and is an important indicator of beef quality. A few miRNAs involved in IMF deposition have been identified in other livestock. However, in the buffalo, the association between miRNA and IMF has not been reported and the miRNA expression profile remains poorly understood. In this study, small RNA sequencing was performed to characterize the miRNA expression pattern in muscle and adipose tissues using the Illumina platform. A total of 108 differentially expressed (DE) miRNAs were identified, including 98 known miRNAs and 10 novel miRNAs. A qRT-PCR experiment confirmed the quality of the DE analysis. Eight DE miRNAs showed high expression in adipose tissue and a considerable expression level in muscle tissue. Functional enrichment indicated that bta-miR-148a, bta-miR-143, bta-miR-10b, bta-let-7i, bta-let-7f, bta-let-7b, bta-miR-30a-5p, and bta-miR-100 were significantly associated with adipogenesis, suggesting these as candidate regulators for IMF deposition in buffalo. However, further functional validation is required. This is the first characterization of the miRNA expression profile in the muscle and adipose tissues of buffalo. These results provide information for the identification of miRNAs with potential effects on IMF deposition in buffalo.

Identification and classification of differentially expressed genes in pyrethroid-resistant Culex pipiens pallens

Culex pipiens pallens is an important vector that transmits Bancroftian filariasis, Japanese encephalitis and other diseases that pose a serious threat to human health. Extensive and improper use of insecticides has caused insecticide resistance in mosquitoes, which has become an important obstacle to the control of mosquito-borne diseases. It is crucial to investigate the underlying mechanism of insecticide resistance. The aims of this study were to identify genes involved in insecticide resistance based on the resistance phenotype, gene expression profile and single-nucleotide polymorphisms (SNPs) and to screen for major genes controlling insecticide resistance. Using a combination of SNP and transcriptome data, gene expression quantitative trait loci (eQTLs) were studied in deltamethrin-resistant mosquitoes. The most differentially expressed pathway in the resistant group was identified, and a regulatory network was built using these SNPs and the differentially expressed genes (DEGs) in this pathway. The major candidate genes involved in the control of insecticide resistance were analyzed by qPCR, siRNA microinjection and CDC bottle bioassays. A total of 85 DEGs that encoded putative detoxification enzymes (including 61 P450s) were identified in this pathway. The resistance regulatory network was built using SNPs, and these metabolic genes, and a major gene, CYP9AL1, were identified. The functional role of CYP9AL1 in insecticide resistance was confirmed by siRNA microinjection and CDC bottle bioassays. Using the eQTL approach, we identified important genes in pyrethroid resistance that may aid in understanding the mechanism underlying insecticide resistance and in targeting new measures for resistance monitoring and management.

Comparative Transcriptome Profiling of mRNA and lncRNA Related to Tail Adipose Tissues of Sheep

The Lanzhou Fat-Tail sheep (LFTS, long fat-tailed sheep) is an endangered sheep breed in China with a fat tail compared to the traditional local varieties, Small Tail Han sheep (STHS, thin-tailed sheep) with a small tail, and Tibetan sheep (TS, short thin-tailed sheep) with a little tail. However, little is known regarding how tail fat deposition is regulated by long noncoding RNA (lncRNA). To evaluate the lncRNA and mRNA associated with tail fat deposition and development among these breeds, high-throughput RNA sequencing of three individuals each of LFTS, STHS, and TS were performed and analyzed in this study. RNA sequencing data from these three groups revealed 10 differentially expressed genes (DEGs) and 37 differentially expressed lncRNAs between the LFTS and STHS groups, 390 DEGs and 59 differentially expressed lncRNAs between the LFTS and TS groups, and 80 DEGs and 16 differentially expressed lncRNAs between the STHS and TS groups (p-value < 0.05 and fold change ≥ 2), respectively. Gene Ontology and pathway analysis of DEGs and target genes of differentially expressed lncRNAs revealed enrichment in fatty acid metabolism and fatty acid elongation-related pathways that contribute to fat deposition. Subsequently, the expression of 14 DEGs and 6 differentially expressed lncRNAs was validated by quantitative real-time PCR. Finally, two co-expression networks of differentially expressed mRNA and lncRNAs were constructed. The results suggested that some differentially expressed lncRNAs (TCONS_00372767, TCONS_00171926, TCONS_00054953, and TCONS_00373007) may play crucial roles as core lncRNAs in tail fat deposition processes. In summary, the present study extends the sheep tail fat lncRNA database and these differentially expressed mRNA and lncRNAs may provide novel candidate regulators for future genetic and molecular studies on tail fat deposition of sheep.

Comparison of transcriptomic landscapes of different lamb muscles using RNA-Seq

Transcriptome deep sequencing is a powerful tool for exploring the genetic architecture of complex traits. Gene expression patterns may explain a high degree of the observed phenotypic differences in histochemical and metabolic parameters related to meat quality among different muscles. In this study, we sequenced by RNA-Seq the whole transcriptome of nine lamb muscles: Semimembranosus (SM), Semitendinosus (ST), Cranial gluteobiceps, Gluteus medius (GM), Rectus femoris, Supraspinatus (SS), Longissimus lumborum (LL), Adductor and Psoas major. Significant gene expression differences were detected between almost all pairwise comparisons, being more pronounced between SS and ST, SM and LL, and ST and GM. These differences can be explained in terms of ATPase and glycolytic activities, muscle fiber typing and oxidative score, clustering muscles as fast glycolytic, intermediate or slow oxidative. ST showed up-regulation of gene pathways related to carbohydrate metabolism, energy generation and protein turnover as expected from a fast white muscle. SS showed myosin isoforms typical of slow muscles and high expression of genes related to calcium homeostasis and vascularization. SM, LL and GM showed in general intermediate gene expression patterns. Several novel transcripts were detected, mostly related to muscle contraction and structure, oxidative metabolism, lipid metabolism and protein phosphorylation. Expression profiles were consistent with previous histochemical and metabolic characterization of these muscles. Up-regulation of ion transport genes may account for significant differences in water holding capacity. High expression of genes related to cell adhesion, cytoskeleton organization, extracellular matrix components and protein phosphorylation may be related to meat yellowness and lower tenderness scores. Differential expression of genes related to glycolytic activity and lactic acid generation among fast, intermediate and slow muscles may explain the detected final meat pH differences. These results reveal new candidate genes associated with lamb meat quality, and give a deeper insight into the genetic architecture of these complex traits.

Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds

The goal of the current study is to analyse the gene expression profile of the ovine skeletal muscle as well as to characterize the genetic variation of transcripts expressed in such tissue. This aim has been achieved by sequencing the longissimus dorsi transcriptomes of 50 sheep distributed in five pools representing the Canaria de Pelo, Roja Mallorquina, Gallega, Xisqueta and Ripollesa Spanish autochthonous breeds. Approximately, 363 million reads per pool have been produced and 71.9–82.9% have been successfully mapped to the ovine genome in a paired-end mode (2 × 75 bp). The 200 most expressed muscle transcripts (≈1% of the total transcript count) account for 51% (Canaria de Pelo) to 67% (Gallega) of the total ovine skeletal muscle mRNA expression. These highly expressed genes play key roles in pathways related with striated muscle contraction, gluconeogenesis, glycolysis, citric acid cycle and respiratory electron transport. RNA-Sequencing of muscle transcripts has also revealed that ~72% of the SNPs detected with this approach are shared by at least two pools, and 10% of them segregate in the five pools under analysis. Most of the substitutions detected by RNA-Seq are synonymous or missense and only a minority are predicted to have consequences on protein function.

Beta-hydroxybutyrate infusion identifies acutely differentially expressed genes related to metabolism and reproduction in the hypothalamus and pituitary of castrated male sheep

To identify molecular pathways that couple metabolic imbalances and reproduction, we randomly assigned 10 castrated male sheep to be centrally injected into the lateral ventricle through intracerebroventricular cannulas with 1 ml of β-hydroxybutyric acid sodium salt solution (BHB; 12,800 µmol/l) or saline solution (CON; 0.9% NaCl). Approximately 2 h postinjection, sheep were humanely euthanized, and hypothalamus and pituitary tissues were harvested for transcriptome characterization by RNA sequencing. RNA was extracted from the hypothalamus and pituitary and sequenced at a high depth (hypothalamus: 468,912,732 reads; pituitary: 515,106,092 reads) with the Illumina Hi-Seq 2500 platform and aligned to Bos taurus and Ovis aries genomes. Of the total raw reads, 87% (hypothalamus) and 90.5% (pituitary) mapped to the reference O. aries genome. Within these read sets, ~56% in hypothalamus and 69% in pituitary mapped to either known or putative protein coding genes. Fragments per kilobase of transcripts per million normalized counts were averaged and ranked to identify the transcript expression level. Gene Ontology analysis (DAVID Bioinformatics Resources) was utilized to identify biological process functions related to genes shared between tissues, as well as functional categories with tissue-specific enrichment. Between CON- and BHB-treated sheep, 11 and 44 genes were differentially expressed (adj. P < 0.05) within the pituitary and hypothalamus, respectively. Functional enrichment analyses revealed BHB altered expression of genes in pathways related to stimulus perception, inflammation, and cell cycle control. The set of genes altered by BHB creates a foundation from which to identify the signaling pathways that impact reproduction during metabolic imbalances.

Identification and characterization of microRNA in the lung tissue of pigs with different susceptibilities to PCV2 infection

Porcine circovirus type 2 (PCV2) is the primary cause of post-weaning multisystemic wasting syndrome (PMWS) and other PCV-associated diseases. According to our previous RNA-sequencing analysis, the differences in the susceptibility to PCV2 infection depended on the genetic differences between the Laiwu (LW) and Yorkshire × Landrace crossbred (YL) pigs, but the cellular microRNA (miRNA) that are differentially expressed between the LW and YL pigs before and after PCV2 infection remain to be determined. In this study, high-throughput sequencing was performed to determine the abundance and differential expression of miRNA in lung tissues from PCV2-infected and PCV2-uninfected LW and YL pigs. In total, 295 known and 95 novel miRNA were identified, and 23 known and 25 novel miRNA were significantly differentially expressed in the PCV2-infected vs. PCV2-uninfected LW pigs and/or the PCV2-infected vs. PCV2-uninfected YL pigs. The expression levels of ssc-miR-122, ssc-miR-192, ssc-miR-451, ssc-miR-486, and ssc-miR-504 were confirmed by quantitative real-time PCR (qRT-PCR). Analysis of the potential targets of the four up-regulated miRNA (i.e., ssc-miR-122, ssc-miR-192, ssc-miR-451 and ssc-miR-486) identified pathways and genes that may be important for disease resistance. Among the up-regulated miRNA, ssc-miR-122 can repress the protein expression and viral DNA replication of PCV2 and down-regulate the expression of the nuclear factor of activated T-cells 5 (NFAT5) and aminopeptidase puromycin sensitive (NPEPPS) by binding to their 3′ untranslated region (3′UTR) in PK15 cells. Therefore, ssc-miR-122 may indirectly suppress PCV2 infection by targeting genes related to the host immune system, such as NFAT5 and NPEPPS.

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.

Transcriptome Analysis of Three Sheep Intestinal Regions reveals Key Pathways and Hub Regulatory Genes of Large Intestinal Lipid Metabolism

The large intestine, also known as the hindgut, is an important part of the animal digestive system. Recent studies on digestive system development in ruminants have focused on the rumen and the small intestine, but the molecular mechanisms underlying sheep large intestine metabolism remain poorly understood. To identify genes related to intestinal metabolism and to reveal molecular regulation mechanisms, we sequenced and compared the transcriptomes of mucosal epithelial tissues among the cecum, proximal colon and duodenum. A total of 4,221 transcripts from 3,254 genes were identified as differentially expressed transcripts. Between the large intestine and duodenum, differentially expressed transcripts were found to be significantly enriched in 6 metabolism-related pathways, among which PPAR signaling was identified as a key pathway. Three genes, CPT1A, LPL and PCK1, were identified as higher expression hub genes in the large intestine. Between the cecum and colon, differentially expressed transcripts were significantly enriched in 5 lipid metabolism related pathways, and CEPT1 and MBOAT1 were identified as hub genes. This study provides important information regarding the molecular mechanisms of intestinal metabolism in sheep and may provide a basis for further study.

Single nucleotide polymorphism discovery in bovine liver using RNA-seq technology

BACKGROUND

RNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits.

RESULTS

The RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds.

CONCLUSIONS

This study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.

Comparative Transcriptome Analysis Reveals that a Ubiquitin-Mediated Proteolysis Pathway Is Important for Primary and Secondary Hair Follicle Development in Cashmere Goats

Background

The fleece of cashmere goats contains two distinct populations of fibers, a short and fine non-medullated insulating cashmere fiber and a long and coarse medullated guard hair. The former is produced by secondary follicles (SFs) and the later by primary follicles (PFs). Evidence suggests that the induction of PFs and SFs may require different signaling pathways. The regulation of BMP2/4 signaling by noggin and Edar signaling via Downless genes are essential for the induction of SFs and PFs, respectively. However, these differently expressed genes of the signaling pathway cannot directly distinguish between the PFs and SFs.

Results

In this study, we selected RNA samples from 11 PFs and 7 SFs that included 145,525 exons. The pathway analysis of 4512 differentially expressed exons revealed that the most statistically significant metabolic pathway was related to the ubiquitin–mediated proteolysis pathway (UMPP) (P<3.32x 10⁻⁷). In addition, the 51 exons of the UMPP that were differentially expressed between the different types of hair follicle (HFs) were compared by cluster analysis. This resulted in the PFs and SFs being divided into two classes. The expression level of two selected exons was analyzed by qRT-PCR, and the results indicated that the expression patterns were consistent with the deep sequencing results obtained by RNA-Seq.

Conclusions

Based on the comparative transcriptome analysis of 18 HFs from cashmere goats, a large number of differentially expressed exons were identified using a high-throughput sequencing approach. This study suggests that UMPP activation is a prominent signaling pathway for distinguishing the PFs and SFs of cashmere goats. It is also a meaningful contribution to the theoretical basis of the biological study of the HFs of cashmere goats and other mammals.

Generation of biallelic knock-out sheep via gene-editing and somatic cell nuclear transfer

Transgenic sheep can be used to achieve genetic improvements in breeds and as an important large-animal model for biomedical research. In this study, we generated a TALEN plasmid specific for ovine MSTN and transfected it into fetal fibroblast cells of STH sheep. MSTN biallelic-KO somatic cells were selected as nuclear donor cells for SCNT. In total, cloned embryos were transferred into 37 recipient gilts, 28 (75.7%) becoming pregnant and 15 delivering, resulting in 23 lambs, 12 of which were alive. Mutations in the lambs were verified via sequencing and T7EI assay, and the gene mutation site was consistent with that in the donor cells. Off-target analysis was performed, and no off-target mutations were detected. MSTN KO affected the mRNA expression of MSTN relative genes. The growth curve for the resulting sheep suggested that MSTN KO caused a remarkable increase in body weight compared with those of wild-type sheep. Histological analyses revealed that MSTN KO resulted in muscle fiber hypertrophy. These findings demonstrate the successful generation of MSTN biallelic-KO STH sheep via gene editing in somatic cells using TALEN technology and SCNT. These MSTN mutant sheep developed and grew normally, and exhibited increased body weight and muscle growth.

Comparative transcriptome profiling of longissimus muscle tissues from Qianhua Mutton Merino and Small Tail Han sheep

The Qianhua Mutton Merino (QHMM) is a new sheep (Ovis aries) variety with better meat performance compared with the traditional local variety Small Tail Han (STH) sheep. We aimed to evaluate the transcriptome regulators associated with muscle growth and development between the QHMM and STH. We used RNA-Seq to obtain the transcriptome profiles of the longissimus muscle from the QHMM and STH. The results showed that 960 genes were differentially expressed (405 were up-regulated and 555 were down-regulated). Among these, 463 differently expressed genes (DEGs) were probably associated with muscle growth and development and were involved in biological processes such as skeletal muscle tissue development and muscle cell differentiation; molecular functions such as catalytic activity and oxidoreductase activity; cellular components such as mitochondrion and sarcoplasmic reticulum; and pathways such as metabolic pathways and citrate cycle. From the potential genes, a gene-act-network and co-expression-network closely related to muscle growth and development were identified and established. Finally, the expressions of nine genes were validated by real-time PCR. The results suggested that some DEGs, including MRFs, GXP1 and STAC3, play crucial roles in muscle growth and development processes. This genome-wide transcriptome analysis of QHMM and STH muscle is reported for the first time.

The developmental transcriptome landscape of bovine skeletal muscle defined by Ribo-Zero ribonucleic acid sequencing

Ribonucleic acid sequencing (RNA-Seq) libraries are normally prepared with oligo(dT) selection of poly(A)+ mRNA, but it depends on intact total RNA samples. Recent studies have described Ribo-Zero technology, a novel method that can capture both poly(A)+ and poly(A)- transcripts from intact or fragmented RNA samples. We report here the first application of Ribo-Zero RNA-Seq for the analysis of the bovine embryonic, neonatal, and adult skeletal muscle whole transcriptome at an unprecedented depth. Overall, 19,893 genes were found to be expressed, with a high correlation of expression levels between the calf and the adult. Hundreds of genes were found to be highly expressed in the embryo and decreased at least 10-fold after birth, indicating their potential roles in embryonic muscle development. In addition, we present for the first time the analysis of global transcript isoform discovery in bovine skeletal muscle and identified 36,694 transcript isoforms. Transcriptomic data were also analyzed to unravel sequence variations; 185,036 putative SNP and 12,428 putative short insertions-deletions (InDel) were detected. Specifically, many stop-gain, stop-loss, and frameshift mutations were identified that probably change the relative protein production and sequentially affect the gene function. Notably, the numbers of stage-specific transcripts, alternative splicing events, SNP, and InDel were greater in the embryo than in the calf and the adult, suggesting that gene expression is most active in the embryo. The resulting view of the transcriptome at a single-base resolution greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during bovine skeletal muscle development.

Identification and Classification of New Transcripts in Dorper and Small-Tailed Han Sheep Skeletal Muscle Transcriptomes

High-throughput mRNA sequencing enables the discovery of new transcripts and additional parts of incompletely annotated transcripts. Compared with the human and cow genomes, the reference annotation level of the sheep genome is still low. An investigation of new transcripts in sheep skeletal muscle will improve our understanding of muscle development. Therefore, applying high-throughput sequencing, two cDNA libraries from the biceps brachii of small-tailed Han sheep and Dorper sheep were constructed, and whole-transcriptome analysis was performed to determine the unknown transcript catalogue of this tissue. In this study, 40,129 transcripts were finally mapped to the sheep genome. Among them, 3,467 transcripts were determined to be unannotated in the current reference sheep genome and were defined as new transcripts. Based on protein-coding capacity prediction and comparative analysis of sequence similarity, 246 transcripts were classified as portions of unannotated genes or incompletely annotated genes. Another 1,520 transcripts were predicted with high confidence to be long non-coding RNAs. Our analysis also revealed 334 new transcripts that displayed specific expression in ruminants and uncovered a number of new transcripts without intergenus homology but with specific expression in sheep skeletal muscle. The results confirmed a complex transcript pattern of coding and non-coding RNA in sheep skeletal muscle. This study provided important information concerning the sheep genome and transcriptome annotation, which could provide a basis for further study.

Characterization and comparative analyses of transcriptomes for in vivo and in vitro produced peri-implantation conceptuses and endometria from sheep

An increasing number of reports indicate that in vitro fertilization (IVF) is highly associated with long‑term side effects on embryonic and postnatal development, and can sometimes result in embryonic implant failure. While high‑throughput gene expression analysis has been used to explore the mechanisms underlying IVF-induced side effects on embryonic development, little is known about the effects of IVF on conceptus–endometrial interactions during the peri-implantation period. Using sheep as a model, we performed a comparative transcriptome analysis between in vivo (IVO; in vivo fertilized followed by further development in the uterus) and in vitro produced (IVP; IVF with further culture in the incubator) conceptuses, and the caruncular and intercaruncular areas of the ovine endometrium. We identified several genes that were differentially expressed between the IVO and IVP groups on day 17, when adhesion between the trophoblast and the uterine luminal epithelium begins in sheep. By performing Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, we found that, in the conceptus, differentially expressed genes (DEGs) were associated mainly with functions relating to cell binding and the cell cycle. In the endometrial caruncular area, DEGs were involved in cell adhesion/migration and apoptosis, and in the intercaruncular area, they were significantly enriched in pathways of signal transduction and transport. Thus, these DEGs are potential candidates for further exploring the mechanism underlying IVF/IVP-induced embryonic implant failure that occurs due to a loss of interaction between the conceptus and endometrium during the peri-implantation period.

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.

Characterization and Comparative Analysis of the Milk Transcriptome in Two Dairy Sheep Breeds using RNA Sequencing

This study presents a dynamic characterization of the sheep milk transcriptome aiming at achieving a better understanding of the sheep lactating mammary gland. Transcriptome sequencing (RNA-seq) was performed on total RNA extracted from milk somatic cells from ewes on days 10, 50, 120 and 150 after lambing. The experiment was performed in Spanish Churra and Assaf breeds, which differ in their milk production traits. Nearly 67% of the annotated genes in the reference genome (Oar_v3.1) were expressed in ovine milk somatic cells. For the two breeds and across the four lactation stages studied, the most highly expressed genes encoded caseins and whey proteins. We detected 573 differentially expressed genes (DEGs) across lactation points, with the largest differences being found, between day 10 and day 150. Upregulated GO terms at late lactation stages were linked mainly to developmental processes linked to extracellular matrix remodeling. A total of 256 annotated DEGs were detected in the Assaf and Churra comparison. Some genes selectively upregulated in the Churra breed grouped under the endopeptidase and channel activity GO terms. These genes could be related to the higher cheese yield of this breed. Overall, this study provides the first integrated overview on sheep milk gene expression.

Genome-wide mRNA-seq profiling reveals predominant down-regulation of lipid metabolic processes in adipose tissues of Small Tail Han than Dorset sheep

Small Tail Han and Dorset sheep are two different sheep with distinguished morphologies in fat depositions. In order to characterize their gene expression profiles, our present study took the advantages of RNA sequencing technology with the aims to identify important genes regulating the metabolisms in adipose tissues of two different sheep. In obtained high quality sequencing reads, 85.9 (Han) and 86.1% (Dorset) were uniquely aligned to Oar v3.1 sheep reference genome, and over 76% of bases in mapped reads corresponded to mRNA. Using R package EBSeq, we identified 602 differentially expressed genes. Using the 602 genes, GO analysis showed that 30 out of 56 significantly enriched biological processes were metabolism related, of which the most significant one was triglyceride biosynthetic process. The KEGG pathway analysis indicated the down-regulation of several fat metabolic pathways. The predominant down-regulation of massive metabolic processes, particularly the lipid metabolism, in adipose tissues of Han sheep could explain, at least in part, the distinguished fat deposition between two different sheep, and our data constitute a basic picture of transcriptomes in these sheep for better understanding of underline biological mechanism in their lipid metabolisms.

Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq

The widespread and improper use of pyrethroid insecticides, such as deltamethrin, has resulted in the evolution of resistance in many mosquito species, including Culex pipiens pallens. With the development of high-throughput sequencing, it is possible to massively screen pyrethroid resistance-associated gene. In this study, we used Illumina-Solexa transcriptome sequencing to identify genes that are expressed differently in deltamethrin-susceptible and -resistant strains of Culex pipiens pallens as a critical knowledge base for further studies. A total of 4,961,197,620 base pairs and 55,124,418 reads were sequenced, mapped to the Culex quinquefasciatus genome and assembled into 17,679 known genes. We recorded 1826 significantly differentially expressed genes (DEGs). Among them, 1078 genes were up-regulated and 748 genes were down-regulated in the deltamethrin-resistant strain compared to -susceptible strain. These DEGs contained cytochrome P450 s, cuticle proteins, UDP-glucuronosyltransferases, lipases, serine proteases, heat shock proteins, esterases and others. Among the 1826 DEGs, we found that the transcriptional levels of CYP6AA9 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Moreover, the expression levels of the CYP6AA9 were significantly higher in the resistant strains than the susceptible strains in three different field populations. We further confirmed the association between the CYP6AA9 gene and deltamethrin resistance in mosquitoes by RNA interfering (RNAi). Altogether, we explored massive potential pyrethroid resistance-associated genes and demonstrated that CYP6AA9 participated in the pyrethroid resistance in mosquitoes.

Global genome and transcriptome analyses of Magnaporthe oryzae epidemic isolate 98-06 uncover novel effectors and pathogenicity-related genes, revealing gene gain and lose dynamics in genome evolution

Genome dynamics of pathogenic organisms are driven by pathogen and host co-evolution, in which pathogen genomes are shaped to overcome stresses imposed by hosts with various genetic backgrounds through generation of a variety of isolates. This same principle applies to the rice blast pathogen Magnaporthe oryzae and the rice host; however, genetic variations among different isolates of M. oryzae remain largely unknown, particularly at genome and transcriptome levels. Here, we applied genomic and transcriptomic analytical tools to investigate M. oryzae isolate 98-06 that is the most aggressive in infection of susceptible rice cultivars. A unique 1.4 Mb of genomic sequences was found in isolate 98-06 in comparison to reference strain 70-15. Genome-wide expression profiling revealed the presence of two critical expression patterns of M. oryzae based on 64 known pathogenicity-related (PaR) genes. In addition, 134 candidate effectors with various segregation patterns were identified. Five tested proteins could suppress BAX-mediated programmed cell death in Nicotiana benthamiana leaves. Characterization of isolate-specific effector candidates Iug6 and Iug9 and PaR candidate Iug18 revealed that they have a role in fungal propagation and pathogenicity. Moreover, Iug6 and Iug9 are located exclusively in the biotrophic interfacial complex (BIC) and their overexpression leads to suppression of defense-related gene expression in rice, suggesting that they might participate in biotrophy by inhibiting the SA and ET pathways within the host. Thus, our studies identify novel effector and PaR proteins involved in pathogenicity of the highly aggressive M. oryzae field isolate 98-06, and reveal molecular and genomic dynamics in the evolution of M. oryzae and rice host interactions.

Genetic variations in pathogens, such as the causal agent of rice blast Magnaporthe oryzae, often lead to circumvention of disease-resistance cultivars. Previous genome-wide analyses of model organisms suggest that pathogen effectors are also rapidly evolving, especially in regions with high genome plasticity. However, genetic variations among different isolates remain largely unknown in M. oryzae, particularly at the genome and transcriptome levels. In this study, we provided a systematic genomic and interaction transcriptome profile for a dominant rice blast field isolate, resulting in identification of 134 candidate effectors. Two effectors, Iug6 and Iug9, and one pathogenicity-related (PaR) gene product, Iug18, were subjected to functional characterization. We found that Iug6 and Iug9 are located in the biotrophic interfacial complex (BIC) and their overexpression leads to suppression of defense-related gene expression in rice, while Iug18 appears to be a novel PaR protein. Our studies support the hypothesis that isolate-unique genes may serve as a source of genetic variability in the M. oryzae population encountering different environments. Our studies also facilitate further understanding of effectors and genomic variations in pathogenicity of M. oryzae.

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

The Qira black sheep and the Hetian sheep are two local breeds in the Northwest of China, which are characterized by high-fecundity and low-fecundity breed respectively. The elucidation of mRNA expression profiles in the ovaries among different sheep breeds representing fecundity extremes will helpful for identification and utilization of major prolificacy genes in sheep. In the present study, we performed RNA-seq technology to compare the difference in ovarian mRNA expression profiles between Qira black sheep and Hetian sheep. From the Qira black sheep and the Hetian sheep libraries, we obtained a total of 11,747,582 and 11,879,968 sequencing reads, respectively. After aligning to the reference sequences, the two libraries included 16,763 and 16,814 genes respectively. A total of 1,252 genes were significantly differentially expressed at Hetian sheep compared with Qira black sheep. Eight differentially expressed genes were randomly selected for validation by real-time RT-PCR. This study provides a basic data for future research of the sheep reproduction.

Cloning and characterization of a ribosomal protein L23a gene from Small Tail Han sheep by screening of a cDNA expression library

As an indispensable component of the eukaryotic ribosome, ribosomal protein L23a plays an important role in protein synthesis, folding and sorting. In this study, the cDNA fragment of ribosomal protein L23a with 471 bp in size was screened from the Small Tail Han sheep ear marginal tissue cDNA expression library, it has 157 amino acids and a molecular weight of 17.69 kDa. The nucleotide sequence of L23a shares a high homology with those of human, mouse, cattle and pig of 91.51%, 88.32%, 96.18% and 93.84%, respectively. L23a is highly basic, containing a combined 45 Arg, Lys, and His residues and only 14 Asp and Glu residues. The expression pattern and intra-cellular distribution of recombinant L23a proteins in Ujumqin sheep fibroblast cells were analyzed after transfected with the plasmid pEGFP-N3-RPL23A, there were green fluorescence signals both in the cytoplasm and nucleolus of transfected cells after 24 h, the number of positive cells was increased with time, and they reached the peak level after 48 h of transfection. The transfection efficiency was 22.8%. Expression patterns of recombinant L23a gene in Escherichia coli were different with induction temperature, inductor concentration and induction time, when the IPTG concentration was 0.1 mmol/L and induction temperature was 37°, L23a protein expression was increased with induction time.

RNA-Seq analysis implicates detoxification pathways in ovine mycotoxin resistance

Mycotoxin induced hepatoxocity has been linked to oxidative stress, resulting from either an increase in levels of reactive oxygen species (ROS) above normal levels and/or the suppression of antioxidant protective pathways. However, few detailed molecular studies of mycotoxicoses in animals have been carried out. This study use current RNA-seq based approaches to investigate the effects of mycotoxin exposure in a ruminant model. Having first assembled a de novo reference transcriptome, we use RNA-Seq technology to define in vivo hepatic gene expression changes resulting from mycotoxin exposure in relationship to pathological effect. As expected, characteristic oxidative stress related gene expression is markedly different in animals exhibiting poorer outcomes. However, expression of multiple genes critical for detoxification, particularly members of the cytochrome P450 gene family, was significantly higher in animals exhibiting mycotoxin tolerance ('resistance'). Further, we present novel evidence for the amplification of Wnt signalling pathway activity in 'resistant' animals, resulting from the marked suppression of multiple key Wnt inhibitor genes. Notably, 'resistance' may be determined primarily by the ability of an individual to detoxify secondary metabolites generated by the metabolism of mycotoxins and the potentiation of Wnt signalling may be pivotal to achieving a favourable outcome upon challenge.

Analysis of differential gene expression and novel transcript units of ovine muscle transcriptomes

In this study, we characterized differentially expressed genes (DEGs) between the muscle transcriptomes of Small-tailed Han sheep and Dorper sheep and predicted novel transcript units using high-throughput RNA sequencing technology. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that 1,300 DEGs were involved in cellular processes, metabolic pathways, and the actin cytoskeleton pathway. Importantly, we identified 34 DEGs related to muscle cell development and differentiation. Additionally, we were able to optimize the gene structure and predict the untranslated regions (UTRs) for some of the DEGs. Among the 123,678 novel predicted transcript units (TUs), 15,015 units were predicted protein sequences. The reliability of the sequencing data was verified through qRT-PCR analysis of 12 genes. These results will provide useful information for functional genetic research in the future.