The CD4 Versus CD8 T Cell Fate Decision: A Multiomics-Informed Perspective.

The choice of developing thymocytes to become CD8+ cytotoxic or CD4+ helper T cells has been intensely studied, but many of the underlying mechanisms remain to be elucidated. Recent multiomics approaches have provided much higher resolution analysis of gene expression in developing thymocytes than was previously achievable, thereby offering a fresh perspective on this question. Focusing on our recent studies using CITE-seq (cellular indexing of transcriptomes and epitopes) analyses of mouse thymocytes, we present a detailed timeline of RNA and protein expression changes during CD8 versus CD4 T cell differentiation. We also revisit our current understanding of the links between T cell receptor signaling and expression of the lineage-defining transcription factors ThPOK and RUNX3. Finally, we propose a sequential selection model to explain the tight linkage between MHC-I versus MHC-II recognition and T cell lineage choice. This model incorporates key aspects of previously proposed kinetic signaling, instructive, and stochastic/selection models.

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data.

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Molecular and Spatial Signatures of Mouse Embryonic Endothelial Cells at Single-Cell Resolution.

BACKGROUND: Mature endothelial cells (ECs) are heterogeneous, with subtypes defined by tissue origin and position within the vascular bed (ie, artery, capillary, vein, and lymphatic). How this heterogeneity is established during the development of the vascular system, especially arteriovenous specification of ECs, remains incompletely characterized. METHODS: We used droplet-based single-cell RNA sequencing and multiplexed error-robust fluorescence in situ hybridization to define EC and EC progenitor subtypes from E9.5, E12.5, and E15.5 mouse embryos. We used trajectory inference to analyze the specification of arterial ECs (aECs) and venous ECs (vECs) from EC progenitors. Network analysis identified candidate transcriptional regulators of arteriovenous differentiation, which we tested by CRISPR (clustered regularly interspaced short palindromic repeats) loss of function in human-induced pluripotent stem cells undergoing directed differentiation to aECs or vECs (human-induced pluripotent stem cell-aECs or human-induced pluripotent stem cell-vECs). RESULTS: From the single-cell transcriptomes of 7682 E9.5 to E15.5 ECs, we identified 19 EC subtypes, including Etv2+Bnip3+ EC progenitors. Spatial transcriptomic analysis of 15 448 ECs provided orthogonal validation of these EC subtypes and established their spatial distribution. Most embryonic ECs were grouped by their vascular-bed types, while ECs from the brain, heart, liver, and lung were grouped by their tissue origins. Arterial (Eln, Dkk2, Vegfc, and Egfl8), venous (Fam174b and Clec14a), and capillary (Kcne3) marker genes were identified. Compared with aECs, embryonic vECs and capillary ECs shared fewer markers than their adult counterparts. Early capillary ECs with venous characteristics functioned as a branch point for differentiation of aEC and vEC lineages. CONCLUSIONS: Our results provide a spatiotemporal map of embryonic EC heterogeneity at single-cell resolution and demonstrate that the diversity of ECs in the embryo arises from both tissue origin and vascular-bed position. Developing aECs and vECs share common venous-featured capillary precursors and are regulated by distinct transcriptional regulatory networks.

mTORC1 Signaling in Brain Endothelial Progenitors Contributes to CCM Pathogenesis.

BACKGROUND: Cerebral vascular malformations (CCMs) are primarily found within the brain, where they result in increased risk for stroke, seizures, and focal neurological deficits. The unique feature of the brain vasculature is the blood-brain barrier formed by the brain neurovascular unit. Recent studies suggest that loss of CCM genes causes disruptions of blood-brain barrier integrity as the inciting events for CCM development. CCM lesions are proposed to be initially derived from a single clonal expansion of a subset of angiogenic venous capillary endothelial cells (ECs) and respective resident endothelial progenitor cells (EPCs). However, the critical signaling events in the subclass of brain ECs/EPCs for CCM lesion initiation and progression are unclear. METHODS: Brain EC-specific CCM3-deficient (Pdcd10BECKO) mice were generated by crossing Pdcd10fl/fl mice with Mfsd2a-CreERT2 mice. Single-cell RNA-sequencing analyses were performed by the chromium single-cell platform (10× genomics). Cell clusters were annotated into EC subtypes based on visual inspection and GO analyses. Cerebral vessels were visualized by 2-photon in vivo imaging and tissue immunofluorescence analyses. Regulation of mTOR (mechanistic target of rapamycin) signaling by CCM3 and Cav1 (caveolin-1) was performed by cell biology and biochemical approaches. RESULTS: Single-cell RNA-sequencing analyses from P10 Pdcd10BECKO mice harboring visible CCM lesions identified upregulated CCM lesion signature and mitotic EC clusters but decreased blood-brain barrier-associated EC clusters. However, a unique EPC cluster with high expression levels of stem cell markers enriched with mTOR signaling was identified from early stages of the P6 Pdcd10BECKO brain. Indeed, mTOR signaling was upregulated in both mouse and human CCM lesions. Genetic deficiency of Raptor (regulatory-associated protein of mTOR), but not of Rictor (rapamycin-insensitive companion of mTOR), prevented CCM lesion formation in the Pdcd10BECKO model. Importantly, the mTORC1 (mTOR complex 1) pharmacological inhibitor rapamycin suppressed EPC proliferation and ameliorated CCM pathogenesis in Pdcd10BECKO mice. Mechanistic studies suggested that Cav1/caveolae increased in CCM3-depleted EPC-mediated intracellular trafficking and complex formation of the mTORC1 signaling proteins. CONCLUSIONS: CCM3 is critical for maintaining blood-brain barrier integrity and CCM3 loss-induced mTORC1 signaling in brain EPCs initiates and facilitates CCM pathogenesis.

Directed Differentiation of Human Induced Pluripotent Stem Cells to Heart Valve Cells.

BACKGROUND: A main obstacle in current valvular heart disease research is the lack of high-quality homogeneous functional heart valve cells. Human induced pluripotent stem cells (hiPSCs)-derived heart valve cells may help with this dilemma. However, there are no well-established protocols to induce hiPSCs to differentiate into functional heart valve cells, and the networks that mediate the differentiation have not been fully elucidated. METHODS: To generate heart valve cells from hiPSCs, we sequentially activated the Wnt, BMP4, VEGF (vascular endothelial growth factor), and NFATc1 signaling pathways using CHIR-99021, BMP4, VEGF-165, and forskolin, respectively. The transcriptional and functional similarity of hiPSC-derived heart valve cells compared with primary heart valve cells were characterized. Longitudinal single-cell RNA sequencing was used to uncover the trajectory, switch genes, pathways, and transcription factors of the differentiation. RESULTS: An efficient protocol was developed to induce hiPSCs to differentiate into functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells. After 6-day differentiation and CD144 magnetic bead sorting, ≈70% CD144+ cells and 30% CD144- cells were obtained. On the basis of single-cell RNA sequencing data, the CD144+ cells and CD144- cells were found to be highly similar to primary heart valve endothelial cells and primary heart valve interstitial cells in gene expression profile. Furthermore, CD144+ cells had the typical function of primary heart valve endothelial cells, including tube formation, uptake of low-density lipoprotein, generation of endothelial nitric oxide synthase, and response to shear stress. Meanwhile, CD144- cells could secret collagen and matrix metalloproteinases, and differentiate into osteogenic or adipogenic lineages like primary heart valve interstitial cells. Therefore, we identified CD144+ cells and CD144- cells as hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells, respectively. Using single-cell RNA sequencing analysis, we demonstrated that the trajectory of heart valve cell differentiation was consistent with embryonic valve development. We identified the main switch genes (NOTCH1, HEY1, and MEF2C), signaling pathways (TGF-ß, Wnt, and NOTCH), and transcription factors (MSX1, SP5, and MECOM) that mediated the differentiation. Finally, we found that hiPSC-derived valve interstitial-like cells might derive from hiPSC-derived valve endothelial-like cells undergoing endocardial-mesenchymal transition. CONCLUSIONS: In summary, this is the first study to report an efficient strategy to generate functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells from hiPSCs, as well as to elucidate the differentiation trajectory and transcriptional dynamics of hiPSCs differentiated into heart valve cells.

Interferon-γ inducible factor 16 (IFI16) restricts adeno-associated virus type 2 (AAV2) transduction in an immune-modulatory independent way.

We determined the transcription profile of adeno-associated virus type 2 (AAV2)-infected primary human fibroblasts. Subsequent analysis revealed that cells respond to AAV infection through changes in several significantly affected pathways, including cell cycle regulation, chromatin modulation, and innate immune responses. Various assays were performed to validate selected differentially expressed genes and to confirm not only the quality but also the robustness of the raw data. One of the genes upregulated in AAV2-infected cells was interferon-γ inducible factor 16 (IFI16). IFI16 is known as a multifunctional cytosolic and nuclear innate immune sensor for double-stranded as well as single-stranded DNA, exerting its effects through various mechanisms, such as interferon response, epigenetic modifications, or transcriptional regulation. IFI16 thereby constitutes a restriction factor for many different viruses among them, as shown here, AAV2 and thereof derived vectors. Indeed, the post-transcriptional silencing of IFI16 significantly increased AAV2 transduction efficiency, independent of the structure of the virus/vector genome. We also show that IFI16 exerts its inhibitory effect on AAV2 transduction in an immune-modulatory independent way by interfering with Sp1-dependent transactivation of wild-type AAV2 and AAV2 vector promoters. IMPORTANCE: Adeno-associated virus (AAV) vectors are among the most frequently used viral vectors for gene therapy. The lack of pathogenicity of the parental virus, the long-term persistence as episomes in non-proliferating cells, and the availability of a variety of AAV serotypes differing in their cellular tropism are advantageous features of this biological nanoparticle. To deepen our understanding of virus-host interactions, especially in terms of antiviral responses, we present here the first transcriptome analysis of AAV serotype 2 (AAV2)-infected human primary fibroblasts. Our findings indicate that interferon-γ inducible factor 16 acts as an antiviral factor in AAV2 infection and AAV2 vector-mediated cell transduction in an immune-modulatory independent way by interrupting the Sp1-dependent gene expression from viral or vector genomes.

UNMF: a unified nonnegative matrix factorization for multi-dimensional omics data.

Factor analysis, ranging from principal component analysis to nonnegative matrix factorization, represents a foremost approach in analyzing multi-dimensional data to extract valuable patterns, and is increasingly being applied in the context of multi-dimensional omics datasets represented in tensor form. However, traditional analytical methods are heavily dependent on the format and structure of the data itself, and if these change even slightly, the analyst must change their data analysis strategy and techniques and spend a considerable amount of time on data preprocessing. Additionally, many traditional methods cannot be applied as-is in the presence of missing values in the data. We present a new statistical framework, unified nonnegative matrix factorization (UNMF), for finding informative patterns in messy biological data sets. UNMF is designed for tidy data format and structure, making data analysis easier and simplifying the development of data analysis tools. UNMF can handle a wide range of data structures and formats, and works seamlessly with tensor data including missing observations and repeated measurements. The usefulness of UNMF is demonstrated through its application to several multi-dimensional omics data, offering user-friendly and unified features for analysis and integration. Its application holds great potential for the life science community. UNMF is implemented with R and is available from GitHub (https://github.com/abikoushi/moltenNMF).

Model-based evaluation of spatiotemporal data reduction methods with unknown ground truth through optimal visualization and interpretability metrics.

Optimizing and benchmarking data reduction methods for dynamic or spatial visualization and interpretation (DSVI) face challenges due to many factors, including data complexity, lack of ground truth, time-dependent metrics, dimensionality bias and different visual mappings of the same data. Current studies often focus on independent static visualization or interpretability metrics that require ground truth. To overcome this limitation, we propose the MIBCOVIS framework, a comprehensive and interpretable benchmarking and computational approach. MIBCOVIS enhances the visualization and interpretability of high-dimensional data without relying on ground truth by integrating five robust metrics, including a novel time-ordered Markov-based structural metric, into a semi-supervised hierarchical Bayesian model. The framework assesses method accuracy and considers interaction effects among metric features. We apply MIBCOVIS using linear and nonlinear dimensionality reduction methods to evaluate optimal DSVI for four distinct dynamic and spatial biological processes captured by three single-cell data modalities: CyTOF, scRNA-seq and CODEX. These data vary in complexity based on feature dimensionality, unknown cell types and dynamic or spatial differences. Unlike traditional single-summary score approaches, MIBCOVIS compares accuracy distributions across methods. Our findings underscore the joint evaluation of visualization and interpretability, rather than relying on separate metrics. We reveal that prioritizing average performance can obscure method feature performance. Additionally, we explore the impact of data complexity on visualization and interpretability. Specifically, we provide optimal parameters and features and recommend methods, like the optimized variational contractive autoencoder, for targeted DSVI for various data complexities. MIBCOVIS shows promise for evaluating dynamic single-cell atlases and spatiotemporal data reduction models.

Division-Independent Differentiation of Muscle Stem Cells During a Growth Stimulus.

Adult muscle stem cells (MuSCs) are known to replicate upon activation before differentiating and fusing to regenerate myofibers. It is unclear whether MuSC differentiation is intrinsically linked to cell division, which has implications for stem cell population maintenance. We use single-cell RNA-sequencing to identify transcriptionally diverse subpopulations of MuSCs after 5 days of a growth stimulus in adult muscle. Trajectory inference in combination with a novel mouse model for tracking MuSC-derived myonuclei and in vivo labeling of DNA replication revealed an MuSC population that exhibited division-independent differentiation and fusion. These findings demonstrate that in response to a growth stimulus in the presence of intact myofibers, MuSC division is not obligatory.

Single-Cell Meta-Analysis Uncovers the Pancreatic Endothelial Cell Transcriptomic Signature and Reveals a Key Role for NKX2-3 in PLVAP Expression.

BACKGROUND: The pancreatic vasculature displays tissue-specific physiological and functional adaptations that support rapid insulin response by ß-cells. However, the digestive enzymes have made it difficult to characterize pancreatic endothelial cells (ECs), resulting in the poor understanding of pancreatic EC specialization. METHODS: Available single-nuclei/single-cell RNA-sequencing data sets were mined to identify pancreatic EC-enriched signature genes and to develop an integrated atlas of human pancreatic ECs. We validated the findings using independent single-nuclei/single-cell RNA-sequencing data, bulk RNA-sequencing data of isolated ECs, spatial transcriptomics data, immunofluorescence, and RNAScope of selected markers. The TF (transcription factor) NKX2-3 was expressed in HUVECs via gene transfection, and the expression of pancreatic EC-enriched signature genes was assessed via RT-qPCR. RESULTS: We defined a pancreatic EC-enriched gene signature conserved across species and developmental stages that included genes involved in ECM (extracellular matrix) composition (COL15A1 and COL4A1), permeability and barrier function (PLVAP, EHD4, CAVIN3, HSPG2, ROBO4, HEG1, and CLEC14A), and key signaling pathways (S1P, TGF-ß [transforming growth factor-ß], RHO-RAC GTPase, PI3k-AKT, and PDGF [platelet-derived growth factor]). The integrated atlas revealed the vascular hierarchy within the pancreas. We identified and validated a specialized islet capillary subpopulation characterized by genes involved in permeability (PLVAP and EHD4), immune-modulation (FABP5, HLA-C, and B2M), ECM composition (SPARC and SPARCL1), IGF (insulin-like growth factor) signaling (IGFBP7), and membrane transport (SLCO2A1, SLC2A3, and CD320). Importantly, we identified NKX2-3 as a key TF enriched in pancreatic ECs. DNA-binding motif analysis found NKX2-3 motifs in ≈40% of the signature genes. Induction of NKX2-3 in HUVECs promoted the expression of the islet capillary EC-enriched genes PLVAP and SPARCL1. CONCLUSIONS: We defined a validated transcriptomic signature of pancreatic ECs and uncovered their intratissue transcriptomic heterogeneity. We showed that NKX2-3 acts upstream of PLVAP and provided a single-cell online resource that can be further explored by the community: https://vasconcelos.shinyapps.io/pancreatic_endothelial/.

Single-Cell Transcriptome Reveals Potential Mechanisms for Coronary Artery Lesions in Kawasaki Disease.

BACKGROUND: Coronary artery lesions (CALs) are the most common and major complication of Kawasaki disease (KD) in developed countries. However, the underlying immunologic mechanisms of CAL development in KD remain unclear. METHODS: Here, we conducted single-cell transcriptome analyses of 212 210 peripheral blood mononuclear cells collected from a cross-sectional cohort of 16 children, including 4 patients with KD with CALs, 5 patients with KD without CALs, 4 healthy controls, and 3 febrile controls. RESULTS: KD altered the proportion of peripheral blood mononuclear cells, including an increasing trend in inflammatory cells (megakaryocytes and monocytes) and a decreasing trend in lymphocytes (eg, CD4+ T, CD8+ T, mucosal-associated invariant T, natural killer, and γδ T cells), highlighting the potential presence of lymphopenia phenomenon in KD. Our data indicated the presence of inflammatory cytokine storm in patients with KD with CALs, caused by systemic upregulation of TNFSF13B (tumor necrosis factor superfamily member 13b), CXCL16 (C-X-C motif chemokine ligand 16), TNFSF10 (tumor necrosis factor superfamily member 10), and IL1RN (interleukin 1 receptor antagonist), mainly produced by monocytes (especially for the Mono_CD14-CD16 cluster) and megakaryocytes. We also found that myeloid cells of patients with KD, particularly in those with CALs, might play a role in vascular injury (eg, increased MMP [matrix metalloproteinase] 9, MMP17, and MMP25) and immune cell recruitment. The immune landscape of patients with KD with CALs was featured by lower exhaustion levels in natural killer cells, a high cytotoxic state in the CD8_Pro cluster, and activation of the complement system in monocytes. Additionally, the activation of B cells was more pronounced in the early stage of KD. CONCLUSIONS: Collectively, this study provides a comprehensive understanding of the roles of various immune cells and inflammatory cytokine storms in the development of CALs in KD and offers a valuable resource for identifying novel therapeutic targets for patients with KD with CALs.

Expression analysis of genes including Zfhx4 in mice and zebrafish reveals a temporospatial conserved molecular basis underlying craniofacial development.

BACKGROUND: Embryonic craniofacial development involves several cellular and molecular events that are evolutionarily conserved among vertebrates. Vertebrate models such as mice and zebrafish have been used to investigate the molecular and cellular etiologies underlying human craniofacial disorders, including orofacial clefts. However, the molecular mechanisms underlying embryonic development in these two species are unknown. Therefore, elucidating the shared mechanisms of craniofacial development between disease models is crucial to understanding the underlying mechanisms of phenotypes in individual species. RESULTS: We selected mice and zebrafish as model organisms to compare various events during embryonic craniofacial development. We identified genes (Sox9, Zfhx3 and 4, Cjun, and Six1) exhibiting similar temporal expression patterns between these species through comprehensive and stage-matched gene expression analyses. Expression analysis revealed similar gene expression in hypothetically corresponding tissues, such as the mice palate and zebrafish ethmoid plate. Furthermore, loss-of-function analysis of Zfhx4/zfhx4, a causative gene of human craniofacial anomalies including orofacial cleft, in both species resulted in deformed skeletal elements such as the palatine and ethmoid plate in mice and zebrafish, respectively. CONCLUSIONS: These results demonstrate that these disease models share common molecular mechanisms, highlighting their usefulness in modeling craniofacial defects in humans.

Unraveling the Immune Signature of Herpes Zoster: Insights Into the Pathophysiology and Human Leukocyte Antigen Risk Profile.

The varicella-zoster virus (VZV) infects >95% of the population. VZV reactivation causes herpes zoster (HZ), known as shingles, primarily affecting the elderly and individuals who are immunocompromised. However, HZ can occur in otherwise healthy individuals. We analyzed the immune signature and risk profile in patients with HZ using a genome-wide association study across different UK Biobank HZ cohorts. Additionally, we conducted one of the largest HZ human leukocyte antigen association studies to date, coupled with transcriptomic analysis of pathways underlying HZ susceptibility. Our findings highlight the significance of the major histocompatibility complex locus for HZ development, identifying 5 protective and 4 risk human leukocyte antigen alleles. This demonstrates that HZ susceptibility is largely governed by variations in the major histocompatibility complex. Furthermore, functional analyses revealed the upregulation of type I interferon and adaptive immune responses. These findings provide fresh molecular insights into the pathophysiology and activation of innate and adaptive immune responses triggered by symptomatic VZV reactivation.

DNA-protein quasi-mapping for rapid differential gene expression analysis in non-model organisms.

BACKGROUND: Conventional differential gene expression analysis pipelines for non-model organisms require computationally expensive transcriptome assembly. We recently proposed an alternative strategy of directly aligning RNA-seq reads to a protein database, and demonstrated drastic improvements in speed, memory usage, and accuracy in identifying differentially expressed genes. RESULT: Here we report a further speed-up by replacing DNA-protein alignment by quasi-mapping, making our pipeline > 1000× faster than assembly-based approach, and still more accurate. We also compare quasi-mapping to other mapping techniques, and show that it is faster but at the cost of sensitivity. CONCLUSION: We provide a quick-and-dirty differential gene expression analysis pipeline for non-model organisms without a reference transcriptome, which directly quasi-maps RNA-seq reads to a reference protein database, avoiding computationally expensive transcriptome assembly.

Targeting Epsins to Inhibit Fibroblast Growth Factor Signaling While Potentiating Transforming Growth Factor-ß Signaling Constrains Endothelial-to-Mesenchymal Transition in Atherosclerosis.

BACKGROUND: Epsin endocytic adaptor proteins are implicated in the progression of atherosclerosis; however, the underlying molecular mechanisms have not yet been fully defined. In this study, we determined how epsins enhance endothelial-to-mesenchymal transition (EndoMT) in atherosclerosis and assessed the efficacy of a therapeutic peptide in a preclinical model of this disease. METHODS: Using single-cell RNA sequencing combined with molecular, cellular, and biochemical analyses, we investigated the role of epsins in stimulating EndoMT using knockout in Apoe-/- and lineage tracing/proprotein convertase subtilisin/kexin type 9 serine protease mutant viral-induced atherosclerotic mouse models. The therapeutic efficacy of a synthetic peptide targeting atherosclerotic plaques was then assessed in Apoe-/- mice. RESULTS: Single-cell RNA sequencing and lineage tracing revealed that epsins 1 and 2 promote EndoMT and that the loss of endothelial epsins inhibits EndoMT marker expression and transforming growth factor-ß signaling in vitro and in atherosclerotic mice, which is associated with smaller lesions in the Apoe-/- mouse model. Mechanistically, the loss of endothelial cell epsins results in increased fibroblast growth factor receptor-1 expression, which inhibits transforming growth factor-ß signaling and EndoMT. Epsins directly bind ubiquitinated fibroblast growth factor receptor-1 through their ubiquitin-interacting motif, which results in endocytosis and degradation of this receptor complex. Consequently, administration of a synthetic ubiquitin-interacting motif-containing peptide atheroma ubiquitin-interacting motif peptide inhibitor significantly attenuates EndoMT and progression of atherosclerosis. CONCLUSIONS: We conclude that epsins potentiate EndoMT during atherogenesis by increasing transforming growth factor-ß signaling through fibroblast growth factor receptor-1 internalization and degradation. Inhibition of EndoMT by reducing epsin-fibroblast growth factor receptor-1 interaction with a therapeutic peptide may represent a novel treatment strategy for atherosclerosis.

Genome-wide identification and expression characterization of the GH3 gene family of tea plant (Camellia sinensis).

To comprehensively understand the characteristics of the GH3 gene family in tea plants (Camellia sinensis), we identified 17 CsGH3 genes and analyzed their physicochemical properties, phylogenetic relationships, gene structures, promoters, and expression patterns in different tissues. The study showed that the 17 CsGH3 genes are distributed on 9 chromosomes, and based on evolutionary analysis, the CsGH3 members were divided into three subgroups. Gene duplication analysis revealed that segmental duplications have a significant impact on the amplification of CsGH3 genes. In addition, we identified and classified cis-elements in the CsGH3 gene promoters and detected elements related to plant hormone responses and non-biotic stress responses. Through expression pattern analysis, we observed tissue-specific expression of CsGH3.3 and CsGH3.10 in flower buds and roots. Moreover, based on predictive analysis of upstream regulatory transcription factors of CsGH3, we identified the potential transcriptional regulatory role of gibberellin response factor CsDELLA in CsGH3.14 and CsGH3.15. In this study, we found that CsGH3 genes are involved in a wide range of activities, such as growth and development, stress response, and transcription. This is the first report on CsGH3 genes and their potential roles in tea plants. In conclusion, these results provide a theoretical basis for elucidating the role of GH3 genes in the development of perennial woody plants and offer new insights into the synergistic effects of multiple hormones on plant growth and development in tea plants.

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.

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.

De novo transcriptome assembly and global analysis of differential gene expression of aphid tolerant wild mustard Rorippa indica (L.) Hiern infested by mustard aphid Lipaphis Erysimi (L.) Kaltenbach.

Rapeseed-mustard, the oleiferous Brassica species are important oilseed crops cultivated all over the globe. Mustard aphid Lipaphis erysimi (L.) Kaltenbach is a major threat to the cultivation of rapeseed-mustard. Wild mustard Rorippa indica (L.) Hiern shows tolerance to mustard aphids as a nonhost and hence is an important source for the bioprospecting of potential resistance genes and defense measures to manage mustard aphids sustainably. We performed mRNA sequencing of the R. indica plant uninfested and infested by the mustard aphids, harvested at 24 hours post-infestation. Following quality control, the high-quality reads were subjected to de novo assembly of the transcriptome. As there is no genomic information available for this potential wild plant, the raw reads will be useful for further bioinformatics analysis and the sequence information of the assembled transcripts will be helpful to design primers for the characterization of specific gene sequences. In this study, we also used the generated resource to comprehensively analyse the global profile of differential gene expression in R. indica in response to infestation by mustard aphids. The functional enrichment analysis of the differentially expressed genes reveals a significant immune response and suggests the possibility of chitin-induced defense signaling.

Gene expression analysis of drought tolerance and cuticular wax biosynthesis in diploid and tetraploid induced wallflowers.

Whole-genome doubling leads to cell reprogramming, upregulation of stress genes, and establishment of new pathways of drought stress responses in plants. This study investigated the molecular mechanisms of drought tolerance and cuticular wax characteristics in diploid and tetraploid-induced Erysimum cheiri. According to real-time PCR analysis, tetraploid induced wallflowers exhibited increased expression of several genes encoding transcription factors (TFs), including AREB1 and AREB3; the stress response genes RD29A and ERD1 under drought stress conditions. Furthermore, two cuticular wax biosynthetic pathway genes, CER1 and SHN1, were upregulated in tetraploid plants under drought conditions. Leaf morphological studies revealed that tetraploid leaves were covered with unique cuticular wax crystalloids, which produced a white fluffy appearance, while the diploid leaves were green and smooth. The greater content of epicuticular wax in tetraploid leaves than in diploid leaves can explain the decrease in cuticle permeability as well as the decrease in water loss and improvement in drought tolerance in wallflowers. GC‒MS analysis revealed that the wax components included alkanes, alcohols, aldehydes, and fatty acids. The most abundant wax compound in this plant was alkanes (50%), the most predominant of which was C29. The relative abundance of these compounds increased significantly in tetraploid plants under drought stress conditions. These findings revealed that tetraploid-induced wallflowers presented upregulation of multiple drought-related and wax biosynthesis genes; therefore, polyploidization has proved useful for improving plant drought tolerance.