Multiregion single-cell sequencing reveals the transcriptional landscape of the immune microenvironment of colorectal cancer.

The tumor microenvironment is a complex ecosystem formed by distinct and interacting cell populations, and its composition is related to cancer prognosis and response to clinical treatment. In this study, we have taken the advantage of two single-cell RNA sequencing technologies (Smart-seq2 and DNBelab C4) to generate an atlas of 15,115 immune and nonimmune cells from primary tumors and hepatic metastases of 18 colorectal cancer (CRC) patients. We observed extensive changes in the proportions and functional states of T cells and B cells in tumor tissues, compared to those of paired non-tumor tissues. Importantly, we found that B cells from early CRC tumor were identified to be pre-B like expressing tumor suppressors, whereas B cells from advanced CRC tumors tended to be developed into plasma cells. We also identified the association of IgA+ IGLC2+ plasma cells with poor CRC prognosis, and demonstrated a significant interaction between B-cell and myeloid-cell signaling, and found CCL8+ cycling B cells/CCR5+ T-cell interactions as a potential antitumoral mechanism in advanced CRC tumors. Our results provide deeper insights into the immune infiltration within CRC, and a new perspective for the future research in immunotherapies for CRC.

Characterization of a pathogenic variant in GBA for Parkinson's disease with mild cognitive impairment patients.

Parkinson's disease (PD) is the second most common neurodegenerative disease, and mild cognitive impairment (MCI) is a well-established risk factor for the development of dementia in PD. A growing body of evidence suggests that low expression of glucocerebrosidase (GBA) promotes the transmission of α-synuclein (α-Syn) interpolymers and the progression of PD. However, how GBA mutations affect the pathogenesis of PD via abnormal aggregation of α-Syn is unclear, and no clinically valid PD-MCI genetic markers have been identified. Here, we first located a GBA eQTL, rs12411216, by analysing DHS, eQTL SNP, and transcription factor binding site data using the UCSC database. Subsequently, we found that rs12411216 was significantly associated with PD-MCI (P < 0.05) in 306 PD patients by genotyping. In exploring the relationship between rs12411216 and GBA expression, the SNP was found to be associated with GBA expression in 50 PD patients through qPCR verification. In a further CRISPR/Cas9-mediated genome editing module, the SNP was identified to cause a decrease in GBA expression, weaken enzymatic activity and enhance the abnormal aggregation of α-Syn in SH-SY5Y cells. Additionally, using an electrophoretic mobility shift assay, we confirmed that the binding efficiency of transcription factor E2F4 was affected by the rs12411216 SNP. In conclusion, our results showed that rs12411216 regulated GBA expression, supporting its potential role as a PD-MCI genetic biomarker and highlighting novel mechanisms underlying Parkinson's disease.

Transcriptome Analysis of the Acid Stress Response of Desulfovibrio vulgaris ATCC 7757.

The application of sulfate-reducing bacteria (SRB) shows great potential in the anaerobic biological treatment of acid mine wastewater; therefore, it has attracted much attention. The low pH in acidic wastewater affects the growth and reducing power of SRB. To uncover the mechanism underlying the reduction efficiency of SRB under acidic conditions, in this study, transcriptomic analysis was performed with Desulfovibrio vulgaris ATCC 7757 under three different pH conditions (pH 4.0, 5.5 and 7.0) and in the initial inoculation, logarithmic growth and plateau phases. Our results showed that ATCC 7757 still had biological activity at pH 4.0 and exhibited gene expression patterns at pH 4.0 that were different from those at pH 5.5 and pH 7. Importantly, the gene expression pattern was similar between pH 5.5 and pH 7. Transcriptomic analysis identified differentially expressed genes that affected the growth of ATCC 7757 under pH 7.0 at 22 h compared to 15 h; 196 of these genes were upregulated and 575 were downregulated. These differentially expressed genes were mainly enriched in genetic information processing and metabolism. Additionally, we identified 57 candidate genes associated with low-pH tolerance. Adaptation to low pH was reflected by an increase in the expression of genes involved in cell membrane structure and proton transport. The expression of genes involved in the reduction process decreased, including the genes DVU0499 and sat, which encode proteins that affect the sulfate reduction process. Both gene activities were validated by qPCR. Our results will contribute to further promoting the reducing power of SRB in acid mine wastewater and the development of successful bioremediation strategies.

Evolution and Comprehensive Analysis of DNaseI Hypersensitive Sites in Regulatory Regions of Primate Brain-Related Genes.

How the human brain differs from those of non-human primates is largely unknown and the complex drivers underlying such differences at the genomic level remain unclear. In this study, we selected 243 brain-related genes, based on Gene Ontology, and identified 184,113 DNaseI hypersensitive sites (DHSs) within their regulatory regions. To performed comprehensive evolutionary analyses, we set strict filtering criteria for alignment quality and filtered 39,132 DHSs for inclusion in the investigation and found that 2,397 (~6%) exhibited evidence of accelerated evolution (aceDHSs), which was a much higher proportion that DHSs genome-wide. Target genes predicted to be regulated by brain-aceDHSs were functionally enriched for brain development and exhibited differential expression between human and chimpanzee. Alignments indicated 61 potential human-specific transcription factor binding sites in brain-aceDHSs, including for CTCF, FOXH1, and FOXQ1. Furthermore, based on GWAS, Hi-C, and eQTL data, 16 GWAS SNPs, and 82 eQTL SNPs were in brain-aceDHSs that regulate genes related to brain development or disease. Among these brain-aceDHSs, we confirmed that one enhanced the expression of GPR133, using CRISPR-Cas9 and western blotting. The GPR133 gene is associated with glioblastoma, indicating that SNPs within DHSs could be related to brain disorders. These findings suggest that brain-related gene regulatory regions are under adaptive evolution and contribute to the differential expression profiles among primates, providing new insights into the genetic basis of brain phenotypes or disorders between humans and other primates.

Polymorphisms in the APOBEC3G gene of Chinese rhesus macaques affect resistance to ubiquitination and degradation mediated by HIV-2 Vif.

Animal cells have multiple innate effector mechanisms that inhibit viral replication. For the pathogenic retrovirus human immunodeficiency virus 1 (HIV-1), there are widely expressed restriction factors, such as APOBEC3 proteins, tetherin/BST2, SAMHD1 and MX2, as well as TRIM5α. We previously found that the TRIM5α gene clearly affects SIVmac or HIV-2 replication, but the major determinant of the combinatorial effect caused by multiple host restriction factors is still not fully clear. APOBEC3G (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3G), a host restriction factor that restricts HIV replication by causing cytosine deamination, can be targeted and degraded by the SIV/HIV-1/HIV-2 accessory protein Vif. Although rhesus macaques are widely used in HIV/AIDS research, little is known regarding the impact of APOBEC3G gene polymorphisms on viral Vif-mediated ubiquitin degradation in Chinese-origin rhesus macaques. In this study, we therefore genotyped APOBEC3G in 35 Chinese rhesus macaques. We identified a novel transcript and 27 APOBEC3G polymorphisms, including 20 non-synonymous variants and 7 synonymous mutation sites, of which 10 were novel. According to the predicted structure of the A3G protein, we predicted that the E88K and G212D mutations, both on the surface of the A3G protein, would have a significant effect on Vif-induced A3G degradation. However, an in vitro overexpression assay showed that these mutations did not influence HIV-2-Vif-mediated degradation of APOBEC3G. Unexpectedly, another polymorphism L71R, conferred resistance to Vif-mediated ubiquitin degradation, strongly suggesting that L71R might play an important role in antiviral defense mechanisms.

Evolution of DNAase I Hypersensitive Sites in MHC Regulatory Regions of Primates.

It has been challenging to determine the disease-causing variant(s) for most major histocompatibility complex (MHC)-associated diseases. However, it is becoming increasingly clear that regulatory variation is pervasive and a fundamentally important mechanism governing phenotypic diversity and disease susceptibility. We gathered DNase I data from 136 human cells to characterize the regulatory landscape of the MHC region, including 4867 DNase I hypersensitive sites (DHSs). We identified thousands of regulatory elements that have been gained or lost in the human or chimpanzee genomes since their evolutionary divergence. We compared alignments of the DHS across six primates and found 149 DHSs with convincing evidence of positive and/or purifying selection. Of these DHSs, compared to neutral sequences, 24 evolved rapidly in the human lineage. We identified 15 instances of transcription-factor-binding motif gains, such as USF, MYC, MAX, MAFK, STAT1, PBX3, etc, and observed 16 GWAS (genome-wide association study) SNPs associated with diseases within these 24 DHSs using FIMO (Find Individual Motif Occurrences) and UCSC (University of California, Santa Cruz) ChIP-seq data. Combining eQTL and Hi-C data, our results indicated that there were five SNPs located in human gains motifs affecting the corresponding gene's expression, two of which closely matched DHS target genes. In addition, a significant SNP, rs7756521, at genome-wide significant level likely affects DDR expression and represents a causal genetic variant for HIV-1 control. These results indicated that species-specific motif gains or losses of rapidly evolving DHSs in the primate genomes might play a role during adaptation evolution and provided some new evidence for a potentially causal role for these GWAS SNPs.

Positive selection of the TRIM family regulatory region in primate genomes.

Viral selection pressure has acted on restriction factors that play an important role in the innate immune system by inhibiting the replication of viruses during primate evolution. Tripartite motif-containing (TRIM) family members are some of these restriction factors. It is becoming increasingly clear that gene expression differences, rather than protein-coding regions changes, could play a vital role in the anti-retroviral immune mechanism. Increasingly, recent studies have created genome-scale catalogues of DNase I hypersensitive sites (DHSs), which demark potentially functional regulatory DNA. To improve our understanding of the molecular evolution mechanism of antiviral differences between species, we leveraged 14 130 DHSs derived from 145 cell types to characterize the regulatory landscape of the TRIM region. Subsequently, we compared the alignments of the DHSs across six primates and found 375 DHSs that are conserved in non-human primates but exhibit significantly accelerated rates of evolution in the human lineage (haDHSs). Furthermore, we discovered 31 human-specific potential transcription factor motifs within haDHSs, including the KROX and SP1, that both interact with HIV-1 Importantly, the corresponding haDHS was correlated with antiviral factor TRIM23 Thus, our results suggested that some viruses may contribute, through regulatory DNA differences, to organismal evolution by mediating TRIM gene expression to escape immune surveillance.

Comprehensive identification of high-frequency and combination MHC-DMA and -DMB alleles in a cohort of Chinese rhesus macaques and cynomolgus macaques of Vietnamese origin.

Rhesus and cynomolgus macaques are currently used as ideal animal models of immune response. Major histocompatibility complex (MHC) molecules play important roles in the susceptibility and/or resistance to many diseases. In this study, MHC-DMA and -DMB were first characterized by sequencing and cloning in 28 unrelated cynomolgus macaques from Vietnam and 34 unrelated Chinese rhesus macaques. A total of 23 novel alleles, including six high frequency alleles, were identified in this study. Our results showed that the alleles with the highest phenotypic frequencies were Mafa-DMA(∗)02:04:03 (57.1%), Mafa-DMB(∗)03:01:02 (76.9%), Mamu-DMA(∗)02:01:04 (88.2%), and Mamu-DMB(∗)03:02:02 (85%), respectively, indicating that distribution and frequencies of alleles had a few differences between Chinese rhesus macaques and cynomolgus macaques from Vietnam. Interestingly, compared to the cynomolgus macaques, we found that the combination of Mamu-DMA(∗)02:01:04-DMB(∗)03:02:02 was detected in 27 (79.4%) of 34 monkeys, suggesting that the combination of the MHC-DMA and -DMB alleles was probably a characteristic feature of the Chinese rhesus macaques. Our results will greatly increase the value of the two species as models for biomedical research.