Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae.

BACKGROUND: Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown. RESULTS: We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions. CONCLUSIONS: This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.

Cellular heterogeneity and transcriptomic profiles during intrahepatic cholangiocarcinoma initiation and progression.

BACKGROUND AND AIMS: Intrahepatic cholangiocarcinoma (ICC) is not fully investigated, and how stromal cells contribute to ICC formation is poorly understood. We aimed to uncover ICC origin, cellular heterogeneity, and critical modulators during ICC initiation/progression, and to decipher how fibroblast and endothelial cells in the stromal compartment favor ICC progression. APPROACH AND RESULTS: We performed single-cell RNA sequencing (scRNA-seq) using AKT/Notch intracellular domain-induced mouse ICC tissues at early, middle, and late stages. We analyzed the transcriptomic landscape, cellular classification and evolution, and intercellular communication during ICC initiation/progression. We confirmed the findings using quantitative real-time PCR, western blotting, immunohistochemistry or immunofluorescence, and gene knockout/knockdown analysis. We identified stress-responding and proliferating subpopulations in late-stage mouse ICC tissues and validated them using human scRNA-seq data sets. By integrating weighted correlation network analysis and protein-protein interaction through least absolute shrinkage and selection operator regression, we identified zinc finger, MIZ-type containing 1 (Zmiz1) and Y box protein 1 (Ybx1) as core transcription factors required by stress-responding and proliferating ICC cells, respectively. Knockout of either one led to the blockade of ICC initiation/progression. Using two other ICC mouse models (YAP/AKT, KRAS/p19) and human ICC scRNA-seq data sets, we confirmed the orchestrating roles of Zmiz1 and Ybx1 in ICC occurrence and development. In addition, hes family bHLH transcription factor 1, cofilin 1, and inhibitor of DNA binding 1 were identified as driver genes for ICC. Moreover, periportal liver sinusoidal endothelial cells could differentiate into tip endothelial cells to promote ICC development, and this was Dll4-Notch4-Efnb2 signaling-dependent. CONCLUSIONS: Stress-responding and ICC proliferating subtypes were identified, and Zmiz1 and Ybx1 were revealed as core transcription factors in these subtypes. Fibroblast-endothelial cell interaction promotes ICC development.

gCAnno: a graph-based single cell type annotation method.

BACKGROUND: Current single cell analysis methods annotate cell types at cluster-level rather than ideally at single cell level. Multiple exchangeable clustering methods and many tunable parameters have a substantial impact on the clustering outcome, often leading to incorrect cluster-level annotation or multiple runs of subsequent clustering steps. To address these limitations, methods based on well-annotated reference atlas has been proposed. However, these methods are currently not robust enough to handle datasets with different noise levels or from different platforms. RESULTS: Here, we present gCAnno, a graph-based Cell type Annotation method. First, gCAnno constructs cell type-gene bipartite graph and adopts graph embedding to obtain cell type specific genes. Then, naïve Bayes (gCAnno-Bayes) and SVM (gCAnno-SVM) classifiers are built for annotation. We compared the performance of gCAnno to other state-of-art methods on multiple single cell datasets, either with various noise levels or from different platforms. The results showed that gCAnno outperforms other state-of-art methods with higher accuracy and robustness. CONCLUSIONS: gCAnno is a robust and accurate cell type annotation tool for single cell RNA analysis. The source code of gCAnno is publicly available at https://github.com/xjtu-omics/gCAnno .

Haplotype-resolved assemblies and variant benchmark of a Chinese Quartet.

BACKGROUND: Recent state-of-the-art sequencing technologies enable the investigation of challenging regions in the human genome and expand the scope of variant benchmarking datasets. Herein, we sequence a Chinese Quartet, comprising two monozygotic twin daughters and their biological parents, using four short and long sequencing platforms (Illumina, BGI, PacBio, and Oxford Nanopore Technology). RESULTS: The long reads from the monozygotic twin daughters are phased into paternal and maternal haplotypes using the parent-child genetic map and for each haplotype. We also use long reads to generate haplotype-resolved whole-genome assemblies with completeness and continuity exceeding that of GRCh38. Using this Quartet, we comprehensively catalogue the human variant landscape, generating a dataset of 3,962,453 SNVs, 886,648 indels (< 50 bp), 9726 large deletions (≥ 50 bp), 15,600 large insertions (≥ 50 bp), 40 inversions, 31 complex structural variants, and 68 de novo mutations which are shared between the monozygotic twin daughters. Variants underrepresented in previous benchmarks owing to their complexity-including those located at long repeat regions, complex structural variants, and de novo mutations-are systematically examined in this study. CONCLUSIONS: In summary, this study provides high-quality haplotype-resolved assemblies and a comprehensive set of benchmarking resources for two Chinese monozygotic twin samples which, relative to existing benchmarks, offers expanded genomic coverage and insight into complex variant categories.

Integrating bulk and single-cell RNA sequencing reveals cellular heterogeneity and immune infiltration in hepatocellular carcinoma.

Efficacy of immunotherapy in hepatocellular carcinoma (HCC) is blocked by its high degree of inter- and intra-tumor heterogeneity and immunosuppressive tumor microenvironment. However, the correlation between tumor heterogeneity and immunosuppressive microenvironment in HCC has not been well addressed. Here, we endeavored to dissect inter- and intra-tumor heterogeneity in HCC and uncover how they contribute to the immunosuppressive microenvironment. We performed consensus molecular subtyping with non-negative matrix factorization (NMF) clustering to stratify the inter-heterogeneity profile of HCC tumors. We grouped HCC tumors from the Cancer Genome Atlas (TCGA) patients into three subtypes (S1, S2 and S3), where S1 was characterized as a 'hot tumor' profile with high expression level of T cell genes and rate of immune scores. S2 was characterized as a 'cold tumor' profile with the highest tumor purity score, and S3 as an 'immunosuppressed tumor' profile with the poorest prognosis and a high expression level of immunosuppressive genes such as cytotoxic T-lymphocyte-associated protein-4, TIGIT, and PDCD1. Moreover, we combined weighted gene co-expression network analysis and single-cell regulatory network inference and clustering (SCENIC) in the single-cell dataset of the S3-like subtype (CS3) and identified a transcription factor, BATF, which could upregulate immunosuppressive genes. Finally, we identified a cell interaction network in which a myeloid-derived suppressor cell-like macrophage subtype could promote the formation of immunosuppressive T-cells.

Integrating cell interaction with transcription factors to obtain a robust gene panel for prognostic prediction and therapies in cholangiocarcinoma.

Objective: The efficacy of immunotherapy for cholangiocarcinoma (CCA) is blocked by a high degree of tumor heterogeneity. Cell communication contributes to heterogeneity in the tumor microenvironment. This study aimed to explore critical cell signaling and biomarkers induced via cell communication during immune exhaustion in CCA. Methods: We constructed empirical Bayes and Markov random field models eLBP to determine transcription factors, interacting genes, and associated signaling pathways involved in cell-cell communication using single-cell RNAseq data. We then analyzed the mechanism of immune exhaustion during CCA progression. Results: We found that VEGFA-positive macrophages with high levels of LGALS9 could interact with HAVCR2 to promote the exhaustion of CD8+ T cells in CCA. Transcription factors SPI1 and IRF1 can upregulate the expression of LGALS9 in VEGFA-positive macrophages. Subsequently, we obtained a panel containing 54 genes through the model, which identified subtype S2 with high expression of immune checkpoint genes that are suitable for immunotherapy. Moreover, we found that patients with subtype S2 with a higher mutation ratio of MUC16 had immune-exhausted genes, such as HAVCR2 and TIGIT. Finally, we constructed a nine-gene eLBP-LASSO-COX risk model, which was designated the tumor microenvironment risk score (TMRS). Conclusion: Cell communication-related genes can be used as important markers for predicting patient prognosis and immunotherapy responses. The TMRS panel is a reliable tool for prognostic prediction and chemotherapeutic decision-making in CCA.

High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads.

Arabidopsis thaliana is an important and long-established model species for plant molecular biology, genetics, epigenetics, and genomics. However, the latest version of reference genome still contains a significant number of missing segments. Here, we reported a high-quality and almost complete Col-0 genome assembly with two gaps (named Col-XJTU) by combining the Oxford Nanopore Technologies ultra-long reads, Pacific Biosciences high-fidelity long reads, and Hi-C data. The total genome assembly size is 133,725,193 bp, introducing 14.6 Mb of novel sequences compared to the TAIR10.1 reference genome. All five chromosomes of the Col-XJTU assembly are highly accurate with consensus quality (QV) scores > 60 (ranging from 62 to 68), which are higher than those of the TAIR10.1 reference (ranging from 45 to 52). We completely resolved chromosome (Chr) 3 and Chr5 in a telomere-to-telomere manner. Chr4 was completely resolved except the nucleolar organizing regions, which comprise long repetitive DNA fragments. The Chr1 centromere (CEN1), reportedly around 9 Mb in length, is particularly challenging to assemble due to the presence of tens of thousands of CEN180 satellite repeats. Using the cutting-edge sequencing data and novel computational approaches, we assembled a 3.8-Mb-long CEN1 and a 3.5-Mb-long CEN2. We also investigated the structure and epigenetics of centromeres. Four clusters of CEN180 monomers were detected, and the centromere-specific histone H3-like protein (CENH3) exhibited a strong preference for CEN180 Cluster 3. Moreover, we observed hypomethylation patterns in CENH3-enriched regions. We believe that this high-quality genome assembly, Col-XJTU, would serve as a valuable reference to better understand the global pattern of centromeric polymorphisms, as well as the genetic and epigenetic features in plants.

A novel CD4+ CTL subtype characterized by chemotaxis and inflammation is involved in the pathogenesis of Graves' orbitopathy.

Graves' orbitopathy (GO), the most severe manifestation of Graves' hyperthyroidism (GH), is an autoimmune-mediated inflammatory disorder, and treatments often exhibit a low efficacy. CD4+ T cells have been reported to play vital roles in GO progression. To explore the pathogenic CD4+ T cell types that drive GO progression, we applied single-cell RNA sequencing (scRNA-Seq), T cell receptor sequencing (TCR-Seq), flow cytometry, immunofluorescence and mixed lymphocyte reaction (MLR) assays to evaluate CD4+ T cells from GO and GH patients. scRNA-Seq revealed the novel GO-specific cell type CD4+ cytotoxic T lymphocytes (CTLs), which are characterized by chemotactic and inflammatory features. The clonal expansion of this CD4+ CTL population, as demonstrated by TCR-Seq, along with their strong cytotoxic response to autoantigens, localization in orbital sites, and potential relationship with disease relapse provide strong evidence for the pathogenic roles of GZMB and IFN-γ-secreting CD4+ CTLs in GO. Therefore, cytotoxic pathways may become potential therapeutic targets for GO.

Whole-genome sequencing of 175 Mongolians uncovers population-specific genetic architecture and gene flow throughout North and East Asia.

The genetic variation in Northern Asian populations is currently undersampled. To address this, we generated a new genetic variation reference panel by whole-genome sequencing of 175 ethnic Mongolians, representing six tribes. The cataloged variation in the panel shows strong population stratification among these tribes, which correlates with the diverse demographic histories in the region. Incorporating our results with the 1000 Genomes Project panel identifies derived alleles shared between Finns and Mongolians/Siberians, suggesting that substantial gene flow between northern Eurasian populations has occurred in the past. Furthermore, we highlight that North, East, and Southeast Asian populations are more aligned with each other than these groups are with South Asian and Oceanian populations.

DNA barcoding and the identification of tree frogs (Amphibia: Anura: Rhacophoridae).

The DNA barcoding gene COI (cytochrome c oxidase subunit I) effectively identifies many species. Herein, we barcoded 172 individuals from 37 species belonging to nine genera in Rhacophoridae to test if the gene serves equally well to identify species of tree frogs. Phenetic neighbor joining and phylogenetic Bayesian inference were used to construct phylogenetic trees, which resolved all nine genera as monophyletic taxa except for Rhacophorus, two new matrilines for Liuixalus, and Polypedates leucomystax species complex. Intraspecific genetic distances ranged from 0.000 to 0.119 and interspecific genetic distances ranged from 0.015 to 0.334. Within Rhacophorus and Kurixalus, the intra- and interspecific genetic distances did not reveal an obvious barcode gap. Notwithstanding, we found that COI sequences unambiguously identified rhacophorid species and helped to discover likely new cryptic species via the synthesis of genealogical relationships and divergence patterns. Our results supported that COI is an effective DNA barcoding marker for Rhacophoridae.