ABSTRACT
In categorical data visualization, appropriate color arrangements can avoid perceptual ambiguity and help perceive underlying data patterns. We introduce a protocol to assign contrastive colors to neighboring categories using both Python and R packages. We describe steps for calculating the interlacement between clusters and generating a proper color palette and calculating color contrast. We then detail procedures for aligning cluster interlacement and color contrast to get an optimized cluster-color assignment, achieving clear categorical visualization. For complete details on the use and execution of this protocol, please refer to Jing et al.1.
Subject(s)
Software , Data Visualization , Color , Cluster AnalysisABSTRACT
Understanding tissue architecture and niche-specific microenvironments in spatially resolved transcriptomics (SRT) requires in situ annotation and labeling of cells. Effective spatial visualization of these data demands appropriate colorization of numerous cell types. However, current colorization frameworks often inadequately account for the spatial relationships between cell types. This results in perceptual ambiguity in neighboring cells of biological distinct types, particularly in complex environments such as brain or tumor. To address this, we introduce Spaco, a potent tool for spatially aware colorization. Spaco utilizes the Degree of Interlacement metric to construct a weighted graph that evaluates the spatial relationships among different cell types, refining color assignments. Furthermore, Spaco incorporates an adaptive palette selection approach to amplify chromatic distinctions. When benchmarked on four diverse datasets, Spaco outperforms existing solutions, capturing complex spatial relationships and boosting visual clarity. Spaco ensures broad accessibility by accommodating color vision deficiency and offering open-accessible code in both Python and R.
ABSTRACT
Onco-fetal reprogramming of the tumor ecosystem induces fetal developmental signatures in the tumor microenvironment, leading to immunosuppressive features. Here, we employed single-cell RNA sequencing, spatial transcriptomics and bulk RNA sequencing to delineate specific cell subsets involved in hepatocellular carcinoma (HCC) relapse and response to immunotherapy. We identified POSTN+ extracellular matrix cancer-associated fibroblasts (EM CAFs) as a prominent onco-fetal interacting hub, promoting tumor progression. Cell-cell communication and spatial transcriptomics analysis revealed crosstalk and co-localization of onco-fetal cells, including POSTN+ CAFs, FOLR2+ macrophages and PLVAP+ endothelial cells. Further analyses suggest an association between onco-fetal reprogramming and epithelial-mesenchymal transition (EMT), tumor cell proliferation and recruitment of Treg cells, ultimately influencing early relapse and response to immunotherapy. In summary, our study identifies POSTN+ CAFs as part of the HCC onco-fetal niche and highlights its potential influence in EMT, relapse and immunotherapy response, paving the way for the use of onco-fetal signatures for therapeutic stratification.
Subject(s)
Carcinoma, Hepatocellular , Folate Receptor 2 , Liver Neoplasms , Humans , Carcinoma, Hepatocellular/therapy , Carcinoma, Hepatocellular/genetics , Liver Neoplasms/therapy , Liver Neoplasms/genetics , Ecosystem , Endothelial Cells , Cell Movement/genetics , Chronic Disease , Recurrence , Immunotherapy , Tumor Microenvironment/geneticsABSTRACT
Dissecting and understanding the cancer ecosystem, especially that around the tumor margins, which have strong implications for tumor cell infiltration and invasion, are essential for exploring the mechanisms of tumor metastasis and developing effective new treatments. Using a novel tumor border scanning and digitization model enabled by nanoscale resolution-SpaTial Enhanced REsolution Omics-sequencing (Stereo-seq), we identified a 500 µm-wide zone centered around the tumor border in patients with liver cancer, referred to as "the invasive zone". We detected strong immunosuppression, metabolic reprogramming, and severely damaged hepatocytes in this zone. We also identified a subpopulation of damaged hepatocytes with increased expression of serum amyloid A1 and A2 (referred to collectively as SAAs) located close to the border on the paratumor side. Overexpression of CXCL6 in adjacent malignant cells could induce activation of the JAK-STAT3 pathway in nearby hepatocytes, which subsequently caused SAAs' overexpression in these hepatocytes. Furthermore, overexpression and secretion of SAAs by hepatocytes in the invasive zone could lead to the recruitment of macrophages and M2 polarization, further promoting local immunosuppression, potentially resulting in tumor progression. Clinical association analysis in additional five independent cohorts of patients with primary and secondary liver cancer (n = 423) showed that patients with overexpression of SAAs in the invasive zone had a worse prognosis. Further in vivo experiments using mouse liver tumor models in situ confirmed that the knockdown of genes encoding SAAs in hepatocytes decreased macrophage accumulation around the tumor border and delayed tumor growth. The identification and characterization of a novel invasive zone in human cancer patients not only add an important layer of understanding regarding the mechanisms of tumor invasion and metastasis, but may also pave the way for developing novel therapeutic strategies for advanced liver cancer and other solid tumors.
Subject(s)
Ecosystem , Liver Neoplasms , Mice , Animals , Humans , Liver Neoplasms/pathology , Hepatocytes/metabolism , Immunosuppression Therapy , Cell Line, TumorABSTRACT
BACKGROUND: Long noncoding RNAs (lncRNAs) are an important class of pervasive genes involved in a variety of biological functions. They are aberrantly expressed in many types of cancers. In this study, we described lncRNAs profiles in 6 pairs of human renal clear cell carcinoma (RCCC) and the corresponding adjacent nontumorous tissues (NT) by microarray. METHODOLOGY/PRINCIPAL FINDINGS: With abundant and varied probes accounting 33,045 LncRNAs in our microarray, the number of lncRNAs that expressed at a certain level could be detected is 17157. From the data we found there were thousands of lncRNAs that differentially expressed (≥ 2 fold-change) in RCCC tissues compared with NT and 916 lncRNAs differentially expressed in five or more of six RCCC samples. Compared with NT, many lncRNAs were significantly up-regulated or down-regulated in RCCC. Our data showed that down-regulated lncRNAs were more common than up-regulated ones. ENST00000456816, X91348, BC029135, NR_024418 were evaluated by qPCR in sixty-three pairs of RCCC and NT samples. The four lncRNAs were aberrantly expressed in RCCC compared with matched histologically normal renal tissues. CONCLUSIONS/SIGNIFICANCE: Our study is the first one to determine genome-wide lncRNAs expression patterns in RCCC by microarray. The results displayed that clusters of lncRNAs were aberrantly expressed in RCCC compared with NT samples, which revealed that lncRNAs differentially expressed in tumor tissues and normal tissues may exert a partial or key role in tumor development. Taken together, this study may provide potential targets for future treatment of RCCC and novel insights into cancer biology.
Subject(s)
Carcinoma, Renal Cell/genetics , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Kidney Neoplasms/genetics , Oligonucleotide Array Sequence Analysis/methods , RNA, Long Noncoding/genetics , Down-Regulation/genetics , Gene Regulatory Networks/genetics , Humans , RNA, Long Noncoding/metabolism , RNA, Messenger/genetics , RNA, Messenger/metabolism , Real-Time Polymerase Chain Reaction , Reproducibility of Results , Up-Regulation/geneticsABSTRACT
The nonhuman primates most commonly used in medical research are from the genus Macaca. To better understand the genetic differences between these animal models, we present high-quality draft genome sequences from two macaque species, the cynomolgus/crab-eating macaque and the Chinese rhesus macaque. Comparison with the previously sequenced Indian rhesus macaque reveals that all three macaques maintain abundant genetic heterogeneity, including millions of single-nucleotide substitutions and many insertions, deletions and gross chromosomal rearrangements. By assessing genetic regions with reduced variability, we identify genes in each macaque species that may have experienced positive selection. Genetic divergence patterns suggest that the cynomolgus macaque genome has been shaped by introgression after hybridization with the Chinese rhesus macaque. Macaque genes display a high degree of sequence similarity with human disease gene orthologs and drug targets. However, we identify several putatively dysfunctional genetic differences between the three macaque species, which may explain functional differences between them previously observed in clinical studies.
Subject(s)
Genome , Macaca fascicularis/genetics , Macaca mulatta/genetics , Animals , Base Sequence , Chromosome Aberrations , Evolution, Molecular , Humans , Models, Animal , Molecular Sequence Data , Polymorphism, Genetic , Sequence Analysis, DNA , Species SpecificityABSTRACT
We report the annotation and analysis of the draft genome sequence of Brassica rapa accession Chiifu-401-42, a Chinese cabbage. We modeled 41,174 protein coding genes in the B. rapa genome, which has undergone genome triplication. We used Arabidopsis thaliana as an outgroup for investigating the consequences of genome triplication, such as structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one of the three copies consistently retaining a disproportionately large fraction of the genes expected to have been present in its ancestor. Variation in the number of members of gene families present in the genome may contribute to the remarkable morphological plasticity of Brassica species. The B. rapa genome sequence provides an important resource for studying the evolution of polyploid genomes and underpins the genetic improvement of Brassica oil and vegetable crops.
Subject(s)
Brassica rapa/genetics , Genome, Plant , Polyploidy , Arabidopsis/genetics , Chromosomes, Artificial, Bacterial/genetics , Chromosomes, Plant/genetics , Contig Mapping , Evolution, Molecular , Gene Duplication , Genes, Plant , Molecular Sequence Annotation , Molecular Sequence Data , Phylogeny , Sequence Analysis, DNAABSTRACT
Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.
Subject(s)
Genome/genetics , Genomics , Ursidae/genetics , Algorithms , Animals , China , Conserved Sequence/genetics , Contig Mapping , Diet/veterinary , Dogs , Evolution, Molecular , Female , Fertility/genetics , Fertility/physiology , Heterozygote , Humans , Multigene Family/genetics , Polymorphism, Single Nucleotide/genetics , Receptors, G-Protein-Coupled/genetics , Sequence Alignment , Sequence Analysis, DNA , Synteny/genetics , Ursidae/classification , Ursidae/physiologyABSTRACT
Cucumber is an economically important crop as well as a model system for sex determination studies and plant vascular biology. Here we report the draft genome sequence of Cucumis sativus var. sativus L., assembled using a novel combination of traditional Sanger and next-generation Illumina GA sequencing technologies to obtain 72.2-fold genome coverage. The absence of recent whole-genome duplication, along with the presence of few tandem duplications, explains the small number of genes in the cucumber. Our study establishes that five of the cucumber's seven chromosomes arose from fusions of ten ancestral chromosomes after divergence from Cucumis melo. The sequenced cucumber genome affords insight into traits such as its sex expression, disease resistance, biosynthesis of cucurbitacin and 'fresh green' odor. We also identify 686 gene clusters related to phloem function. The cucumber genome provides a valuable resource for developing elite cultivars and for studying the evolution and function of the plant vascular system.