RESUMEN

Background: Yak is the main livestock species in the plateau area, and its reproductive performance is low, usually two years or three years. A very few of yaks recover within a certain period of time after delivery and smoothly enter the next estrous cycle, while most of them enter the postpartum anestrus and show no estrus performance. However, the key biological factors and influencing mechanisms that cause postpartum anestrus in yaks are not clear. Objective: To study the expression of differential transcripts in ovaries of yak during pregnancy and postpartum anestrus. Methods: Each three yaks in pregnancy and anestrus under natural grazing conditions in Haiyan County, Qinghai Province were selected and slaughtered, and their ovaries were collected and sent to Biomarker Technologies. Oxford Nanopore Technologies single-molecule real-time electrical signal sequencing technology was used to perform full-length transcriptome sequencing. Astalavista software was used to identify the types of alternative splicing events in yak estrus and pregnancy, and TAPIS pipeline was used to identify alternative polyadenylation. Results: The results showed that there were 1751 differentially expressed transcripts (DETs) between pregnancy and anestrus in yak, of which 808 were upregulated and 943 were downregulated. GO analysis showed that the biological processes of DETs were mainly reproductive, reproductive and rhythmic processes. KEGG analysis showed that the DET cell junction-related adhesion junction protein (ß-catenin) and amino terminal kinase (JNK) were involved in FAs (local adhesion). Phosphatidylinositol-3-kinase (PI3K) is involved in the PI3K/AKT/mTOR signaling pathway. Circadian rhythm output cycle failure (Clock) and brain and muscle tissue aromatic hydrocarbon receptor nuclear transporter-like protein 1 (Bmal1) are involved in circadian rhythm signaling pathway. Conclusion: This study found that ß-catenin, JNK, PI3K, Clock and Bmal1 were closely related to postpartum anestrus in yak.

Asunto(s)

RESUMEN

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

RESUMEN

Bovine leukemia virus (BLV) infection in cattle is omnipresent, which causes significantly economical losses worldwide. The objective of this study was to determine microRNA (miRNA) and transcript profiles and to establish their relationship in response to exposure to the virus. Small noncoding and messenger RNA were extracted and sequenced from serum and white blood cells (WBCs) derived from seven BLV seropositive and seven seronegative cows. Transcriptomic profiles were generated by sequencing RNA libraries from WBC. Bta-miR-206 and bta-miR-133a-3p were differentially expressed in serum (P < 0.05). In WBC, bta-miR-335-3p, bta-miR-375, and bta-novel-miR76-3p were differentially expressed (P < 0.03). There were 64 differentially expressed transcripts (DETs). Gene ontology (GO) analysis of the DETs overexpressed in the seropositive group with GOs of response to stimulus and immune system process predicted that the DETs could potentially negatively regulate viral life cycle and viral entry or release from host cells. In addition, the DETs depleted in the seropositive group could play a role in the downregulation of antigen processing and presentation of endogenous peptide antigen via MHC class I. The differentially expressed miRNAs targeted 17 DETs, among which the expressions of bta-miR-133a-3p and bta-miR-335-3p were significantly negatively correlated with the expressions of ENSBTAT00000079143 and ENSBTAT00000066733, respectively. Under high prediction criteria, 90 targets of the differentially expressed miRNAs were all non-DETs. The most enriched biological process GO term of the targets was the RNA-dependent DNA biosynthetic process, which could be associated with virus replication. These results suggested that the differentially expressed miRNAs fine-tune most of the target genes in responding to BLV exposure. In addition, Bta-miR-206 interacted with BLV regulatory genes rex and tax by targeting their coding regions. A further study of the miRNAs and the genes may reveal the molecular mechanisms of BLV infection and uncover possible ways to prevent the infection.

RESUMEN

The commercial value of red tilapia is hampered by variations in skin color during overwintering. In this study, three types of skin of red tilapia, including the skin remained pink color during and after overwintering (P), the skin changed from pink color to black color during overwintering and remained black color after overwintering (P-B), and the skin changed from pink color to black color during overwintering but recovered to pink color when the temperature rose after overwintering (P-B-P), were used to analyze their molecular mechanisms of color variation. The transcriptome results revealed that the P, P-B, and P-B-P libraries had 43, 42, and 43 million clean reads, respectively. The top 10 abundance mRNAs and specific mRNAs (specificity measure SPM > 0.9) were screened. After comparing intergroup gene expression levels, there were 2528, 1924, and 1939 differentially expressed genes (DEGs) between P-B-P and P-B, P-B-P and P, and P-B and P, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of color-related mRNAs showed that a number of DEGs, including tyrp1, tyr, pmel, mitf, mc1r, asip, tat, hpdb, and foxd3, might play a potential role in pigmentation. Additionally, the co-expression patterns of genes were detected within the pigment-related pathways by the PPI network from P-B vs. P group. Furthermore, DEGs from the apoptosis and autophagy pathways, such as baxα, beclin1, and atg7, might be involved in the fading of red tilapia melanocytes. The findings will aid in understanding the molecular mechanism underlying skin color variation in red tilapia during and after overwintering as well as lay a foundation for future research aimed at improving red tilapia skin color characteristics.

Asunto(s)

RESUMEN

BACKGROUND: Gossypium hirsutum seeds are rich in gossypol. In addition to its diverse beneficial properties, it is a known anti-fertility inducing agent in humans. Oxycarenus laetus feeds on the cottonseeds and yet its courtship, mating and reproduction is unaffected. METHODS AND RESULTS: In this study, we performed a transcriptome profiling of O. laetus fed on Abutilon indicum (AB-no gossypol), G. hirsutum (GH-natural gossypol) and 1400 ppm commercial gossypol-soaked GH seeds (GHGO). Illumina NextSeq-500 paired-end 75 bp reads were generated and de novo assembled (48,214 genes) to identify the differentially expressed transcripts (DET) between the samples. Gene enrichment, KEGG pathway and cluster profiling of the DETs resulted in the identification of vital genes involved in the detoxification, pheromone biosynthesis, cuticle protein in the GHGO sample. Cyp4C1, Cyp6a13, Cyp6a14, Cyp4g15, Cyp4em8, Cyp303a1 were the detoxification related genes identified. Similarly, SDR dehydrogenase family 11 and fatty acid synthase in pheromone biosynthesis and cuticle proteins (RR1 and RR2) coding transcripts were found to be differentially expressed. CONCLUSION: This is the first study to report the expression of genes induced by gossypol in O. laetus. Based on the findings from the DET analysis, we conclude that the detoxification related genes of gossypol treated samples were affected.

Asunto(s)

RESUMEN

BACKGROUND: The yak (Bos grunniens) is an important livestock species that can survive the extremely cold, harsh, and oxygen-poor conditions of the Qinghai-Tibetan Plateau and provide meat, milk, and transportation for the Tibetans living there. However, the regulatory network that drive this hypoxic adaptation remain elusive. RESULTS: The heart tissues from LeiRoqi (LWQY) yak and their related cattle (Bos Taurus) breeds, which are two native cattle breeds located in high altitude (HAC) and low altitude (LAC) regions, respectively, were collected for RNA sequencing. A total of 178 co-differentially expressed protein-coding transcripts (co-DETs) were discovered in each of the LAC-vs-LWQY and LAC-vs-HAC comparison groups, including NFATC2, NFATC1, ENPP2, ACSL4, BAD, and many other genes whose functions were reported to be associated with the immune-system, endocrine-system, and lipid metabolism. Two and 230 lncRNA transcripts were differentially expressed in the LAC-vs-LWQY and LAC-vs-HAC comparisons' respectively, but no lncRNA transcripts that were co-differentially expressed. Among the 58 miRNAs that were co-differentially expressed, 18 were up-regulated and 40 were down-regulated. In addition, 640 (501 up-regulated and 139 down-regulated) and 152 (152 up-regulated and one down-regulated) circRNAs showed differential expression in LAC-vs-LWQY and LAC-vs-HAC comparison groups, respectively, and 53 up-regulated co-differentially expressed circRNAs were shared. Multiple co-DETs, which are the targets of miRNAs/lncRNAs, are significantly enriched in high-altitude adaptation related processes, such as, T cell receptor signaling, VEGF signaling, and cAMP signaling. A competing endogenous RNA (ceRNA) network was constructed by integrating the competing relationships among co-differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs. Furthermore, the hypoxic adaptation related ceRNA network was constructed, and the six mRNAs (MAPKAPK3, PXN, NFATC2, ATP7A, DIAPH1, and F2R), the eight miRNAs (including miR-195), and 15 circRNAs (including novel-circ-017096 and novel-circ-018073) are proposed as novel and promising candidates for regulation of hypoxic adaptation in the heart. CONCLUSION: In conclusion, the data recorded in the present study provides new insights into the molecular network of high-altitude adaptation along with more detailed information of protein-coding transcripts and non-coding transcripts involved in this physiological process, the detailed mechanisms behind how these transcripts "crosstalk" with each other during the plateau adaptation are worthy of future research efforts.

RESUMEN

Genus Gnetum, of which the majority species are pantropical liana, have broad industrial uses including for string, nets, and paper production. Although numerous studies have investigated anatomical structures during stem development, the underlying molecular mechanisms that regulate this developmental trajectory in Gnetum species remain poorly understood. A total of 12 full-length transcriptomes were generated from four stem developmental stages of an arborescent representative of this genus, Gnetum luofuense, using Oxford Nanopore Technologies. The results of this analysis reveal a total of 24,151 alternative splicing (AS) and 134,391 alternative polyadenylation events. A remarkably dynamic pattern of AS events, especially in the case of intron retentions, was found across the four developmental stages while no dynamic pattern was found among transcript numbers with varied poly(A) sites. A total of 728 long non-coding RNAs were also detected; the number of cis-regulated target genes dramatically increased while no changes were found among trans-regulated target genes. In addition, a K-means clustering analysis of all full-length transcripts revealed that primary growth is associated with carbohydrate metabolism and fungi defense, while secondary growth is closely linked with photosynthesis, nitrogen transportation, and leaf ontogenesis. The use of weighted gene co-expression network analysis as well as differentially expressed transcripts reveals that bHLH, GRF, and MYB-related transcription factors are involved in primary growth, while AP2/ERF, MYB, NAC, PLAZ, and bZIP participate in G. luofuense stem secondary growth. The results of this study provide further evidence that Nanopore sequencing technology provides a cost-effective method for generating full-length transcriptome data as well as for investigating seed plant organ development.

RESUMEN

Deciphering sperm transcriptome is the key to understanding the molecular mechanisms governing peri-fertilization, embryonic development, and pregnancy establishment. This study aimed to profile sperm transcriptome to identify signature transcripts regulating male fertility. Semen samples were collected from 47 bulls with varied fertility rates. The sperm total RNA was isolated (n = 8) and subjected to transcriptome sequencing. Based on the expression pattern obtained from RNA profiling, the bulls were grouped (p = 0.03) into high-fertile and sub-fertile, and signature transcripts controlling sperm functions and fertility were identified. The results were validated using the OMIM database, qPCR, and sperm function tests. The sperm contains 1100 to 1700 intact transcripts, of which BCL2L11 and CAPZA3 were abundant and associated (p < 0.05) with spermatogenesis and post-embryonic organ morphogenesis. The upregulated genes in the acrosome integrity and functional membrane integrity groups had a close association with the fertility rate. The biological functions of these upregulated genes (p < 0.05) in the high-fertile bulls were associated with spermatogenesis (AFF4 and BRIP1), sperm motility (AK6 and ATP6V1G3), capacitation and zona binding (AGFG1), embryo development (TCF7 and AKIRIN2), and placental development (KRT19). The transcripts involved in pathways regulating embryonic development such as translation (EEF1B2 and MTIF3, p = 8.87E-05) and nonsense-mediated decay (RPL23 and RPL7A, p = 5.01E-27) were upregulated in high-fertile bulls. The identified transcripts may significantly impact oocyte function, embryogenesis, trophectoderm development, and pregnancy establishment. In addition, the study also reveals that the genes governing sperm functional membrane integrity and acrosome integrity have a prospective effect on male fertility.

Asunto(s)

RESUMEN

BACKGROUND: The pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms. RESULTS: Herein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment. CONCLUSIONS: Our results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.

Asunto(s)

RESUMEN

BACKGROUND: The Intramuscular fat (IMF) content in meat products, which is positively correlated with meat quality, is an important trait considered by consumers. The regulation of IMF deposition is species specific. However, the IMF-deposition-related mRNA and non-coding RNA and their regulatory network in yak (Bos grunniens) remain unknown. High-throughput sequencing technology provides a powerful approach for analyzing the association between transcriptome-related differences and specific traits in animals. Thus, the whole transcriptomes of yak muscle and adipose tissues were screened and analyzed to elucidate the IMF deposition-related genes. The muscle tissues were used for IMF content measurements. RESULTS: Significant differences were observed between the 0.5- and 2.5-year-old yaks. Several mRNAs, miRNAs, lncRNAs and circRNAs were generally expressed in both muscle and adipose tissues. Between the 0.5- and 2.5-year-old yaks, 149 mRNAs, 62 miRNAs, 4 lncRNAs, and 223 circRNAs were differentially expressed in muscle tissue, and 72 mRNAs, 15 miRNAs, 9 lncRNAs, and 211 circRNAs were differentially expressed in adipose tissue. KEGG annotation revelved that these differentially expressed genes were related to pathways that maintain normal biological functions of muscle and adipose tissues. Moreover, 16 mRNAs, 5 miRNAs, 3 lncRNAs, and 5 circRNAs were co-differentially expressed in both types of tissue. We suspected that these co-differentially expressed genes were involved in IMF-deposition in the yak. Additionally, LPL, ACADL, SCD, and FASN, which were previously shown to be associated with the IMF content, were identified in the competing endogenous RNA (ceRNA) regulatory network that was constructed on the basis of the IMF deposition-related genes. Three ceRNA subnetworks also revealed that TCONS-00016416 and its target SIRT1 "talk" to each other through the same miR-381-y and miR-208 response elements, whereas TCONS-00061798 and its target PRKCA, and TCONS-00084092 and its target LPL "talk" to each other through miR-122-x and miR-499-y response elements, respectively. CONCLUSION: Taken together, our results reveal the potential mRNA and noncoding RNAs involved in IMF deposition in the yak, providing a useful resource for further research on IMF deposition in this animal species.

Asunto(s)

RESUMEN

Differential Display (DD) is a technique widely used in studies of differential expression. Most of these analyses, especially those involving fish species, are restricted to species from North America and Europe or to commercial species, as salmonids. Studies related to South American fish species are underexplored. Thus, the present work aimed to describe DD technique modifications in order to improve outcomes related to the isolation of DETs (Differentially Expressed Transcripts), using Leporinus macrocephalus, a large commercially exploited South American species, as a fish design. Different DDRT-PCR approaches were applied to brain samples and the products of the reactions were analyzed on 6% polyacrylamide gels stained with 0.17% Silver Nitrate (AgNO3). The use of PCR reactions under high stringency conditions and longer oligonucleotides based on VNTR (Variable Number of Tandem Repeats) core sequences led to better results when compared to low stringency PCR conditions and the use of decamer oligonucleotides. The improved approach led to the isolation of differentially expressed transcripts on adult males and females of L. macrocephalus. This study indicates that some modifications on the DDRT-PCR method can ensure isolation of DETs from different fish tissues and the development of robust data related to this approach.

Display Diferencial (DD) é uma técnica amplamente utilizada em estudos de expressão diferencial. A maioria desses estudos envolvendo espécies de peixes está restrita a espécies da América do Norte e Europa ou a espécies comerciais, como os salmoniformes. Estudos relacionados a peixes da América do Sul são ainda pouco explorados. Desse modo, o presente trabalho teve como objetivo descrever modificações na técnica de DD, a fim de melhorar os resultados relacionados ao isolamento de DETs (Transcritos Diferencialmente Expressos), utilizando Leporinus macrocephalus, peixe explorado comercialmente na América do Sul, como espécie para tal delineamento. Diferentes abordagens de DDRT-PCR foram desenvolvidas a partir de amostras de tecido cerebral e os produtos das reações foram analisados em gel de poliacrilamida 6% corados com 0,17% de nitrato de Prata (AgNO3). A utilização de reações de PCR sob condições de elevada estringência e oligonucleotídeos mais longos, com base em sequências cerne de VNTR (Número Variável de Repetições em Tandem), mostrou melhores resultados quando comparada a condições de baixa estringência e ao uso de oligonucleotídeos decâmeros. A estratégia empregada permitiu o isolamento de transcritos diferencialmente expressos em machos e fêmeas adultos de L. macrocephalus. Este estudo evidencia que modificações no método de DDRT-PCR garantem o melhor isolamento de DETs a partir de diferentes tecidos de peixes e asseguram a obtenção de dados mais sólidos relacionados a essa abordagem.

Asunto(s)