Comparison of Single Cell Transcriptome Sequencing Methods: Of Mice and Men.

Single cell RNAseq has been a big leap in many areas of biology. Rather than investigating gene expression on a whole organism level, this technology enables scientists to get a detailed look at rare single cells or within their cell population of interest. The field is growing, and many new methods appear each year. We compared methods utilized in our core facility: Smart-seq3, PlexWell, FLASH-seq, VASA-seq, SORT-seq, 10X, Evercode, and HIVE. We characterized the equipment requirements for each method. We evaluated the performances of these methods based on detected features, transcriptome diversity, mitochondrial RNA abundance and multiplets, among others and benchmarked them against bulk RNA sequencing. Here, we show that bulk transcriptome detects more unique transcripts than any single cell method. While most methods are comparable in many regards, FLASH-seq and VASA-seq yielded the best metrics, e.g., in number of features. If no equipment for automation is available or many cells are desired, then HIVE or 10X yield good results. In general, more recently developed methods perform better. This also leads to the conclusion that older methods should be phased out, and that the development of single cell RNAseq methods is still progressing considerably.

CABA-V7: a prospective biomarker selected trial of cabazitaxel treatment in AR-V7 positive prostate cancer patients.

BACKGROUND: Metastatic castration-resistant prostate cancer (mCRPC) patients with positive AR-V7 expression in their circulating tumour cells (CTCs) rarely derive benefit from abiraterone and enzalutamide. DESIGN: We performed a prospective, multicenter, single arm phase II clinical trial (CABA-V7) in mCRPC patients previously treated with docetaxel and androgen deprivation therapy. OBJECTIVE: In this trial, we investigated whether cabazitaxel treatment resulted in clinically meaningful PSA response rates in patients with positive CTC-based AR-V7 expression and collected liquid biopsies for genomic profiling. RESULTS: Cabazitaxel was found to be modestly effective, with only 12% of these patients obtaining a PSA response. Genomic profiling revealed that CTC-based AR-V7 expression was not associated with other known mCRPC-associated alterations. CTC-based AR-V7 status and dichotomised CTC counts were observed as independent prognostic markers at baseline. CONCLUSIONS: AR-V7 positivity predicted poor overall survival (OS). However, cabazitaxel-treated AR-V7 positive patients and those lacking AR-V7 positivity, who received cabazitaxel as standard of care, appeared to have similar OS. Therefore, despite the low response rate, cabazitaxel may still be an effective treatment in this poor prognosis, AR-V7 positive patient population.

CFTR Correctors and Antioxidants Partially Normalize Lipid Imbalance but not Abnormal Basal Inflammatory Cytokine Profile in CF Bronchial Epithelial Cells.

A deficiency in cystic fibrosis transmembrane conductance regulator (CFTR) function in CF leads to chronic lung disease. CF is associated with abnormalities in fatty acids, ceramides, and cholesterol, their relationship with CF lung pathology is not completely understood. Therefore, we examined the impact of CFTR deficiency on lipid metabolism and pro-inflammatory signaling in airway epithelium using mass spectrometric, protein array. We observed a striking imbalance in fatty acid and ceramide metabolism, associated with chronic oxidative stress under basal conditions in CF mouse lung and well-differentiated bronchial epithelial cell cultures of CFTR knock out pig and CF patients. Cell-autonomous features of all three CF models included high ratios of ω-6- to ω-3-polyunsaturated fatty acids and of long- to very long-chain ceramide species (LCC/VLCC), reduced levels of total ceramides and ceramide precursors. In addition to the retinoic acid analog fenretinide, the anti-oxidants glutathione (GSH) and deferoxamine partially corrected the lipid profile indicating that oxidative stress may promote the lipid abnormalities. CFTR-targeted modulators reduced the lipid imbalance and oxidative stress, confirming the CFTR dependence of lipid ratios. However, despite functional correction of CF cells up to 60% of non-CF in Ussing chamber experiments, a 72-h triple compound treatment (elexacaftor/tezacaftor/ivacaftor surrogate) did not completely normalize lipid imbalance or oxidative stress. Protein array analysis revealed differential expression and shedding of cytokines and growth factors from CF epithelial cells compared to non-CF cells, consistent with sterile inflammation and tissue remodeling under basal conditions, including enhanced secretion of the neutrophil activator CXCL5, and the T-cell activator CCL17. However, treatment with antioxidants or CFTR modulators that mimic the approved combination therapies, ivacaftor/lumacaftor and ivacaftor/tezacaftor/elexacaftor, did not effectively suppress the inflammatory phenotype. We propose that CFTR deficiency causes oxidative stress in CF airway epithelium, affecting multiple bioactive lipid metabolic pathways, which likely play a role in CF lung disease progression. A combination of anti-oxidant, anti-inflammatory and CFTR targeted therapeutics may be required for full correction of the CF phenotype.

Redundant and specific roles of cohesin STAG subunits in chromatin looping and transcriptional control.

Cohesin is a ring-shaped multiprotein complex that is crucial for 3D genome organization and transcriptional regulation during differentiation and development. It also confers sister chromatid cohesion and facilitates DNA damage repair. Besides its core subunits SMC3, SMC1A, and RAD21, cohesin in somatic cells contains one of two orthologous STAG subunits, STAG1 or STAG2. How these variable subunits affect the function of the cohesin complex is still unclear. STAG1- and STAG2-cohesin were initially proposed to organize cohesion at telomeres and centromeres, respectively. Here, we uncover redundant and specific roles of STAG1 and STAG2 in gene regulation and chromatin looping using HCT116 cells with an auxin-inducible degron (AID) tag fused to either STAG1 or STAG2. Following rapid depletion of either subunit, we perform high-resolution Hi-C, gene expression, and sequential ChIP studies to show that STAG1 and STAG2 do not co-occupy individual binding sites and have distinct ways by which they affect looping and gene expression. These findings are further supported by single-molecule localizations via direct stochastic optical reconstruction microscopy (dSTORM) super-resolution imaging. Since somatic and congenital mutations of the STAG subunits are associated with cancer (STAG2) and intellectual disability syndromes with congenital abnormalities (STAG1 and STAG2), we verified STAG1-/STAG2-dependencies using human neural stem cells, hence highlighting their importance in particular disease contexts.

DOC1-Dependent Recruitment of NURD Reveals Antagonism with SWI/SNF during Epithelial-Mesenchymal Transition in Oral Cancer Cells.

The Nucleosome Remodeling and Deacetylase (NURD) complex is a key regulator of cell differentiation that has also been implicated in tumorigenesis. Loss of the NURD subunit Deleted in Oral Cancer 1 (DOC1) is associated with human oral squamous cell carcinomas (OSCCs). Here, we show that restoration of DOC1 expression in OSCC cells leads to a reversal of epithelial-mesenchymal transition (EMT). This is caused by the DOC1-dependent targeting of NURD to repress key transcriptional regulators of EMT. NURD recruitment drives extensive epigenetic reprogramming, including eviction of the SWI/SNF remodeler, formation of inaccessible chromatin, H3K27 deacetylation, and binding of PRC2 and KDM1A, followed by H3K27 methylation and H3K4 demethylation. Strikingly, depletion of SWI/SNF mimics the effects of DOC1 re-expression. Our results suggest that SWI/SNF and NURD function antagonistically to control chromatin state and transcription. We propose that disturbance of this dynamic equilibrium may lead to defects in gene expression that promote oncogenesis.

Overexpression of LMO2 causes aberrant human T-Cell development in vivo by three potentially distinct cellular mechanisms.

Overexpression of LMO2 is known to be one of the causes of T-cell acute lymphoblastic leukemia (T-ALL) development; however, the mechanisms behind its oncogenic activity are incompletely understood. LMO2-overexpressing transgenic mouse models suggest an accumulation of immature T-cell progenitors in the thymus as the main preleukemic event. The effects of LMO2 overexpression on human T-cell development in vivo are unknown. Here, we report studies of a humanized mouse model transplanted with LMO2-transduced human hematopoietic stem/progenitor cells. The effects of LMO2 overexpression were confined to the T-cell lineage; however, initially, multipotent cells were transduced. Three effects of LMO2 on human T-cell development were observed: (1) a block at the double-negative/immature single-positive stage, (2) an accumulation of CD4(+)CD8(+) double-positive CD3(-) cells, and (3) an altered CD8/CD4 ratio with enhanced peripheral T lymphocytes. Microarray analysis of sorted double-positive cells overexpressing LMO2 led to the identification of an LMO2 gene set that clustered with human T-ALL patient samples of the described "proliferative" cluster. In this article, we demonstrate previously unrecognized mechanisms by which LMO2 alters human T-cell development in vivo; these mechanisms correlate with human T-ALL leukemogenesis.

Genetic variants in RBFOX3 are associated with sleep latency.

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10(-08), 6.59 × 10(-)(08) and 9.17 × 10(-)(08)). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10(-)(02), 7.0 × 10(-)(03) and 2.5 × 10(-)(03); combined meta-analysis P-values=5.5 × 10(-07), 5.4 × 10(-07) and 1.0 × 10(-07)). A functional prediction of RBFOX3 based on co-expression with other genes shows that this gene is predominantly expressed in brain (P-value=1.4 × 10(-316)) and the central nervous system (P-value=7.5 × 10(-)(321)). The predicted function of RBFOX3 based on co-expression analysis with other genes shows that this gene is significantly involved in the release cycle of neurotransmitters including gamma-aminobutyric acid and various monoamines (P-values<2.9 × 10(-11)) that are crucial in triggering the onset of sleep. To conclude, in this first large-scale GWAS of sleep latency we report a novel association of variants in RBFOX3 gene. Further, a functional prediction of RBFOX3 supports the involvement of RBFOX3 with sleep latency.

The dystrophin gene and cognitive function in the general population.

The aim of our study is to investigate whether single-nucleotide dystrophin gene (DMD) variants associate with variability in cognitive functions in healthy populations. The study included 1240 participants from the Erasmus Rucphen family (ERF) study and 1464 individuals from the Rotterdam Study (RS). The participants whose exomes were sequenced and who were assessed for various cognitive traits were included in the analysis. To determine the association between DMD variants and cognitive ability, linear (mixed) modeling with adjustment for age, sex and education was used. Moreover, Sequence Kernel Association Test (SKAT) was used to test the overall association of the rare genetic variants present in the DMD with cognitive traits. Although no DMD variant surpassed the prespecified significance threshold (P<1 × 10(-4)), rs147546024:A>G showed strong association (ß = 1.786, P-value = 2.56 × 10(-4)) with block-design test in the ERF study, while another variant rs1800273:G>A showed suggestive association (ß = -0.465, P-value = 0.002) with Mini-Mental State Examination test in the RS. Both variants are highly conserved, although rs147546024:A>G is an intronic variant, whereas rs1800273:G>A is a missense variant in the DMD which has a predicted damaging effect on the protein. Further gene-based analysis of DMD revealed suggestive association (P-values = 0.087 and 0.074) with general cognitive ability in both cohorts. In conclusion, both single variant and gene-based analyses suggest the existence of variants in the DMD which may affect cognitive functioning in the general populations.

Gene expression analysis reveals inhibition of radiation-induced TGFß-signaling by hyperbaric oxygen therapy in mouse salivary glands.

A side effect of radiation therapy in the head and neck region is injury to surrounding healthy tissues such as irreversible impaired function of the salivary glands. Hyperbaric oxygen therapy (HBOT) is clinically used to treat radiation-induced damage but its mechanism of action is largely unknown. In this study, we investigated the molecular pathways that are affected by HBOT in mouse salivary glands two weeks after radiation therapy by microarray analysis. Interestingly, HBOT led to significant attenuation of the radiation-induced expression of a set of genes and upstream regulators that are involved in processes such as fibrosis and tissue regeneration. Our data suggest that the TGFß-pathway, which is involved in radiation-induced fibrosis and chronic loss of function after radiation therapy, is affected by HBOT. On the longer term, HBOT reduced the expression of the fibrosis-associated factor α-smooth muscle actin in irradiated salivary glands. This study highlights the potential of HBOT to inhibit the TGFß-pathway in irradiated salivary glands and to restrain consequential radiation induced tissue injury.