Prognostic significance of circulating tumor cells and tumor related transcripts in small cell lung cancer: A step further to clinical implementation

Small-cell lung cancer (SCLC) is a fatal disease with limited treatment options. Circulating tumor cells (CTCs) in liquid biopsy samples may serve as predictive and prognostic biomarkers; but the analysis of CTCs is still challenging. By using microfluidic or density gradient CTC enrichment in combination with immunofluorescent (IF) staining or qPCR of CTC-related transcripts, we achieved a 60.8% to 88.0% positivity in SCLC blood samples. Epithelial and neuroendocrine transcripts including the druggable target DLL3 were associated with shorter overall survival (OS), indicating the clinical value of these markers in terms of differential diagnosis and treatment decisions. High CTC counts and the presence of CTC duplets detected by IF staining were prognostic for OS, and thus may serve as indicators of disease progression or therapy failure. In patient samples with high CTC load detected by IF staining, a concordance of the transcripts positivity in circulating free plasma RNA and CTCs was observed. Our data emphasize the role of CTCs and CTC-related transcripts and underline the clinical value of liquid biopsy analysis in SCLC.

Characterization of PANoptosis-related genes and the immune landscape in moyamoya disease

Moyamoya disease (MMD) is a cerebrovascular narrowing and occlusive condition characterized by progressive stenosis of the terminal portion of the internal carotid artery and the formation of an abnormal network of dilated, fragile perforators at the base of the brain. However, the role of PANoptosis, an apoptotic mechanism associated with vascular disease, has not been elucidated in MMD. In our study, a total of 40 patients' genetic data were included, and a total of 815 MMD-related differential genes were screened, including 215 upregulated genes and 600 downregulated genes. Among them, DNAJA3, ESR1, H19, KRT18 and STK3 were five key genes. These five key genes were associated with a variety of immune cells and immune factors. Moreover, GSEA (gene set enrichment analysis) and GSVA (gene set variation analysis) showed that the different expression levels of the five key genes affected multiple signaling pathways associated with MMD. In addition, they were associated with the expression of MMD-related genes. Then, based on the five key genes, a transcription factor regulatory network was constructed. In addition, targeted therapeutic drugs against MMD-related genes were obtained by the Cmap drug prediction method: MST-312, bisacodyl, indirubin, and tropanyl-3,5-dimethylbenzoate. These results suggest that the PANoptosis-related genes may contribute to the pathogenesis of MMD through multiple mechanisms.

JAK-STAT signaling maintains homeostasis in T cells and macrophages

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Transcriptomic and Metatranscriptomic Analyses Provide New Insights into the Response of the Pea Aphid Acyrthosiphon pisum (Hemiptera: Aphididae) to Acetamiprid

Acetamiprid is a broad-spectrum neonicotinoid insecticide used in agriculture to control aphids. While recent studies have documented resistance to acetamiprid in several aphid species, the underlying mechanisms are still not fully understood. In this study, we analyzed the transcriptome and metatranscriptome of a laboratory strain of the pea aphid, Acyrthosiphon pisum (Harris, 1776), with reduced susceptibility to acetamiprid after nine generations of exposure to identify candidate genes and the microbiome involved in the adaptation process. Sequencing of the transcriptome of both selected (RS) and non-selected (SS) strains allowed the identification of 14,858 genes and 4938 new transcripts. Most of the differentially expressed genes were associated with catalytic activities and metabolic pathways involving carbon and fatty acids. Specifically, alcohol-forming fatty acyl-CoA reductase (FAR) and acyl-CoA synthetase (ACSF2), both involved in the synthesis of epidermal wax layer components, were significantly upregulated in RS, suggesting that adaptation to acetamiprid involves the synthesis of a thicker protective layer. Metatranscriptomic analyses revealed subtle shifts in the microbiome of RS. These results contribute to a deeper understanding of acetamiprid adaptation by the pea aphid and provide new insights for aphid control strategies.

Protein restriction during pregnancy alters Cdkn1c silencing, dopamine circuitry and offspring behaviour without changing expression of key neuronal marker genes

We tracked the consequences of in utero protein restriction in mice throughout their development and life course using a luciferase-based allelic reporter of imprinted Cdkn1c. Exposure to gestational low-protein diet (LPD) results in the inappropriate expression of paternally inherited Cdkn1c in the brains of embryonic and juvenile mice. These animals were characterised by a developmental delay in motor skills, and by behavioural alterations indicative of reduced anxiety. Exposure to LPD in utero resulted in significantly more tyrosine hydroxylase positive (dopaminergic) neurons in the midbrain of adult offspring as compared to age-matched, control-diet equivalents. Positron emission tomography (PET) imaging revealed an increase in striatal dopamine synthesis capacity in LPD-exposed offspring, where elevated levels of dopamine correlated with an enhanced sensitivity to cocaine. These data highlight a profound sensitivity of the developing epigenome to gestational protein restriction. Our data also suggest that loss of Cdkn1c imprinting and p57KIP2 upregulation alters the cellular composition of the developing midbrain, compromises dopamine circuitry, and thereby provokes behavioural abnormalities in early postnatal life. Molecular analyses showed that despite this phenotype, exposure to LPD solely during pregnancy did not significantly change the expression of key neuronal- or dopamine-associated marker genes in adult offspring.

Age-dependent loss of Crls1 causes myopathy and skeletal muscle regeneration failure

Skeletal muscle aging results in the gradual suppression of myogenesis, leading to muscle mass loss. However, the specific role of cardiolipin in myogenesis has not been determined. This study investigated the crucial role of mitochondrial cardiolipin and cardiolipin synthase 1 (Crls1) in age-related muscle deterioration and myogenesis. Our findings demonstrated that cardiolipin and Crls1 are downregulated in aged skeletal muscle. Moreover, the knockdown of Crls1 in myoblasts reduced mitochondrial mass, activity, and OXPHOS complex IV expression and disrupted the structure of the mitochondrial cristae. AAV9-shCrls1-mediated downregulation of Crls1 impaired muscle regeneration in a mouse model of cardiotoxin (CTX)-induced muscle damage, whereas AAV9-mCrls1-mediated Crls1 overexpression improved regeneration. Overall, our results highlight that the age-dependent decrease in CRLS1 expression contributes to muscle loss by diminishing mitochondrial quality in skeletal muscle myoblasts. Hence, modulating CRLS1 expression is a promising therapeutic strategy for mitigating muscle deterioration associated with aging, suggesting potential avenues for developing interventions to improve overall muscle health and quality of life in elderly individuals.

An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

TM4SF19-mediated control of lysosomal activity in macrophages contributes to obesity-induced inflammation and metabolic dysfunction

Adipose tissue (AT) adapts to overnutrition in a complex process, wherein specialized immune cells remove and replace dysfunctional and stressed adipocytes with new fat cells. Among immune cells recruited to AT, lipid-associated macrophages (LAMs) have emerged as key players in obesity and in diseases involving lipid stress and inflammation. Here, we show that LAMs selectively express transmembrane 4 L six family member 19 (TM4SF19), a lysosomal protein that represses acidification through its interaction with Vacuolar-ATPase. Inactivation of TM4SF19 elevates lysosomal acidification and accelerates the clearance of dying/dead adipocytes in vitro and in vivo. TM4SF19 deletion reduces the LAM accumulation and increases the proportion of restorative macrophages in AT of male mice fed a high-fat diet. Importantly, male mice lacking TM4SF19 adapt to high-fat feeding through adipocyte hyperplasia, rather than hypertrophy. This adaptation significantly improves local and systemic insulin sensitivity, and energy expenditure, offering a potential avenue to combat obesity-related metabolic dysfunction.

Comparative Genomic Analysis of PEBP Genes in Cucurbits Explores the Interactors of Cucumber CsPEBPs Related to Flowering Time

The family of phosphatidylethanolamine-binding proteins (PEBPs) participates in various plant biological processes, mainly flowering regulation and seed germination. In cucurbit crops, several PEBP genes have been recognized to be responsible for flowering time. However, the investigation of PEBP family members across the genomes of cucurbit species has not been reported, and their conservation and divergence in structure and function remain largely unclear. Herein, PEBP genes were identified from seven cucurbit crops and were used to perform a comparative genomics analysis. The cucurbit PEBP proteins could be classified into MFT, FT, TFL, and PEBP clades, and further, the TFL clade was divided into BFT-like, CEN-like, and TFL1-like subclades. The MFT-like, FT-like, and TFL-like proteins were clearly distinguished by a critical amino acid residue at the 85th position of the Arabidopsis FT protein. In gene expression analysis, CsaPEBP1 was highly expressed in flowers, and its expression levels in females and males were 70.5 and 89.2 times higher, respectively, than those in leaves. CsaPEBP5, CsaPEBP6, and CsaPEBP7 were specifically expressed in male flowers, with expression levels 58.1, 17.3, and 15.7 times higher, respectively, than those of leaves. At least five CsaPEBP genes exhibited the highest expression during the later stages of corolla opening. Through clustering of time-series-based RNA-seq data, several potential transcription factors (TFs) interacting with four CsaPEBPs were identified during cucumber corolla opening. Because of the tandem repeats of binding sites in promoters, NF-YB (Csa4G037610) and GATA (Csa7G64580) TFs appeared to be better able to regulate the CsaPEBP2 and CsaPEBP5 genes, respectively. This study would provide helpful information for further investigating the roles of PEBP genes and their interacting TFs in growth and development processes, such as flowering time regulation in cucurbit crops.

Implementable assay for monitoring minimum residual disease after radical treatment for colorectal cancer

Considering the cost and invasiveness of monitoring postoperative minimal residual disease (MRD) of colorectal cancer (CRC) after adjuvant chemoradiotherapy (ACT), we developed a favorable approach based on methylated circulating tumor DNA to detect MRD after radical resection. Analyzing the public database, we identified the methylated promoter regions of the genes FGD5, GPC6, and MSC. Using digital polymerase chain reaction (dPCR), we termed the "amplicon of methylated sites using a specific enzyme" assay as "AMUSE." We examined 180 and 114 pre- and postoperative serial plasma samples from 28 recurrent and 19 recurrence-free pathological stage III CRC patients, respectively. The results showed 22 AMUSE-positive of 28 recurrent patients (sensitivity, 78.6%) and 17 AMUSE-negative of 19 recurrence-free patients (specificity, 89.5%). AMUSE predicted recurrence 208 days before conventional diagnosis using radiological imaging. Regarding ACT evaluation by the reactive response, 19 AMUSE-positive patients during their second or third blood samples showed a significantly poorer prognosis than the other patients (p = 9E-04). The AMUSE assay stratified four groups by the altered patterns of tumor burden postoperatively. Interestingly, only 34.8% of cases tested AMUSE-negative during ACT treatment, indicating eligibility for ACT. The AMUSE assay addresses the clinical need for accurate MRD monitoring with universal applicability, minimal invasiveness, and cost-effectiveness, thereby enabling the timely detection of recurrences. This assay can effectively evaluate the efficacy of ACT in patients with stage III CRC following curative resection. Our study strongly recommends reevaluating the clinical application of ACT using the AMUSE assay.

Evaluation of Collection and Processing Conditions for Gene Expression Analysis Using Human Myeloid Cells

Background: The population of blast cells among peripheral blood mononuclear cells (PBMCs) obtained from patients is a desirable specimen for analyzing gene expression in diseases including acute myeloid leukemia. Although the enrichment of blast cells often needs to be performed at a central laboratory, acceptable conditions for sample transport from clinical sites remain to be established. Methods: We evaluated storage temperature, duration, and tube type before initiating sample processing for the analysis of cluster of differentiation (CD)33+ myeloid cells among PBMCs as an alternative to CD34+/CD33+ blast cells. Results: CD33+ myeloid cells were successfully purified by MACS. The cell viability and the RNA integrity were sustained during storage up to 48 hours before sample processing. Storage at 4°C had minimal effects on gene expression, whereas storage at room temperature induced the senescence pathway, characterized by the expression of stress-inducible genes. A CPT tube was also better than an ethylenediaminetetraacetic acid tube for minimizing gene expression change. Conclusions: Our study provided important clues for establishing a sample handling approach for gene expression analysis with purified cell fractions from human PBMCs. To keep the variation of gene expression to a minimum, samples should be delivered at 4°C within 48 hours before processing.

Pregnenolone sulfate induces transcriptional and immunoregulatory effects on T cells

Pregnenolone sulfate is a steroid metabolite of the steroidogenesis precursor, pregnenolone, with similar functional properties, including immunosuppression. We recently reported an elevation in serum levels of pregnenolone sulfate in children with malaria, contributing to an immunosuppressed state. Yet, the molecular mechanisms in which this steroid exerts its immunoregulatory functions are lacking. In this study, we examined the effects of pregnenolone sulfate on T cell viability, proliferation and transcriptome. We observed a pregnenolone sulfate dose-dependent induction of T cell death and reduction in proliferation. RNA sequencing analysis of pregnenolone sulfate-treated T cells for 2 and 24 h revealed the downregulation of pro-inflammatory genes and the upregulation of the steroid nuclear receptor superfamily, NR4A, as early-response genes. We also report a strong activation of the integrated stress response mediated by the upregulation of EIF2AK3. These results contribute to the knowledge on transcriptional regulation driving the immunoregulatory effects of pregnenolone sulfate on T cells.

Local effect of allopregnanolone in rat ovarian steroidogenesis, follicular and corpora lutea development

Allopregnanolone (ALLO) is a known neurosteroid and a progesterone metabolite synthesized in the ovary, CNS, PNS, adrenals and placenta. Its role in the neuroendocrine control of ovarian physiology has been studied, but its in situ ovarian effects are still largely unknown. The aims of this work were to characterize the effects of intrabursal ALLO administration on different ovarian parameters, and the probable mechanism of action. ALLO administration increased serum progesterone concentration and ovarian 3β-HSD2 while decreasing 20α-HSD mRNA expression. ALLO increased the number of atretic follicles and the number of positive TUNEL granulosa and theca cells, while decreasing positive PCNA immunostaining. On the other hand, there was an increase in corpora lutea diameter and PCNA immunostaining, whereas the count of TUNEL-positive luteal cells decreased. Ovarian angiogenesis and the immunohistochemical expression of GABAA receptor increased after ALLO treatment. To evaluate if the ovarian GABAA receptor was involved in these effects, we conducted a functional experiment with a specific antagonist, bicuculline. The administration of bicuculline restored the number of atretic follicles and the diameter of corpora lutea to normal values. These results show the actions of ALLO on the ovarian physiology of the female rat during the follicular phase, some of them through the GABAA receptor. Intrabursal ALLO administration alters several processes of the ovarian morpho-physiology of the female rat, related to fertility and oocyte quality.

Determination of internal controls for quantitative gene expression of Spodoptera litura under microbial pesticide stress

Quantitative real-time polymerase chain reaction (qRT-PCR) has become a commonly used method for the quantification of gene expression. However, accurate qRT-PCR analysis requires a valid internal reference for data normalization. To determine the valid reference characterized with low expression variability among Spodoptera litura samples after microbial pesticide treatments, nine housekeeping genes, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), arginine kinase, ubiquitin C, actin-5C (ACT5C), actin, ribosomal protein S13 (RPS13), tubulin, acidic ribosomal protein P0 (RPLP0) and ubiquinol-cytochrome c reductase, were evaluated for their suitability using geNorm, Normfinder, BestKeeper, RefFinder and the comparative delta CT methods in this study. S. litura larvae after direct treatment (larvae were immersed in biopesticides), indirect treatment (larvae were fed with biopesticide immersed artificial diets) and comprehensive treatment (larvae were treated with the first two treatments in sequence), respectively with Metarhizium anisopliae, Empedobacter brevis and Bacillus thuringiensis, were investigated. The results indicated that the best sets of internal references were as follows: RPLP0 and ACT5C for direct treatment conditions; RPLP0 and RPS13 for indirect treatment conditions; RPS13 and GAPDH for comprehensive treatment conditions; RPS13 and RPLP0 for all the samples. These results provide valuable bases for further genetic researches in S. litura.

Kdm1a safeguards the topological boundaries of PRC2-repressed genes and prevents aging-related euchromatinization in neurons

Kdm1a is a histone demethylase linked to intellectual disability with essential roles during gastrulation and the terminal differentiation of specialized cell types, including neurons, that remains highly expressed in the adult brain. To explore Kdm1a's function in adult neurons, we develop inducible and forebrain-restricted Kdm1a knockouts. By applying multi-omic transcriptome, epigenome and chromatin conformation data, combined with super-resolution microscopy, we find that Kdm1a elimination causes the neuronal activation of nonneuronal genes that are silenced by the polycomb repressor complex and interspersed with active genes. Functional assays demonstrate that the N-terminus of Kdm1a contains an intrinsically disordered region that is essential to segregate Kdm1a-repressed genes from the neighboring active chromatin environment. Finally, we show that the segregation of Kdm1a-target genes is weakened in neurons during natural aging, underscoring the role of Kdm1a safeguarding neuronal genome organization and gene silencing throughout life.

Deciphering Breast Cancer Metastasis Cascade: A Systems Biology Approach Integrating Transcriptome and Interactome Insights for Target Discovery

Breast cancer is the lead cause of cancer-related deaths among women globally. Breast cancer metastasis is a complex and still inadequately understood process and a key dimension of mortality attendant to breast cancer. This study reports dysregulated genes across metastatic stages and tissues, shedding light on their molecular interplay in disease pathogenesis and new possibilities for drug discovery. Comprehensive analyses of gene expression data from primary breast tumor, circulating tumor cells, and distant metastatic sites in the brain, lung, liver, and bone were conducted. Genes dysregulated across multiple stages and tissues were identified as metastatic cascade genes, and are further classified based on functional associations with metastasis-related mechanisms. Their interactions with HUB genes in interactome networks were scrutinized, followed by pathway enrichment analysis. Validation for their potential as targets included assessments for survival, druggability, prognostic marker status, secretome annotation, protein expression, and cell type marker association. Results displayed critical genes in the metastatic cascade and those specific to metastatic sites, revealing the involvement of the collagen degradation and assembly of collagen fibrils and other multimeric structure pathways in driving metastasis. Notably, pivotal cascade genes FABP4, CXCL12, APOD, and IGF1 emerged with high metastatic potential, linked to significant druggability and survival scores, establishing them as potential molecular targets. The significance of this research lies in its potential to uncover novel biomarkers for early detection, therapeutic targets, and a deeper understanding of the molecular mechanisms underpinning the metastatic cascade in breast cancer, and with an eye to precision/personalized medicine.

A single cell atlas of frozen shoulder capsule identifies features associated with inflammatory fibrosis resolution

Frozen shoulder is a spontaneously self-resolving chronic inflammatory fibrotic human disease, which distinguishes the condition from most fibrotic diseases that are progressive and irreversible. Using single-cell analysis, we identify pro-inflammatory MERTKlowCD48+ macrophages and MERTK + LYVE1 + MRC1+ macrophages enriched for negative regulators of inflammation which co-exist in frozen shoulder capsule tissues. Micro-cultures of patient-derived cells identify integrin-mediated cell-matrix interactions between MERTK+ macrophages and pro-resolving DKK3+ and POSTN+ fibroblasts, suggesting that matrix remodelling plays a role in frozen shoulder resolution. Cross-tissue analysis reveals a shared gene expression cassette between shoulder capsule MERTK+ macrophages and a respective population enriched in synovial tissues of rheumatoid arthritis patients in disease remission, supporting the concept that MERTK+ macrophages mediate resolution of inflammation and fibrosis. Single-cell transcriptomic profiling and spatial analysis of human foetal shoulder tissues identify MERTK + LYVE1 + MRC1+ macrophages and DKK3+ and POSTN+ fibroblast populations analogous to those in frozen shoulder, suggesting that the template to resolve fibrosis is established during shoulder development. Crosstalk between MerTK+ macrophages and pro-resolving DKK3+ and POSTN+ fibroblasts could facilitate resolution of frozen shoulder, providing a basis for potential therapeutic resolution of persistent fibrotic diseases.

GAiN: An integrative tool utilizing generative adversarial neural networks for augmented gene expression analysis

Big genomic data and artificial intelligence (AI) are ushering in an era of precision medicine, providing opportunities to study previously under-represented subtypes and rare diseases rather than categorize them as variances. However, clinical researchers face challenges in accessing such novel technologies as well as reliable methods to study small datasets or subcohorts with unique phenotypes. To address this need, we developed an integrative approach, GAiN, to capture patterns of gene expression from small datasets on the basis of an ensemble of generative adversarial networks (GANs) while leveraging big population data. Where conventional biostatistical methods fail, GAiN reliably discovers differentially expressed genes (DEGs) and enriched pathways between two cohorts with limited numbers of samples (n = 10) when benchmarked against a gold standard. GAiN is freely available at GitHub. Thus, GAiN may serve as a crucial tool for gene expression analysis in scenarios with limited samples, as in the context of rare diseases, under-represented populations, or limited investigator resources.

SPP1, a potential therapeutic target and biomarker for lung cancer: functional insights through computational studies

NIH reported 128 different types of cancer of which lung cancer is the leading cause of mortality. Globally, it is estimated that on average one in every seventeen hospitalized patients was deceased. There are plenty of studies that have been reported on lung cancer draggability and therapeutics, but yet a protein that plays a central specific to cure the disease remains unclear. So, this study is designed to identify the possible therapeutic targets and biomarkers that can be used for the potential treatment of lung cancers. In order to identify differentially expressed genes, 39 microarray datasets of lung cancer patients were obtained from various demographic regions of the GEO database available at NCBI. After annotating statistically, 6229 up-regulated genes and 10324 down-regulated genes were found. Out of 17 up-regulated genes and significant genes, we selected SPP1 (osteopontin) through virtual screening studies. We found functional interactions with the other cancer-associated genes such as VEGF, FGA, JUN, EGFR, and TGFB1. For the virtual screening studies,198 biological compounds were retrieved from the ACNPD database and docked with SPP1 protein (PDBID: 3DSF). In the results, two highly potential compounds secoisolariciresinol diglucoside (-12.9 kcal/mol), and Hesperidin (-12.0 kcal/mol) showed the highest binding affinity. The stability of the complex was accessed by 100 ns simulation in an SPC water model. From the functional insights obtained through these computational studies, we report that SPP1 could be a potential biomarker and successive therapeutic protein target for lung cancer treatment.

Implications of stress-induced gene expression for hematopoietic stem cell aging studies

A decline in hematopoietic stem cell (HSC) function is believed to underlie hematological shortcomings with age; however, a comprehensive molecular understanding of these changes is currently lacking. Here we provide evidence that a transcriptional signature reported in several previous studies on HSC aging is linked to stress-induced changes in gene expression rather than aging. Our findings have strong implications for the design and interpretation of HSC aging studies.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

Anti-correlated feature selection prevents false discovery of subpopulations in scRNAseq

While sub-clustering cell-populations has become popular in single cell-omics, negative controls for this process are lacking. Popular feature-selection/clustering algorithms fail the null-dataset problem, allowing erroneous subdivisions of homogenous clusters until nearly each cell is called its own cluster. Using real and synthetic datasets, we find that anti-correlated gene selection reduces or eliminates erroneous subdivisions, increases marker-gene selection efficacy, and efficiently scales to millions of cells.

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle

Breeding bulls are well suited to investigate inherited variation in male fertility because they are genotyped and their reproductive success is monitored through semen analyses and thousands of artificial inseminations. However, functional data from relevant tissues are lacking in cattle, which prevents fine-mapping fertility-associated genomic regions. Here, we characterize gene expression and splicing variation in testis, epididymis, and vas deferens transcriptomes of 118 mature bulls and conduct association tests between 414,667 molecular phenotypes and 21,501,032 genome-wide variants to identify 41,156 regulatory loci. We show broad consensus in tissue-specific and tissue-enriched gene expression between the three bovine tissues and their human and murine counterparts. Expression- and splicing-mediating variants are more than three times as frequent in testis than epididymis and vas deferens, highlighting the transcriptional complexity of testis. Finally, we identify genes (WDR19, SPATA16, KCTD19, ZDHHC1) and molecular phenotypes that are associated with quantitative variation in male fertility through transcriptome-wide association and colocalization analyses.

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

LncRNA SOX9-AS1 triggers a transcriptional program involved in lipid metabolic reprogramming, cell migration and invasion in triple-negative breast cancer

At the molecular level, triple-negative breast cancer (TNBC) is frequently categorized as PAM50 basal-like subtype, but despite the advances in molecular analyses, the clinical outcome for these subtypes is uncertain. Long non-coding RNAs (lncRNAs) are master regulators of genes involved in hallmarks of cancer, which makes them suitable biomarkers for breast cancer (BRCA) diagnosis and prognosis. Here, we evaluated the regulatory role of lncRNA SOX9-AS1 in these subtypes. Using the BRCA-TCGA cohort, we observed that SOX9-AS1 was significantly overexpressed in basal-like and TNBC in comparison with other BRCA subtypes. Survival analyzes showed that SOX9-AS1 overexpression was associated with a favorable prognosis in TNBC and basal-like patients. To study the functions of SOX9-AS1, we determined the expression levels in a panel of nine BRCA cell lines finding increased levels in MDA-MB-468 and HCC1187 TNBC. Using subcellular fractionation in these cell lines, we ascertained that SOX9-AS1 was located in the cytoplasmic compartment. In addition, we performed SOX9-AS1 gene silencing using two short-harping constructs, which were transfected in both cell models and performed a genome-wide RNA-seq analysis. Data showed that 351 lncRNAs and 740 mRNAs were differentially expressed in MDA-MB-468 while 56 lncRNAs and 100 mRNAs were modulated in HCC1187 cells (Log2FC < - 1.5 and > 1.5, p.adj value < 0.05). Pathway analysis revealed that the protein-encoding genes potentially regulate lipid metabolic reprogramming, and epithelial-mesenchymal transition (EMT). Expression of lipid metabolic-related genes LIPE, REEP6, GABRE, FBP1, SCD1, UGT2B11, APOC1 was confirmed by RT-qPCR. Functional analysis demonstrated that the knockdown of SOX9-AS1 increases the triglyceride synthesis, cell migration and invasion in both two TNBC cell lines. In conclusion, high SOX9-AS1 expression predicts an improved clinical course in patients, while the loss of SOX9-AS1 expression enhances the aggressiveness of TNBC cells.

Associations between low serum levels of ANRIL and some common gene SNPs in Iranian patients with premature coronary artery disease

Coronary artery disease (CAD) is the major cause of mortality in the world. Premature development of CAD can be attributed to women under 55 and men under 45. Many genetic factors play a part in premature CAD. Among them, ANRIL, a long noncoding RNA is located at the 9p21 risk locus, and its expression seems to be correlated with CAD. In the current study, premature CAD and control blood samples, with and without Type 2 Diabetes (T2D), were genotyped for six SNPs at the 9p21 locus. Additionally, ANRIL serum expression was assessed in both groups using real-time PCR. It was performed using different primers targeting exons 1, 5-6, and 19. The χ2 test for association, along with t-tests and ANOVA, was employed for statistical analysis. In this study, we did not find any significant correlation between premature coronary artery disease and rs10757274, rs2383206, rs2383207, rs496892, rs10757278 and rs10738605. However, a lower ANRIL expression was correlated with each SNP risk genotype. Despite the correlation between lower ANRIL expression and CAD, Type 2 diabetes was associated with higher ANRIL expression. Altogether, the correlation between ANRIL expression and the genotypes of the studied SNPs indicated that genetic variants, even those in intronic regions, affect long noncoding RNA expression levels. In conclusion, we recommend combining genetic variants with expression analysis when developing screening strategies for families with premature CAD. To prevent the devastating outcomes of CAD in young adults, it is crucial to discover noninvasive genetic-based screening tests.

Expression of Inflammatory Genes in Murine Lungs in a Model of Experimental Pulmonary Hypertension: Effects of an Antibody-Based Targeted Delivery of Interleukin-9

BACKGROUND

Pathogenesis of pulmonary hypertension (PH) is a multifactorial process driven by inflammation and pulmonary vascular remodeling. To target these two aspects of PH, we recently tested a novel treatment: Interleukin-9 (IL9) fused to F8, an antibody that binds to the extra-domain A of fibronectin (EDA+ Fn). As EDA+ Fn is not found in healthy adult tissue but is expressed during PH, IL9 is delivered specifically to the tissue affected by PH. We found that F8IL9 reduced pulmonary vascular remodeling and attenuated PH compared with sham-treated mice.

PURPOSE

To evaluate possible F8IL9 effects on PH-associated inflammatory processes, we analysed the expression of genes involved in pulmonary immune responses.

METHODS

We applied the monocrotaline (MCT) model of PH in mice (n = 44). Animals were divided into five experimental groups: sham-induced animals without PH (control, n = 4), MCT-induced PH without treatment (PH, n = 8), dual endothelin receptor antagonist treatment (dual ERA, n = 8), F8IL9 treatment (n = 12, 2 formats with n = 6 each), or with KSFIL9 treatment (KSFIL9, n = 12, 2 formats with n = 6 each, KSF: control antibody with irrelevant antigen specificity). After 28 days, a RT-PCR gene expression analysis of inflammatory response (84 genes) was performed in the lung.

RESULTS

Compared with the controls, 19 genes exhibited relevant (+2.5-fold) upregulation in the PH group without treatment. Gene expression levels in F8IL9-treated lung tissue were reduced compared to the PH group without treatment. This was the case especially for CCL20, CXCL5, C-reactive protein, pentraxin related (CRPPR), and Kininogen-1 (KNG1).

CONCLUSION

In accordance with the hypothesis stated above, F8IL9 treatment diminished the upregulation of some genes associated with inflammation in a PH animal model. Therefore, we hypothesize that IL9-based immunocytokine treatment will likely modulate various inflammatory pathways.

Modelling bronchial epithelial-fibroblast cross-talk in idiopathic pulmonary fibrosis (IPF) using a human-derived in vitro air liquid interface (ALI) culture

Idiopathic Pulmonary Fibrosis (IPF) is a devastating form of respiratory disease with a life expectancy of 3-4 years. Inflammation, epithelial injury and myofibroblast proliferation have been implicated in disease initiation and, recently, epithelial-fibroblastic crosstalk has been identified as a central driver. However, the ability to interrogate this crosstalk is limited due to the absence of in vitro models that mimic physiological conditions. To investigate IPF dysregulated cross-talk, primary normal human bronchial epithelial (NHBE) cells and primary normal human lung fibroblasts (NHLF) or diseased human lung fibroblasts (DHLF) from IPF patients, were co-cultured in direct contact at the air-liquid interface (ALI). Intercellular crosstalk was assessed by comparing cellular phenotypes of co-cultures to respective monocultures, through optical, biomolecular and electrical methods. A co-culture-dependent decrease in epithelium thickness, basal cell mRNA (P63, KRT5) and an increase in transepithelial electrical resistance (TEER) was observed. This effect was significantly enhanced in DHLF co-cultures and lead to the induction of epithelial to mesenchymal transition (EMT) and increased mRNA expression of TGFβ-2, ZO-1 and DN12. When stimulated with exogenous TGFβ, NHBE and NHLF monocultures showed a significant upregulation of EMT (COL1A1, FN1, VIM, ASMA) and senescence (P21) markers, respectively. In contrast, direct NHLF/NHBE co-culture indicated a protective role of epithelial-fibroblastic cross-talk against TGFβ-induced EMT, fibroblast-to-myofibroblast transition (FMT) and inflammatory cytokine release (IL-6, IL-8, IL-13, IL-1β, TNF-α). DHLF co-cultures showed no significant phenotypic transition upon stimulation, likely due to the constitutively high expression of TGFβ isoforms prior to any exogenous stimulation. The model developed provides an alternative method to generate IPF-related bronchial epithelial phenotypes in vitro, through the direct co-culture of human lung fibroblasts with NHBEs. These findings highlight the importance of fibroblast TGFβ signaling in EMT but that monocultures give rise to differential responses compared to co-cultures, when exposed to this pro-inflammatory stimulus. This holds implications for any translation conclusions drawn from monoculture studies and is an important step in development of more biomimetic models of IPF. In summary, we believe this in vitro system to study fibroblast-epithelial crosstalk, within the context of IPF, provides a platform which will aid in the identification and validation of novel targets.

MicroRNAs modulate immunological and inflammatory responses in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

MicroRNAs (miRNAs) regulate the post-transcriptional expression of genes by binding to their target mRNAs. In this study, whole miRNA sequencing was used to compare the expression of miRNAs in ileocecal valve (ICV) and peripheral blood (PB) samples of cows with focal or diffuse paratuberculosis (PTB)-associated lesions in gut tissues versus (vs) control cows without lesions. Among the eight miRNAs differentially expressed in the PB samples from cows with diffuse lesions vs controls, three (miR-19a, miR-144, miR32) were also down-regulated in cows with diffuse vs focal lesions. In the ICV samples, we identified a total of 4, 5, and 18 miRNAs differentially expressed in cows with focal lesions vs controls, diffuse lesions vs controls, and diffuse vs focal lesions, respectively. The differential expression of five microRNAs (miR-19a, miR-144, miR-2425-3p, miR-139, miR-101) was confirmed by RT-qPCR. Next, mRNA target prediction was performed for each differentially expressed miRNA. A functional analysis using the predicted gene targets revealed a significant enrichment of the RNA polymerase and MAPK signaling pathways in the comparison of cows with focal vs no lesions and with diffuse vs focal lesions, respectively. The identified miRNAs could be used for the development of novel diagnostic and therapeutical tools for PTB control.

Slide-tags enables single-nucleus barcoding for multimodal spatial genomics

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed1-6. However, missing from these measurements is the ability to routinely and easily spatially localize these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are tagged with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as an input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 μm spatial resolution and delivered whole-transcriptome data that are indistinguishable in quality from ordinary single-nucleus RNA-sequencing data. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualized receptor-ligand interactions driving B cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to almost any single-cell measurement technology. As a proof of principle, we performed multiomic measurements of open chromatin, RNA and T cell receptor (TCR) sequences in the same cells from metastatic melanoma, identifying transcription factor motifs driving cancer cell state transitions in spatially distinct microenvironments. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

Identification and expression analysis of MAPK cascade gene family in foxtail millet (Setaria italica)

The mitogen-activated protein kinase (MAPK) cascade pathway is a highly conserved plant cell signaling pathway that plays an important role in plant growth and development and stress response. Currently, MAPK cascade genes have been identified and reported in a variety of plants including Arabidopsis thaliana, Oryza sativa, and Triticum aestivum, but have not been identified in foxtail millet (Setaria italica). In this study, a total of 93 MAPK cascade genes, including 15 SiMAPKs, 10 SiMAPKKs and 68 SiMAPKKKs genes, were identified by genome-wide analysis of foxtail millet, and these genes were distributed on nine chromosomes of foxtail millet. Using phylogenetic analysis, we divided the SiMAPKs and SiMAPKKs into four subgroups, respectively, and the SiMAPKKKs into three subgroups (Raf, ZIK, and MEKK). Whole-genome duplication analysis revealed that there are 14 duplication pairs in the MAPK cascade family in foxtail millet, and they are expanded by segmental replication events. Results from quantitative real-time PCR (qRT-PCR) revealed that the expression levels of most SiMAPKs and SiMAPKKs were changed under both exogenous hormone and abiotic stress treatments, with SiMAPK3 and SiMAPKK4-2 being induced under almost all treatments, while the expression of SiMAPKK5 was repressed. In a nutshell, this study will shed some light on the evolution of MAPK cascade genes and the functional mechanisms underlying MAPK cascade genes in response to hormonal and abiotic stress signaling pathways in foxtail millet (Setaria italica).

Analysis of Sporulation in Bacillus cereus Biovar anthracis Which Contains an Insertion in the Gene for the Sporulation Factor σK

Bacillus cereus biovar anthracis (Bcbva) is an untypical pathogen causing a fatal anthrax-like disease in a variety of wildlife species in African rainforest areas. In contrast to Bacillus anthracis and most species of the B. cereus group, all strains of the Bcbva cluster contain a 22 kb insertion in the sigK gene which encodes the essential late sporulation sigma factor σK. This insertion is excised during sporulation in a site-specific recombination process resulting in an intact sigK gene and a circular molecule. The sporulation kinetics of two strains each of Bcbva and B. anthracis were compared by the expression analysis of eight sporulation-associated genes, including sigK, using reverse transcriptase quantitative real-time PCR. In addition, morphological sporulation stages were analyzed and quantified by electron microscopy. Our results indicated that the necessary excision of the insertion in Bcbva neither delayed nor inhibited its sporulation. In two spontaneous mutants of Bcbva, the excision of the sigK insertion and sporulation were impeded due to mutations in the spo0A and spoVG regulator genes, respectively. The spo0A frameshift mutation was overcome by intragenic suppression in a revertant which was able to sporulate normally, despite an M171S amino acid exchange in the global regulator Spo0A. A screening of the NCBI database identified further strains of the B. cereus group which possess unrelated insertions in the sigK gene, and two strains containing almost identical insertions at the same gene position. Some of the sigK insertions encode putative prophages, whereas the Bcbva insertion encoded a type I restriction-modification system. The function of these insertions and if they are possibly essential for sporulation remains to be assessed.

Identification of ITPR1 gene as a novel target for hsa-miR-34b-5p in non-obstructive azoospermia: a Ca2+/apoptosis pathway cross-talk

MiR-34b-5p has been reported as a non-invasive diagnostic biomarker for infertility. However, no gene targets regulating the mechanism of cation of this miRNA are known. In this study, using gene set enrichment analysis the Inositol 1,4,5-Trisphosphate Receptor Type 1 (ITPR1) gene was identified as the sole target for hsa-miR-34b-5p, and found significantly overexpressed in non-obstructive azoospermia (NOA) patients. This finding was confirmed by qRT-PCR on fresh testicular tissues from NOA patients. Then, pathway enrichment analysis as well as the diagnostic value analysis of hsa-miR-34b-5p/ITPR1 indicated ITPR1 as a hub gene in the calcium (Ca2+)-apoptosis pathway, and a valuable predictive biomarker for NOA. Moreover, gene expression and histological assays showed the association of the effects of ITPR1's increased expression on spermatogenesis failure through induction of apoptosis in NOA patients. These data suggested that the hsa-miR-34b-5p/ITPR1 axis could serve as a potential regulatory predictive biomarker for human spermatogenesis through the Ca2+-apoptosis pathway cross-talk.

A simplified and rapid in situ hybridization protocol for planarians

Whole-mount in situ hybridization is a critical technique for analyzing gene expression in planarians. While robust in situ protocols have been developed, these protocols are laborious, making them challenging to incorporate in an academic setting, reducing throughput and increasing time to results. Here, the authors systematically tested modifications to all phases of the protocol with the goal of eliminating steps and reducing time without impacting quality. This modified protocol allows for whole-mount colorimetric in situ hybridization and multicolor fluorescence in situ hybridization to be completed in two days with a significant reduction in steps and hands-on processing time.

KRAS4A and KRAS4B in liquid biopsy of metastatic lung adenocarcinoma patients treated with Pembrolizumab or chemotherapy plus Pembrolizumab

KRAS is involved in the stability and expression of PD-L1. We investigated the expression of circulating mRNA (cmRNA) of KRAS4A and KRAS4B and the possible impact on progression-free survival (PFS) of patients with metastatic lung adenocarcinoma treated with immunotherapy. Patients without driver mutations undergoing Pembrolizumab (P) or P plus chemotherapy (PC) were prospectively accrued for liquid biopsy analysis of KRAS4A, KRAS4B, and PD-L1 cmRNA. Both KRAS isoforms were also studied for association with PD-L1 cmRNA. Of 56 patients, 28 received P and 28 PC. Patients with high levels of both KRAS isoforms showed significantly better PFS. The median PFS for KRAS4A was 29 months (95% CI 22-29 months) and KRAS4B 24 months (95% CI 13-29 months), respectively. The median PFS of patients with low levels of both isoforms was 12 months (95% CI 6-15 months for KRAS4A and 95% CI 5-20 months for KRAS4B). High KRAS4A retained a significant positive association with PFS in the multivariate model. An exploratory analysis in treatment subgroups found a positive association between high KRAS4A and KRAS4B with PFS in patients treated with P. PD-L1 cmRNA was significantly higher in patients with high KRAS isoforms levels and this effect was pronounced for high KRAS4A carriers. KRAS4A deserves further investigation as a potential marker for defining patients who may benefit the most from immune checkpoint inhibitors therapy and improving personalized cancer immunotherapeutic strategies.

A comparative analysis of genes differentially expressed between rete testis cells and Sertoli cells of the mouse testis

The rete testis (RT) is a region of the mammalian testis that plays an important role in testicular physiology. The RT epithelium consists of cells sharing some well-known gene markers with supporting Sertoli cells (SCs). However, little is known about the differences in gene expression between these two cell populations. Here, we used fluorescence-activated cell sorting (FACS) to obtain pure cultures of neonatal RT cells and SCs and identified differentially expressed genes (DEGs) between these cell types using RNA sequencing (RNA-seq). We then compared our data with the RNA-seq data of other studies that examined RT cells and SCs of mice of different ages and generated a list of DEGs permanently upregulated in RT cells throughout testis development and in culture, which included 86 genes, and a list of 79 DEGs permanently upregulated in SCs. The analysis of studies on DMRT1 function revealed that nearly half of the permanent DEGs could be regulated by this SC upregulated transcription factor. We suggest that useful cell lineage markers and candidate genes for the specification of both RT cells and SCs may be present among these permanent DEGs.

Cochlear organoids reveal transcriptional programs of postnatal hair cell differentiation from supporting cells

We explore the changes in chromatin accessibility and transcriptional programs for cochlear hair cell differentiation from postmitotic supporting cells using organoids from postnatal cochlea. The organoids contain cells with transcriptional signatures of differentiating vestibular and cochlear hair cells. Construction of trajectories identifies Lgr5+ cells as progenitors for hair cells, and the genomic data reveal gene regulatory networks leading to hair cells. We validate these networks, demonstrating dynamic changes both in expression and predicted binding sites of transcription factors (TFs) during organoid differentiation. We identify known regulators of hair cell development, Atoh1, Pou4f3, and Gfi1, and the analysis predicts the regulatory factors Tcf4, an E-protein and heterodimerization partner of Atoh1, and Ddit3, a CCAAT/enhancer-binding protein (C/EBP) that represses Hes1 and activates transcription of Wnt-signaling-related genes. Deciphering the signals for hair cell regeneration from mammalian cochlear supporting cells reveals candidates for hair cell (HC) regeneration, which is limited in the adult.

Fragment-sequencing unveils local tissue microenvironments at single-cell resolution

Cells collectively determine biological functions by communicating with each other-both through direct physical contact and secreted factors. Consequently, the local microenvironment of a cell influences its behavior, gene expression, and cellular crosstalk. Disruption of this microenvironment causes reciprocal changes in those features, which can lead to the development and progression of diseases. Hence, assessing the cellular transcriptome while simultaneously capturing the spatial relationships of cells within a tissue provides highly valuable insights into how cells communicate in health and disease. Yet, methods to probe the transcriptome often fail to preserve native spatial relationships, lack single-cell resolution, or are highly limited in throughput, i.e. lack the capacity to assess multiple environments simultaneously. Here, we introduce fragment-sequencing (fragment-seq), a method that enables the characterization of single-cell transcriptomes within multiple spatially distinct tissue microenvironments. We apply fragment-seq to a murine model of the metastatic liver to study liver zonation and the metastatic niche. This analysis reveals zonated genes and ligand-receptor interactions enriched in specific hepatic microenvironments. Finally, we apply fragment-seq to other tissues and species, demonstrating the adaptability of our method.

A pan-cancer analysis of the role of argininosuccinate synthase 1 in human tumors

Aim

There is accumulating evidence indicating that ASS1 is closely related to tumors. No pan-cancer analysis of ASS1 was available.

Methods

Here we explored the gene expression and survival analysis of ASS1 across thirty-three tumors based on the datasets of the TCGA (Cancer Genome Atlas), the GEO (Gene Expression Omnibus), and the GEPIA2 (Gene Expression Profiling Interactive Analysis, version 2).

Results

ASS1 is highly expressed in most normal tissues and is related to the progression of some tumors. We also report ASS1 genetic alteration and their association with tumor prognosis and report differences in ASS1 phosphorylation sites between tumors and control normal tissues. ASS1 expression was associated with the infiltration of cancer-associated fibroblasts (CAFs) for the TCGA tumors of BRCA (Breast invasive carcinoma), CESC (Cervical squamous cell carcinoma and endocervical adenocarcinoma), COAD (Colon adenocarcinoma), ESCA (Esophageal carcinoma), SKCM (Skin cutaneous melanoma), SKCM-Metastasis, TGCT (Testicular germ cell tumors), and endothelial cell for the tumors of BRCA, BRCA-Basal, CESC, ESCA, KIRC (Kidney renal clear cell carcinoma), LUAD (Lung adenocarcinoma), LUSC (Lung squamous cell carcinoma), SKCM, SKCM-Metastasis, SKCM-Primary, STAD (Stomach adenocarcinoma), and TGCT. The KEGG and GO analysis were used to analyze ASS1-related signaling pathways. Finally, we used Huh7 cell line to verify the function of ASS1 in vitro. After ASS1 knockdown using small interfering RNA (siRNA), the proliferation and invasion of Huh7 were enhanced, cyclin D1 was up-regulated, and anti-apoptotic protein bax was down-regulated, suggesting that ASS1 is a tumor suppressor gene in hepatocellular carcinoma.

Conclusion

Our first pan-cancer study offers a relatively comprehensive understanding of the roles of ASS1 in different tumors.

Determining the biocontrol capacities of Trichoderma spp. originating from Turkey on Fusarium culmorum by transcriptional and antagonistic analyses

In this study aiming to investigate potential fungal biocontrol agents for Fusarium culmorum, several isolates of Trichoderma spp. were evaluated for their antagonistic effects by means of transcriptional analyses. At first, 21 monosporic Trichoderma spp. isolates were obtained from natural wood debris and wood area soils in Manisa, Turkey. Trichoderma spp. Isolates were identified as belonging to four different species (T. atroviride, T. harzianum, T. koningii, and T. brevicompactum) by tef1-α sequencing. Then, the linear growth rate (LGR) of each species was calculated and determined to be in a range between 13.22 ± 0.71 mm/day (T. atroviride TR2) and 25.06 ± 1.45 mm/day (T. harzianum K30). Inter-simple sequence repeat (ISSR) genotyping validated the tef1-α sequencing results by presenting two sub-clusters in the dendrogram. We determined the genetically most similar (TR1 & TR2; 97.77%) and dissimilar (K9 & K17; 40.40%) individuals belonging to the same and different species, respectively. Dual sandwich culture tests (which are useful for antagonism studies) revealed that T. harzianum K21 (the least suppressive) and T. brevicompactum K26 (the most suppressive) isolates suppressed F. culmorum with growth rates of 3% and 46%, respectively. Expressions of genes previously associated with mycoparasitism-plant protection-secondary metabolism (nag1, tgf-1, and tmk-1) were tested by quantitative real-time polymerase chain reaction (qRT-PCR) in both those isolates. While there were no significant differences (p>0.05) in expression that were present in the K21 isolate, those three genes were upregulated with fold change values of 2.69 ± 0.26 (p<0.001), 2.23 ± 0.16 (p<0.001), and 5.38 ± 2.01 (p<0.05) in K26, meaning that the presence of significant alteration in the physiological processes of the fungus. Also, its mycoparasitism potential was tested on Triticum aestivum L. cv Basribey in planta, which was infected with the F. culmorum FcUK99 strain. Results of the trials, including specific plant growth parameters (weight or length of plantlets), confirmed the mycoparasitic potential of the isolate. It can be concluded that (i) nag1, tgf-1, and tmk-1 genes could be approved as reliable markers for evaluation of BCA capacities of Trichoderma spp. and (ii) the T. brevicompactum K26 strain can be suggested as a promising candidate for combating in F. culmorum diseases following the necessary procedures to ensure it is non-hazardous and safe.

Transcriptomic analysis reveals that RasGEF1b deletion alters basal and LPS-induced expression of genes involved in chemotaxis and cytokine responses in macrophages

Ras guanine nucleotide exchange factor member 1b (RasGEF1b) of the RasGEF/CDC25 domain-containing family is preferentially expressed by macrophages. However, information is lacking about its role in macrophage function. In this study, we generated mice with ubiquitous deletion of Rasgef1b and used RNA-seq-based transcriptomics to compare the global gene expression in wild-type and knock-out primary bone-marrow-derived macrophages under basal conditions and after lipopolysaccharide (LPS) treatment. Transcriptional filtering identified several genes with significantly different transcript levels between wild-type and knock-out macrophages. In total, 49 and 37 differentially expressed genes were identified at baseline and in LPS-activated macrophages, respectively. Distinct biological processes were significantly linked to down-regulated genes at the basal condition only, and largely included chemotaxis, response to cytokines, and positive regulation of GTPase activity. Importantly, validation by RT-qPCR revealed that the expression of genes identified as down-regulated after LPS stimulation was also decreased in the knock-out cells under basal conditions. We used a luciferase-based reporter assay to showcase the capability of RasGEF1b in activating the Serpinb2 promoter. Notably, knockdown of RasGEF1b in RAW264.7 macrophages resulted in impaired transcriptional activation of the Serpinb2 promoter, both in constitutive and LPS-stimulated conditions. This study provides a small collection of genes that shows relative expression changes effected by the absence of RasGEF1b in macrophages. Thus, we present the first evidence that RasGEF1b mediates the regulation of both steady-state and signal-dependent expression of genes and propose that this GEF plays a role in the maintenance of the basal transcriptional level in macrophages.

Characterization of the FLA Gene Family in Tomato (Solanum lycopersicum L.) and the Expression Analysis of SlFLAs in Response to Hormone and Abiotic Stresses

Fasciclin-like arabinogalactan proteins (FLAs), a subclass of arabinogalactan proteins (AGPs), participate in mediating plant growth, development, and response to abiotic stress. However, the characterization and function of FLAs in tomato are currently unknown. In this study, members of the tomato FLA family are characterized and analyzed in relation to their response to phytohormonal and abiotic stresses. The results show that a total of 24 FLA members were characterized in tomato. The structural domain analysis showed that these members have a high protein similarity. The expression profiles of different tissues indicated that the genes of most members of the tomato FLA gene family are highly expressed in roots, but to a lower extent in fruits. qRT-PCR analysis revealed that all 24 tomato FLA genes are responsive to ABA and MeJA. SlFLAs showed a positive response to salt and cold stress. SlFLA1, SlFLA12, and SlFLA14 are significantly induced under darkness. SlFLA1 and SlFLA3 are significantly induced under drought stress. This study provides a basis for a further understanding of the role of tomato FLA homologous genes in plant response to abiotic stress and lays the foundation for further research on the function of FLAs in tomato.

Comprehensive and evolutionary analysis of Spodoptera litura-inducible Cytochrome P450 monooxygenase gene family in Glycine max elucidate their role in defense

Plants being sessile organisms and lacking both circulating phagocytic cells and somatic adaptive immune response, have thrived on various defense mechanisms to fend off insect pests and invasion of pathogens. CYP450s are the versatile enzymes, which thwart plants against insect pests by ubiquitous biosynthesis of phytohormones, antioxidants, and secondary metabolites, utilizing them as feeding deterrents and direct toxins. Therefore, a comprehensive analysis of biotic stress-responsive CYPs from Glycine max was performed to ascertain their function against S. litura-infestation. Phylogenetic analysis and evolutionary studies on conserved domains and motifs disclosed the evolutionary correspondence of these GmCYPs with already characterized members of the CYP450 superfamily and close relatedness to Medicago truncatula. These GmCYPs were mapped on 13 chromosomes; they possess 1-8 exons; they have evolved due to duplication and are localized in endoplasmic reticulumn. Further, identification of methyl-jasmonate, salicylic acid, defense responsive and flavonoid biosynthesis regulating cis-acting elements, their interaction with biotic stress regulating proteins and their differential expression in diverse types of tissues, and during herbivory, depicted their responsiveness to biotic stress. Three-dimensional homology modelling of GmCYPs, docking with heme cofactor required for their catalytic activity and enzyme-substrate interactions were performed to understand the functional mechanism of their action. Moreover, to gain insight into their involvement in plant defense, gene expression analysis was evaluated, which revealed differential expression of 11 GmCYPs upon S. litura-infestation, 12 GmCYPs on wounding while foliar spray of ethylene, methyl-jasmonate and salicylic acid differentially regulated 11 GmCYPs, 6 GmCYPs, and 10 GmCYPs respectively. Our study comprehensively analysed the underlying mechanism of GmCYPs function during S. litura-infestation, which can be further utilized for functional characterization to develop new strategies for enhancing soybean resistance to insect pests.

Clinical features and prognostic outcomes of angioimmunoblastic T cell lymphoma in an Asian multicenter study

In our Asian multicenter retrospective study, we investigated the clinical prognostic factors affecting the outcomes of AITL patients and identified a novel prognostic index relevant in the Asian context. In our 174-patient cohort, the median PFS and OS was 1.8 years and 5.6 years respectively. Age > 60, bone marrow involvement, total white cell count >12 × 109/L and raised serum lactate dehydrogenase were associated with poorer PFS and OS in multivariate analyses. This allowed for a prognostic index (AITL-PI) differentiating patients into low (0-1 factors, n = 64), moderate (2 factors, n = 59) and high-risk (3-4 factors, n = 49) subgroups with 5-year OS of 84.0%, 44.0% and 28.0% respectively (p < 0.0001). POD24 proved to be strongly prognostic (5-year OS 24% vs 89%, p < 0.0001). Exploratory gene expression studies were performed and disparate immune cell profiles and cell signaling signatures were seen in the low risk group as compared to the intermediate and high risk groups.

Genomic disturbance of vitellogenin 2 (vtg2) leads to vitellin membrane deficiencies and significant mortalities at early stages of embryonic development in zebrafish (Danio rerio)

The specific functions and essentiality of type II vitellogenin (Vtg2) in early zebrafish development were investigated in this study. A vtg2-mutant zebrafish line was produced and effects of genomic disturbance were observed in F2 females and F3 offspring. No change in vtg2 transcript has been detected, however, Vtg2 abundance in F2 female liver was 5×, and in 1 hpf F3 vtg2-mutant embryos was 3.8× less than Wt (p < 0.05). Fecundity was unaffected while fertilization rate was more than halved in F2 vtg2-mutant females (p < 0.05). Hatching rate was significantly higher in F3 vtg2-mutant embryos in comparison to Wt embryos. Survival rate declined drastically to 29% and 18% at 24 hpf and 20 dpf, respectively, in F3 vtg2-mutant embryos. The introduced mutation caused vitelline membrane deficiencies, significant mortalities at early embryonic stages, and morphological abnormalities in the surviving F3 vtg2-mutant larvae. Overrepresentation of histones, zona pellucida proteins, lectins, and protein degradation related proteins in F3 vtg2-mutant embryos provide evidence to impaired mechanisms involved in vitellin membrane formation. Overall findings imply a potential function of Vtg2 in acquisition of vitellin membrane integrity, among other reproductive functions, and therefore, its essentiality in early zebrafish embryo development.

Microfluidics-free single-cell genomics with templated emulsification

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.