RNA-Seq-based transcriptomic profiling of primary interstitial cells of Cajal in response to bovine viral diarrhea virus infection

Transcriptomic Analysis of Metformin's Effect on Bovine Viral Diarrhea Virus Infection

(1) Background: Bovine viral diarrhea virus (BVDV) causes calf diarrhea, bovine respiratory syndrome, and cow abortion, resulting in substantial economic losses in the cattle industry. Owing to its persistent infection mechanism, BVDV is a major challenge in the treatment of cattle. (2) Methods: To determine how metformin (Met) inhibits the interaction between BVDV and host cells, we treated BVDV-infected cells with Met. We then performed an RNA sequencing (RNA-seq) analysis of Met-treated cells infected with BVDV to identify differentially expressed genes (DEGs). Consequently, the RNA-seq results were validated through real-time quantitative PCR (qPCR). (3) Results: Our analysis revealed 3169 DEGs in the Met-treated cells (Met group) vs. the negative controls (NC group) and 2510 DEGs in the BVDV-infected cells after pretreatment with Met (MetBVDV group) vs. the BVDV-infected cells (BVDV group). The DEGs were involved in MDBK interactions during BVDV infection, as indicated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The potential interactions of the DEGs were confirmed via a protein-protein interaction (PPI) network. Met treatment induced autophagy signaling activity and the expression of the autophagy-related genes ATG2A, ATG4B, ATG10, and ATG12 in BVDV-infected Met-pretreated cells. (4) Conclusions: We found that the host transcriptomic profile was affected by BVDV infection and Met pretreatment. These findings offer valuable new insights and provide support for future studies on the inhibition of BVDV replication by Met.

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.

Transcriptome-based comparison reveals key genes regulating allometry growth of forelimb and hindlimb bone in duck embryos

Allometric growth of the forelimb and hindlimb is a widespread phenomenon observed in vertebrates. As a typical precocial bird, ducks exhibit more advanced development of their hindlimbs compared to their forelimbs, enabling them to walk shortly after hatching. This phenomenon is closely associated with the development of long bones in the embryonic stage. However, the molecular mechanism governing the allometric growth of duck forelimb and hindlimb bones is remains elusive. In this study, we employed phenotypic, histological, and gene expression analyses to investigate developmental differences between the humerus (forelimb bone) and tibia/femur (hindlimb bones) in duck embryos. Our results revealed a gradual increase in weight and length disparity between the tibia and humerus from E12 to E28 (embryo age). At E12, endochondral ossification was observed solely in the tibia but not in the humerus. The number of differentially expressed genes (DEGs) gradually increased at H12 vs. T12, H20 vs. T20, and H28 vs. T28 stages consistent with phenotypic variations. A total of 38 DEGs were found across all 3 stages. Protein-protein interaction network analysis demonstrated strong interactions among members of HOXD gene family (HOXD3/8/9/10/11/12), HOXB gene family (HOXB8/9), TBX gene family (TBX4/5/20), HOXA11, SHOX2, and MEIS2. Gene expression profiling indicated higher expression levels for all HOXD genes in the humerus compared to tibia while opposite trends were observed for HOXA/HOXB genes with low or no expression detected in the humerus. These findings suggest distinct roles played by different clusters within HOX gene family during skeletal development regulation of duck embryo's forelimbs versus hind limbs. Notably, TBX4 exhibited high expression levels specifically in tibia whereas TBX5 showed similar patterns exclusively within humerus as seen previously across other species' studies. In summary, this study identified key regulatory genes involved in allometric growth of duck forelimb and hindlimb bones during embryonic development. Skeletal development is a complex physiological process, and further research is needed to elucidate the regulatory role of candidate genes in endochondral ossification.

Effects of testis testosterone deficiency on gene expression in the adrenal gland and skeletal muscle of ducks

1. Testosterone has an anabolic effect on skeletal muscle. The testes produce most of the testosterone in vivo, while the adrenal glands contribute smaller amounts. When testis testosterone is deficient the adrenal gland increases steroid hormone synthesis, which is referred to as compensatory testicular adaptation (CTA).2. To reveal the effects of testis testosterone deficiency on adrenal steroid hormones synthesis and skeletal muscle development, gene expression related to adrenal steroid hormones synthesis and skeletal muscle development were determined by RNA-seq.3. The results showed that castrating male ducks had significant effects on their body weight but no significant impact on cross-sectional area (CSA) or density of pectoral muscle fibres. In skeletal muscle protein metabolism, expression levels of the catabolic gene atrogin1/MAFbx and the anabolic gene eEF2 were significantly higher, with concomitant increases after castration. The adrenal glands' alteration of the steroid hormone 11β-hydroxylase (CYP11B1) was significantly lower following castration.4. Expression pattern analysis showed that the adrenal glands' glucocorticoid receptor (NR3C1/GR) had a potential regulatory relationship with the skeletal muscle-related genes (Pax7, mTOR, FBXO32, FOXO3, and FOXO4).5. The data showed that castration affected muscle protein metabolism, adrenal steroid and testosterone synthesis. In addition, it was speculated that, after castration, steroid hormones produced by the adrenal gland could have a compensatory effect, which might mediate the changes in skeletal muscle protein metabolism and development.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Transcriptomic Analysis of MDBK Cells Infected with Cytopathic and Non-Cytopathic Strains of Bovine Viral Diarrhea Virus (BVDV)

Bovine viral diarrhea virus (BVDV) belongs to the Flaviviridae family and the Pestivirus genus. Infection with BVDV causes a disease with a wide spectrum of clinical symptoms, most often mild, although infections with this virus constitute a serious economic problem all over the world. The virus is characterized by a high genetic variability, while the accumulation of single mutations leads to the formation of its new variants. The aim of this study was to better understand the complicated pathogenesis of this disease at the molecular level via the analysis of the transcriptome of cells infected with this virus. The bovine kidney cell line (MDBK), the cytopathic (cp) reference strain, and two non-cytopathic (ncp) BVD virus field strains were used in transcriptomic studies. The cell transcriptome was tested 24 and 72 h after infection. The results of the microarray analysis revealed changes in the expression levels of numerous genes. Genes with changed expression as a result of infection with the cp strain caused changes in the expression levels of a large number of genes and enriched a number of pathways. Genes with increased expression levels were enriched among other pathways involved in the cell cycle, while genes with reduced expression levels enriched pathways mostly related to metabolism. Genes with increased expression levels as a result of infection with ncp strains enriched a much smaller number of pathways, among them, pathways related to signaling activity 24 h post-infection and serine biosynthetic pathways both 24 and 72 h post-infection. Pathways enriched by genes with reduced expression levels were related to the innate immune response (72 h post-infection) or metabolism (24 and 72 h post-infection). The results of microarray studies can help us to better understand the host’s response to BVDV infection.

Transcriptome analysis of infected Crandell Rees Feline Kidney (CRFK) cells by canine parvovirus type 2c Laotian isolates

The advent of RNA sequencing technology provides insight into the dynamic nature of tremendous transcripts within Crandell-Reese feline kidney (CRFK) cells in response to canine parvovirus (CPV-2c) infection. A total of 1,603 genes displayed differentially expressed genes (DEGs), including 789 up-regulated genes and 814 downregulated genes in the infected cells. Gene expression profiles have shown a subtle pattern of defense mechanism and immune response to CPV through significant DEGs when extensively examined via Gene Ontology (GO) and pathway analysis. Prospective GO analysis was performed and identified several enriched GO biological process terms with significant participating roles in the immune system process and defense response to virus pathway. A Gene network was constructed using the 22 most significantly enriched genes of particular interests in defense response to virus pathways to illustrate the key pathways. Eleven genes (C1QBP, CD40, HYAL2, IFNB1, IFNG, IL12B, IL6, IRF3, LSM14A, MAVS, NLRC5) were identified, which are directly related to the defense response to the virus. Results of transcriptome profiling permit us to understand the heterogeneity of DEGs during in vitro experimental study of CPV infection, reflecting a unique transcriptome signature for the CPV virus. Our findings also demonstrate a distinct scenario of enhanced CPV responses in CRFK cells for viral clearance that involved multistep and perplexity of biological processes. Collectively, our data have given a fundamental role in anti-viral immunity as our highlights of this study, thus providing outlooks on future research priorities to be important in studying CPV.

Developmental Transcriptome Profiling of the Tibial Reveals the Underlying Molecular Basis for Why Newly Hatched Quails Can Walk While Newly Hatched Pigeons Cannot

Birds can be classified into altricial and precocial species. The hatchlings of altricial birds cannot stand, whereas precocial birds can walk and run soon after hatching. It might be owing to the development of the hindlimb bones in the embryo stage, but the molecular regulatory basis underlying the divergence is unclear. To address this issue, we chose the altricial pigeon and the precocial Japanese quail as model animals. The data of tibia weight rate, embryonic skeletal staining, and tibia tissues paraffin section during the embryonic stage showed that the Japanese quail and pigeon have similar skeletal development patterns, but the former had a faster calcification rate. We utilized the comparative transcriptome approach to screen the genes and pathways related to this heterochronism. We separately analyzed the gene expression of tibia tissues of quail and pigeon at two consecutive time points from an inability to stand to be able to stand. There were 2910 differentially expressed genes (DEGs) of quail, and 1635 DEGs of pigeon, respectively. A total of 409 DEGs in common in the quail and pigeon. On the other hand, we compared the gene expression profiles of pigeons and quails at four time points, and screened out eight pairs of expression profiles with similar expression trends but delayed expression in pigeons. By screening the common genes in each pair of expression profiles, we obtained a gene set consisting of 152 genes. A total of 79 genes were shared by the 409 DEGs and the 152 genes. Gene Ontology analysis of these common genes showed that 21 genes including the COL gene family (COL11A1, COL9A3, COL9A1), IHH, MSX2, SFRP1, ATP6V1B1, SRGN, CTHRC1, NOG, and GDF5 involved in the process of endochondral ossification. These genes were the candidate genes for the difference of tibial development between pigeon and quail. This is the first known study on the embryo skeletal staining in pigeon. It provides some new insights for studying skeletal development mechanisms and locomotor ability of altricial and precocial bird species.

HPV16-LINC00393 Integration Alters Local 3D Genome Architecture in Cervical Cancer Cells

High-risk human papillomavirus (hrHPV) infection and integration were considered as essential onset factors for the development of cervical cancer. However, the mechanism on how hrHPV integration influences the host genome structure remains not fully understood. In this study, we performed in situ high-throughput chromosome conformation capture (Hi-C) sequencing, chromatin immunoprecipitation and sequencing (ChIP-seq), and RNA-sequencing (RNA-seq) in two cervical cells, 1) NHEK normal human epidermal keratinocyte; and 2) HPV16-integrated SiHa tumorigenic cervical cancer cells. Our results reveal that the HPV-LINC00393 integrated chromosome 13 exhibited significant genomic variation and differential gene expression, which was verified by calibrated CTCF and H3K27ac ChIP-Seq chromatin restructuring. Importantly, HPV16 integration led to differential responses in topologically associated domain (TAD) boundaries, with a decrease in the tumor suppressor KLF12 expression downstream of LINC00393. Overall, this study provides significant insight into the understanding of HPV16 integration induced 3D structural changes and their contributions on tumorigenesis, which supplements the theory basis for the cervical carcinogenic mechanism of HPV16 integration.

Transcriptome reveals genes involving in black skin color formation of ducks

BACKGROUND

Skin color is colorful for birds, which has been reported to be associated with multi-biological functions, such as crypsis, camouflage, social signaling and mate choice, but little is known about its underlying molecular mechanism.

OBJECTIVE

Studies on the major genes affecting the black skin color of ducks.

METHODS

For this purpose, Silver ammonia staining and RNA-seq analysis were carried out to identify the differences in tissue morphology and gene expressions between black and yellow skin ducks.

RESULTS

The silver ammonia dyes slice results showed that in the development of black duck, the content of melanin in black skin gradually increased and then decreased, and the content of melanin in yellow and black skin was significantly different. Through transcriptome, a total of 102 and 84 differentially expressed genes (DEGs) were identified in beak skin and web skin, respectively. These DEGs were enriched in melanin biosynthesis and play a critical role in melanogenesis pathway. Co-expression analysis showed that EDNRB2 was the only gene associated with black skin color in DEGs, which was also consistent with qRT-PCR.

CONCLUSIONS

The melanin synthesis pathway dominated by EDNRB2 up-regulated the amount of melanin synthesis, leading to the formation of black skin in ducks.

Transcriptome Reveals Multi Pigmentation Genes Affecting Dorsoventral Pattern in Avian Body

Certain animals exhibit a special dorsoventral pattern with a lighter ventral side compared to the dorsal one and this phenomenon was preserved in the long-term evolution process. Birds also retain this trait. Recently, Inaba et al. (2019) found that ASIP (agouti signal protein) regulated interconversion between different melanocyte types leads to dorsal stripe pattern, which may partly explain the birds' dorsoventral plumage color difference. In this study, we used the embryo samples of LBM (light brown mottling) ducks (Anas platyrhynchos) with white ventral and dark dorsal body parts to investigate the mechanism of dorsoventral color variation. Firstly, melanin deposition process of duck embryos was investigated. The result indicated that E13 and E16 were the active stages of melanin synthesis. Moreover, the melanin deposition on the dorsum of LBM ducks was higher than that on the ventral side throughout. Then, RNA-seq was conducted for the dorsal and ventral skin tissues from E7 (early), E13 (middle) and E19 (late) of LBM ducks. Expression pattern analysis showed that the mRNA expression of most melanin synthesis related genes were at the highest level at E13, which was consistent with the section analysis. A correlation was found between melanogenesis pathway and dorsoventral color difference by co-expression analysis. In the DEG (differentially expressed gene) analysis, we added the dorsal skin transcriptome of embryonic white and black duck of same subspecies (Anas platyrhynchos domestica) for horizontal comparison. The results showed that 8 melanogenesis related genes (TYR, TYRP1, MLANA, RAB38, OCA2, TSPAN10, MC1R, and MSLN) were the common DEGs (Differential expressed genes) in the comparisons of body parts and breeds suggesting that the underlying molecular regulatory mechanism of dorsoventral plumage color difference may be similar to that of albino and melanic duck, which were caused by the different expression of multiple genes in melanin synthesis pathway. In addition, the molecular regulation of melanin synthesis pathway in the dorsal and ventral side of LBM ducks was analyzed. In this pathway, ASIP, MC1R, TYR, and TYRP1 have differential mRNA expression. ASIP, as an upstream gene in this pathway, was likely to play a decisive role in determining the dorsoventral plumage pattern.

Detection and characterization of microRNA expression profiling and its target genes in response to canine parvovirus in Crandell Reese Feline Kidney cells

Background

MicroRNAs (miRNAs) play an essential role in gene regulators in many biological and molecular phenomena. Unraveling the involvement of miRNA as a key cellular factor during in vitro canine parvovirus (CPV) infection may facilitate the discovery of potential intervention candidates. However, the examination of miRNA expression profiles in CPV in tissue culture systems has not been fully elucidated.

Method

In the present study, we utilized high-throughput small RNA-seq (sRNA-seq) technology to investigate the altered miRNA profiling in miRNA libraries from uninfected (Control) and CPV-2c infected Crandell Reese Feline Kidney cells.

Results

We identified five of known miRNAs (miR-222-5p, miR-365-2-5p, miR-1247-3p, miR-322-5p and miR-361-3p) and three novel miRNAs (Novel 137, Novel 141 and Novel 102) by sRNA-seq with differentially expressed genes in the miRNA repertoire of CPV-infected cells over control. We further predicted the potential target genes of the aforementioned miRNAs using sequence homology algorithms. Notably, the targets of miR-1247-3p exhibited a potential function associated with cellular defense and humoral response to CPV. To extend the probing scheme for gene targets of miR-1247-3p, we explored and performed Gene Ontology (GO) enrichment analysis of its target genes. We discovered 229 putative targets from a total of 38 enriched GO terms. The top over-represented GO enrichment in biological process were lymphocyte activation and differentiation, marginal zone B cell differentiation, negative regulation of cytokine production, negative regulation of programed cell death, and negative regulation of signaling. We next constructed a GO biological process network composed of 28 target genes of miR-1247-3p, of which, some genes, namely BCL6, DLL1, GATA3, IL6, LEF1, LFNG and WNT1 were among the genes with obviously intersected in multiple GO terms.

Conclusion

The miRNA-1247-3p and its cognate target genes suggested their great potential as novel therapeutic targets or diagnostic biomarkers of CPV or other related viruses.