Single-cell analysis of isoform switching and transposable element expression during preimplantation embryonic development.

Alternative splicing is an essential regulatory mechanism for development and pathogenesis. Through alternative splicing one gene can encode multiple isoforms and be translated into proteins with different functions. Therefore, this diversity is an important dimension to understand the molecular mechanism governing embryo development. Isoform expression in preimplantation embryos has been extensively investigated, leading to the discovery of new isoforms. However, the dynamics of isoform switching of different types of transcripts throughout the development remains unexplored. Here, using single-cell direct isoform sequencing in over 100 single blastomeres from the mouse oocyte to blastocyst stage, we quantified isoform expression and found that 3-prime partial transcripts lacking stop codons are highly accumulated in oocytes and zygotes. These transcripts are not transcription by-products and might play a role in maternal to zygote transition (MZT) process. Long-read sequencing also enabled us to determine the expression of transposable elements (TEs) at specific loci. In this way, we identified 3,894 TE loci that exhibited dynamic changes along the preimplantation development, likely regulating the expression of adjacent genes. Our work provides novel insights into the transcriptional regulation of early embryo development.

Intraocular mRNA delivery with endogenous MmPEG10-based virus-like particles.

Virus-like particles (VLP) are a promising tool for intracellular gene delivery, yet their potential in ocular gene therapy remains underexplored. In this study, we bridged this knowledge gap by demonstrating the successful generation and application of vesicular stomatitis virus glycoprotein (VSVG)-pseudotyped mouse PEG10 (MmPEG10)-VLP for intraocular mRNA delivery. Our findings revealed that PEG10-VLP can efficiently deliver GFP mRNA to adult retinal pigment epithelial cell line-19 (ARPE-19) cells, leading to transient expression. Moreover, we showed that MmPEG10-VLP can transfer SMAD7 to inhibit epithelial-mesenchymal transition (EMT) in RPE cells effectively. In vivo experiments further substantiated the potential of these vectors, as subretinal delivery into adult mice resulted in efficient transduction of retinal pigment epithelial (RPE) cells and GFP reporter gene expression without significant immune response. However, intravitreal injection did not yield efficient ocular expression. We also evaluated the transduction characteristics of MmPEG10-VLP following intracameral delivery, revealing transient GFP protein expression in corneal endothelial cells without significant immunotoxicities. In summary, our study established that VSVG pseudotyped MmPEG10-based VLP can transduce mitotically inactive RPE cells and corneal endothelial cells in vivo without triggering an inflammatory response, underscoring their potential utility in ocular gene therapy.

CD169+ classical monocyte as an important participant in Graves' ophthalmopathy through CXCL12-CXCR4 axis.

Patients with Graves' disease (GD) can develop Graves' ophthalmopathy (GO), but the underlying pathological mechanisms driving this development remain unclear. In our study, which included patients with GD and GO, we utilized single-cell RNA sequencing (scRNA-seq) and multiplatform analyses to investigate CD169+ classical monocytes, which secrete proinflammatory cytokines and are expanded through activated interferon signaling. We found that CD169+ clas_mono was clinically significant in predicting GO progression and prognosis, and differentiated into CD169+ macrophages that promote inflammation, adipogenesis, and fibrosis. Our murine model of early-stage GO showed that CD169+ classical monocytes accumulated in orbital tissue via the Cxcl12-Cxcr4 axis. Further studies are needed to investigate whether targeting circulating monocytes and the Cxcl12-Cxcr4 axis could alleviate GO progression.

Dynamic Transition of Regulatory T Cells to Cytotoxic Phenotype Amid Systemic Inflammation in Graves' Ophthalmopathy.

Graves' disease (GD) is an autoimmune condition that can progress to Graves' Ophthalmopathy (GO), leading to irreversible damage to orbital tissues and potential blindness. The pathogenic mechanism is not fully understood. In this study, we conducted single-cell multi-omics analyses on healthy individuals, GD patients without GO, newly diagnosed GO patients, and treated GO patients. Our findings revealed gradual systemic inflammation during GO progression, marked by overactivation of cytotoxic effector T cell subsets, and expansion of specific T cell receptor clones. Importantly, we observed a decline in the immunosuppressive function of activated regulatory T cells (aTreg) accompanied by a cytotoxic phenotypic transition. In vitro experiments revealed that dysfunction and transition of GO-autoreactive Treg were regulated by the yinyang1 (YY1) upon secondary stimulation of thyroid stimulating hormone receptor (TSHR) under inflammatory conditions. Furthermore, adoptive transfer experiments of GO mouse model confirmed infiltration of these cytotoxic Treg into the orbital lesion tissues. Notably, these cells were found to upregulate inflammation and promote pathogenic fibrosis of orbital fibroblasts (OFs). Our results revealed the dynamic changes in immune landscape during GO progression and provided novel insights into the instability and phenotypic transition of Treg, offering potential targets for therapeutic intervention and prevention of autoimmune diseases.

Single-Cell Transcriptomic Analysis Reveals Dynamic Cellular Processes in Corneal Epithelium During Wound Healing in Cynomolgus Monkeys.

Purpose: Corneal wounding healing is critical for maintaining clear vision, however, a complete understanding of its dynamic regulatory mechanisms remains elusive. Here, we used single-cell RNA sequencing (scRNA-seq) to analyze the cellular activities and transcriptional changes of corneal limbal epithelial cells at different stages after wound healing in cynomolgus monkeys, which exhibit a closer transcriptomic similarity to humans. Methods: Corneal limbal tissues were collected during uninjured, 1-day and 3-day healing stages, dissociated into single cells, and subjected to scRNA-seq using the 10× Genomics platform. Cell types were clustered by graph-based visualization methods and unbiased computational analysis. Additionally, cell migration assays and immunofluorescent staining were performed on cultured human corneal epithelial cells. Results: We characterized nine cell clusters by scRNA-seq analysis of the cynomolgus monkey corneal epithelium. By comparing heterogeneous transcriptional changes in major cell types during corneal healing, we highlighted the importance of limbal epithelial cells (LEPCs) and basal epithelial cells (BEPCs) in extracellular matrix (ECM) formation and wound healing, as well as suprabasal epithelial cells (SEPCs) in epithelial differentiation during the healing processes. We further identified five different sub-clusters in LEPC, including the transit amplifying cell (TAC) sub-cluster that promotes early healing through the activation of thrombospondin-1 (THBS1) expression. Conclusions: Our study represents the first comprehensive exploration of the detailed transcriptome profile of individual corneal cells during the wound healing process in nonhuman primates. We demonstrate the intricate mechanisms involved in corneal healing and provide a promising avenue for potential therapies in corneal wound healing.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.).

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

Study of LZ-word distribution and its application for sequence comparison.

Lempel-Ziv complexity has been widely used for sequence comparison and achieved promising results, but until now components' distribution in exhaustive history has not been studied. This paper investigated the whole distribution of LZ-words and presented a novel statistical method for sequence comparison. With the components' length in mind, we revised Lempel-Ziv complexity and obtained various sets of LZ-words. Instead of calculating the LZ-words' contents, we defined a series of set operations on LZ-word set to compare biological sequences. In order to assess the effectiveness of the proposed method, we performed two sets of experiments and compared it with alignment-based methods.

Single-cell chromatin accessibility and transcriptomic characterization of Behcet's disease.

Behect's disease is a chronic vasculitis characterized by complex multi-organ immune aberrations. However, a comprehensive understanding of the gene-regulatory profile of peripheral autoimmunity and the diverse immune responses across distinct cell types in Behcet's disease (BD) is still lacking. Here, we present a multi-omic single-cell study of 424,817 cells in BD patients and non-BD individuals. This study maps chromatin accessibility and gene expression in the same biological samples, unraveling vast cellular heterogeneity. We identify widespread cell-type-specific, disease-associated active and pro-inflammatory immunity in both transcript and epigenomic aspects. Notably, integrative multi-omic analysis reveals putative TF regulators that might contribute to chromatin accessibility and gene expression in BD. Moreover, we predicted gene-regulatory networks within nominated TF activators, including AP-1, NF-kB, and ETS transcript factor families, which may regulate cellular interaction and govern inflammation. Our study illustrates the epigenetic and transcriptional landscape in BD peripheral blood and expands understanding of potential epigenomic immunopathology in this disease.

Assessing the effects of aging on the liver endothelial cell landscape using single-cell RNA sequencing.

Endothelial cell (EC) function declines with age and contributes to the development of many vascular-related disease processes. Currently, the effects of aging on the molecular regulatory mechanisms of liver ECs have not been fully elucidated. Here, we employed single-cell RNA sequencing to map the transcriptome of ECs and analyzed their relationship with aging. We identified 8 different EC subtypes, interestingly, 2 of which were specially expressed in aged mice ECs namely aged capillary ECs (Aged ECs) and pro-inflammation capillary ECs (Proinfla.ECs). Double immunostaining for an EC marker (Cd31) and a marker of these specialized EC phenotypes confirmed the single-cell RNA sequencing data. Gene ontology analysis revealed that Aged ECs and Proinfla.ECs were associated with inflammatory response. Then we found that liver proliferating capillary ECs (Prolife.ECs) were most affected by senescence. Single-cell transcript analysis suggests that Prolife.ECs and angiogenic capillary ECs may form a poor microenvironment that promotes angiogenesis and tumorigenesis. Pseudo-temporal trajectories revealed that Prolife.ECs have different differentiation pathways in young and aged mice. In aged mice, Prolife.ECs could specifically differentiate into an unstable state, which was mainly composed of angiogenic capillary ECs. Intercellular communication revealed inflammatory activation in old group. Overall, this work compared the single-cell RNA profiles of liver ECs in young and aged mice. These findings provide a new insight into liver aging and its molecular mechanisms, and further exploration of Aged ECs and Proinfla.ECs may help to elucidate the molecular mechanisms associated with senescence.

Insights gained from single-cell RNA analysis of murine endothelial cells in aging hearts.

Aging is the strongest risk factor for cardiovascular disease, with progressive decline in the function of vascular endothelial cells (ECs) with age. Systematic analyses of the effects of aging on different cardiac EC types remain limited. Here, we constructed a scRNA atlas of EC transcriptomes in young and old mouse hearts. We identified 10 EC subclusters. The multidimensionally differential genes (DEGs) analysis across different EC clusters shows molecular changes with aging, showing the increase in the overall inflammatory microenvironment and the decrease in angiogenesis and cytoskeletal support capacity of aged ECs. And we performed an in-depth analysis of 3 special ECs, Immunology, Proliferating and Angiogenic. The Immunology EC seems highly associated with some immune regulatory functions, which decline with aging at different degrees. Analysis of two types of neovascular ECs, Proliferating, Angiogenic, implied that Angiogenic ECs can differentiate into multiple EC directions after initially originating from proliferating ECs. And aging leads to a decrease in the ability of vascular angiogenesis and differentiation. Finally, we summarized the effects of aging on cell signaling communication between different EC clusters. This cardiac EC atlas offers comprehensive insights into the molecular regulations of cardiovascular aging, and provides new directions for the prevention and treatment of age-related cardiovascular disease.

Integrated Transcriptome Analysis of Long Noncoding RNA and mRNA in Developing and Aging Mouse Retina.

Mice have emerged as a widely employed model for investigating various retinal diseases. However, the availability of comprehensive datasets capturing the entire developmental and aging stages of the mouse retina, particularly during the elderly period, encompassing integrated lncRNA and mRNA expression profiles, is limited. In this study, we assembled a total of 18 retina samples from mice across 6 distinct stages of development and aging (5 days, 3 weeks, 6 weeks, 10 weeks, 6 months, and 15 months) to conduct integrated lncRNA and mRNA sequencing analysis. This invaluable dataset offers a comprehensive transcriptomic resource of mRNA and lncRNA expression profiles during the natural progression of retinal development and aging. The discoveries stemming from this investigation will significantly contribute to the elucidation of the underlying molecular mechanisms associated with various retinal diseases, such as congenital retinal dysplasia and retinal degenerative diseases.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Mouse Corneal Epithelial Cells.

Purpose: Corneal epithelial homeostasis is maintained by coordinated gene expression across distinct cell populations, but the gene regulatory programs underlying this cellular diversity remain to be characterized. Here we applied single-cell multi-omics analysis to delineate the gene regulatory profile of mouse corneal epithelial cells under normal homeostasis. Methods: Single cells isolated from the cornea epithelium (with marginal conjunctiva) of adult mice were subjected to scRNA-seq and scATAC-seq using the 10×Genomics platform. Cell types were clustered by the graph-based visualization method uniform manifold approximation and projection and unbiased computational informatics analysis. The scRNA-seq and scATAC-seq datasets were integrated following the integration pipeline described in ArchR and Seurat. Results: We characterized diverse corneal epithelial cell types based on gene expression signatures and chromatin accessibility. We found that cell type-specific accessibility regions were mainly located at distal regions, suggesting essential roles of distal regulatory elements in determining corneal epithelial cell diversity. Trajectory analyses revealed a continuum of cell state transition and higher coordination between transcription factor (TF) motif accessibility and gene expression during corneal epithelial cell differentiation. By integrating transcriptomic and chromatin accessibility analysis, we identified cell type-specific and shared gene regulation programs. We also uncovered critical TFs driving corneal epithelial cell differentiation, such as nuclear factor I (NFI) family members, Rarg, Elf3. We found that nuclear factor-κB (NF-κB) family members were positive TFs in limbal cells and some superficial cells, but they were involved in regulating distinct biological processes. Conclusions: Our study presents a comprehensive gene regulatory landscape of mouse cornea epithelial cells, and provides valuable foundations for future investigation of corneal epithelial homeostasis in the context of cornea pathologies and regenerative medicine.

High-throughput and high-accuracy single-cell RNA isoform analysis using PacBio circular consensus sequencing.

Although long-read single-cell RNA isoform sequencing (scISO-Seq) can reveal alternative RNA splicing in individual cells, it suffers from a low read throughput. Here, we introduce HIT-scISOseq, a method that removes most artifact cDNAs and concatenates multiple cDNAs for PacBio circular consensus sequencing (CCS) to achieve high-throughput and high-accuracy single-cell RNA isoform sequencing. HIT-scISOseq can yield >10 million high-accuracy long-reads in a single PacBio Sequel II SMRT Cell 8M. We also report the development of scISA-Tools that demultiplex HIT-scISOseq concatenated reads into single-cell cDNA reads with >99.99% accuracy and specificity. We apply HIT-scISOseq to characterize the transcriptomes of 3375 corneal limbus cells and reveal cell-type-specific isoform expression in them. HIT-scISOseq is a high-throughput, high-accuracy, technically accessible method and it can accelerate the burgeoning field of long-read single-cell transcriptomics.

Chromatin accessibility analysis reveals regulatory dynamics and therapeutic relevance of Vogt-Koyanagi-Harada disease.

The barrier to curing Vogt-Koyanagi-Harada disease (VKH) is thought to reside in a lack of understanding in the roles and regulations of peripheral inflammatory immune cells. Here we perform a single-cell multi-omic study of 166,149 cells in peripheral blood mononuclear cells from patients with VKH, profile the chromatin accessibility and gene expression in the same blood samples, and uncover prominent cellular heterogeneity. Immune cells in VKH blood are highly activated and pro-inflammatory. Notably, we describe an enrichment of transcription targets for nuclear factor kappa B in conventional dendritic cells (cDCs) that governed inflammation. Integrative analysis of transcriptomic and chromatin maps shows that the RELA in cDCs is related to disease complications and poor prognosis. Ligand-receptor interaction pairs also identify cDC as an important predictor that regulated multiple immune subsets. Our results reveal epigenetic and transcriptional dynamics in auto-inflammation, especially the cDC subtype that might lead to therapeutic strategies in VKH.

Genome assembly of a tropical maize inbred line provides insights into structural variation and crop improvement.

Maize is one of the most important crops globally, and it shows remarkable genetic diversity. Knowledge of this diversity could help in crop improvement; however, gold-standard genomes have been elucidated only for modern temperate varieties. Here, we present a high-quality reference genome (contig N50 of 15.78 megabases) of the maize small-kernel inbred line, which is derived from a tropical landrace. Using haplotype maps derived from B73, Mo17 and SK, we identified 80,614 polymorphic structural variants across 521 diverse lines. Approximately 22% of these variants could not be detected by traditional single-nucleotide-polymorphism-based approaches, and some of them could affect gene expression and trait performance. To illustrate the utility of the diverse SK line, we used it to perform map-based cloning of a major effect quantitative trait locus controlling kernel weight-a key trait selected during maize improvement. The underlying candidate gene ZmBARELY ANY MERISTEM1d provides a target for increasing crop yields.

Subcellular localization prediction of apoptosis proteins based on evolutionary information and support vector machine.

OBJECTIVES: In this paper, a high-quality sequence encoding scheme is proposed for predicting subcellular location of apoptosis proteins. METHODS: In the proposed methodology, the novel evolutionary-conservative information is introduced to represent protein sequences. Meanwhile, based on the proportion of golden section in mathematics, position-specific scoring matrix (PSSM) is divided into several blocks. Then, these features are predicted by support vector machine (SVM) and the predictive capability of proposed method is implemented by jackknife test RESULTS: The results show that the golden section method is better than no segmentation method. The overall accuracy for ZD98 and CL317 is 98.98% and 91.11%, respectively, which indicates that our method can play a complimentary role to the existing methods in the relevant areas. CONCLUSIONS: The proposed feature representation is powerful and the prediction accuracy will be improved greatly, which denotes our method provides the state-of-the-art performance for predicting subcellular location of apoptosis proteins.

A protein structural classes prediction method based on predicted secondary structure and PSI-BLAST profile.

Knowledge of protein secondary structural classes plays an important role in understanding protein folding patterns. In this paper, 25 features based on position-specific scoring matrices are selected to reflect evolutionary information. In combination with other 11 rational features based on predicted protein secondary structure sequences proposed by the previous researchers, a 36-dimensional representation feature vector is presented to predict protein secondary structural classes for low-similarity sequences. ASTRALtraining dataset is used to train and design our method, other three low-similarity datasets ASTRALtest, 25PDB and 1189 are used to test the proposed method. Comparisons with other methods show that our method is effective to predict protein secondary structural classes. Stand alone version of the proposed method (PSSS-PSSM) is written in MATLAB language and it can be downloaded from http://letsgob.com/bioinfo_PSSS_PSSM/.