RNA-seq screening and gene function analysis uncover GPR183 as a key mediator for methionine to stimulate milk synthesis in mouse mammary epithelial cells

Methionine (Met) can activate mTOR to promote milk synthesis in mammary epithelial cells (MECs). However, it is largely unknown which G protein-coupled receptor (GPCR) can mediate the stimulation of Met on mTOR activation. In this study, we employed transcriptome sequencing to analyze which GPCRs were associated with the role of Met, and further used gene function study approaches to explore the role of GPR183 in Met stimulation on mTOR activation in HC11 cells. We identified 9 GPCRs including GPR183 which expression levels were upregulated by Met treatment through RNA-seq and subsequent RT-qPCR analysis. Using GPR183 knockdown and overexpression technology, we demonstrate that GPR183 is a positive regulator of milk protein and fat synthesis and proliferation of HC11 cells. Met affected GPR183 expression in a dose-dependent manner, and GPR183 mediated the stimulation of Met (0.6 mM) on milk protein and fat synthesis, cell proliferation, and mTOR phosphorylation and mRNA expression. The inhibition of PI3K blocked the phosphorylation of mTOR and AKT stimulated by GPR183 activation. In summary, through RNA-seq and gene function study, we uncover that GPR183 is a key mediator for Met to activate the PI3K-mTOR signaling and milk synthesis in mouse MECs.

RNA-clique: a method for computing genetic distances from RNA-seq data

BACKGROUND

Although RNA-seq data are traditionally used for quantifying gene expression levels, the same data could be useful in an integrated approach to compute genetic distances as well. Challenges to using mRNA sequences for computing genetic distances include the relatively high conservation of coding sequences and the presence of paralogous and, in some species, homeologous genes.

RESULTS

We developed a new computational method, RNA-clique, for calculating genetic distances using assembled RNA-seq data and assessed the efficacy of the method using biological and simulated data. The method employs reciprocal BLASTn followed by graph-based filtering to ensure that only orthologous genes are compared. Each vertex in the graph constructed for filtering represents a gene in a specific sample under comparison, and an edge connects a pair of vertices if the genes they represent are best matches for each other in their respective samples. The distance computation is a function of the BLAST alignment statistics and the constructed graph and incorporates only those genes that are present in some complete connected component of this graph. As a biological testbed we used RNA-seq data of tall fescue (Lolium arundinaceum), an allohexaploid plant ( 2 n = 14 Gb ), and bluehead wrasse (Thalassoma bifasciatum), a teleost fish. RNA-clique reliably distinguished individual tall fescue plants by genotype and distinguished bluehead wrasse RNA-seq samples by individual. In tests with simulated RNA-seq data, the ground truth phylogeny was accurately recovered from the computed distances. Moreover, tests of the algorithm parameters indicated that, even with stringent filtering for orthologs, sufficient sequence data were retained for the distance computations. Although comparisons with an alternative method revealed that RNA-clique has relatively high time and memory requirements, the comparisons also showed that RNA-clique's results were at least as reliable as the alternative's for tall fescue data and were much more reliable for the bluehead wrasse data.

CONCLUSION

Results of this work indicate that RNA-clique works well as a way of deriving genetic distances from RNA-seq data, thus providing a methodological integration of functional and genetic diversity studies.

Identification of genes associated with sperm storage capacity in hens at different times after insemination by RNA-seq and Ribo-seq

BACKGROUND

Sperm storage capacity (SSC) determines the duration of fertility in hens and is an important reproduction trait that cannot be ignored in production. Currently, the genetic mechanism of SSC is still unclear in hens. Therefore, to explore the genetic basis of SSC, we analyzed the uterus-vagina junction (UVJ) of hens with different SSC at different times after insemination by RNA-seq and Ribo-seq.

RESULTS

Our results showed that 589, 596, and 527 differentially expressed genes (DEGs), 730, 783, and 324 differentially translated genes (DTGs), and 804, 625, and 467 differential translation efficiency genes (DTEGs) were detected on the 5th, 10th, and 15th days after insemination, respectively. In transcription levels, we found that the differences of SSC at different times after insemination were mainly reflected in the transmission of information between cells, the composition of intercellular adhesion complexes, the regulation of ion channels, the regulation of cellular physiological activities, the composition of cells, and the composition of cell membranes. In translation efficiency (TE) levels, the differences of SSC were mainly related to the physiological and metabolic activities in the cell, the composition of the organelle membrane, the physiological activities of oxidation, cell components, and cell growth processes. According to pathway analysis, SSC was related to neuroactive ligand-receptor interaction, histidine metabolism, and PPAR signaling pathway at the transcriptional level and glutathione metabolism, oxidative phosphorylation, calcium signaling pathway, cell adhesion molecules, galactose metabolism, and Wnt signaling pathway at the TE level. We screened candidate genes affecting SSC at transcriptional levels (COL4A4, MUC6, MCHR2, TACR1, AVPR1A, COL1A1, HK2, RB1, VIPR2, HMGCS2) and TE levels(COL4A4, MUC6, CYCS, NDUFA13, CYTB, RRM2, CAMK4, HRH2, LCT, GCK, GALT). Among them, COL4A4 and MUC6 were the key candidate genes differing in transcription, translation, and translation efficiency.

CONCLUSIONS

Our study used the combined analysis of RNA-seq and Ribo-seq for the first time to investigate the SSC and reveal the physiological processes associated with SSC. The key candidate genes affecting SSC were screened, and the theoretical basis was provided for the analysis of the molecular regulation mechanism of SSC.

HBeAg induces neutrophils activation impairing NK cells function in patients with chronic hepatitis B

BACKGROUND

The role of neutrophils in hepatitis B virus (HBV) infection has been a subject of debate due to their involvement in antiviral responses and immune regulation. This study aimed to elucidate the neutrophil characteristics in patients with chronic hepatitis B (CHB).

METHODS

Through flow cytometry and ribonucleic acid-sequencing analysis, the phenotypes and counts of neutrophils were analyzed in patients with CHB. Moreover, the effects of HBeAg on neutrophils and the corresponding pattern recognition receptors were identified. Simultaneously, the cross-talk between neutrophils and natural killer (NK) cells was investigated.

RESULTS

Neutrophils were activated in patients with CHB, characterized by higher expression levels of programmed death-ligand 1 (PD-L1), cluster of differentiation 86, and interleukin-8, and lower levels of CXC motif chemokine receptor (CXCR) 1 and CXCR2. Hepatitis B e antigen (HBeAg) partially induces neutrophil activation through the Toll-like receptor 2 (TLR2). A consistent upregulation of the TLR2 and HBeAg expression was observed in patients with CHB. Notably, the genes encoding molecules pivotal for NK-cell function upon NK receptor engagement enriched in neutrophils after HBeAg activation. The HBeAg-activated neutrophils demonstrated the ability to decrease the production of interferon-gamma (IFN-γ) and tumor necrosis factor-alpha (TNF-α) in NK cells, while the PD-1 and PD-L1 pathways partially mediated the immunosuppression.

CONCLUSIONS

The immunosuppression of neutrophils induced by HBeAg suggests a novel pathogenic mechanism contributing to immune tolerance in patients with CHB.

Competing endogenous RNA networks were associated with fat accumulation in skeletal muscle of aged male mice

Muscle aging contributed to morbidity and mortality in the elderly adults by leading to severe outcomes such as frailty, falls and fractures. Post-transcriptional regulation especially competing endogenous RNA (ceRNA) mechanism may modulate the process of skeletal muscle aging. RNA-seq was performed in quadriceps of 6-month-old (adult) and 22-month-old (aged) male mice to identify differentially expressed ncRNAs and mRNAs and further construct ceRNA networks. Decreased quadriceps-body weight ratio and muscle fiber cross-sectional area as well as histological characteristics of aging were observed in the aged mice. Besides, there were higher expressions of atrogin-1 and MuRF-1 and lower expression of Myog, Myf4 and Myod1 in the quadriceps of aged mice relative to that of adult mice. The expression of 85 lncRNAs, 52 circRNAs, 10 miRNAs and 277 mRNAs were significantly dysregulated in quadriceps between the two groups, among which two ceRNA networks lncRNA 2700081O15Rik/circRNA_0000820-miR-673-3p-Tmem120b were constructed. Level of triglycerides and expression of PPARγ, C/EBPα, FASN and leptin were elevated and the expression of adiponectin was reduced in quadriceps of aged mice compared with that of adult mice. LncRNA 2700081O15Rik/circRNA_0000820-miR-673-3p-Tmem120b were possibly associated with the adipogenesis and fat accumulation in skeletal muscle of age male mice.

Differential Impact of Smoking on Methylome and Transcriptome in Crohn's Disease and Ulcerative Colitis

BACKGROUND

Smoking is an environmental factor that differentially impacts Crohn's disease (CD) and ulcerative colitis (UC). The mechanism of impact of smoking on disease risk and clinical outcomes remains to be established.

METHODS

This study used a prospective cohort of patients with CD or UC. Self-reported smoking status was validated using serum cotinine measurement. We profiled methylation changes in peripheral blood using the Illumina Methylation BeadChip. Transcriptomic profiling was performed on ileal and colonic tissue using an Illumina TruSeq platform. We compared the methylation and transcriptional changes in current, former, and never smokers stratified by disease type.

RESULTS

Our cohort included 200 patients with CD or UC with methylation profiles and 160 with transcriptomic data. The mean serum cotinine level was higher in current compared with former or never smokers. Epigenetic changes common to both CD and UC included hypomethylation at AHRR. Smoking-associated MGAT3 hypomethylation was associated with severe disease course only in UC, while IER3 hypomethylation was associated with worse course only in CD. Smoking downregulated several inflammatory pathways in UC. Current smoking in CD but not in UC was associated with upregulation of several genes mediating Paneth cell function. Genes with opposite direction of effects in CD and UC include HSD3B2 and GSTA1.

CONCLUSIONS

Our findings suggest both common and differential effects of cigarette smoking on CD and UC. Paneth cell dysfunction may mediate adverse impact of smoking on CD. Bile acid and oxidative stress pathways may be relevant for the differential effect of smoking on CD and UC.

Microglia contribute to methamphetamine reinforcement and reflect persistent transcriptional and morphological adaptations to the drug

Methamphetamine use disorder (MUD) is a chronic, relapsing disease that is characterized by repeated drug use despite negative consequences and for which there are currently no FDA-approved cessation therapeutics. Repeated methamphetamine (METH) use induces long-term gene expression changes in brain regions associated with reward processing and drug-seeking behavior, and recent evidence suggests that methamphetamine-induced neuroinflammation may also shape behavioral and molecular responses to the drug. Microglia, the resident immune cells in the brain, are principal drivers of neuroinflammatory responses and contribute to the pathophysiology of substance use disorders. Here, we investigated transcriptional and morphological changes in dorsal striatal microglia in response to methamphetamine-taking and during methamphetamine abstinence, as well as their functional contribution to drug-taking behavior. We show that methamphetamine self-administration induces transcriptional changes associated with protein folding, mRNA processing, immune signaling, and neurotransmission in dorsal striatal microglia. Importantly, many of these transcriptional changes persist through abstinence, a finding supported by morphological analyses. Functionally, we report that microglial ablation increases methamphetamine-taking, possibly involving neuroimmune and neurotransmitter regulation. In contrast, microglial depletion during abstinence does not alter methamphetamine-seeking. Taken together, these results suggest that methamphetamine induces both short and long-term changes in dorsal striatal microglia that contribute to altered drug-taking behavior and may provide valuable insights into the pathophysiology of MUD.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development

STUDY QUESTION

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development?

SUMMARY ANSWER

Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce.

STUDY DESIGN, SIZE, DURATION

This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome.

MAIN RESULTS AND THE ROLE OF CHANCE

A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation.

LARGE SCALE DATA

All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501.

LIMITATIONS, REASONS FOR CAUTION

Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells.

WIDER IMPLICATIONS OF THE FINDINGS

The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

Cutin-Derived Oligomers Induce Hallmark Plant Immune Responses

The cuticle constitutes the outermost defensive barrier of most land plants. It comprises a polymeric matrix - cutin, surrounded by soluble waxes. Moreover, the cuticle constitutes the first line of defense against pathogen invasion, while also protecting the plant from many abiotic stresses. Aliphatic monomers in cutin have been suggested to act as immune elicitors in plants. This study analyses the potential of cutin oligomers to activate rapid signaling outputs reminiscent of pattern-triggered immunity (PTI) in the model plant Arabidopsis. Cutin oligomeric mixtures led to Ca2+ influx and MAPK activation. Comparable responses were measured for cutin, which was also able to induce a reactive oxygen species (ROS) burst. Furthermore, cutin oligomer treatment resulted in a unique transcriptional reprogramming profile, having many archetypal features of PTI. Targeted spectroscopic and spectrometric analyses of the cutin oligomers suggest that the elicitors compounds consist mostly of two up to three 10,16-dihydroxyhexadecanoic acid monomers linked together through ester bonds. This study demonstrates that cutin breakdown products can act as inducers of early plant immune responses, which underlying mechanisms of perception and potential use in agriculture warrant further investigation.

Transcriptomic sequence dataset of a potential new model species for studying biomineralization Onchidoris muricata (Nudibranchia, Gastropoda, Mollusca)

Onchidoris muricata is a widespread shell-less species of nudibranch molluscs, which has unique for Gastropoda skeletal elements - subepidermal calcite spicules. The general and fine morphology of the spicules, as well as their maturation process in ontogenesis, have been studied in detail by authors. The uniqueness of spicules lies in their intracellular formation and location under the ectodermal epithelium, which is more typical for deuterostomes. We present O. muricata as a potentially new model species for studying calcification of intracellular protein structure. A total of 96 individuals were collected in the Kandalaksha Bay of the White Sea, both manually and by scuba diving. All individuals were divided into three groups based on morphological characteristics such as specimens' size, spicule condition etc. This division suggests the existence of three stages in postembryonic ontogenesis of O. muricata reflecting the maturation of the spicule complex. Total RNA samples were isolated from three size groups of molluscs in three biological replicates. Libraries were prepared from the polyadenylated RNA fraction and sequenced at NovaSeq6000 (Illumina), yielding a total of 112.8 Gb of 150 bp paired-end reads, corresponding to almost 1,000-fold coverage of the transcriptome. Representative transcriptome assembled de novo with Trinity. In addition to obtaining the transcriptome sequences of O. muricata, differential expression analysis was also performed for these three size groups. This allows us to trace the dynamics of molecular and biological processes during the life of a mollusc. The obtained data can then be used as a reference transcriptome for closely related species, to study specific expressed genes, to identify various unique sequences, including protein-coding ones, to understand biological processes, including biomineralization and much more.

Transcriptome dataset for comparative analysis of differentially expressed genes between wild type and transgenic potato plants overexpressing Nuclear Factor Y subunit A7

Nuclear Factor Y (NF-Y) is divided into three different types of subunits, A, B, and C. NF-Ys play crucial roles in plants for controlling gene expression associated with various developmental processes and abiotic stresses, but it is mostly unknown the downstream genes regulated by NF-Ys in plant. One of the potato NF-Y genes, StNF-YA7, increased potato's drought tolerance when overexpressed under the control of constitutive CaMV 35S promoter. Therefore, it was of interest what genes are regulated by the increased expression level of StNF-YA7. To investigate the downstream genes of StNF-YA7, the transcriptome sequencing was carried out for four potato lines, including Solanum tuberosum L 'Superior' as wild type (WT), empty vector control (VC), and two StNF-YA7 overexpressor lines (designated to StNF-YA7 #19 & #26). The RNA sequencing data was produced by the Illumina NovaSeq 6000 sequencing system. The number of total raw reads obtained from the RNA sequencing was 36.7 million for WT, 36.2 for VC, 29.3 for StNF-YA7 #19, and 29.5 million for StNF-YA7 #26, respectively. The length of total raw reads for each sample was between 5.92 Gb (StNF-YA7 #19) and 7.42 Gb (WT), and after filtering raw quality reads, the total length was between 5.81 Gb (StNF-YA7 #19) and 7.29 Gb (WT). Each filtered clear read set of four transcriptomes was mapped on the potato reference genome, SolTub_3.0, and the percentage of mapped reads ranged from 89.8 % (VC) to 90.3 % (WT). GC contents range between 43.01 % (StNF-YA7 #19) and 42.44 % (StNF-YA7 #26). Q20 quality score ranges between 98.63 % (StNF-YA7 #26) and 98.74 % (VC).

Multiomics analyses of Jining Grey goat and Boer goat reveal genomic regions associated with fatty acid and amino acid metabolism and muscle development

OBJECTIVE

Jining Grey goat is a local Chinese goat breed that is well known for its high fertility and excellent meat quality but shows low meat production performance. Numerous studies have focused on revealing the genetic mechanism of its high fertility, but its highlighting meat quality and muscle growth mechanism still need to be studied.

METHODS

In this research, an integrative analysis of the genomics and transcriptomics of Jining Grey goats compared with Boer goats was performed to identify candidate genes and pathways related to the mechanisms of meat quality and muscle development.

RESULTS

Our results overlap among five genes (ABHD2, FN1, PGM2L1, PRKAG3, RAVER2) and detected a set of candidate genes associated with fatty acid metabolism (PRKAG3, HADHB, FASN, ACADM), amino acid metabolism (KMT2C, PLOD3, NSD2, SETDB1, STT3B, MAN1A2, BCKDHB, NAT8L, P4HA3) and muscle development (MSTN, PPARGC1A, ANKRD2). Several pathways have also been detected, such as the FoxO signaling pathway and Apelin signaling pathway that play roles in lipid metabolism, lysine degradation, N-glycan biosynthesis, valine, leucine and isoleucine degradation that involving with amino acid metabolism.

CONCLUSION

The comparative genomic and transcriptomic analysis of Jining Grey goat and Boer goat revealed the mechanisms underlying the meat quality and meat productive performance of goats. These results provide valuable information for future breeding of goats.

PvGeneExpDB: An integrative gene expression database for in-depth understanding on the Pacific white shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei) is a high-valued economic farming species. With the development of high-throughput sequencing technology, cumulative large-scale transcriptomic studies have been revealing molecular landscape of various biological conditions including genetic selection, breeding, evolution, disease landscape, etc. However, no single experiment or databases allow thorough investigations of transcriptomic dynamics for these progressions. Meanwhile, the available datasets are often scattered and lack management. Here, we have established PvGeneExpDB, the first gene expression database for L. vannamei (www.bio-marine-scau.com/pv_ex/), which encompasses gene expression profiles, differential expression, and co-expression analyses under various biological conditions. Based on the analyses of 7 datasets, which include 53 samples with accurate and detailed records, PvGeneExpDB identifies 20,599 novel transcripts, shows expression profiles of a total of 20,817 genes, and implements Gene Ontology (GO) reconstruction of 76.7 % of these genes. Besides, 26 co-expressed groups were first identified by large-scale, cross-sample Weighted Gene Co-expression Network Analysis (WGCNA). By integrating the gene expression data in the database, our goal is to deepen the biological understanding of L. vannamei.

Disturbed flow in the juxta-anastomotic area of an arteriovenous fistula correlates with endothelial loss, acute thrombus formation, and neointimal hyperplasia

Clinical failure of arteriovenous neointimal hyperplasia (NIH) fistulae (AVF) is frequently due to juxta-anastomotic NIH (JANIH). Although the mouse AVF model recapitulates human AVF maturation, previous studies focused on the outflow vein distal to the anastomosis. We hypothesized that the juxta-anastomotic area (JAA) has increased NIH compared with the outflow vein. AVF was created in C57BL/6 mice without or with chronic kidney disease (CKD). Temporal and spatial changes of the JAA were examined using histology and immunofluorescence. Computational techniques were used to model the AVF. RNA-seq and bioinformatic analyses were performed to compare the JAA with the outflow vein. The jugular vein to carotid artery AVF model was created in Wistar rats. The neointima in the JAA shows increased volume compared with the outflow vein. Computational modeling shows an increased volume of disturbed flow at the JAA compared with the outflow vein. Endothelial cells are immediately lost from the wall contralateral to the fistula exit, followed by thrombus formation and JANIH. Gene Ontology (GO) enrichment analysis of the 1,862 differentially expressed genes (DEG) between the JANIH and the outflow vein identified 525 overexpressed genes. The rat jugular vein to carotid artery AVF showed changes similar to the mouse AVF. Disturbed flow through the JAA correlates with rapid endothelial cell loss, thrombus formation, and JANIH; late endothelialization of the JAA channel correlates with late AVF patency. Early thrombus formation in the JAA may influence the later development of JANIH.NEW & NOTEWORTHY Disturbed flow and focal endothelial cell loss in the juxta-anastomotic area of the mouse AVF colocalizes with acute thrombus formation followed by late neointimal hyperplasia. Differential flow patterns between the juxta-anastomotic area and the outflow vein correlate with differential expression of genes regulating coagulation, proliferation, collagen metabolism, and the immune response. The rat jugular vein to carotid artery AVF model shows changes similar to the mouse AVF model.

The hypoxia response pathway in the Antarctic fish Notothenia coriiceps is functional despite a poly Q/E insertion mutation in HIF-1α

Antarctic notothenioid fishes, inhabiting the oxygen-rich Southern Ocean, possess a polyglutamine and glutamic acid (poly Q/E) insertion mutation in the master transcriptional regulator of oxygen homeostasis, hypoxia- inducible factor-1α (HIF-1α). To determine if this mutation impairs the ability of HIF-1 to regulate gene expression in response to hypoxia, we exposed Notothenia coriiceps, with a poly Q/E insertion mutation in HIF-1α that is 9 amino acids long, to hypoxia (2.3 mg L-1 O2) or normoxia (10 mg L -1 O2) for 12 h. Heart ventricles, brain, liver, and gill tissue were harvested and changes in gene expression quantified using RNA sequencing. Levels of glycogen and lactate were also quantified to determine if anaerobic metabolism increases in response to hypoxia. Exposure to hypoxia resulted in 818 unique differentially expressed genes (DEGs) in liver tissue of N. coriiceps. Many hypoxic genes were induced, including ones involved in the MAP kinase and FoxO pathways, glycolytic metabolism, and vascular remodeling. In contrast, there were fewer than 104 unique DEGs in each of the other tissues sampled. Lactate levels significantly increased in liver in response to hypoxia, indicating that anaerobic metabolism increases in response to hypoxia in this tissue. Overall, our results indicate that the hypoxia response pathway is functional in N. coriiceps despite a poly Q/E mutation in HIF-1α, and confirm that Antarctic fishes are capable of altering gene expression in response to hypoxia.

Methylosome protein 50 is necessary for oogenesis in medaka

Methylosome protein 50 (Mep50) functions as a partner to protein arginine methyltransferase 5. MEP50 serves as a coactivator for both the androgen receptor and estrogen receptor in humans. Mep50 plays a crucial role in the development of germ cells in Drosophila. The precise role of Mep50 in oogenesis remains unclear in vertebrates. The objective of this study was to investigate the role of Mep50 in oogenesis in medaka fish. Disruption of Mep50 resulted in impaired oogenesis and the formation of multiple oocyte follicles in medaka. RNA-seq analysis revealed significant differential gene expression in the mutant ovary, with 4542 genes up-regulated and 1264 genes down-regulated. The regulated genes were found to be enriched in cellular matrices and ECM-receptor interaction, the Notch signaling pathway, the PI3K-Akt signaling pathway, the MAPK signaling pathway, the Hippo signaling pathway, and the Jak-Stat pathway, among others. In addition, the genes related to the hypothalamus-pituitary-gonad axis, steroid metabolism, and IGF system were impacted. Furthermore, the mutation of mep50 caused significant alterations in alternative splicing of pre-mRNA in ovarian cells. Quantitative RT-PCR results validated the findings from RNA-seq analysis in the specific genes, including akt2, map3k5, yap1, fshr, cyp17a, igf1, ythdc2, cdk6, and col1, among others. The findings of this study demonstrate that Mep50 plays a crucial role in oogenesis, participating in a diverse range of biological processes such as steroid metabolism, cell matrix regulation, and signal pathways. This may be achieved through the regulation of gene expression via mRNA splicing in medaka ovarian cells.

Multi-omic profiling of simultaneous ductal carcinoma in situ and invasive breast cancer

PURPOSE

The progression of ductal carcinoma in situ (DCIS) to invasive breast carcinoma (IBC) in humans is highly variable. To better understand the relationship between them, we performed a multi-omic characterization of co-occurring DCIS and IBC lesions in a cohort of individuals.

METHODS

Formalin-fixed paraffin-embedded tissue samples from 50 patients with co-occurring DCIS and IBC lesions were subjected to DNA-seq and whole transcriptome RNA-seq. Paired DCIS and IBC multi-omics profiles were then interrogated for DNA mutations, gene expression profiles and pathway analysis.

RESULTS

Most small variants and copy number variations were shared between co-occurring DCIS and IBC lesions, with IBC exhibiting on average a higher degree of additional mutations. However, 36% of co-occurring lesions shared no common mutations and 49% shared no common copy number variations. The most frequent genomic variants in both DCIS and IBC were PIK3CA, TP53, KMT2C, MAP3K1, GATA3 and SF3B1, with KMT2C being more frequent in DCIS and TP53 and MAP3K1 more frequent in IBC, though the numbers are too small for definitive conclusions. The most frequent copy number variations were seen in MCL1, CKSB1 and ERBB2. ERBB2 changes were not seen in IBC unless present in the corresponding DCIS. Transcriptional profiles were highly distinct between DCIS and IBC, with DCIS exhibiting upregulation of immune-related signatures, while IBC showed significant overexpression in genes and pathways associated with cell division and proliferation. Interestingly, DCIS and IBC exhibited significant differential expression of different components of extracellular matrix (ECM) formation and regulation, with DCIS showing overexpression of ECM-membrane interaction components while IBC showed upregulation of genes associated with fibronectin and invadopodia.

CONCLUSION

While most co-occurring DCIS and IBC were mutationally similar and suggestive of a common clonal progenitor, transcriptionally the lesions are highly distinct, with IBC expressing key pathways that facilitate invasion and proliferation. These results are suggestive of additional levels of regulation, epigenetic or other, that facilitate the acquisition of invasive properties during tumor evolution.

Modulation of macrophage polarization by secondary cross-linked hyaluronan-dopamine hydrogels

The inflammatory response plays a critical role in standard tissue repair processes, wherein active modulation of macrophage polarization is necessary for wound healing. Dopamine, a mussel-inspired bioactive material, is widely involved in wound healing, neural/bone/myocardial regeneration, and more. Recent studies indicated that dopamine-modified biomaterials can potentially alter macrophages polarization towards a pro-healing phenotype, thereby enhancing tissue regeneration. Nevertheless the immunoregulatory activity of dopamine on macrophage polarization remains unclear. This study introduces a novel interpenetrating hydrogel to bridge this research gap. The hydrogel, combining varying concentrations of oxidized dopamine with hyaluronic acid hydrogel, allows precise regulation of mechanical properties, antioxidant bioactivity, and biocompatibility. Surprisingly, both in vivo and in vitro outcomes demonstrated that dopamine concentration modulates macrophage polarization, but not linearly. Lower concentration (2 mg/mL) potentially decrease inflammation and facilitate M2 type macrophage polarization. In contrast, higher concentration (10 mg/mL) exhibited a pro-inflammatory tendency in the late stages of implantation. RNA-seq analysis revealed that lower dopamine concentrations induced the M1/M2 transition of macrophages by modulating the NF-κB signaling pathway. Collectively, this research offers valuable insights into the immunoregulation effects of dopamine-integrated biomaterials in tissue repair and regeneration.

Generation and characterization of cytochrome P450 3A74 CRISPR/Cas9 knockout bovine foetal hepatocyte cell line (BFH12)

In human, the cytochrome P450 3A (CYP3A) subfamily of drug-metabolizing enzymes (DMEs) is responsible for a significant number of phase I reactions, with the CYP3A4 isoform superintending the hepatic and intestinal metabolism of diverse endobiotic and xenobiotic compounds. The CYP3A4-dependent bioactivation of chemicals may result in hepatotoxicity and trigger carcinogenesis. In cattle, four CYP3A genes (CYP3A74, CYP3A76, CYP3A28 and CYP3A24) have been identified. Despite cattle being daily exposed to xenobiotics (e.g., mycotoxins, food additives, drugs and pesticides), the existing knowledge about the contribution of CYP3A in bovine hepatic metabolism is still incomplete. Nowadays, CRISPR/Cas9 mediated knockout (KO) is a valuable method to generate in vivo and in vitro models for studying the metabolism of xenobiotics. In the present study, we successfully performed CRISPR/Cas9-mediated KO of bovine CYP3A74, human CYP3A4-like, in a bovine foetal hepatocyte cell line (BFH12). After clonal expansion and selection, CYP3A74 ablation was confirmed at the DNA, mRNA, and protein level. The subsequent characterization of the CYP3A74 KO clone highlighted significant transcriptomic changes (RNA-sequencing) associated with the regulation of cell cycle and proliferation, immune and inflammatory response, as well as metabolic processes. Overall, this study successfully developed a new CYP3A74 KO in vitro model by using CRISPR/Cas9 technology, which represents a novel resource for xenobiotic metabolism studies in cattle. Furthermore, the transcriptomic analysis suggests a key role of CYP3A74 in bovine hepatocyte cell cycle regulation and metabolic homeostasis.

Bifenthrin induces changes in clinical poisoning symptoms, oxidative stress, DNA damage, histological characteristics, and transcriptome in Chinese giant salamander (Andrias davidianus) larvae

Bifenthrin (BF) is a broad-spectrum insecticide that has gained widespread use due to its high effectiveness. However, there is limited research on the potential toxic effects of bifenthrin pollution on amphibians. This study aimed to investigate the 50 % lethal concentration (LC50) and safety concentration of Chinese giant salamanders (CGS) exposed to BF (at 0, 6.25,12.5,25 and 50 μg/L BF) for 96 h. Subsequently, CGS were exposed to BF (at 0, 0.04, and 4 μg/L BF) for one week to investigate its toxic effects. Clinical poisoning symptoms, liver pathology, oxidative stress factors, DNA damage, and transcriptome differences were observed and analyzed. The results indicate that exposure to BF at 4 μg/L significantly decreased the adenosine-triphosphate (ATP), superoxide dismutase (SOD), glutathione (GSH), and catalase (CAT) contents in the brain, liver, and kidney of CGS. Additionally, the study found that the malondialdehyde (MDA), reactive oxygen species (ROS), and 8-hydroxydeoxyguanosine (8-OHdG) contents were increased. The liver tissue exhibited significant inflammatory reactions and structural malformations. RNA-seq analysis of the liver showed that BF caused abnormal antioxidant indices of CGS. This affected molecular function genes such as catalytic activity, ATP-dependent activity, metabolic processes, signaling and immune system processes, behavior, and detoxification, which were significantly upregulated, resulting in the differential genes significantly enriched in the calcium signaling pathway, PPARα signaling pathway and NF-kB signaling pathway. The results suggest that BF induces the abnormal production of free radicals, which overwhelms the body's self-defense system, leading to varying degrees of oxidative stress. This can result in oxidative damage, DNA damage, abnormal lipid metabolism, autoimmune diseases, clinical poisoning symptoms, and tissue inflammation. This work provides a theoretical basis for the rational application of bifenthrin and environmental risk assessment, as well as scientific guidance for the conservation of amphibian populations.

Transcriptomic analysis identifies dysregulated pathways and therapeutic targets in PMM2-CDG

PMM2-CDG (MIM # 212065), the most common congenital disorder of glycosylation, is caused by the deficiency of phosphomannomutase 2 (PMM2). It is a multisystemic disease of variable severity that particularly affects the nervous system; however, its molecular pathophysiology remains poorly understood. Currently, there is no effective treatment. We performed an RNA-seq based transcriptomic study using patient-derived fibroblasts to gain insight into the mechanisms underlying the clinical symptomatology and to identify druggable targets. Systems biology methods were used to identify cellular pathways potentially affected by PMM2 deficiency, including Senescence, Bone regulation, Cell adhesion and Extracellular Matrix (ECM) and Response to cytokines. Functional validation assays using patients' fibroblasts revealed defects related to cell proliferation, cell cycle, the composition of the ECM and cell migration, and showed a potential role of the inflammatory response in the pathophysiology of the disease. Furthermore, treatment with a previously described pharmacological chaperone reverted the differential expression of some of the dysregulated genes. The results presented from transcriptomic data might serve as a platform for identifying therapeutic targets for PMM2-CDG, as well as for monitoring the effectiveness of therapeutic strategies, including pharmacological candidates and mannose-1-P, drug repurposing.

Hippocampal transcriptomic analyses reveal the potential antiapoptotic mechanism of a novel anticonvulsant agent Q808 on pentylenetetrazol-induced epilepsy in rats

Brain apoptosis is one of the main causes of epileptogenesis. The antiapoptotic effect and potential mechanism of Q808, an innovative anticonvulsant chemical, have never been reported. In this study, the seizure stage and latency to reach stage 2 of pentylenetetrazol (PTZ) seizure rat model treated with Q808 were investigated. The morphological change and neuronal apoptosis in the hippocampus were detected by hematoxylin and eosin (HE) and terminal deoxynucleotidyl transferase-mediated dUTP nick end-labeling (TUNEL) staining, respectively. The hippocampal transcriptomic changes were observed using RNA sequencing (RNA-seq). The expression levels of hub genes were verified by quantitative reverse-transcription PCR (qRT-PCR). Results revealed that Q808 could allay the seizure score and prolong the stage 2 latency in seizure rats. The morphological changes of neurons and the number of apoptotic cells in the DG area were diminished by Q808 treatment. RNA-seq analysis revealed eight hub genes, including Map2k3, Nfs1, Chchd4, Hdac6, Siglec5, Slc35d3, Entpd1, and LOC103690108, and nine hub pathways among the control, PTZ, and Q808 groups. Hub gene Nfs1 was involved in the hub pathway sulfur relay system, and Map2k3 was involved in the eight remaining hub pathways, including Amyotrophic lateral sclerosis, Cellular senescence, Fc epsilon RI signaling pathway, GnRH signaling pathway, Influenza A, Rap1 signaling pathway, TNF signaling pathway, and Toll-like receptor signaling pathway. qRT-PCR confirmed that the mRNA levels of these hub genes were consistent with the RNA-seq results. Our findings might contribute to further studies exploring the new apoptosis mechanism and actions of Q808.

Molecular, behavioral, and growth responses of juvenile yellow catfish (Tachysurus fulvidraco) exposed to carbamazepine

Carbamazepine (CBZ) is an anticonvulsant medication used to treat epilepsy and bipolar disorder. Due to its persistence and low removal rate in wastewater treatment plants, it is frequently detected in the environment, raising concerns regarding its potential adverse effects on aquatic organisms and ecosystems. In this study, we aimed to assess the impact of CBZ on the behavior and growth of juvenile yellow catfish Tachysurus fulvidraco, a native and economically important species in China. Fish were exposed to CBZ at three concentrations of 1, 10, or 100 µg/L for 14 days. The fish exposed to 10 and 100 μg/L of CBZ exhibited decreased feeding, and a significant increase in cannibalistic tendencies was observed in fish exposed to 100 μg/L CBZ. Acetylcholinesterase activity was increased in the brain of fish exposed to 100 μg/L CBZ. CBZ also inhibited the growth of yellow catfish. To better elucidate mechanisms of toxicity, transcriptomics was conducted in both the brain and liver. In the brain, gene networks associated with neurotransmitter dysfunction were altered by CBZ, as well as networks associated with mitochondrial dysfunction and metabolism. In the liver, gene networks associated with the immune system were altered by CBZ. The current study improves comprehension of the sub-lethal effects of CBZ and reveals novel insight into molecular and biochemical pathways disrupted by CBZ, identifying putative key events associated with reduced growth and altered behavior. This study emphasizes the necessity for improved comprehension of the effects of pharmaceutical contaminants on fish at environmentally relevant levels.

Genes involved in the cholecystokinin receptor signaling map were differentially expressed in the jejunum of steers with variation in residual feed intake

The jejunum is a critical site for nutrient digestion and absorption, and variation in its ability to take up nutrients within the jejunum is likely to affect feed efficiency. The purpose of this study was to determine differences in gene expression in the jejunum of beef steers divergent for residual feed intake (RFI) in one cohort of steers (Year 1), and to validate those genes in animals from a second study (Year 2). Steers from Year 1 (n = 16) were selected for high and low RFI. Jejunum mucosal tissue was obtained for RNA-seq. Thirty-two genes were differentially expressed (PFDR≤0.15), and five were over-represented in pathways including inflammatory mediator, cholecystokinin receptor (CCKR) signaling, and p38 MAPK pathways. Several differentially expressed genes (ALOX12, ALPI, FABP6, FABP7, FLT1, GSTA2, MEF2B, PDK4, SPP1, and TTF2) have been previously associated with RFI in other studies. Real-time qPCR was used to validate nine differentially expressed genes in the Year 1 steers used for RNA-seq, and in the Year 2 validation cohort. Six genes were validated as differentially expressed (P < 0.1) using RT-qPCR in the Year 1 population. In the Year 2 population, five genes displayed the same direction of expression as the Year 1 population and 3 were differentially expressed (P < 0.1). The CCKR pathway is involved in digestion, appetite control, and regulation of body weight making it a compelling candidate for feed efficiency in cattle, and the validation of these genes in a second population of cattle is suggestive of a role in feed efficiency.

Characterization of transcriptome changes in saline stress adaptation on Leuciscus merzbacheri using PacBio Iso-Seq and RNA-Seq

Leuciscus merzbacheri is a native fish species found exclusively in the Junggar Basin in Xinjiang. It exhibits remarkable adaptability, thriving in varying water conditions such as the saline waters, the semi-saline water, and the freshwater. Despite its significant economic and ecological value, the underlying mechanisms of its remarkable salinity tolerance remain elusive. Our study marks the first time the full-length transcriptome of L. merzbacheri has been reported, utilizing RNA-Seq and PacBio Iso-Seq technologies. We found that the average length of the full-length transcriptome is 1,780 bp, with an N50 length of 2,358 bp. We collected RNA-Seq data from gill, liver, and kidney tissues of L. merzbacheri from both saline water and freshwater environments and conducted comparative analyses across these tissues. Further analysis revealed significant enrichment in several key functional gene categories and signalling pathways related to stress response and environmental adaptation. The findings provide a valuable genetic resource for further investigation into saline-responsive candidate genes, which will deepen our understanding of teleost adaptation to extreme environmental stress. This knowledge is crucial for the future breeding and conservation of native fish species.