13C metabolic flux analysis clarifies distinct metabolic phenotypes of cancer cell spheroid mimicking tumor hypoxia

Cancer cells adapt their intracellular energy metabolism to the oxygen-deprived tumor microenvironment (TME) to ensure tumor progression. This adaptive mechanism has focused attention on the metabolic phenotypes of tumor cells under hypoxic TME for developing novel cancer therapies. Although widely used monolayer (2D) culture does not fully reflect in vivo hypoxic TME, spheroid (3D) culture can produce a milieu similar to the TME in vivo. However, how different metabolic phenotypes are expressed in 3D cultures mimicking tumor hypoxia compared with 2D cultures under hypoxia remains unclear. To address this issue, we investigated the metabolic phenotypes of 2D- and 3D-cultured cancer cells by ¹³C-metabolic flux analysis (¹³C-MFA). Principal component analysis of ¹³C mass isotopomer distributions clearly demonstrated distinct metabolic phenotypes of 3D-cultured cells. ¹³C-MFA clarified that 3D culture significantly upregulated pyruvate carboxylase flux in line with the pyruvate carboxylase protein expression level. On the other hand, 3D culture downregulated glutaminolytic flux. Consistent with our findings, 3D-cultured cells are more resistant to a glutaminase inhibitor than 2D-cultured cells. This study suggests the importance of considering the metabolic characteristics of the particular in vitro model used for research on cancer metabolism.     

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.     

Cyclic mechanical stretch contributes to network development of osteocyte-like cells with morphological change and autophagy promotion but without preferential cell alignment in rat

Osteocytes play important roles in controlling bone quality as well as preferential alignment of biological apatite c-axis/collagen fibers. However, the relationship between osteocytes and mechanical stress remains unclear due to the difficulty of three-dimensional (3D) culture of osteocytes in vitro. The aim of this study was to investigate the effect of cyclic mechanical stretch on 3D-cultured osteocyte-like cells. Osteocyte-like cells were established using rat calvarial osteoblasts cultured in a 3D culture system. Cyclic mechanical stretch (8% amplitude at a rate of 2 cycles min^?1) was applied for 24, 48 and 96 consecutive hours. Morphology, cell number and preferential cell alignment were evaluated. Apoptosis- and autophagy-related gene expression levels were measured using quantitative PCR. 3D-cultured osteoblasts became osteocyte-like cells that expressed osteocyte-specific genes such as Dmp1, Cx43, Sost, Fgf23 and RANKL, with morphological changes similar to osteocytes. Cell number was significantly decreased in a time-dependent manner under non-loaded conditions, whereas cyclic mechanical stretch significantly prevented decreased cell numbers with increased expression of anti-apoptosis-related genes. Moreover, cyclic mechanical stretch significantly decreased cell size and ellipticity with increased expression of autophagy-related genes, LC3b and atg7. Interestingly, preferential cell alignment did not occur, irrespective of mechanical stretch. These findings suggest that an anti-apoptotic effect contributes to network development of osteocyte-like cells under loaded condition. Spherical change of osteocyte-like cells induced by mechanical stretch may be associated with autophagy upregulation. Preferential alignment of osteocytes induced by mechanical load in vivo may be partially predetermined before osteoblasts differentiate into osteocytes and embed into bone matrix.     

Relaxed 3D genome conformation facilitates the pluripotent to totipotent-like state transition in embryonic stem cells

The 3D genome organization is crucial for gene regulation. Although recent studies have revealed a uniquely relaxed genome conformation in totipotent early blastomeres of both fertilized and cloned embryos, how weakened higher-order chromatin structure is functionally linked to totipotency acquisition remains elusive. Using low-input Hi-C, ATAC-seq and ChIP-seq, we systematically examined the dynamics of 3D genome and epigenome during pluripotent to totipotent-like state transition in mouse embryonic stem cells (ESCs). The spontaneously converted 2-cell-embryo-like cells (2CLCs) exhibited more relaxed chromatin architecture compared to ESCs, including global weakening of both enhancer-promoter interactions and TAD insulation. While the former correlated with inactivation of ESC enhancers and down-regulation of pluripotent genes, the latter might facilitate contacts between the putative new enhancers arising in 2CLCs and neighboring 2C genes. Importantly, disruption of chromatin organization by depleting CTCF or the cohesin complex promoted the ESC to 2CLC transition. Our results thus establish a critical role of 3D genome organization in totipotency acquisition.     

Tn5-FISH,a novel cytogenetic method to image chromatin interactions with sub-kilobase resolution

There is an increasing interest in understanding how three-dimensional (3D) organization of the genome is regulated. Different strategies have been employed to identify genome-wide chromatin interactions. However, due to current limitations in resolving genomic contacts, visualization and validation of these genomic loci with sub-kilobase resolution remain unsolved to date. Here, we describe Tn5 transposase-based Fluorescencein situhybridization (Tn5-FISH), a PCR-based, cost-effective imaging method, which can co-localize the genomic loci with sub-kilobase resolution, dissect genome architecture, and verify chromatin interactions detected by chromatin configuration capture (3C)-derived methods. To validate this method, short-range interactions in keratin-encoding gene (KRT) locus in topologically associated domain (TAD) were imaged by triple-color Tn5-FISH, indicating that Tn5-FISH is very useful to verify short-range chromatin interactions inside the contact domain and TAD. Therefore, Tn5-FISH can be a powerful molecular tool for the clinical detection of cytogenetic changes in numerous genetic diseases such as cancers.     

Three-dimensional differentiation of adipose-derived mesenchymal stem cells into insulin-producing cells

The aim of this study is to evaluate the collagen/hyaluronic acid (Col/HA) scaffold effect on the differentiation of insulin-producing cells (IPCs) from adipose-derived mesenchymal stem cells (ASCs). In this experimental study, ASCs were cultured and seeded in a Col/HA scaffold (3D culture) and then treated with induction media. After induction, the presence of IPCs was evaluated using gene expression (PDX-1, GLUT-2 and insulin) analysis and immunocytochemistry, while functional maturity was determined by measuring insulin release in response to low- and high-glucose media. The induced IPCs were morphologically similar to pancreatic islet-like cells. Expression of the islet-associated genes PDX-1, GLUT-2 and insulin genes in 3D-cultured cells was markedly higher than the 2D-cultured cells exposure differentiation media. Compared to the 2D culture of ASCs-derived IPCs, the insulin release from 3D ASCs-derived IPCs showed a nearly 4-fold (p?<?0.05) increase when exposed to a high glucose (25 mmol) medium. The percentage of insulin-positive cells in the 3D experimental group showed an approximately 4-fold increase compared to the 2D experimental culture cells. The results of this study demonstrated that the COL/HA scaffold can enhance the differentiation of IPCs from rat ASCs.     

Epithelial cancer cells exhibit different electrical properties when cultured in 2D and 3D environments

Background

Many drug development and toxicology studies are performed using cells grown in monolayers in well-plates and flasks, despite the fact that these are widely held to be different to cells found in the native environment. 3D, tissue engineered, organotypical tissue culture systems have been developed to be more representative of the native tissue environment than standard monolayer cultures. Whilst the biochemical differences between cells grown in 2D and 3D culture have been explored, the changes on the electrophysiological properties of the cells have not.

Methods

We compared the electrophysiological properties of primary normal oral keratinocytes (nOK) and cancerous abnormal oral keratinocytes (aOK), cultured in standard monolayer and reconstituted 3D organotypical tissue cultures. The electrophysiological properties of populations of the cells were analysed using dielectrophoresis. The intracellular conductivity of aOK was significantly increased when grown in organotypical cultures compared to counterpart cells grown in monolayer cultures.

Results

3D cultured aOK showed almost identical intracellular conductivity to nOK also grown in organotypical cultures, but significantly different to aOK grown in monolayers. The effective membrane capacitance of aOK grown in 3D was found to be significantly higher than nOK, but there was no significant difference between the electrophysiological properties of nOK grown in 2D and 3D cultures.

General significance

This work suggests that factors such as cell shape and cytoplasmic trafficking between cells play an important role in their electrophysiology, and highlights the need to use in vitro models more representative of native tissue when studying cell electrophysiological properties.     

SATB1 Mediates Long-Range Chromatin Interactions: A Dual Regulator of Anti-Apoptotic BCL2 and Pro-Apoptotic NOXA Genes

Aberrant expression of special AT-rich binding protein 1 (SATB1), a global genomic organizer, has been associated with various cancers, which raises the question of how higher-order chromatin structure contributes to carcinogenesis. Disruption of apoptosis is one of the hallmarks of cancer. We previously demonstrated that SATB1 mediated specific long-range chromosomal interactions between the mbr enhancer located within 3’-UTR of the BCL2 gene and the promoter to regulate BCL2 expression during early apoptosis. In the present study, we used chromosome conformation capture (3C) assays and molecular analyses to further investigate the function of the SATB1-mediated higher-order chromatin structure in co-regulation of the anti-apoptotic BCL2 gene and the pro-apoptotic NOXA gene located 3.4Mb downstream on Chromosome 18. We demonstrated that the mbr enhancer spatially juxtaposed the promoters of BCL2 and NOXA genes through SATB1-mediated chromatin-loop in Jurkat cells. Decreased SATB1 levels switched the mbr-BCL2 loop to mbr-NOXA loop, and thus changed expression of these two genes. The SATB1-mediated dynamic switch of the chromatin loop structures was essential for the cooperative expression of the BCL2 and NOXA genes in apoptosis. Notably, the role of SATB1 was specific, since inhibition of SATB1 degradation by caspase-6 inhibitor or caspase-6-resistant SATB1 mutant reversed expression of BCL-2 and NOXA in response to apoptotic stimulation. This study reveals the critical role of SATB1-organized higher-order chromatin structure in regulating the dynamic equilibrium of apoptosis-controlling genes with antagonistic functions and suggests that aberrant SATB1 expression might contribute to cancer development by disrupting the co-regulated genes in apoptosis pathways.     

G9a/GLP-sensitivity of H3K9me2 Demarcates Two Types of Genomic Compartments

In the nucleus, chromatin is folded into hierarchical architecture that is tightly linked to various nuclear functions. However, the underlying molecular mechanisms that confer these architectures remain incompletely understood. Here, we investigated the functional roles of H3 lysine 9 dimethylation (H3K9me2), one of the abundant histone modifications, in three-dimensional (3D) genome organization. Unlike in mouse embryonic stem cells, inhibition of methyltransferases G9a and GLP in differentiated cells eliminated H3K9me2 predominantly at A-type (active) genomic compartments, and the level of residual H3K9me2 modifications was strongly associated with B-type (inactive) genomic compartments. Furthermore, chemical inhibition of G9a/GLP in mouse hepatocytes led to decreased chromatin-nuclear lamina interactions mainly at G9a/GLP-sensitive regions, increased degree of genomic compartmentalization, and up-regulation of hundreds of genes that were associated with alterations of the 3D chromatin. Collectively, our data demonstrated essential roles of H3K9me2 in 3D genome organization.     

Involvement of RhoA/ROCK in insulin secretion of pancreatic β-cells in 3D culture

Cell-cell contacts and interactions between pancreatic β-cells and/or other cell populations within islets are essential for cell survival, insulin secretion, and functional synchronization. Three-dimensional (3D) culture systems supply the ideal microenvironment for islet-like cluster formation and functional maintenance. However, the underlying mechanisms remain unclear. In this study, mouse insulinoma 6 (MIN6) cells were cultured in a rotating 3D culture system to form islet-like aggregates. Glucose-stimulated insulin secretion (GSIS) and the RhoA/ROCK pathway were investigated. In the 3D-cultured MIN6 cells, more endocrine-specific genes were up-regulated, and GSIS was increased to a greater extent than in cells grown in monolayers. RhoA/ROCK inactivation led to F-actin remodeling in the MIN6 cell aggregates and greater insulin exocytosis. The gap junction protein, connexin 36 (Cx36), was up-regulated in MIN6 cell aggregates and RhoA/ROCK-inactivated monolayer cells. GSIS dramatically decreased when Cx36 was knocked down by short interfering RNA and could not be reversed by RhoA/ROCK inactivation. Thus, the RhoA/ROCK signaling pathway is involved in insulin release through the up-regulation of Cx36 expression in 3D-cultured MIN6 cells.     

Variation in tissue-specific gene expression among natural populations

Background

Variation in gene expression is extensive among tissues, individuals, strains, populations and species. The interactions among these sources of variation are relevant for physiological studies such as disease or toxic stress; for example, it is common for pathologies such as cancer, heart failure and metabolic disease to be associated with changes in tissue-specific gene expression or changes in metabolic gene expression. But how conserved these differences are among outbred individuals and among populations has not been well documented. To address this we examined the expression of a selected suite of 192 metabolic genes in brain, heart and liver in three populations of the teleost fish Fundulus heteroclitus using a highly replicated experimental design.

Results

Half of the genes (48%) were differentially expressed among individuals within a population-tissue group and 76% were differentially expressed among tissues. Differences among tissues reflected well established tissue-specific metabolic requirements, suggesting that these measures of gene expression accurately reflect changes in proteins and their phenotypic effects. Remarkably, only a small subset (31%) of tissue-specific differences was consistent in all three populations.

Conclusions

These data indicate that many tissue-specific differences in gene expression are unique to one population and thus are unlikely to contribute to fundamental differences between tissue types. We suggest that those subsets of treatment-specific gene expression patterns that are conserved between taxa are most likely to be functionally related to the physiological state in question.     

Integrative Analysis of Genome, 3D Genome,and Transcriptome Alterations of Clinical Lung Cancer Samples

Genomic studies of cancer cell alterations, such as mutations, copy number variations (CNVs), and translocations, greatly promote our understanding of the genesis and development of cancers. However, the 3D genome architecture of cancers remains less studied due to the complexity of cancer genomes and technical difficulties. To explore the 3D genome structure in clinical lung cancer, we performed Hi-C experiments using paired normal and tumor cells harvested from patients with lung cancer, combining with RNA sequenceing analysis. We demonstrated the feasibility of studying 3D genome of clinical lung cancer samples with a small number of cells (1 × 10⁴), compared the genome architecture between clinical samples and cell lines of lung cancer, and identified conserved and changed spatial chromatin structures between normal and cancer samples. We also showed that Hi-C data can be used to infer CNVs and point mutations in cancer. By integrating those different types of cancer alterations, we showed significant associations between CNVs, 3D genome, and gene expression. We propose that 3D genome mediates the effects of cancer genomic alterations on gene expression through altering regulatory chromatin structures. Our study highlights the importance of analyzing 3D genomes of clinical cancer samples in addition to cancer cell lines and provides an integrative genomic analysis pipeline for future larger-scale studies in lung cancer and other cancers.