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.