Integrating spatial and single-cell transcriptomics reveals tumor heterogeneity and intercellular networks in colorectal cancer.

Single cell RNA sequencing (scRNA-seq), a powerful tool for studying the tumor microenvironment (TME), does not preserve/provide spatial information on tissue morphology and cellular interactions. To understand the crosstalk between diverse cellular components in proximity in the TME, we performed scRNA-seq coupled with spatial transcriptomic (ST) assay to profile 41,700 cells from three colorectal cancer (CRC) tumor-normal-blood pairs. Standalone scRNA-seq analyses revealed eight major cell populations, including B cells, T cells, Monocytes, NK cells, Epithelial cells, Fibroblasts, Mast cells, Endothelial cells. After the identification of malignant cells from epithelial cells, we observed seven subtypes of malignant cells that reflect heterogeneous status in tumor, including tumor_CAV1, tumor_ATF3_JUN | FOS, tumor_ZEB2, tumor_VIM, tumor_WSB1, tumor_LXN, and tumor_PGM1. By transferring the cellular annotations obtained by scRNA-seq to ST spots, we annotated four regions in a cryosection from CRC patients, including tumor, stroma, immune infiltration, and colon epithelium regions. Furthermore, we observed intensive intercellular interactions between stroma and tumor regions which were extremely proximal in the cryosection. In particular, one pair of ligands and receptors (C5AR1 and RPS19) was inferred to play key roles in the crosstalk of stroma and tumor regions. For the tumor region, a typical feature of TMSB4X-high expression was identified, which could be a potential marker of CRC. The stroma region was found to be characterized by VIM-high expression, suggesting it fostered a stromal niche in the TME. Collectively, single cell and spatial analysis in our study reveal the tumor heterogeneity and molecular interactions in CRC TME, which provides insights into the mechanisms underlying CRC progression and may contribute to the development of anticancer therapies targeting on non-tumor components, such as the extracellular matrix (ECM) in CRC. The typical genes we identified may facilitate to new molecular subtypes of CRC.

singleCellBase: a high-quality manually curated database of cell markers for single cell annotation across multiple species.

Annotating cells in the analysis of single-cell RNA-seq (scRNA-seq) data is one of the most challenging tasks that researchers are actively addressing. Manual cell annotation is generally considered the gold standard method, although it is labor intensive and independent of prior knowledge. At present, the relationship between high-quality, known marker genes and cell types is very limited, especially for a variety of species other than humans and mice. The singleCellBase is a manually curated resource of high-quality cell types and gene markers associations across multiple species. In details, it offers 9,158 entries spanning a total of 1,221 cell types and linking with 8,740 genes (cell markers), covering 464 diseases/status, and 165 types of tissues across 31 species. The singleCellBase provides a user-friendly interface to the scientific community to browse, search, download and submit records of marker genes and cell types. The resource providing ineluctable prior knowledge required by manual cell annotation, which is valuable to interpret scRNA-seq data and elucidate what cell type or cell state that a cell population represents.

Single-cell analysis of peripheral blood from high-altitude pulmonary hypertension patients identifies a distinct monocyte phenotype.

Immune and inflammatory responses have an important function in the pathophysiology of pulmonary hypertension (PH). However, little is known about the immune landscape in peripheral circulation in patients with high-altitude pulmonary hypertension (HAPH). We apply single-cell transcriptomics to characterize the monocytes that are significantly enriched in the peripheral blood mononuclear cells (PBMC) of HAPH patients. We discover an increase in C1 (non-classical) and C2 (intermediate) monocytes in PBMCs and a decrease in hypoxia-inducible transcription factor-1α (HIF-1α) in all monocyte subsets associated with HAPH. In addition, we demonstrate that similar immune adaptations may exist in HAPH and PH. Overall, we characterize an immune cell atlas of the peripheral blood in HAPH patients. Our data provide evidence that specific monocyte subsets and HIF-1α downregulation might be implicated in the pathogenesis of HAPH.

Prognostic Pathways Guide Drug Indications in Pan-Cancers.

Pathway-level analysis is a powerful approach enabling the interpretation of post-genomic data at a higher level than that of individual molecules. Molecular-targeted therapy focusing on cascade signaling pathways has become a new paradigm in anticancer therapy, instead of a single protein. However, the approaches to narrowing down the long list of biological pathways are limited. Here, we proposed a strategy for in silico Drug Prescription on biological pathways across pan-Cancers (CDP), by connecting drugs to candidate pathways. Applying on a list of 120 traditional Chinese medicines (TCM), we especially identified the "TCM-pathways-cancers" triplet and constructed it into a heterogeneous network across pan-cancers. Applying them into TCMs, the computational prescribing methods deepened the understanding of the efficacy of TCM at the molecular level. Further applying them into Western medicines, CDP could promote drug reposition avoiding time-consuming developments of new drugs.

Development of a maize 55 K SNP array with improved genome coverage for molecular breeding.

With the decrease of cost in genotyping, single nucleotide polymorphisms (SNPs) have gained wide acceptance because of their abundance, even distribution throughout the maize (Zea mays L.) genome, and suitability for high-throughput analysis. In this study, a maize 55 K SNP array with improved genome coverage for molecular breeding was developed on an Affymetrix® Axiom® platform with 55,229 SNPs evenly distributed across the genome, including 22,278 exonic and 19,425 intronic SNPs. This array contains 451 markers that are associated with 368 known genes and two traits of agronomic importance (drought tolerance and kernel oil biosynthesis), 4067 markers that are not covered by the current reference genome, 734 markers that are differentiated significantly between heterotic groups, and 132 markers that are tags for important transgenic events. To evaluate the performance of 55 K array, we genotyped 593 inbred lines with diverse genetic backgrounds. Compared with the widely-used Illumina® MaizeSNP50 BeadChip, our 55 K array has lower missing and heterozygous rates and more SNPs with lower minor allele frequency (MAF) in tropical maize, facilitating in-depth dissection of rare but possibly valuable variation in tropical germplasm resources. Population structure and genetic diversity analysis revealed that this 55 K array is also quite efficient in resolving heterotic groups and performing fine fingerprinting of germplasm. Therefore, this maize 55 K SNP array is a potentially powerful tool for germplasm evaluation (including germplasm fingerprinting, genetic diversity analysis, and heterotic grouping), marker-assisted breeding, and primary quantitative trait loci (QTL) mapping and genome-wide association study (GWAS) for both tropical and temperate maize.

Genome-wide association study identifies two new susceptibility loci for colorectal cancer at 5q23.3 and 17q12 in Han Chinese.

Genome-wide association studies (GWAS) have reported a number of loci harboring common variants that influence risk of colorectal cancer (CRC) in European descent. But all the SNPs identified explained a small fraction of total heritability. To identify more genetic factors that modify the risk of CRC, especially Chinese Han specific, we conducted a three-stage GWAS including a screening stage (932 CRC cases and 966 controls) and two independent validations (Stage 2: 1,759 CRC cases and 1,875 controls; Stage 3: 943 CRC cases and 1,838 controls). In the combined analyses, we discovered two novel loci associated with CRC: rs12522693 at 5q23.3 (CDC42SE2-CHSY3, OR = 1.31, P = 2.08 × 10-8) and rs17836917 at 17q12 (ASIC2-CCL2, OR = 0.75, P = 4.55 × 10-8). Additionally, we confirmed two previously reported risk loci, rs6983267 at 8q24.21 (OR = 1.17, P = 7.17 × 10-7) and rs10795668 at 10p14 (OR = 0.86, P = 2.96 × 10-6) in our cohorts. These results bring further insights into the CRC susceptibility and advance our understanding on etiology of CRC.

A fully sealed plastic chip for multiplex PCR and its application in bacteria identification.

Multiplex PCR is an effective tool for simultaneous multiple target detection but is limited by the intrinsic interference and competition among primer pairs when it is performed in one reaction tube. Dividing a multiplex PCR into many single PCRs is a simple strategy to overcome this issue. Here, we constructed a plastic, easy-to-use, fully sealed multiplex PCR chip based on reversible centrifugation for the simultaneous detection of 63 target DNA sequences. The structure of the chip is quite simple, which contains sine-shaped infusing channels and a number of reaction chambers connecting to one side of these channels. Primer pairs for multiplex PCR were sequentially preloaded in the different reaction chambers, and the chip was enclosed with PCR-compatible adhesive tape. For usage, the PCR master mix containing a DNA template is pipetted into the infusing channels and centrifuged into the reaction chambers, leaving the infusing channels filled with air to avoid cross-contamination of the different chambers. Then, the chip is sealed and placed on a flat thermal cycler for PCR. Finally, amplification products can be detected in situ using a fluorescence scanner or recovered by reverse centrifugation for further analyses. Therefore, our chip possesses two functions: 1) it can be used for multi-target detection based on end-point in situ fluorescence detection; and 2) it can work as a sample preparation unit for analyses that need multiplex PCR such as hybridization and target sequencing. The performance of this chip was carefully examined and further illustrated in the identification of 8 pathogenic bacterial genomic DNA samples and 13 drug-resistance genes. Due to simplicity of its structure and operation, accuracy and generality, high-throughput capacity, and versatile functions (i.e., for in situ detection and sample preparation), our multiplex PCR chip has great potential in clinical diagnostics and nucleic acid-based point-of-care testing.

Identification of differentially-expressed genes in intestinal gastric cancer by microarray analysis.

Gastric cancer (GC) is one of the most frequent malignant tumors. In order to systematically characterize the cellular and molecular mechanisms of intestinal GC development, in this study, we used 22K oligonucleotide microarrays and bioinformatics analysis to evaluate the gene expression profiles of GC in 45 tissue samples, including 20 intestinal GC tissue samples, 20 normal appearing tissues (NATs) adjacent to tumors and 5 noncancerous gastric mucosa tissue samples. These profiles allowed us to explore the transcriptional characteristics of GC and determine the change patterns in gene expression that may be of clinical significance. 1519 and 1255 differentially-expressed genes (DEGs) were identified in intestinal GC tissues and NATs, respectively, as determined by Bayesian analysis (P<0.001). These genes were associated with diverse functions such as mucosa secretion, metabolism, proliferation, signaling and development, which occur at different stages of GC development.

Construction of a novel oligonucleotide array-based transcription factor interaction assay platform and its uses for profiling STAT1 cofactors in mouse fibroblast cells.

Here, we describe a novel oligonucleotide array-based transcription factor (TF) interaction assay platform that can directly identify cointeracting TF complexes following binding to their regulatory DNA elements. This platform that combines immuno-coprecipitation technology with our previously reported oligonucleotide array-based TF assay (OATFA), is named targeted immuno-coprecipitation OATFA (TIC-OATFA). We illustrate use of the system to identify interaction partners of STAT1 (signal transducer and activator of transcription proteins 1) in mouse fibroblasts. Several previously known partners of STAT1, as well as new partners, were identified by TIC-OATFA, including the upstream stimulatory factors 1 and 2 (USF1, USF2), nuclear factor of activated T cells, TATA box-binding protein, nuclear factor erythroid-derived 2, nuclear factor-kappa B, and nuclear factor 1. Both USF1 and nuclear factor-kappa B are well known to interact with STAT1, but the other five TFs are previously unreported STAT1 interaction partners. We examined interactions between one new TF, USF2, and STAT1 in detail. USF2 belongs to the group of bHLH-zip transcription factors, which in a number of diseases including cancers, has enhanced activity. In summary, a novel oligonucleotide array-based assay platform was developed and used to study interactions between STAT1 and functional TF binding partners, revealing that USF2 and potentially four other new TFs are partners of STAT1 in an IFN-γ stimulated mouse fibroblast cell line.

Association between paraoxonase gene and stroke in the Han Chinese population.

BACKGROUND: The human paraoxonase (PON) gene family has three isoforms: PON1, PON2 and PON3. These genes are implicated as potential risk factors of cerebrovascular disease and can prevent oxidative modification of low-density lipoproteins and atherosclerosis. This study evaluated the association between the genetic variants of all three PON genes and the risks of total stroke, ischemic stroke and hemorrhagic stroke in the Han Chinese population. METHODS: A total of 1016 subjects were recruited, including 508 healthy controls and 498 patients (328 with ischemic stroke and 170 with hemorrhagic stroke). A total of 11 single nucleotide polymorphisms (SNPs) covering the PON genes were genotyped for statistical analysis. Two of the 11 SNPs (rs662 and rs854560) were contextualized in a meta-analysis of ischemic stroke. RESULTS: The presence of rs705381 (-162) in the promoter region of PON1 was significantly associated with total stroke (P(adjusted) = 0.0007, OR = 0.57 [95% CI = 0.41-0.79]) and ischemic stroke (P(adjusted) = 0.0017, OR = 0.54 [95% CI = 0.37-0.79]) when analyzed using a dominant model, but was not associated with hemorrhagic stroke. There was also a nominal association between rs854571 (-824) and total stroke. Meta-analysis demonstrated a significant nominal association between rs662 and ischemic stroke, but there was no evidence of an association between rs662 and ischemic stroke risk in a single site association study. CONCLUSIONS: These findings indicate that polymorphisms of PON1 gene may be a risk factor of stroke.

[The development of fire-needle and the improvement of its application].

Based on the literature research on the development and application of fire-needle, the materials, the manufacturing process and application in clinic practices were studied. The major type, application and the representative experts of fire-needle in the contemporary era were summarized in this paper. The application fire-needle is more and more expanding in clinic practices, and the modern type of fire-needle should be promoted at present.

Transcriptome analysis of Hpa1Xoo transformed cotton revealed constitutive expression of genes in multiple signalling pathways related to disease resistance.

The transcriptome profile in leaves and roots of the transgenic cotton line T-34 expressing hpa1(Xoo) from Xanthomonas oryzae pv. oryzae was analysed using a customized 12k cotton cDNA microarray. A total of 530 cDNA transcripts involved in 34 pathways were differentially expressed in the transgenic line T-34, in which 123 differentially expressed genes were related to the cotton defence responses including the hypersensitive reaction, defence responses associated with the recognition of pathogen-derived elicitors, and defence signalling pathways mediated by salicylic acid, jasmonic acid, ethylene, auxin, abscicic acid, and Ca(2+). Furthermore, transcripts encoding various leucine-rich protein kinases and mitogen-activated protein kinases were up-regulated in the transgenic line T-34 and expression of transcripts related to the energy producing and consuming pathway was also increased, which suggested that the enhanced metabolism related to the host defence response in the transgenic line T-34 imposed an increased energy demand on the transgenic plant.

Genetic diversity of toxigenic and nontoxigenic Vibrio cholerae serogroups O1 and O139 revealed by array-based comparative genomic hybridization.

Toxigenic serogroups O1 and O139 of Vibrio cholerae may cause cholera epidemics or pandemics. Nontoxigenic strains within these serogroups also exist in the environment, and also some may cause sporadic cases of disease. Herein, we investigate the genomic diversity among toxigenic and nontoxigenic O1 and O139 strains by comparative genomic microarray hybridization with the genome of El Tor strain N16961 as a base. Conservation of the toxigenic O1 El Tor and O139 strains is found as previously reported, whereas accumulation of genome changes was documented in toxigenic El Tor strains isolated within the 40 years of the seventh pandemic. High phylogenetic diversity in nontoxigenic O1 and O139 strains is observed, and most of the genes absent from nontoxigenic strains are clustered together in the N16961 genome. By comparing these toxigenic and nontoxigenic strains, we observed that the small chromosome of V. cholerae is quite conservative and stable, outside of the superintegron region. In contrast to the general stability of the genome, the superintegron demonstrates pronounced divergence among toxigenic and nontoxigenic strains. Additionally, sequence variation in virulence-related genes is found in nontoxigenic El Tor strains, and we speculate that these intermediate strains may have pathogenic potential should they acquire CTX prophage alleles and other gene clusters. This genome-wide comparison of toxigenic and nontoxigenic V. cholerae strains may promote understanding of clonal differentiation of V. cholerae and contribute to an understanding of the origins and clonal selection of epidemic strains.

Genomewide oligonucleotide microarray analysis on placentae of pre-eclamptic pregnancies.

OBJECTIVE: Human placentae from normal and pre-eclamptic pregnancies were evaluated for possible changes in gene expression by microarray analysis to uncover new clues for the research of the etiology of pre-eclampsia. METHODS: Placentae from five normal pregnancies and five pregnancies complicated by pre-eclampsia were collected. mRNA levels of five pre-eclamptic placentae were examined using genome-wide 70-mer oligonucleotide microarrays (CapitalBio, Beijing, China) in comparison with the pooled control consisting of total RNA from five normotensive placentae. RESULTS: Ninety-six genes were found consistently down- or up-regulated in at least four pre-eclamptic samples. Most of them were related to an imbalance of reactive oxygen metabolites in placenta, abnormal trophoblast invasion, disorders of lipoprotein metabolism and signal transduction, or some that have been reported to have close correlation to the pathology of pre-eclampsia. The microarray data were also confirmed by the measurement of real-time PCR. CONCLUSION: DNA microarray is a high throughput and time-saving method to monitor altered gene expression. The results could provide interesting clues to the etiology of pre-eclampsia and lead to further studies in a more targeted fashion.