The cynomolgus and Chinese rhesus macaques are used as animal models in biomedical research. Yan et al. sequence their genomes and compare the sequences to that of the Indian rhesus macaque, providing a genetic foundation for interpreting research results. The nonhuman primates most commonly used in medical research are from the genus Macaca1. 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.

The extracellular matrix (ECM) spatiotemporally controls cell fate; however, dysregulation of ECM remodeling can lead to tumorigenesis and cancer development by providing favorable conditions for tumor cells. Proteoglycans (PGs) and glycosaminoglycans (GAGs) are the major macromolecules composing ECM. They influence both cell behavior and matrix properties through direct and indirect interactions with various cytokines, growth factors, cell surface receptors, adhesion molecules, enzymes, and glycoproteins within the ECM. The classical features of PGs/GAGs play well-known roles in cancer angiogenesis, proliferation, invasion, and metastasis. Several lines of evidence suggest that PGs/GAGs critically affect broader aspects in cancer initiation and the progression process, including regulation of cell metabolism, serving as a sensor of ECM’s mechanical properties, affecting immune supervision, and participating in therapeutic resistance to various forms of treatment. These functions may be implemented through the characteristics of PGs/GAGs as molecular bridges linking ECM and cells in cell-specific and context-specific manners within the tumor microenvironment (TME). In this review, we intend to present a comprehensive illustration of the ways in which PGs/GAGs participate in and regulate several aspects of tumorigenesis; we put forward a perspective regarding their effects as biomarkers or targets for diagnoses and therapeutic interventions.

Altered metabolism is a hallmark of cancer and glycolysis is one of the important factors promoting tumor development. There is however still a lack of molecular characterization glycolysis and comprehensive studies related to tumor glycolysis in the pan-cancer landscape. Here, we applied a gene expression signature to quantify glycolysis in 9229 tumors across 25 cancer types and 7875 human lung cancer single cells and verified the robustness of signature using defined glycolysis samples from previous studies. We classified tumors and cells into glycolysis score-high and -low groups, demonstrated their prognostic associations, and identified genome and transcriptome molecular features associated with glycolysis activity. We observed that glycolysis score-high tumors were associated with worse prognosis across cancer types. High glycolysis tumors exhibited specific driver genes altered by copy number aberrations (CNAs) in most cancer types. Tricarboxylic acid (TCA) cycle, DNA replication, tumor proliferation and other cancer hallmarks were more active in glycolysis-high tumors. Glycolysis signature was strongly correlated with hypoxia signature in all 25 cancer tissues (r > 0.7) and cancer single cells (r > 0.8). In addition, HSPA8 and P4HA1 were screened out as the potential modulating factors to glycolysis as their expression were highly correlated with glycolysis score and glycolysis genes, which enables future efforts for therapeutic options to block the glycolysis and control tumor progression. Our study provides a comprehensive molecular-level understanding of glycolysis with a large sample data and demonstrates the hypoxia pressure, growth signals, oncogene mutation and other potential signals could activate glycolysis, thereby to regulate cell cycle, energy material synthesis, cell proliferation and cancer progression.

Hypoxia is a characteristic of tumor microenvironment (TME) and is a major contributor to tumor progression. Yet, subtype identification of tumor-associated non-malignant cells at single-cell resolution and how they influence cancer progression under hypoxia TME remain largely unexplored. Here, we used RNA-seq data of 424,194 single cells from 108 patients to identify the subtypes of cancer cells, stromal cells, and immune cells; to evaluate their hypoxia score; and also to uncover potential interaction signals between these cells in vivo across six cancer types. We identified SPP1+ tumor-associated macrophage (TAM) subpopulation potentially enhanced epithelial–mesenchymal transition (EMT) by interaction with cancer cells through paracrine pattern. We prioritized SPP1 as a TAM-secreted factor to act on cancer cells and found a significant enhanced migration phenotype and invasion ability in A549 lung cancer cells induced by recombinant protein SPP1. Besides, prognostic analysis indicated that a higher expression of SPP1 was found to be related to worse clinical outcome in six cancer types. SPP1 expression was higher in hypoxia-high macrophages based on single-cell data, which was further validated by an in vitro experiment that SPP1 was upregulated in macrophages under hypoxia-cultured compared with normoxic conditions. Additionally, a differential analysis demonstrated that hypoxia potentially influences extracellular matrix remodeling, glycolysis, and interleukin-10 signal activation in various cancer types. Our work illuminates the clearer underlying mechanism in the intricate interaction between different cell subtypes within hypoxia TME and proposes the guidelines for the development of therapeutic targets specifically for patients with high proportion of SPP1+ TAMs in hypoxic lesions.

The development and innovation of next generation sequencing (NGS) and the subsequent analysis tools have gain popularity in scientific researches and clinical diagnostic applications. Hence, a systematic comparison of the sequencing platforms and variant calling pipelines could provide significant guidance to NGS-based scientific and clinical genomics. In this study, we compared the performance, concordance and operating efficiency of 27 combinations of sequencing platforms and variant calling pipelines, testing three variant calling pipelines—Genome Analysis Tool Kit HaplotypeCaller, Strelka2 and Samtools-Varscan2 for nine data sets for the NA12878 genome sequenced by different platforms including BGISEQ500, MGISEQ2000, HiSeq4000, NovaSeq and HiSeq Xten. For the variants calling performance of 12 combinations in WES datasets, all combinations displayed good performance in calling SNPs, with their F-scores entirely higher than 0.96, and their performance in calling INDELs varies from 0.75 to 0.91. And all 15 combinations in WGS datasets also manifested good performance, with F-scores in calling SNPs were entirely higher than 0.975 and their performance in calling INDELs varies from 0.71 to 0.93. All of these combinations manifested high concordance in variant identification, while the divergence of variants identification in WGS datasets were larger than that in WES datasets. We also down-sampled the original WES and WGS datasets at a series of gradient coverage across multiple platforms, then the variants calling period consumed by the three pipelines at each coverage were counted, respectively. For the GIAB datasets on both BGI and Illumina platforms, Strelka2 manifested its ultra-performance in detecting accuracy and processing efficiency compared with other two pipelines on each sequencing platform, which was recommended in the further promotion and application of next generation sequencing technology. The results of our researches will provide useful and comprehensive guidelines for personal or organizational researchers in reliable and consistent variants identification.