Transcriptome sequencing analysis of bovine mammary epithelial cells induced by lipopolysaccharide

Mastitis in cows is caused by the inflammation of the mammary glands due to an infection by external pathogenic bacteria. Mammary gland epithelial cells, which are in direct contact with the external environment, are responsible for the first line of defense of the mammary gland against pathogenic bacteria, playing an essential role in immune defense. To investigate the mechanism of bovine mammary epithelial cells in the inflammatory process, we treated the cells with LPS for 12 hours and analyzed the changes in mRNA by transcriptome sequencing. The results showed that compared to the control group, the LPS treatment group had 121 up-regulated genes and 18 down-regulated genes. GO and KEGG enrichment analysis revealed that these differential genes were mainly enriched in the IL-17 signaling pathway, Legionellosis, Cytokine-cytokine receptor interaction, NF-kappa B signaling pathway, and other signaling pathways. Furthermore, the expression of GRO1 and CXCL3 mRNAs increased significantly after LPS treatment. These findings provide new insights for the treatment of mastitis in cows in the future.

Microbial-transcriptome integrative analysis of heat stress effects on amino acid metabolism and lipid peroxidation in poultry jejunum

Despite the significant threat of heat stress to livestock animals, only a few studies have considered the potential relationship between broiler chickens and their microbiota. Therefore, this study examined microbial modifications, transcriptional changes and host-microbiome interactions using a predicted metabolome data-based approach to understand the impact of heat stress on poultry. After the analysis, the host functional enrichment analysis revealed that pathways related to lipid and protein metabolism were elevated under heat stress conditions. In contrast, pathways related to the cell cycle were suppressed under normal environmental temperatures. In line with the transcriptome analysis, the microbial analysis results indicate that taxonomic changes affect lipid degradation. Heat stress engendered statistically significant difference in the abundance of 11 microorganisms, including Bacteroides and Peptostreptococcacea. Together, integrative approach analysis suggests that microbiota-induced metabolites affect host fatty acid peroxidation metabolism, which is correlated with the gene families of Acyl-CoA dehydrogenase long chain (ACADL), Acyl-CoA Oxidase (ACOX) and Acetyl-CoA Acyltransferase (ACAA). This integrated approach provides novel insights into heat stress problems and identifies potential biomarkers associated with heat stress.

A comprehensive study on the longissius dorsi muscle of Ashdan yaks under different feeding regimes based on transcriptomic and metabolomic analyses

Yak is an important dominant livestock species at high altitude, and the growth performance of yak has obvious differences under different feeding methods. This experiment was conducted to compare the effects of different feeding practices on growth performance and meat quality of yaks through combined transcriptomic and metabolomic analyses. In terms of yak growth performance, compared with traditional grazing, in-house feeding can significantly improve the average daily weight gain, carcass weight and net meat weight of yaks; in terms of yak meat quality, in-house feeding can effectively improve the quality of yak meat. A combined transcriptomic and metabolomic analysis revealed 31 co-enriched pathways, among which arginine metabolism, proline metabolism and glycerophospholipid metabolism may be involved in the development of the longissimus dorsi muscle of yak and the regulation of meat quality-related traits. The experimental results increased our understanding of yak meat quality and provided data materials for subsequent deep excavation of the mechanism of yak meat quality.

Comparative transcriptome profile analysis of granulosa cells from porcine ovarian follicles during early atresia

At various stages of ovarian follicular development, more than 99% of follicles will be eliminated through a degenerative process called atresia. The regulatory mechanisms of atresia have been elucidated to some extent, involving hormones, growth factors, cytokines, and other factors. However, the stimuli initiating atresia in follicular granulosa cells remain unknown. In this study, we isolated the granulosa cells from porcine ovarian follicles (3-5 mm diameter) divided into healthy follicles (HFs) and early atretic follicles (EAFs). We applied high-throughput RNA sequencing to identify and compare differentially expressed genes (DEGs) between HFs and EAFs. A total of 31,694 genes were detected, of which 21,806 were co-expressed in six samples, and 243 genes (p < 0.05; FDR < 0.05) were differentially expressed (DEGs), including 123 downregulated and 120 upregulated in EAFs. GO analysis highlighted hormone metabolism, plasma membrane localization, and transporter activity. The pathway analysis indicated that 51 DEGs, involved in steroidogenesis, cell adhesion molecules, and TGF-beta signaling pathways, were highly related to atresia. Additionally, the interaction network of DEGs (p < 0.01; FDR < 0.05) using STRING highlighted LHR, ACACB, and CXCR4 as central nodes. In summary, this transcriptome analysis enriched our knowledge of the shifted mechanisms in granulosa cells during early atresia and provided novel perspectives into the atresia initiation.

Complete genome sequence and transcriptome response to vitamin C supplementation of Novacetimonas hansenii SI1 - producer of highly-stretchable cellulose

Novacetimonas hansenii SI1, previously known as Komagataeibacter hansenii, produces bacterial nanocellulose (BNC) with unique ability to stretch. The addition of vitamin C in the culture medium increases the porosity of the membranes and their stretchability making them highly moldable. To better understand the genetic background of this strain, we obtained its complete genome sequence using a hybrid sequencing and assembly strategy. We described the functional regions in the genome which are important for the synthesis of BNC and acetan-like II polymer. We next investigated the effect of 1% vitamin C supplementation on the global gene expression profile using RNA sequencing. Our transcriptomic readouts imply that vitamin C functions mainly as a reducing agent. We found that the changes in cellular redox status are balanced by strong repression of the sulfur assimilation pathway. Moreover, in the reduced conditions, glucose oxidation is decreased and alternative pathways for energy generation, such as acetate accumulation, are activated. The presence of vitamin C negatively influences acetan-like II polymer biosynthesis, which may explain the lowered yield and changed mechanical properties of BNC. The results of this study enrich the functional characteristics of the genomes of the efficient producers of the N. hansenii species. Improved understanding of the adaptation to the presence of vitamin C at the molecular level has important guiding significance for influencing the biosynthesis of BNC and its morphology.

Elucidating the hepatoprotective mechanisms of cholic acid against CCl4-Induced acute liver injury: A transcriptomic and metabolomic study

ETHNOPHARMACOLOGICAL RELEVANCE

Cholic acid (CA) is one of the main active ingredients in Calculus Bovis, a traditional Chinese medicine, which helps to regulate the heart and liver meridians, clearing the heart, opening the mouth, cooling the liver and calming the wind. However, the molecular mechanism of its liver protective effect is still unclear.

AIM OF THE STUDY

Growing attention has been directed towards traditional Chinese medicine (TCM), particularly Calculus Bovis, as a potential solution for liver protection. Despite this interest, a comprehensive understanding of its hepatoprotective mechanisms remains lacking. This research seeks to explore the potential protective properties of cholic acid (CA) against CCl4-induced acute liver injury (ALI) in mice, while also examining the mechanisms involved.

MATERIALS AND METHODS

In the experiment, a mouse model was employed to ALI using CCl4, and the potential therapeutic effects of orally administered CA at varying doses (15, 30, and 60 mg/kg) were assessed. The study employed a multi-faceted approach, integrating liver transcriptomics with serum metabolomics, and conducting thorough analyses of serum biochemical markers and liver histopathological sections.

RESULTS

Oral CA administration markedly reduced the organ indices of the liver, spleen, and thymus in comparison with the model group. It also elevated the expression of superoxide dismutase (SOD) in serum while diminishing the concentrations of ALT, AST, MDA, IL-6, and TNF-α. Moreover, CA ameliorated the pathological damage induced by CCl4. Integrated metabolomic and transcriptomic analyses indicated that the hepatoprotective action of CA on ALI is mediated through the modulation of lipid metabolic pathways-specifically, metabolisms of glycerophospholipid, arachidonic acid, as well as linoleic acid-and by altering the expression of genes such as Ptgr1, PLpp1, Tbxas1, and Cyp2c37.

CONCLUSIONS

The current investigation offers insights into the hepatoprotective mechanisms by which CA mitigates ALI caused by CCl4 exposure, thus supporting the further evaluation and development of CA-based therapeutics for ALI.

RNA-seq transcriptome and pathway analysis of the medicinal mushroom Lignosus tigris (Polyporaceae) offer insights into its bioactive compounds with anticancer and antioxidant potential

ETHNOPHARMACOLOGICAL RELEVANCE

Medicinal mushrooms belonging to the Lignosus spp., colloquially known as Tiger Milk mushrooms (TMMs), are used as traditional medicine by communities across various regions of China and Southeast Asia to enhance immunity and to treat various diseases. At present, three Lignosus species have been identified in Malaysia: L. rhinocerus, L. tigris, and L. cameronensis. Similarities in their macroscopic morphologies and the nearly indistinguishable appearance of their sclerotia often lead to interchangeability between them. Hence, substantiation of their traditional applications via identification of their individual bioactive properties is imperative in ensuring that they are safe for consumption. L. tigris was first identified in 2013. Thus far, studies on L. tigris cultivar sclerotia (Ligno TG-K) have shown that it possesses significant antioxidant activities and has greater antiproliferative action against selected cancer cells in vitro compared to its sister species, L. rhinocerus TM02®. Our previous genomics study also revealed significant genetic dissimilarities between them. Further omics investigations on Ligno TG-K hold immense potential in facilitating the identification of its bioactive compounds and their associated bioactivities.

AIM OF STUDY

The overall aim of this study was to investigate the gene expression profile of Ligno TG-K via de novo RNA-seq and pathway analysis. We also aimed to identify highly expressed genes encoding compounds that contribute to its cytotoxic and antioxidant properties, as well as perform a comparative transcriptomics analysis between Ligno TG-K and its sister species, L. rhinocerus TM02®.

MATERIALS AND METHODS

Total RNA from fresh 3-month-old cultivated L. tigris sclerotia (Ligno TG-K) was extracted and analyzed via de novo RNA sequencing. Expressed genes were analyzed using InterPro and NCBI-Nr databases for domain identification and homology search. Functional categorization based on gene functions and pathways was performed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Clusters of Orthologous Genes (COG) databases. Selected genes were subsequently subjected to phylogenetic analysis.

RESULTS

Our transcriptomics analysis of Ligno TG-K revealed that 68.06% of its genes are expressed in the sclerotium; 80.38% of these were coding transcripts. Our analysis identified highly expressed transcripts encoding proteins with prospective medicinal properties. These included serine proteases (FPKM = 7356.68), deoxyribonucleases (FPKM = 3777.98), lectins (FPKM = 3690.87), and fungal immunomodulatory proteins (FPKM = 2337.84), all of which have known associations with anticancer activities. Transcripts linked to proteins with antioxidant activities, such as superoxide dismutase (FPKM = 1161.69) and catalase (FPKM = 1905.83), were also highly expressed. Results of our sequence alignments revealed that these genes and their orthologs can be found in other mushrooms. They exhibit significant sequence similarities, suggesting possible parallels in their anticancer and antioxidant bioactivities.

CONCLUSION

This study is the first to provide a reference transcriptome profile of genes expressed in the sclerotia of L. tigris. The current study also presents distinct COG profiles of highly expressed genes in Ligno TG-K and L. rhinocerus TM02®, highlighting that any distinctions uncovered may be attributed to their interspecies variations and inherent characteristics that are unique to each species. Our findings suggest that Ligno TG-K contains bioactive compounds with prospective medicinal properties that warrant further investigations.

CLASSIFICATION

Systems biology and omics.

AtERF13 and AtERF6 double knockout fine-tunes growth and the transcriptome to promote cadmium tolerance in Arabidopsis

The toxic heavy metal cadmium (Cd) restricts plant growth. However, how plants fine-tune their growth to modulate Cd resistance has not been determined. Ethylene response factors (ERFs) are key regulators of Cd stress, and Arabidopsis thaliana ERF13 and ERF6 (AtERF13 and AtERF6) negatively regulate growth. We previously demonstrated that AtERF13 is a transcriptional activator that binds a Cd-responsive element. Herein, we report that Arabidopsis plants improve Cd tolerance by repressing AtERF13 and AtERF6. We found that AtERF13 and AtERF6 were strongly downregulated by Cd stress and that AtERF6 weakly bound Cd-responsive elements. Moreover, AtERF13 physically interacted with AtERF6. Importantly, AtERF13 and AtERF6 double knockout mutants, but not single mutants or overexpression lines, grew better, tolerated more Cd and had higher Cd contents than did the wild type. Comparative transcriptome analysis revealed that the double mutants regulate the defense response to cope with Cd toxicity. Accordingly, we propose that, upon Cd stress, Arabidopsis plants repress AtERF13 and AtERF6 to relieve their growth inhibition effects and adjust the transcriptome to adapt to Cd stress, leading to increased Cd tolerance. Our findings thereby provide deep mechanical insights into how dual-function transcription factors fine-tune growth and the transcriptome to promote Cd tolerance in plants.

Transcriptome analysis of the synergistic mechanisms between two strains of potato virus Y in Solanum tuberosum L

Many viruses employ a process known as superinfection exclusion (SIE) to block subsequent entry or replication of the same or closely related viruses in the cells they occupy. SIE is also referred to as Cross-protection refers to the situation where a host plant infected by a mild strain of a virus or viroid gains immunity against a more severe strain closely related to the initial infectant. The mechanisms underlying cross-protection are not fully understood. In this study, we performed a comparative transcriptomic analysis of potato (Solanum tuberosum L.) leaves. The strains PVYN-Wi-HLJ-BDH-2 and PVYNTN-NW-INM-W-369-12 are henceforth designated as BDH and 369, respectively. In total, 806 differentially expressed genes (DEGs) were detected between the Control and JZ (preinfected with BDH and challenge with 369) treatment. Gene Ontology (GO) analysis showed that the response to external biological stimulation, signal transduction, kinase, immunity, redox pathways were significantly enriched. Among these pathways, we identified numerous differentially expressed metabolites related to virus infection. Moreover, our data also identified a small set of genes that likely play important roles in the establishment of cross-protection. Specifically, we observed significant differential expression of the A1-II gamma-like gene, elongation factor 1-alpha-like gene, and subtilisin-like protease StSBT1.7 gene, with StSBT1.7 being the most significant in our transcriptome data. These genes can stimulate the expression of defense plant genes, induce plant chemical defense, and participate in the induction of trauma and pathogenic bacteria. Our findings provided insights into the mechanisms underlying the ability of mild viruses to protect host plants against subsequent closely related virus infection in Solanum tuberosum L.

Early statin exposure influences cardiac and skeletal development with implications for ion channel transcriptomes in zebrafish

Statins, widely prescribed for cholesterol management by inhibiting HMG-CoA reductase in the cholesterol biosynthesis pathway, may also influence vertebrate development. In this study, we investigated the developmental effects of two widely used statins, atorvastatin (ATO) and pravastatin (PRA), on zebrafish offspring. For ATO, we administered doses classified as low (1 μM), medium (5 μM), and high (10 μM), while for PRA, the corresponding concentrations were set at low (18 μM), medium (180 μM), and high (270 μM). Our results showed significant reductions in birth and hatching rates, along with decreased body length in offspring at all ATO concentrations and medium to high PRA concentrations. A notable increase in malformation rates, especially in the spine and heart, was observed across all ATO treatments and in medium and high PRA groups. Additionally, we observed reduced heart contraction rates, decreased heart size, lower bone volumes, and diminished expression of mRNA osteogenic markers. Elevated venous sinus-artery bulb (SV-BA) ratios, increased thoracic area, and abnormal cartilage development were also prominent in all ATO-treated groups. Transcriptome analysis revealed alterations in genes predominantly associated with ion channels. These findings provide insights into the potential impacts of specific concentrations of statins on offspring development and highlight potential gene interactions with statins.

Translatability of mouse muscle-aging for humans: the role of sex

Muscle-aging drives sarcopenia and is a major public health issue. Mice are frequently used as a model for human muscle-aging, however, research investigating their translational value is limited. In addition, mechanisms underlying muscle-aging may have sex-specific features in humans, but it is not yet assessed whether these are recapitulated in mice. Here, we studied the effects of aging on a functional, histological and transcriptional level at multiple timepoints in male and female mice (4, 17, 21 and 25 months), with particular emphasis on sex-differences. The effects of natural aging on the transcriptome of quadriceps muscle were compared to humans on pathway level. Significant loss of muscle mass occurred late, at 25 months, in both male (-17%, quadriceps) and female mice (-10%, quadriceps) compared to young control mice. Concomitantly, we found in female, but not male mice, a slower movement speed in the aged groups compared to the young mice (P < 0.001). Consistently, weighted gene co-expression network analysis revealed a stronger association between the aging-related reduction of movement and aging-related changes in muscle transcriptome of female compared to male mice (P < 0.001). In male, but not female mice, major distinctive aging-related changes occurred in the last age group (25 months), which highlights the necessity for careful selection of age using mice as a muscle-aging model. Furthermore, contrasting to humans, more aging-related changes were found in the muscle transcriptome of male mice compared to female mice (4090 vs. 2285 differentially expressed genes at 25 months, respectively). Subsequently, male mice recapitulated more muscle-aging related pathways characteristic for both male and female humans. In conclusion, our data show that sex has a critical effect on the mouse muscle-aging trajectory, although these do not necessarily reflect sex differences observed in the human muscle-aging trajectory.

Transcriptomic analysis to identify genes associated with hypothalamus vulnerability in aging mice with cognitive decline

The normal aging process is accompanied by cognitive decline, and previous studies have indicated the crucial role of the hypothalamus in regulating both aging and cognition. However, the precise molecular mechanism underlying this relationship remains unclear. Therefore, this present study aimed to identify potential predictors of cognitive decline associated with aging specifically within the hypothalamus. To achieve this, we employed Morris water maze (MWM) testing to assess learning and memory differences between young and aged mice. Additionally, transcriptome sequencing was conducted on the hypothalamus of young and aged mice to identify potential genes. Subsequently, GO and KEGG analyses were performed to investigate the functions of differentially expressed genes (DEGs) and their associated biological pathways. Finally, the results obtained from sequencing analysis were further validated using qRT-PCR. Notably, MWM testing revealed a significant decrease in spatial learning and memory ability among aged mice. According to KEGG analysis, the DEGs primarily encompassed various biochemical signaling pathways related to immune system (e.g., C3; C4b; Ccl2; Ccl7; Cebpb; Clec7a; Col3a1; Cxcl10; Cxcl2; Fosb; Fosl1; Gbp5; H2-Ab1; Hspa1a; Hspa1b; Icam1; Il1b; Itga5; Itgax; Lilrb4a; Plaur; Ptprc; Serpine1; Tnfrsf10b; Tnfsf10), neurodegenerative disease (e.g., Atp2a1; Creb5; Fzd10; Hspa1a; Hspa1b; Il1b; Kcnj10; Nxf3; Slc6a3; Tubb6; Uba1y; Wnt9b), nervous system function (e.g., Chrna4; Chrna6; Creb5; Slc6a3),and aging (e.g., Creb5; Hspa1a; Hspa1b) among others. These identified genes may serve as potential predictors for cognitive function in elderly individuals and will provide a crucial foundation for further exploration into the underlying molecular mechanisms.

Integrating high-throughput phenotypic profiling and transcriptomic analyses to predict the hepatosteatosis effects induced by per- and polyfluoroalkyl substances

Per- and polyfluoroalkyl substances (PFAS) is a large compound class (n > 12,000) that is extensively present in food, drinking water, and aquatic environments. Reduced serum triglycerides and hepatosteatosis appear to be the common phenotypes for different PFAS chemicals. However, the hepatosteatosis potential of most PFAS chemicals remains largely unknown. This study aims to investigate PFAS-induced hepatosteatosis using in vitro high-throughput phenotype profiling (HTPP) and high-throughput transcriptomic (HTTr) data. We quantified the in vitro hepatosteatosis effects and mitochondrial damage using high-content imaging, curated the transcriptomic data from the Gene Expression Omnibus (GEO) database, and then calculated the point of departure (POD) values for HTPP phenotypes or HTTr transcripts, using the Bayesian benchmark dose modeling approach. Our results indicated that PFAS compounds with fully saturated C-F bonds, sulfur- and nitrogen-containing functional groups, and a fluorinated carbon chain length greater than 8 have the potential to produce biological effects consistent with hepatosteatosis. PFAS primarily induced hepatosteatosis via disturbance in lipid transport and storage. The potency rankings of PFAS compounds are highly concordant among in vitro HTPP, HTTr, and in vivo hepatosteatosis phenotypes (ρ = 0.60-0.73). In conclusion, integrating the information from in vitro HTPP and HTTr analyses can accurately project in vivo hepatosteatosis effects induced by PFAS compounds.

Transcriptomic analysis of Vibrio alginolyticus challenged by Rhizoma coptidis reveals mechanisms of virulence genes

Rhizoma coptidis, a Chinese herbal medicine widely used to treat various bacterial infections, has the potential to develop antibiotic substitutes to overcome the drug resistance of Vibrio alginolyticus. To study the inhibitory effect of R. coptidis on V. alginolyticus, we sequenced the transcriptomes of three groups of samples of wild-type V. alginolyticus (CK) and V. alginolyticus, which were stressed by 5 mg/mL R. coptidis for 2 h (RC_2 h) and 4 h (RC_4 h). CK was compared with RC_2 h and RC_4 h, respectively, and a total of 1565 differentially expressed genes (DEGs) (988 up-regulated and 577 down-regulated) and 1737 DEGs (1152 up-regulated and 585 down-regulated) were identified. Comparing RC_2 h with RC_4 h, 156 DEGs (114 up-regulated and 42 down-regulated) were identified. The ability of biofilm formation and motility of V. alginolyticus altered upon with different concentrations of R. coptidis. Interestingly, relative expression patterns of virulence genes appeared statistically significantly varied, upon different concentrations of R. coptidis extract. DEGs were annotated to the Gene Ontology (GO) database for function enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the results showed that the main enriched pathways, was those related to the virulence of V. alginolyticus. This study provides a new perspective for understanding the complex pathogenic mechanism of V. alginolyticus. R. coptidis could potnetially be used as alternative or complimnetary to antibiotics to treat infections after further research.

Transcriptome sequencing analysis reveals miR-30c-5p promotes ferroptosis in cervical cancer and inhibits growth and metastasis of cervical cancer xenografts by targeting the METTL3/KRAS axis

Cervical cancer is the most common malignant tumor in the female reproductive system worldwide, and its molecular mechanisms remain complex and poorly understood. Various techniques, including transcriptome sequencing, RT-qPCR, ELISA, immunofluorescence, Western blot, CCK-8 assay, Transwell assay, and xenograft models, were employed to investigate gene/miRNA expression, cellular proliferation, migration, and the interactions between miR-30c-5p, METTL3, and KRAS. Our transcriptome sequencing results demonstrated a significant downregulation of miR-30c-5p in cervical cancer cells. Further investigations using RNA pull-down, dual-luciferase reporter assay, Me-RIP, and PAR-CLIP confirmed METTL3 as one of the downstream targets of miR-30c-5p, while KRAS was identified as an iron-death suppressor gene susceptible to m6A modification. Notably, our Me-RIP analysis demonstrated the involvement of METTL3 in m6A modification of KRAS. In vitro experiments revealed that miR-30c-5p facilitated ferroptosis in cervical cancer cells by inhibiting the METTL3/KRAS axis, thus suppressing proliferation and migration. Additionally, in vivo studies demonstrated that miR-30c-5p repressed the growth and metastasis of cervical cancer xenografts through the inhibition of the METTL3/KRAS axis. Overall, this study highlights the critical role of miR-30c-5p in modulating cervical cancer progression by targeting the METTL3/KRAS axis, providing new insights into the molecular mechanisms underlying cervical cancer growth and metastasis.

Liver transcriptome analysis reveal the metabolic and apoptotic responses of Trachinotus ovatus under acute cold stress

Trachinotus ovatus is an economically important fish and has been recommended as a high-quality aquaculture fish breed for the high-quality development of sea ranches in the South China Sea. However, T. ovatus shows intolerance to low temperature, greatly limiting the extension of farming scale, reducing production efficiency in winter, and increasing farming risks. In this study, liver transcriptome analysis was investigated in T. ovatus under acute low temperature conditions (20 and 15 °C) using RNA sequencing (RNA-Seq) technology. Inter-groups differential expression analysis and trend analysis screened 1219 DEGs and four significant profiles (profiles 0, 3, 4, and 7), respectively. GO enrichment analysis showed that these DEGs were mainly related to metabolic process and cell growth and death process. KEGG enrichment analysis found that DEGs were mainly associated with lipid metabolism, carbohydrate metabolism, and cell growth and death, such as gluconeogenesis, glycolysis, fatty acid oxidation, cholesterol biosynthesis, p53 signaling pathway, cell cycle arrest, and apoptotic cell death. Moreover, protein-protein interaction networks identified two hub genes (FOS and JUNB) and some important genes related to metabolic process and cell growth and death process, that corresponding to enrichment analysis. Overall, gluconeogenesis, lipid mobilization, and fatty acid oxidation in metabolic process and cell cycle arrest and apoptotic cell death in cell growth and death process were enhanced, while glycolysis, liver glycogen synthesis and cholesterol biosynthesis in metabolic process were inhibited. The enhancement or attenuatment of metabolic process and cell growth and death process is conducive to maintain energy balance, normal fluidity of cell membrane, normal physiological functions of liver cell, enhancing the tolerance of T. ovatus to cold stress. These results suggested that metabolic process and cell growth and death process play important roles in response to acute cold stress in the liver of T. ovatus. Gene expreesion level analysis showed that acute cold stress at 15 °C was identified as a critical temperature point for T. ovatus in term of cellular metabolism alteration and apoptosis inducement, and rewarming intervention should be timely implemented above 15 °C. Our study can provide theoretical support for breeding cold-tolerant cultivars of T. ovatus, which is contributed to high-quality productions fish production.

Transcriptomic contributions to a modern cytoarchitectonic parcellation of the human cerebral cortex

Transcriptomic contributions to the anatomical, functional, and network layout of the human cerebral cortex (HCC) have become a major interest in cognitive and systems neuroscience. Here, we tested if transcriptomic differences support a modern, algorithmic cytoarchitectonic parcellation of HCC. Using a data-driven approach, we identified a sparse subset of genes that differentially contributed to the cytoarchitectonic parcellation of HCC. A combined metric of cortical thickness and myelination (CT/M ratio), as well as cell density, correlated with gene expression. Enrichment analyses showed that genes specific to the cytoarchitectonic parcellation of the HCC were related to molecular functions such as transmembrane transport and ion channel activity. Together, the relationship between transcriptomics and cytoarchitecture bridges the gap among (i) gradients at the macro-scale (including thickness and myelination), (ii) areas at the meso-scale, and (iii) cell density at the microscale, as well as supports the recently proposed cortical spectrum theory and structural model.

yunnanense King)

BACKGROUND

Yunnan hulled wheat grains (YHWs) have abundant phenolic compounds (PCs). However, a systematic elucidation of the phenolic characteristics and molecular basis in YHWs is currently lacking. The aim of the study, for the first time, was to conduct metabolomic and transcriptomic analyses of YHWs at different developmental stages.

RESULTS

A total of five phenolic metabolite classes (phenolic acids, flavonoids, quinones, lignans and coumarins, and tannins) and 361 PCs were identified, with flavonoids and phenolic acids being the most abundant components. The relative abundance of the identified PCs showed a dynamic decreasing pattern with grain development, and the most significant differences in accumulation were between the enlargement and mature stage, which is consistent with the gene regulation patterns of the corresponding phenolic biosynthesis pathway. Through co-expression and co-network analysis, PAL, HCT, CCR, F3H, CHS, CHI and bZIP were identified and predicted as candidate key enzymes and transcription factors.

CONCLUSION

The results broaden our understanding of PC accumulation in wheat whole grains, especially the differential transfer between immature and mature grains. The identified PCs and potential regulatory factors provide important information for future in-depth research on the biosynthesis of PCs and the improvement of wheat nutritional quality. © 2024 Society of Chemical Industry.

Whole transcriptome mapping reveals the lncRNA regulatory network of TFP5 treatment in diabetic nephropathy

BACKGROUND

TFP5 is a Cdk5 inhibitor peptide, which could restore insulin production. However, the role of TFP5 in diabetic nephropathy (DN) is still unclear.

OBJECTIVE

This study aims to characterize the transcriptome profiles of mRNA and lncRNA in TFP5-treated DN mice to mine key lncRNAs associated with TFP5 efficacy.

METHODS

We evaluated the role of TFP5 in DN pathology and performed RNA sequencing in C57BL/6J control mice, C57BL/6J db/db model mice, and TFP5 treatment C57BL/6J db/db model mice. The differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) were analyzed. WGCNA was used to screen hub-gene of TFP5 in treatment of DN.

RESULTS

Our results showed that TFP5 therapy ameliorated renal tubular injury in DN mice. In addition, compared with the control group, the expression profile of lncRNAs in the model group was significantly disordered, while TFP5 alleviated the abnormal expression of lncRNAs. A total of 67 DElncRNAs shared among the three groups, 39 DElncRNAs showed a trend of increasing in the DN group and decreasing after TFP treatment, while the remaining 28 showed the opposite trend. DElncRNAs were enriched in glycosphingolipid biosynthesis signaling pathways, NF-κB signaling pathways, and complement activation signaling pathways. There were 1028 up-regulated and 1117 down-regulated DEmRNAs in the model group compared to control group, and 123 up-regulated and 153 down-regulated DEmRNAs in the TFP5 group compared to the model group. The DEmRNAs were involved in PPAR and MAPK signaling pathway. We confirmed that MSTRG.28304.1 is a key DElncRNA for TFP5 treatment of DN. TFP5 ameliorated DN maybe by inhibiting MSTRG.28304.1 through regulating the insulin resistance and PPAR signaling pathway. The qRT-PCR results confirmed the reliability of the sequencing data through verifying the expression of ENSMUST00000211209, MSTRG.31814.5, MSTRG.28304.1, and MSTRG.45642.14.

CONCLUSION

Overall, the present study provides novel insights into molecular mechanisms of TFP5 treatment in DN.

Attenuated cell cycle and DNA damage response transcriptome signatures and overrepresented cell adhesion processes imply accelerated progression in patients with lower-risk myelodysplastic neoplasms

Patients with myelodysplastic neoplasms (MDS) are classified according to the risk of acute myeloid leukemia transformation. Some lower-risk MDS patients (LR-MDS) progress rapidly despite expected good prognosis. Using diagnostic samples, we aimed to uncover the mechanisms of this accelerated progression at the transcriptome level. RNAseq was performed on CD34+ ribodepleted RNA samples from 53 LR-MDS patients without accelerated progression (stMDS) and 8 who progressed within 20 months (prMDS); 845 genes were differentially expressed (ІlogFCІ > 1, FDR < 0.01) between these groups. stMDS CD34+ cells exhibited transcriptional signatures of actively cycling, megakaryocyte/erythrocyte lineage-primed progenitors, with upregulation of cell cycle checkpoints and stress pathways, which presumably form a tumor-suppressing barrier. Conversely, cell cycle, DNA damage response (DDR) and energy metabolism-related pathways were downregulated in prMDS samples, whereas cell adhesion processes were upregulated. Also, prMDS samples showed high levels of aberrant splicing and global lncRNA expression that may contribute to the attenuation of DDR pathways. We observed overexpression of multiple oncogenes and diminished differentiation in prMDS; the expression of ZEB1 and NEK3, genes not previously associated with MDS prognosis, might serve as potential biomarkers for LR-MDS progression. Our 19-gene DDR signature showed a significant predictive power for LR-MDS progression. In validation samples (stMDS = 3, prMDS = 4), the key markers and signatures retained their significance. Collectively, accelerated progression of LR-MDS appears to be associated with transcriptome patterns of a quiescent-like cell state, reduced lineage differentiation and suppressed DDR, inherent to CD34+ cells. The attenuation of DDR-related gene-expression signature may refine risk assessment in LR-MDS patients.

Host gene expression signatures to identify infection type and organ dysfunction in children evaluated for sepsis: a multicentre cohort study

BACKGROUND

Sepsis is defined as dysregulated host response to infection that leads to life-threatening organ dysfunction. Biomarkers characterising the dysregulated host response in sepsis are lacking. We aimed to develop host gene expression signatures to predict organ dysfunction in children with bacterial or viral infection.

METHODS

This cohort study was done in emergency departments and intensive care units of four hospitals in Queensland, Australia, and recruited children aged 1 month to 17 years who, upon admission, underwent a diagnostic test, including blood cultures, for suspected sepsis. Whole-blood RNA sequencing of blood was performed with Illumina NovaSeq (San Diego, CA, USA). Samples with completed phenotyping, monitoring, and RNA extraction by March 31, 2020, were included in the discovery cohort; samples collected or completed thereafter and by Oct 27, 2021, constituted the Rapid Paediatric Infection Diagnosis in Sepsis (RAPIDS) internal validation cohort. An external validation cohort was assembled from RNA sequencing gene expression count data from the observational European Childhood Life-threatening Infectious Disease Study (EUCLIDS), which recruited children with severe infection in nine European countries between 2012 and 2016. Feature selection approaches were applied to derive novel gene signatures for disease class (bacterial vs viral infection) and disease severity (presence vs absence of organ dysfunction 24 h post-sampling). The primary endpoint was the presence of organ dysfunction 24 h after blood sampling in the presence of confirmed bacterial versus viral infection. Gene signature performance is reported as area under the receiver operating characteristic curves (AUCs) and 95% CI.

FINDINGS

Between Sept 25, 2017, and Oct 27, 2021, 907 patients were enrolled. Blood samples from 595 patients were included in the discovery cohort, and samples from 312 children were included in the RAPIDS validation cohort. We derived a ten-gene disease class signature that achieved an AUC of 94·1% (95% CI 90·6-97·7) in distinguishing bacterial from viral infections in the RAPIDS validation cohort. A ten-gene disease severity signature achieved an AUC of 82·2% (95% CI 76·3-88·1) in predicting organ dysfunction within 24 h of sampling in the RAPIDS validation cohort. Used in tandem, the disease class and disease severity signatures predicted organ dysfunction within 24 h of sampling with an AUC of 90·5% (95% CI 83·3-97·6) for patients with predicted bacterial infection and 94·7% (87·8-100·0) for patients with predicted viral infection. In the external EUCLIDS validation dataset (n=362), the disease class and disease severity predicted organ dysfunction at time of sampling with an AUC of 70·1% (95% CI 44·1-96·2) for patients with predicted bacterial infection and 69·6% (53·1-86·0) for patients with predicted viral infection.

INTERPRETATION

In children evaluated for sepsis, novel host transcriptomic signatures specific for bacterial and viral infection can identify dysregulated host response leading to organ dysfunction.

FUNDING

Australian Government Medical Research Future Fund Genomic Health Futures Mission, Children's Hospital Foundation Queensland, Brisbane Diamantina Health Partners, Emergency Medicine Foundation, Gold Coast Hospital Foundation, Far North Queensland Foundation, Townsville Hospital and Health Services SERTA Grant, and Australian Infectious Diseases Research Centre.

Transcriptomic and Proteomic Characterization of the Immune Response to Elective Spinal Reconstructive Surgery: Impact of Aging and Comparison with Traumatic Injury Response

BACKGROUND

Major surgery triggers trauma-like stress responses linked to age, surgery duration, and blood loss, resembling polytrauma. This similarity suggests elective surgery as a surrogate model for studying polytrauma immune responses. We investigated stress responses across age groups and compared them with those of polytrauma patients.

STUDY DESIGN

Patients undergoing major spinal reconstruction surgery were divided into older (age >65 years, n = 5) and young (age 18 to 39 years, n = 6) groups. A comparison group consisted of matched trauma patients (n = 8). Blood samples were collected before, during, and after surgery. Bone marrow mononuclear cells and peripheral blood mononuclear cells were analyzed using cellular indexing of transcriptomes and epitopes sequencing or single-cell RNA sequencing. Plasma was subjected to dual-platform proteomic analysis (SomaLogic and O-link).

RESULTS

Response to polytrauma was highest within 4 hours. By comparison, the response to surgery was highest at 24 hours. Both insults triggered significant changes in cluster of differentiation 14 monocytes, with increased inflammation and lower major histocompatibility complex-class 2 expression. Older patient's cluster of differentiation 14 monocytes displayed higher inflammation and less major histocompatibility complex-class 2 suppression; a trend was also seen in bone marrow mononuclear cells. Although natural killer cells were markedly activated after polytrauma, they were suppressed after surgery, especially in older patients. In plasma, innate immunity proteins dominated at 24 hours, shifting to adaptive immunity proteins by 6 weeks with heightened inflammation in older patients. Senescence-associated secretory phenotype proteins were higher in older patients at baseline and further elevated during and after surgery.

CONCLUSIONS

Although both major surgery and polytrauma initiate immune and stress responses, substantial differences exist in timing and cellular profiles, suggesting major elective surgery is not a suitable surrogate for the polytrauma response. Nonetheless, distinct responses in young vs older patients highlight the utility of elective spinal in studying patient-specific factors affecting outcomes after major elective surgery.

Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

The immune cell transcriptome is modulated by vitamin D3 supplementation in people with a first demyelinating event participating in a randomized placebo-controlled trial

Vitamin D deficiency is a risk factor for developing multiple sclerosis. The PrevANZ trial was conducted to determine if vitamin D3 supplementation can prevent recurrent disease activity in people with a first demyelinating event. As a sub-study of this trial, we investigated the effect of supplementation on peripheral immune cell gene expression. Participants were randomized to 1000, 5000 or 10,000 international units daily of vitamin D3 or placebo. Peripheral blood was collected at baseline and 12 weeks and sent for ribonucleic acid sequencing. Datasets from 55 participants were included. Gene expression was modulated by high dose supplementation. Antigen presentation and viral response pathways were upregulated. Oxidative phosphorylation and immune signaling pathways, including tumor necrosis factor-alpha and interleukin-17 signaling, were downregulated. Overall, vitamin D3 supplementation for 12 weeks modulated the peripheral immune cell transcriptome with induction of anti-inflammatory gene expression profiles. Our results support a dose-dependent effect of vitamin D3 supplementation on immune gene expression.

PANoptosis-related signature in melanoma: Transcriptomic mapping and clinical prognostication

Programmed cell death plays a pivotal role in maintaining tissue homeostasis, and recent advancements in cell biology have uncovered PANoptosis-a novel paradigm integrating pyroptosis, apoptosis, and necroptosis. This study investigates the implications of PANoptosis in melanoma, a formidable skin cancer known for its metastatic potential and resistance to conventional therapies. Leveraging bulk and single-cell transcriptome analyses, machine learning modeling, and immune correlation assessments, we unveil the molecular intricacies of PANoptosis in melanoma. Single-cell sequencing identifies diverse cell types involved in PANoptosis, while bulk transcriptome analysis reveals key gene sets correlated with PANoptosis. Machine learning algorithms construct a robust prognostic model, demonstrating consistent predictive power across diverse cohorts. Patients with different cohorts can be divided into high-risk and low-risk groups according to this PANoptosis score, with the high-risk group having a significantly worse prognosis. Immune correlation analyses unveil a link between PANoptosis and immunotherapy response, with potential therapeutic implications. Mutation analysis and enrichment studies provide insights into the mutational landscape associated with PANoptosis. Finally, we used cell experiments to verify the expression and function of key gene PARVA, showing that PARVA was highly expressed in melanoma cell lines, and after PARVA is knocked down, cell invasion, migration, and colony formation ability were significantly decreased. This study advances our understanding of PANoptosis in melanoma, offering a comprehensive framework for targeted therapeutic interventions and personalized medicine strategies in combating this aggressive malignancy.

BioProtIS: Streamlining protein-ligand interaction pipeline for analysis in genomic and transcriptomic exploration

The identification of protein-ligand interactions plays a pivotal role in elucidating biological processes and discovering potential bioproducts. Harnessing the capabilities of computational methods in drug discovery, we introduce an innovative Inverted Virtual Screening (IVS) pipeline. This pipeline Integrated molecular dynamics and docking analyses to ensure that protein structures are not only energetically favorable but also representative of stable conformations. The primary objective of this pipeline is to automate and streamline the analysis of protein-ligand interactions at both genomic and transcriptomic scales. In the contemporary post-genomic era, high-throughput computational screening for bioproducts, biological systems, and therapeutic drugs has become a cornerstone practice. This approach offers the promise of cost-effectiveness, time efficiency, and optimization of laboratory work. Nevertheless, a notable deficiency persists in the availability of efficient pipelines capable of automating the virtual screening process, seamlessly integrating input and output, and leveraging the full potential of open-source tools. To bridge this critical gap, we have developed a versatile pipeline known as BioProtIS. This tool seamlessly integrates a suite of state-of-the-art tools, including Modeller, AlphaFold, Gromacs, FPOCKET, and AutoDock Vina, thus facilitating the streamlined docking of ligands with an expansive repertoire of proteins sourced from genomes and transcriptomes, and substrates. To assess the pipeline's performance, we employed the transcriptomes of Cereus jamacaru (a cactus species) and Aspisoma lineatum (firefly), along with the genome of Homo sapiens. This integration not only improves the accuracy of ligand-protein interactions by minimizing replicability deviations but also optimizes the discovery process by enabling the simultaneous evaluation of multiple substrates. Furthermore, our pipeline accommodates distinct testing scenarios, such as blind docking or site-specific targeting, which are invaluable in applications ranging from drug repositioning to the exploration of new allosteric binding sites and toxicity assessments. BioProtIS has been designed with modularity at its core. This inherent flexibility empowers users to make custom modifications directly within the source code, tailoring the pipeline to their specific research needs. Moreover, it lays the foundation for seamless integration of diverse docking algorithms in future iterations, promising ongoing advancements in the field of computational biology. This pipeline is available for free distribution and can be download at: https://github.com/BBMDO/BioProtIS.

Cellular and transcriptomic response to pathogenic and non-pathogenic Vibrio parahaemolyticus strains causing acute hepatopancreatic necrosis disease (AHPND) in Litopenaeus vannamei

The shrimp industry has historically been affected by viral and bacterial diseases. One of the most recent emerging diseases is Acute Hepatopancreatic Necrosis Disease (AHPND), which causes severe mortality. Despite its significance to sanitation and economics, little is known about the molecular response of shrimp to this disease. Here, we present the cellular and transcriptomic responses of Litopenaeus vannamei exposed to two Vibrio parahaemolyticus strains for 98 h, wherein one is non-pathogenic (VpN) and the other causes AHPND (VpP). Exposure to the VpN strain resulted in minor alterations in hepatopancreas morphology, including reductions in the size of R and B cells and detachments of small epithelial cells from 72 h onwards. On the other hand, exposure to the VpP strain is characterized by acute detachment of epithelial cells from the hepatopancreatic tubules and infiltration of hemocytes in the inter-tubular spaces. At the end of exposure, RNA-Seq analysis revealed functional enrichment in biological processes, such as the toll3 receptor signaling pathway, apoptotic processes, and production of molecular mediators involved in the inflammatory response of shrimp exposed to VpN treatment. The biological processes identified in the VpP treatment include superoxide anion metabolism, innate immune response, antimicrobial humoral response, and toll3 receptor signaling pathway. Furthermore, KEGG enrichment analysis revealed metabolic pathways associated with survival, cell adhesion, and reactive oxygen species, among others, for shrimp exposed to VpP. Our study proves the differential immune responses to two strains of V. parahaemolyticus, one pathogenic and the other nonpathogenic, enlarges our knowledge on the evolution of AHPND in L. vannamei, and uncovers unique perspectives on establishing genomic resources that may function as a groundwork for detecting probable molecular markers linked to the immune system in shrimp.

Unraveling the transcriptome profile of pulsed electromagnetic field stimulation in bone regeneration using a bioreactor-based investigation platform

INTRODUCTION

Human mesenchymal stem cells (hMSCs) sense and respond to biomechanical and biophysical stimuli, yet the involved signaling pathways are not fully identified. The clinical application of biophysical stimulation including pulsed electromagnetic field (PEMF) has gained momentum in musculoskeletal disorders and bone tissue engineering.

METHODOLOGY

We herein aim to explore the role of PEMF stimulation in bone regeneration by developing trabecular bone-like tissues, and then, culturing them under bone-like mechanical stimulation in an automated perfusion bioreactor combined with a custom-made PEMF stimulator. After selecting the optimal cell seeding and culture conditions for inspecting the effects of PEMF on hMSCs, transcriptomic studies were performed on cells cultured under direct perfusion with and without PEMF stimulation.

RESULTS

We were able to identify a set of signaling pathways and upstream regulators associated with PEMF stimulation and to distinguish those linked to bone regeneration. Our findings suggest that PEMF induces the immune potential of hMSCs by activating and inhibiting various immune-related pathways, such as macrophage classical activation and MSP-RON signaling in macrophages, respectively, while promoting angiogenesis and osteogenesis, which mimics the dynamic interplay of biological processes during bone healing.

CONCLUSIONS

Overall, the adopted bioreactor-based investigation platform can be used to investigate the impact of PEMF stimulation on bone regeneration.

Transcriptomic profile of human iPSC-derived podocyte-like cells exposed to a panel of xenobiotics

Podocytes play a critical role in the formation and maintenance of the glomerular filtration barrier and injury to these cells can lead to a breakdown of the glomerular barrier causing permanent damage leading to progressive chronic kidney disease. Matured podocytes have little proliferative potential, which makes them critical cells from a health perspective, but also challenging cells to maintain in vitro. Differentiating podocyte-like cells from induced pluripotent stem cells (iPSC) provides a novel and continuous source of cells. Here, we investigated the effect of a 24-h exposure to eight compounds, including the known glomerular toxins doxorubicin and pamidronate, on transcriptomic alterations in iPSC derived podocytes. Doxorubicin (50 nM), pamidronate (50 μM), sodium arsenite (10 μM), and cyclosporine A (15 μM) had a strong impact on the transcriptome, gentamicin (450 μg/ml), lead chloride (15 μM) and valproic acid (500 μM) had a mild impact and busulfan (50 μM) exhibited no impact. Gene alterations and pathways analysis provided mechanistic insight for example, doxorubicin exposure affected the p53 pathway and dedifferentiation, pamidronate activated several pathways including HIF1alpha and sodium arsenite up-regulated oxidative stress and metal responses. The results demonstrate the applicability of iPSC derived podocytes for toxicological and mechanistic investigations.

Spatial metabolomics and its application in the liver

Hepatocytes work in highly structured, repetitive hepatic lobules. Blood flow across the radial axis of the lobule generates oxygen, nutrient, and hormone gradients, which result in zoned spatial variability and functional diversity. This large heterogeneity suggests that hepatocytes in different lobule zones may have distinct gene expression profiles, metabolic features, regenerative capacity, and susceptibility to damage. Here, we describe the principles of liver zonation, introduce metabolomic approaches to study the spatial heterogeneity of the liver, and highlight the possibility of exploring the spatial metabolic profile, leading to a deeper understanding of the tissue metabolic organization. Spatial metabolomics can also reveal intercellular heterogeneity and its contribution to liver disease. These approaches facilitate the global characterization of liver metabolic function with high spatial resolution along physiological and pathological time scales. This review summarizes the state of the art for spatially resolved metabolomic analysis and the challenges that hinder the achievement of metabolome coverage at the single-cell level. We also discuss several major contributions to the understanding of liver spatial metabolism and conclude with our opinion on the future developments and applications of these exciting new technologies.

Transcriptome analysis provides new insights into the immune response of Ruditapes philippinarum infected with Vibrio alginolyticus

Manila clam (Ruditapes philippinarum) is a bivalve species with commercial value, but it is easily infected by pathogenic microorganisms in aquaculture, which restricts the shellfish industry. Notably, the impact of Vibrio alginolyticus on clam culture is obvious. In this study, RNA-seq was performed to analyze clam hepatopancreas tissue in 48 h (challenge group, G48h) and 96 h (challenge group, G96h) after infection with V. alginolyticus and 0 h after injection of PBS (control group, C). The results showed that a total of 1670 differentially expressed genes were detected in the G48h vs C group, and 1427 differentially expressed genes were detected in the G96h vs C group. In addition, KEGG analysis showed that DEGs were significantly enriched in pathways such as Lysosome and Mitophagy. Moreover, 15 immune related DEGs were selected for qRT-PCR analysis to verify the accuracy of RNA-seq, and the results showed that the expression level of DEGs was consistent with that of RNA-seq. Therefore, the results obtained in this study provides a preliminary understanding of the immune defense of R. philippinarum and molecular insights for genetic breeding of V. alginolyticus resistance in Manila clam.

The landscape of immune dysregulation in pediatric sepsis at a single-cell resolution

Recognizing immune dysregulation as a hallmark of sepsis pathophysiology, leukocytes have attracted major attention of investigation. While adult and pediatric sepsis are clinically distinct, their immunological delineation remains limited. Single cell technologies facilitated the characterization of immune signatures. We tackled to delineate immunological profiles of pediatric sepsis at a single-cell level by analyzing blood samples from six septic children, at both acute and recovery phases, and four healthy children. 16 single-cell transcriptomic datasets were analyzed and compared to adult sepsis dataset. We showed a unique shift in neutrophil subpopulations and functions between acute and recovery phases, along with the regulatory role of resistin. Neutrophil signatures were comparable between adult and pediatric sepsis. Innate-like CD4 T cells were predominantly and uniquely observed in acute phase of pediatric sepsis. Our study serves as a rich source of information about the phenotypic diversity and trajectory of circulating immune cells during pediatric sepsis.

Transcriptome meta-analysis-based identification of hub transcription factors and RNA-binding proteins potentially orchestrating gene regulatory cascades and crosstalk in response to abiotic stresses in Arabidopsis thaliana

Deteriorating climatic conditions and increasing human population necessitate the development of robust plant varieties resistant to harsh environments. Manipulation of regulatory proteins such as transcription factors (TFs) and RNA-binding proteins (RBPs) would be a beneficial strategy in this regard. Further, understanding the complex interconnections between different classes of regulatory molecules would be essential for the identification of candidate genes/proteins for trait improvement. Most studies to date have analysed the roles of TFs or RBPs individually, in conferring stress resilience. However, it would be important to identify dominant/upstream TFs and RBPs inducing widespread transcriptomic alterations through other regulators (i.e., other TFs/RBPs targeted by the upstream regulators). To this end, the present study employed a transcriptome meta-analysis and computational approaches to obtain a comprehensive overview of regulatory interactions. This work identified dominant TFs and RBPs potentially influencing stress-mediated differential expression of other regulators, which could in turn influence gene expression, and consequently, physiological responses. Twenty transcriptomic studies [related to (i) UV radiation, (ii) wounding, (iii) salinity, (iv) cold, and (v) drought stresses in Arabidopsis thaliana] were analysed for differential gene expression, followed by the identification of differentially expressed TFs and RBPs. Subsequently, other TFs and RBPs which could be influencing these regulators were identified, and their interaction networks and hub nodes were analysed. As a result, an interacting module of Basic Leucine Zipper (bZIP) family TFs as well as Heterogeneous nuclear ribonucleoproteins (hnRNP) and Glycine-rich protein (GRP) family RBPs (among other TFs and RBPs) were shown to potentially influence the stress-induced differential expression of other TFs and RBPs under all the considered stress conditions. Some of the identified hub TFs and RBPs are known to be of major importance in orchestrating stress-induced transcriptomic changes influencing a variety of physiological processes from seed germination to senescence. This study highlighted the gene/protein candidates that could be considered for multiplexed genetic manipulation - a promising approach to develop robust, multi-stress-resilient plant varieties.

Identifying transcriptomic indicators of tertiary treated municipal effluent in longnose dace (Rhinichthys cataractae) caged under semi-controlled conditions in experimental raceways

To evaluate effects of tertiary treated wastewater treatment plant effluent (MWWE) on transcriptomic responses in longnose dace (Rhinichthys cataractae; LND) we conducted a semi-controlled study in experimental raceways (Advancing Canadian Water Assets facility) imbedded in the Pine Creek treatment plant (Calgary, AB). LND collected from a reference site in the Bow River (REF) were caged in raceways containing either 5 % Pine Creek effluent (PC) or Bow River water (BR; control) over 28 d. Liver transcriptomes were analyzed in males and females sampled on days 7, 14 and 28 from BR and PC, and compared to REF fish on day 0. Concurrent with the caging, selected environmental substances of concern were analyzed in the BR and PC. Significantly different unigenes (SDUs) in females (vs males) within both BR and PC raceways increased over time and compared to REF fish. Moreover, SDUs in females and males within the same treatment (i.e., BR, PC) showed a temporal increase as well as compared to REF fish. Time was the dominant factor affecting SDUs, whereas sex and treatment had less of an impact on the transcriptome profiling. Gene Set Enrichment Analysis of BR vs PC over time revealed effects on genes involved in growth, metabolism of carbohydrates and lipids, and immune system on day 7; however, by day 28, 80-100 % of the transcripts localized to enriched biomarkers were associated with tissue immune responses in both sexes. Exposure to 5 % effluent had significant effects on female liver somatic index but no effects were observed on other phenotypic health indices in either sex. BR was used as the source of reference water, but analyses showed trace amounts of ESOCs. Analyses did not point towards definitive response patterns that could be used in field-based ecotoxicogenomic studies on the impacts of well-treated MWWE but suggested compromised adaptive immune responses.

Transcriptome analysis of primary adult B-cell lineage acute lymphoblastic leukemia identifies pathogenic variants and gene fusions, and predicts subtypes for in depth molecular diagnosis

BACKGROUND

B-cell acute lymphoblastic leukemia (B-ALL) is classified into subgroups based on known driver oncogenes and molecular lesions, including translocations and recurrent mutations. However, the current diagnostic tests do not identify subtypes or oncogenic lesions for all B-ALL samples, creating a heterogeneous B-ALL group of unknown subtypes.

METHODS

We sorted primary adult B-ALL cells and performed transcriptome analysis by bulk RNA sequencing (RNA-seq).

RESULTS

Transcriptomic analysis of an adult B-ALL cohort allowed the classification of four patient samples with subtypes that were not previously revealed by standard gene panels. The leukemia of two patients were of the DUX4 subtype and two were CRLF2+ Ph-like B-ALL. Furthermore, single nucleotide variant analysis detected the oncogenic NRAS-G12D, KRAS-G12D, and KRAS-G13D mutations in three of the patient samples, presenting targetable mutations. Additional oncogenic variants and gene fusions were uncovered, as well as multiple variants in the PDE4DIP gene across five of the patient samples.

CONCLUSION

We demonstrate that RNA-seq is an effective tool for precision medicine in B-ALL by providing comprehensive molecular profiling of leukemia cells, identifying subtype and oncogenic lesions, and stratifying patients for appropriate therapy.

Integrative transcriptomic and proteomic profile revealed inhibition of oxidative phosphorylation and peroxisomes during renal interstitial fibrosis

Effective therapies of chronic kidney disease (CKD) are lacking due to the unclear molecular pathogenesis. Previous single omics-studies have described potential molecular regulation mechanism of CKD only at the level of transcription or translation. Therefore, this study generated an integrated transcriptomic and proteomic profile to provide deep insights into the continuous transcription-translation process during CKD. The comprehensive datasets identified 14,948 transcripts and 6423 proteins, 233 up-regulated and 364 down-regulated common differentially expressed genes of transcriptome and proteome were selected to further combined bioinformatics analysis. The obtained results revealed reactive oxygen species (ROS) metabolism and antioxidant system due to imbalance of mitochondria and peroxisomes were significantly repressed in CKD. Overall, this study presents a valuable multi-omics analysis that sheds light on the molecular mechanisms underlying CKD. SIGNIFICANCE: Chronic kidney disease (CKD) is a progressive and irreversible condition that results in abnormal kidney function and structure, and is ranked 18th among the leading causes of death globally, leading to a significant societal burden. Hence, there is an urgent need for research to detect new, sensitive, and specific biomarkers. Omics-based studies offer great potential to identify underlying disease mechanisms, aid in clinical diagnosis, and develop novel treatment strategies for CKD. Previous studies have mainly focused on the regulation of gene expression or protein synthesis in CKD, thereby compelling us to conduct a meticulous analysis of transcriptomic and proteomic data from the UUO mouse model. Here, we have performed a unified analysis of CKD model by integrating transcriptomes and protein suites for the first time. Our study contributes to a deeper understanding of the pathogenesis of CKD and provides a basis for subsequent disease management and drug development.

PCTA, a pan-cancer cell line transcriptome atlas

A substantial volume of RNA sequencing data have been generated from cancer cell lines. However, it requires specific bioinformatics skills to compare gene expression levels across cell lines. This has hindered non-bioinformaticians from fully utilizing these valuable datasets in their research. To bridge this gap, we established a curated Pan-cancer Cell Line Transcriptome Atlas (PCTA) dataset. This resource aims to provide a user-friendly platform, allowing researchers without extensive bioinformatics expertise to access and leverage the wealth of information within the dataset for their studies. The PCTA dataset encompasses the expression matrix of 24,965 genes, featuring data from 84,385 samples derived from 5677 studies. This comprehensive compilation spans 535 cell lines, representing a spectrum of 114 cancer types originating from 30 diverse tissue types. On UMAP plots, cell lines originating from the same type of tissue tend to cluster together, illustrating the dataset's ability to capture biological relationships. Additionally, an interactive and user-friendly web application (https://pcatools.shinyapps.io/PCTA_app/) was developed for researchers to explore the PCTA dataset. This platform allows users to examine the expression of their genes of interest across a diverse array of samples.

Single-cell transcriptomics dissects the transcriptome alterations of hematopoietic stem cells in myelodysplastic neoplasms

BACKGROUND

Myelodysplastic neoplasms (MDS) are myeloid neoplasms characterized by disordered differentiation of hematopoietic stem cells and a predisposition to acute myeloid leukemia (AML). The underline pathogenesis remains unclear.

METHODS

In this study, the trajectory of differentiation and mechanisms of leukemic transformation were explored through bioinformatics analysis of single-cell RNA-Seq data from hematopoietic stem and progenitor cells (HSPCs) in MDS patients.

RESULTS

Among the HSPC clusters, the proportion of common myeloid progenitor (CMP) was the main cell cluster in the patients with excess blasts (EB)/ secondary AML. Cell cycle analysis indicated the CMP of MDS patients were in an active proliferative state. The genes involved in the cell proliferation, such as MAML3 and PLCB1, were up-regulated in MDS CMP. Further validation analysis indicated that the expression levels of MAML3 and PLCB1 in patients with MDS-EB were significantly higher than those without EB. Patients with high expression of PLCB1 had a higher risk of transformation to AML. PLCB1 inhibitor can suppress proliferation, induce cell cycle arrest, and activate apoptosis of leukemic cells in vitro.

CONCLUSION

This study revealed the transcriptomic change of HSPCs in MDS patients along the pseudotime and indicated that PLCB1 plays a key role in the transformation of MDS into leukemia.

DNA-based molecular classifiers for the profiling of gene expression signatures

Although gene expression signatures offer tremendous potential in diseases diagnostic and prognostic, but massive gene expression signatures caused challenges for experimental detection and computational analysis in clinical setting. Here, we introduce a universal DNA-based molecular classifier for profiling gene expression signatures and generating immediate diagnostic outcomes. The molecular classifier begins with feature transformation, a modular and programmable strategy was used to capture relative relationships of low-concentration RNAs and convert them to general coding inputs. Then, competitive inhibition of the DNA catalytic reaction enables strict weight assignment for different inputs according to their importance, followed by summation, annihilation and reporting to accurately implement the mathematical model of the classifier. We validated the entire workflow by utilizing miRNA expression levels for the diagnosis of hepatocellular carcinoma (HCC) in clinical samples with an accuracy 85.7%. The results demonstrate the molecular classifier provides a universal solution to explore the correlation between gene expression patterns and disease diagnostics, monitoring, and prognosis, and supports personalized healthcare in primary care.

Prominent transcriptomic changes in Mycobacterium intracellulare under acidic and oxidative stress

BACKGROUND

Mycobacterium avium complex (MAC), including Mycobacterium intracellulare is a member of slow-growing mycobacteria and contributes to a substantial proportion of nontuberculous mycobacterial lung disease in humans affecting immunocompromised and elderly populations. Adaptation of pathogens in hostile environments is crucial in establishing infection and persistence within the host. However, the sophisticated cellular and molecular mechanisms of stress response in M. intracellulare still need to be fully explored. We aimed to elucidate the transcriptional response of M. intracellulare under acidic and oxidative stress conditions.

RESULTS

At the transcriptome level, 80 genes were shown [FC] ≥ 2.0 and p < 0.05 under oxidative stress with 10 mM hydrogen peroxide. Specifically, 77 genes were upregulated, while 3 genes were downregulated. In functional analysis, oxidative stress conditions activate DNA replication, nucleotide excision repair, mismatch repair, homologous recombination, and tuberculosis pathways. Additionally, our results demonstrate that DNA replication and repair system genes, such as dnaB, dinG, urvB, uvrD2, and recA, are indispensable for resistance to oxidative stress. On the contrary, 878 genes were shown [FC] ≥ 2.0 and p < 0.05 under acidic stress with pH 4.5. Among these genes, 339 were upregulated, while 539 were downregulated. Functional analysis highlighted nitrogen and sulfur metabolism pathways as the primary responses to acidic stress. Our findings provide evidence of the critical role played by nitrogen and sulfur metabolism genes in the response to acidic stress, including narGHIJ, nirBD, narU, narK3, cysND, cysC, cysH, ferredoxin 1 and 2, and formate dehydrogenase.

CONCLUSION

Our results suggest the activation of several pathways potentially critical for the survival of M. intracellulare under a hostile microenvironment within the host. This study indicates the importance of stress responses in M. intracellulare infection and identifies promising therapeutic targets.

Global changes in gene expression during compatible and incompatible interactions of faba bean (Vicia faba L.) during Orobanche foetida parasitism

Orobanche foetida Poiret is the main constraint facing faba bean crop in Tunisia. Indeed, in heavily infested fields with this parasitic plant, yield losses may reach 90%, and the recent estimation of the infested area is around 80,000 ha. Identifying genes involved in the Vicia faba/O. foetida interaction is crucial for the development of effective faba bean breeding programs. However, there is currently no available information on the transcriptome of faba bean responding to O. foetida parasitism. In this study, we employed RNA sequencing to explore the global gene expression changes associated with compatible and incompatible V. faba/O. foetida interactions. In this perspective, two faba bean varieties (susceptible and resistant) were examined at the root level across three stages of O. foetida development (Before Germination (BG), After Germination (AG) and Tubercule Stage (TS)). Our analyses presented an exploration of the transcriptomic profile, including comprehensive assessments of differential gene expression and Gene Ontology (GO) enrichment analyses. Specifically, we investigated key pathways revealing the complexity of molecular responses to O. foetida attack. In this study, we detected differential gene expression of pathways associated with secondary metabolites: flavonoids, auxin, thiamine, and jasmonic acid. To enhance our understanding of the global changes in V. faba response to O. foetida, we specifically examined WRKY genes known to play a role in plant host-parasitic plant interactions. Furthermore, considering the pivotal role of parasitic plant seed germination in this interaction, we investigated genes involved in the orobanchol biosynthesis pathway. Interestingly, we detected the gene expression of VuCYP722C homolog, coding for a key enzyme involved in orobanchol biosynthesis, exclusively in the susceptible host. Clearly, this study enriches our understanding of the V. faba/O. foetida interaction, shedding light on the main differences between susceptible and resistant faba bean varieties during O. foetida infestation at the gene expression level.

New treatment alternatives for primary and metastatic colorectal cancer by an integrated transcriptome and network analyses

Metastatic colorectal cancer (CRC) is still in need of effective treatments. This study applies a holistic approach to propose new targets for treatment of primary and liver metastatic CRC and investigates their therapeutic potential in-vitro. An integrative analysis of primary and metastatic CRC samples was implemented for alternative target and treatment proposals. Integrated microarray samples were grouped based on a co-expression network analysis. Significant gene modules correlated with primary CRC and metastatic phenotypes were identified. Network clustering and pathway enrichments were applied to gene modules to prioritize potential targets, which were shortlisted by independent validation. Finally, drug-target interaction search led to three agents for primary and liver metastatic CRC phenotypes. Hesperadin and BAY-1217389 suppress colony formation over a 14-day period, with Hesperadin showing additional efficacy in reducing cell viability within 48 h. As both candidates target the G2/M phase proteins NEK2 or TTK, we confirmed their anti-proliferative properties by Ki-67 staining. Hesperadinin particular arrested the cell cycle at the G2/M phase. IL-29A treatment reduced migration and invasion capacities of TGF-β induced metastatic cell lines. In addition, this anti-metastatic treatment attenuated TGF-β dependent mesenchymal transition. Network analysis suggests IL-29A induces the JAK/STAT pathway in a preventive manner.

pekinensis)

BACKGROUND

Long non-coding RNAs (lncRNAs) play a crucial role in regulating gene expression vital for the growth and development of plants. Despite this, the role of lncRNAs in Chinese cabbage (Brassica rapa L. ssp. pekinensis) pollen development and male fertility remains poorly understood.

RESULTS

In this study, we characterized a recessive genic male sterile mutant (366-2 S), where the delayed degradation of tapetum and the failure of tetrad separation primarily led to the inability to form single microspores, resulting in male sterility. To analyze the role of lncRNAs in pollen development, we conducted a comparative lncRNA sequencing using anthers from the male sterile mutant line (366-2 S) and the wild-type male fertile line (366-2 F). We identified 385 differentially expressed lncRNAs between the 366-2 F and 366-2 S lines, with 172 of them potentially associated with target genes. To further understand the alterations in mRNA expression and explore potential lncRNA-target genes (mRNAs), we performed comparative mRNA transcriptome analysis in the anthers of 366-2 S and 366-2 F at two stages. We identified 1,176 differentially expressed mRNAs. Remarkably, GO analysis revealed significant enrichment in five GO terms, most notably involving mRNAs annotated as pectinesterase and polygalacturonase, which play roles in cell wall degradation. The considerable downregulation of these genes might contribute to the delayed degradation of tapetum in 366-2 S. Furthermore, we identified 15 lncRNA-mRNA modules through Venn diagram analysis. Among them, MSTRG.9997-BraA04g004630.3 C (β-1,3-glucanase) is associated with callose degradation and tetrad separation. Additionally, MSTRG.5212-BraA02g040020.3 C (pectinesterase) and MSTRG.13,532-BraA05g030320.3 C (pectinesterase) are associated with cell wall degradation of the tapetum, indicating that these three candidate lncRNA-mRNA modules potentially regulate pollen development.

CONCLUSION

This study lays the foundation for understanding the roles of lncRNAs in pollen development and for elucidating their molecular mechanisms in regulating male sterility in Chinese cabbage.

Principles, functions, and biological implications of m6A in plants

Over the past decade, N 6-methyladenosine (m6A) has emerged as a prevalent and dynamically regulated modification across the transcriptome; it has been reversibly installed, removed, and interpreted by specific binding proteins, and has played crucial roles in molecular and biological processes. Within this scope, we consolidate recent advancements of m6A research in plants regarding gene expression regulation, diverse physiologic and pathogenic processes, as well as crop trial implications, to guide discussions on challenges associated with and leveraging epitranscriptome editing for crop improvement.

Integrated analysis of metabolome, transcriptome, and bioclimatic factors of Acer truncatum seeds reveals key candidate genes related to unsaturated fatty acid biosynthesis, and potentially optimal production area

BACKGROUND

Lipids found in plant seeds are essential for controlling seed dormancy, dispersal, and defenses against biotic and abiotic stress. Additionally, these lipids provide nutrition and energy and are therefore important to the human diet as edible oils. Acer truncatum, which belongs to the Aceaceae family, is widely cultivated around the world for its ornamental value. Further because its seed oil is rich in unsaturated fatty acids (UFAs)- i.e. α-linolenic acid (ALA) and nervonic acid (NA)- and because it has been validated as a new food resource in China, the importance of A. truncatum has greatly risen. However, it remains unknown how UFAs are biosynthesized during the growth season, to what extent environmental factors impact their content, and what areas are potentially optimal for their production.

RESULTS

In this study, transcriptome and metabolome of A. truncatum seeds at three representative developmental stages was used to find the accumulation patterns of all major FAs. Cumulatively, 966 metabolites and 87,343 unigenes were detected; the differential expressed unigenes and metabolites were compared between stages as follows: stage 1 vs. 2, stage 1 vs. 3, and stage 2 vs. 3 seeds, respectively. Moreover, 13 fatty acid desaturases (FADs) and 20 β-ketoacyl-CoA synthases (KCSs) were identified, among which the expression level of FAD3 (Cluster-7222.41455) and KCS20 (Cluster-7222.40643) were consistent with the metabolic results of ALA and NA, respectively. Upon analysis of the geographical origin-affected diversity from 17 various locations, we found significant variation in phenotypes and UFA content. Notably, in this study we found that 7 bioclimatic variables showed considerable influence on FAs contents in A. truncatum seeds oil, suggesting their significance as critical environmental parameters. Ultimately, we developed a model for potentially ecological suitable regions in China.

CONCLUSION

This study provides a comprehensive understanding of the relationship between metabolome and transcriptome in A. truncatum at various developmental stages of seeds and a new strategy to enhance seed FA content, especially ALA and NA. This is particularly significant in meeting the increasing demands for high-quality edible oil for human consumption. The study offers a scientific basis for A. truncatum's novel utilization as a woody vegetable oil rather than an ornamental plant, potentially expanding its cultivation worldwide.

New directions for Ψ and m1A decoding in mRNA: deciphering the stoichiometry and function

Over the past decade, advancements in epitranscriptomics have significantly enhanced our understanding of mRNA metabolism and its role in human development and diseases. This period has witnessed breakthroughs in sequencing technologies and the identification of key proteins involved in RNA modification processes. Alongside the well-studied m6A, Ψ and m1A have emerged as key epitranscriptomic markers. Initially identified through transcriptome-wide profiling, these modifications are now recognized for their broad impact on RNA metabolism and gene expression. In this Perspective, we focus on the detections and functions of Ψ and m1A modifications in mRNA and discuss previous discrepancies and future challenges. We summarize recent advances and highlight the latest sequencing technologies for stoichiometric detection and their mechanistic investigations for functional unveiling in mRNA as the new research directions.

Opto-combinatorial indexing enables high-content transcriptomics by linking cell images and transcriptomes

We introduce a simple integrated analysis method that links cellular phenotypic behaviour with single-cell RNA sequencing (scRNA-seq) by utilizing a combination of optical indices from cells and hydrogel beads. With our method, the combinations, referred to as joint colour codes, enable the link via matching the optical combinations measured by conventional epi-fluorescence microscopy with the concatenated DNA molecular barcodes created by cell-hydrogel bead pairs and sequenced by next-generation sequencing. We validated our approach by demonstrating an accurate link between the cell image and scRNA-seq with mixed species experiments, longitudinal cell tagging by electroporation and lipofection, and gene expression analysis. Furthermore, we extended our approach to multiplexed chemical transcriptomics, which enabled us to identify distinct phenotypic behaviours in HeLa cells treated with various concentrations of paclitaxel, and determine the corresponding gene regulation associated with the formation of a multipolar spindle.

Human genetic associations of the airway microbiome in chronic obstructive pulmonary disease

Little is known about the relationships between human genetics and the airway microbiome. Deeply sequenced airway metagenomics, by simultaneously characterizing the microbiome and host genetics, provide a unique opportunity to assess the microbiome-host genetic associations. Here we performed a co-profiling of microbiome and host genetics with the identification of over 5 million single nucleotide polymorphisms (SNPs) through deep metagenomic sequencing in sputum of 99 chronic obstructive pulmonary disease (COPD) and 36 healthy individuals. Host genetic variation was the most significant factor associated with the microbiome except for geography and disease status, with its top 5 principal components accounting for 12.11% of the microbiome variability. Within COPD individuals, 113 SNPs mapped to candidate genes reported as genetically associated with COPD exhibited associations with 29 microbial species and 48 functional modules (P < 1 × 10-5), where Streptococcus salivarius exhibits the strongest association to SNP rs6917641 in TBC1D32 (P = 9.54 × 10-8). Integration of concurrent host transcriptomic data identified correlations between the expression of host genes and their genetically-linked microbiome features, including NUDT1, MAD1L1 and Veillonella parvula, TTLL9 and Stenotrophomonas maltophilia, and LTA4H and Haemophilus influenzae. Mendelian randomization analyses revealed a potential causal link between PARK7 expression and microbial type III secretion system, and a genetically-mediated association between COPD and increased relative abundance of airway Streptococcus intermedius. These results suggest a previously underappreciated role of host genetics in shaping the airway microbiome and provide fresh hypotheses for genetic-based host-microbiome interactions in COPD.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease

BACKGROUND

X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique.

RESULTS

In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis.

CONCLUSIONS

This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Epigenetic regulation of H3K27me3 in laying hens with fatty liver hemorrhagic syndrome induced by high-energy and low-protein diets

BACKGROUND

Fatty liver hemorrhagic syndrome (FLHS) in the modern poultry industry is primarily caused by nutrition. Despite encouraging progress on FLHS, the mechanism through which nutrition influences susceptibility to FLHS is still lacking in terms of epigenetics.

RESULTS

In this study, we analyzed the genome-wide patterns of trimethylated lysine residue 27 of histone H3 (H3K27me3) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes in healthy and FLHS hens. The study results indicated that H3K27me3 levels were increased in the FLHS hens on a genome-wide scale. Additionally, H3K27me3 was found to occupy the entire gene and the distant intergenic region, which may function as silencer-like regulatory elements. The analysis of transcription factor (TF) motifs in hypermethylated peaks has demonstrated that 23 TFs are involved in the regulation of liver metabolism and development. Transcriptomic analysis indicated that differentially expressed genes (DEGs) were enriched in fatty acid metabolism, amino acid, and carbohydrate metabolism. The hub gene identified from PPI network is fatty acid synthase (FASN). Combined ChIP-seq and transcriptome analysis revealed that the increased H3K27me3 and down-regulated genes have significant enrichment in the ECM-receptor interaction, tight junction, cell adhesion molecules, adherens junction, and TGF-beta signaling pathways.

CONCLUSIONS

Overall, the trimethylation modification of H3K27 has been shown to have significant regulatory function in FLHS, mediating the expression of crucial genes associated with the ECM-receptor interaction pathway. This highlights the epigenetic mechanisms of H3K27me3 and provides insights into exploring core regulatory targets and nutritional regulation strategies in FLHS.