Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

HiCAR is a robust and sensitive method to analyze open-chromatin-associated genome organization.

The long-range interactions of cis-regulatory elements (cREs) play a central role in gene regulation. cREs can be characterized as accessible chromatin sequences. However, it remains technically challenging to comprehensively identify their spatial interactions. Here, we report a new method HiCAR (Hi-C on accessible regulatory DNA), which utilizes Tn5 transposase and chromatin proximity ligation, for the analysis of open-chromatin-anchored interactions with low-input cells. By applying HiCAR in human embryonic stem cells and lymphoblastoid cells, we demonstrate that HiCAR identifies high-resolution chromatin contacts with an efficiency comparable with that of in situ Hi-C over all distance ranges. Interestingly, we found that the "poised" gene promoters exhibit silencer-like function to repress the expression of distal genes via promoter-promoter interactions. Lastly, we applied HiCAR to 30,000 primary human muscle stem cells and demonstrated that HiCAR is capable of analyzing chromatin accessibility and looping using low-input primary cells and clinical samples.

Cohesin-mediated loop anchors confine the locations of human replication origins.

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

Human fetal cerebellar cell atlas informs medulloblastoma origin and oncogenesis.

Medulloblastoma (MB) is the most common malignant childhood brain tumour1,2, yet the origin of the most aggressive subgroup-3 form remains elusive, impeding development of effective targeted treatments. Previous analyses of mouse cerebella3-5 have not fully defined the compositional heterogeneity of MBs. Here we undertook single-cell profiling of freshly isolated human fetal cerebella to establish a reference map delineating hierarchical cellular states in MBs. We identified a unique transitional cerebellar progenitor connecting neural stem cells to neuronal lineages in developing fetal cerebella. Intersectional analysis revealed that the transitional progenitors were enriched in aggressive MB subgroups, including group 3 and metastatic tumours. Single-cell multi-omics revealed underlying regulatory networks in the transitional progenitor populations, including transcriptional determinants HNRNPH1 and SOX11, which are correlated with clinical prognosis in group 3 MBs. Genomic and Hi-C profiling identified de novo long-range chromatin loops juxtaposing HNRNPH1/SOX11-targeted super-enhancers to cis-regulatory elements of MYC, an oncogenic driver for group 3 MBs. Targeting the transitional progenitor regulators inhibited MYC expression and MYC-driven group 3 MB growth. Our integrated single-cell atlases of human fetal cerebella and MBs show potential cell populations predisposed to transformation and regulatory circuitries underlying tumour cell states and oncogenesis, highlighting hitherto unrecognized transitional progenitor intermediates predictive of disease prognosis and potential therapeutic vulnerabilities.

Subtype-specific 3D genome alteration in acute myeloid leukaemia.

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome.

Differential methylation of the two parental genomes in placental mammals is essential for genomic imprinting and embryogenesis. To systematically study this epigenetic process, we have generated a base-resolution, allele-specific DNA methylation (ASM) map in the mouse genome. We find parent-of-origin dependent (imprinted) ASM at 1,952 CG dinucleotides. These imprinted CGs form 55 discrete clusters including virtually all known germline differentially methylated regions (DMRs) and 23 previously unknown DMRs, with some occurring at microRNA genes. We also identify sequence-dependent ASM at 131,765 CGs. Interestingly, methylation at these sites exhibits a strong dependence on the immediate adjacent bases, allowing us to define a conserved sequence preference for the mammalian DNA methylation machinery. Finally, we report a surprising presence of non-CG methylation in the adult mouse brain, with some showing evidence of imprinting. Our results provide a resource for understanding the mechanisms of imprinting and allele-specific gene expression in mammalian cells.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding.

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

DDX5 inhibits inflammation by modulating m6A levels of TLR2/4 transcripts during bacterial infection.

DExD/H-box helicases are crucial regulators of RNA metabolism and antiviral innate immune responses; however, their role in bacteria-induced inflammation remains unclear. Here, we report that DDX5 interacts with METTL3 and METTL14 to form an m6A writing complex, which adds N6-methyladenosine to transcripts of toll-like receptor (TLR) 2 and TLR4, promoting their decay via YTHDF2-mediated RNA degradation, resulting in reduced expression of TLR2/4. Upon bacterial infection, DDX5 is recruited to Hrd1 at the endoplasmic reticulum in an MyD88-dependent manner and is degraded by the ubiquitin-proteasome pathway. This process disrupts the DDX5 m6A writing complex and halts m6A modification as well as degradation of TLR2/4 mRNAs, thereby promoting the expression of TLR2 and TLR4 and downstream NF-κB activation. The role of DDX5 in regulating inflammation is also validated in vivo, as DDX5- and METTL3-KO mice exhibit enhanced expression of inflammatory cytokines. Our findings show that DDX5 acts as a molecular switch to regulate inflammation during bacterial infection and shed light on mechanisms of quiescent inflammation during homeostasis.

A map of cis-regulatory elements and 3D genome structures in zebrafish.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species1. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Epigenomic analysis reveals prevalent contribution of transposable elements to cis-regulatory elements, tissue-specific expression, and alternative promoters in zebrafish.

Transposable elements (TEs) encode regulatory elements that impact gene expression in multiple species, yet a comprehensive analysis of zebrafish TEs in the context of gene regulation is lacking. Here, we systematically investigate the epigenomic and transcriptomic landscape of TEs across 11 adult zebrafish tissues using multidimensional sequencing data. We find that TEs contribute substantially to a diverse array of regulatory elements in the zebrafish genome and that 37% of TEs are positioned in active regulatory states in adult zebrafish tissues. We identify TE subfamilies enriched in highly specific regulatory elements among different tissues. We use transcript assembly to discover TE-derived transcriptional units expressed across tissues. Finally, we show that novel TE-derived promoters can initiate tissue-specific transcription of alternate gene isoforms. This work provides a comprehensive profile of TE activity across normal zebrafish tissues, shedding light on mechanisms underlying the regulation of gene expression in this widely used model organism.

Genomic analyses of wild argali, domestic sheep, and their hybrids provide insights into chromosome evolution, phenotypic variation, and germplasm innovation.

Understanding the genetic mechanisms of phenotypic variation in hybrids between domestic animals and their wild relatives may aid germplasm innovation. Here, we report the high-quality genome assemblies of a male Pamir argali (O ammon polii, 2n = 56), a female Tibetan sheep (O aries, 2n = 54), and a male hybrid of Pamir argali and domestic sheep, and the high-throughput sequencing of 425 ovine animals, including the hybrids of argali and domestic sheep. We detected genomic synteny between Chromosome 2 of sheep and two acrocentric chromosomes of argali. We revealed consistent satellite repeats around the chromosome breakpoints, which could have resulted in chromosome fusion. We observed many more hybrids with karyotype 2n = 54 than with 2n = 55, which could be explained by the selfish centromeres, the possible decreased rate of normal/balanced sperm, and the increased incidence of early pregnancy loss in the aneuploid ewes or rams. We identified genes and variants associated with important morphological and production traits (e.g., body weight, cannon circumference, hip height, and tail length) that show significant variations. We revealed a strong selective signature at the mutation (c.334C > A, p.G112W) in TBXT and confirmed its association with tail length among sheep populations of wide geographic and genetic origins. We produced an intercross population of 110 F2 offspring with varied number of vertebrae and validated the causal mutation by whole-genome association analysis. We verified its function using CRISPR-Cas9 genome editing. Our results provide insights into chromosomal speciation and phenotypic evolution and a foundation of genetic variants for the breeding of sheep and other animals.

Accurate identification of structural variations from cancer samples.

Structural variations (SVs) are commonly found in cancer genomes. They can cause gene amplification, deletion and fusion, among other functional consequences. With an average read length of hundreds of kilobases, nano-channel-based optical DNA mapping is powerful in detecting large SVs. However, existing SV calling methods are not tailored for cancer samples, which have special properties such as mixed cell types and sub-clones. Here we propose the Cancer Optical Mapping for detecting Structural Variations (COMSV) method that is specifically designed for cancer samples. It shows high sensitivity and specificity in benchmark comparisons. Applying to cancer cell lines and patient samples, COMSV identifies hundreds of novel SVs per sample.

Single-Nuclei RNA-Sequencing of the Gastrocnemius Muscle in Peripheral Artery Disease.

BACKGROUND: Lower extremity peripheral artery disease (PAD) is a growing epidemic with limited effective treatment options. Here, we provide a single-nuclei atlas of PAD limb muscle to facilitate a better understanding of the composition of cells and transcriptional differences that comprise the diseased limb muscle. METHODS: We obtained gastrocnemius muscle specimens from 20 patients with PAD and 12 non-PAD controls. Nuclei were isolated and single-nuclei RNA-sequencing was performed. The composition of nuclei was characterized by iterative clustering via principal component analysis, differential expression analysis, and the use of known marker genes. Bioinformatics analysis was performed to determine differences in gene expression between PAD and non-PAD nuclei, as well as subsequent analysis of intercellular signaling networks. Additional histological analyses of muscle specimens accompany the single-nuclei RNA-sequencing atlas. RESULTS: Single-nuclei RNA-sequencing analysis indicated a fiber type shift with patients with PAD having fewer type I (slow/oxidative) and more type II (fast/glycolytic) myonuclei compared with non-PAD, which was confirmed using immunostaining of muscle specimens. Myonuclei from PAD displayed global upregulation of genes involved in stress response, autophagy, hypoxia, and atrophy. Subclustering of myonuclei also identified populations that were unique to PAD muscle characterized by metabolic dysregulation. PAD muscles also displayed unique transcriptional profiles and increased diversity of transcriptomes in muscle stem cells, regenerating myonuclei, and fibro-adipogenic progenitor cells. Analysis of intercellular communication networks revealed fibro-adipogenic progenitors as a major signaling hub in PAD muscle, as well as deficiencies in angiogenic and bone morphogenetic protein signaling which may contribute to poor limb function in PAD. CONCLUSIONS: This reference single-nuclei RNA-sequencing atlas provides a comprehensive analysis of the cell composition, transcriptional signature, and intercellular communication pathways that are altered in the PAD condition.

Activation of the Aryl Hydrocarbon Receptor in Muscle Exacerbates Ischemic Pathology in Chronic Kidney Disease.

BACKGROUND: Chronic kidney disease (CKD) accelerates the development of atherosclerosis, decreases muscle function, and increases the risk of amputation or death in patients with peripheral artery disease (PAD). However, the mechanisms underlying this pathobiology are ill-defined. Recent work has indicated that tryptophan-derived uremic solutes, which are ligands for AHR (aryl hydrocarbon receptor), are associated with limb amputation in PAD. Herein, we examined the role of AHR activation in the myopathy of PAD and CKD. METHODS: AHR-related gene expression was evaluated in skeletal muscle obtained from mice and human PAD patients with and without CKD. AHRmKO (skeletal muscle-specific AHR knockout) mice with and without CKD were subjected to femoral artery ligation, and a battery of assessments were performed to evaluate vascular, muscle, and mitochondrial health. Single-nuclei RNA sequencing was performed to explore intercellular communication. Expression of the constitutively active AHR was used to isolate the role of AHR in mice without CKD. RESULTS: PAD patients and mice with CKD displayed significantly higher mRNA expression of classical AHR-dependent genes (Cyp1a1, Cyp1b1, and Aldh3a1) when compared with either muscle from the PAD condition with normal renal function (P<0.05 for all 3 genes) or nonischemic controls. AHRmKO significantly improved limb perfusion recovery and arteriogenesis, preserved vasculogenic paracrine signaling from myofibers, increased muscle mass and strength, as well as enhanced mitochondrial function in an experimental model of PAD/CKD. Moreover, viral-mediated skeletal muscle-specific expression of a constitutively active AHR in mice with normal kidney function exacerbated the ischemic myopathy evidenced by smaller muscle masses, reduced contractile function, histopathology, altered vasculogenic signaling, and lower mitochondrial respiratory function. CONCLUSIONS: These findings establish AHR activation in muscle as a pivotal regulator of the ischemic limb pathology in CKD. Further, the totality of the results provides support for testing of clinical interventions that diminish AHR signaling in these conditions.

PRMT5 links lipid metabolism to contractile function of skeletal muscles.

Skeletal muscle plays a key role in systemic energy homeostasis besides its contractile function, but what links these functions is poorly defined. Protein Arginine Methyl Transferase 5 (PRMT5) is a well-known oncoprotein but also expressed in healthy tissues with unclear physiological functions. As adult muscles express high levels of Prmt5, we generated skeletal muscle-specific Prmt5 knockout (Prmt5MKO ) mice. We observe reduced muscle mass, oxidative capacity, force production, and exercise performance in Prmt5MKO mice. The motor deficiency is associated with scarce lipid droplets in myofibers due to defects in lipid biosynthesis and accelerated degradation. Specifically, PRMT5 deletion reduces dimethylation and stability of Sterol Regulatory Element-Binding Transcription Factor 1a (SREBP1a), a master regulator of de novo lipogenesis. Moreover, Prmt5MKO impairs the repressive H4R3 symmetric dimethylation at the Pnpla2 promoter, elevating the level of its encoded protein ATGL, the rate-limiting enzyme catalyzing lipolysis. Accordingly, skeletal muscle-specific double knockout of Pnpla2 and Prmt5 normalizes muscle mass and function. Together, our findings delineate a physiological function of PRMT5 in linking lipid metabolism to contractile function of myofibers.

Genome-wide detection of enhancer-hijacking events from chromatin interaction data in rearranged genomes.

Recent efforts have shown that structural variations (SVs) can disrupt three-dimensional genome organization and induce enhancer hijacking, yet no computational tools exist to identify such events from chromatin interaction data. Here, we develop NeoLoopFinder, a computational framework to identify the chromatin interactions induced by SVs, including interchromosomal translocations, large deletions and inversions. Our framework can automatically resolve complex SVs, reconstruct local Hi-C maps surrounding the breakpoints, normalize copy number variation and allele effects and predict chromatin loops induced by SVs. We applied NeoLoopFinder in Hi-C data from 50 cancer cell lines and primary tumors and identified tens of recurrent genes associated with enhancer hijacking. To experimentally validate NeoLoopFinder, we deleted the hijacked enhancers in prostate adenocarcinoma cells using CRISPR-Cas9, which significantly reduced expression of the target oncogene. In summary, NeoLoopFinder enables identification of critical oncogenic regulatory elements that can potentially reveal therapeutic targets.

HiCLift: a fast and efficient tool for converting chromatin interaction data between genome assemblies.

MOTIVATION: With the continuous effort to improve the quality of human reference genome and the generation of more and more personal genomes, the conversion of genomic coordinates between genome assemblies is critical in many integrative and comparative studies. While tools have been developed for such task for linear genome signals such as ChIP-Seq, no tool exists to convert genome assemblies for chromatin interaction data, despite the importance of three-dimensional genome organization in gene regulation and disease. RESULTS: Here, we present HiCLift, a fast and efficient tool that can convert the genomic coordinates of chromatin contacts such as Hi-C and Micro-C from one assembly to another, including the latest T2T-CHM13 genome. Comparing with the strategy of directly remapping raw reads to a different genome, HiCLift runs on average 42 times faster (hours vs. days), while outputs nearly identical contact matrices. More importantly, as HiCLift does not need to remap the raw reads, it can directly convert human patient sample data, where the raw sequencing reads are sometimes hard to acquire or not available. AVAILABILITY AND IMPLEMENTATION: HiCLift is publicly available at https://github.com/XiaoTaoWang/HiCLift.

Loss of FoxO1 activates an alternate mechanism of mitochondrial quality control for healthy adipose browning.

Browning of white adipose tissue is hallmarked by increased mitochondrial density and metabolic improvements. However, it remains largely unknown how mitochondrial turnover and quality control are regulated during adipose browning. In the present study, we found that mice lacking adipocyte FoxO1, a transcription factor that regulates autophagy, adopted an alternate mechanism of mitophagy to maintain mitochondrial turnover and quality control during adipose browning. Post-developmental deletion of adipocyte FoxO1 (adO1KO) suppressed Bnip3 but activated Fundc1/Drp1/OPA1 cascade, concurrent with up-regulation of Atg7 and CTSL. In addition, mitochondrial biogenesis was stimulated via the Pgc1α/Tfam pathway in adO1KO mice. These changes were associated with enhanced mitochondrial homeostasis and metabolic health (e.g., improved glucose tolerance and insulin sensitivity). By contrast, silencing Fundc1 or Pgc1α reversed the changes induced by silencing FoxO1, which impaired mitochondrial quality control and function. Ablation of Atg7 suppressed mitochondrial turnover and function, causing metabolic disorder (e.g., impaired glucose tolerance and insulin sensitivity), regardless of elevated markers of adipose browning. Consistently, suppression of autophagy via CTSL by high-fat diet was associated with a reversal of adO1KO-induced benefits. Our data reveal a unique role of FoxO1 in coordinating mitophagy receptors (Bnip3 and Fundc1) for a fine-tuned mitochondrial turnover and quality control, underscoring autophagic clearance of mitochondria as a prerequisite for healthy browning of adipose tissue.

Combined untargeted and targeted lipidomics approaches reveal potential biomarkers in type 2 diabetes mellitus cynomolgus monkeys.

BACKGROUND: The significantly increasing incidence of type 2 diabetes mellitus (T2DM) over the last few decades triggers the demands of T2DM animal models to explore the pathogenesis, prevention, and therapy of the disease. The altered lipid metabolism may play an important role in the pathogenesis and progression of T2DM. However, the characterization of molecular lipid species in fasting serum related to T2DM cynomolgus monkeys is still underrecognized. METHODS: Untargeted and targeted LC-mass spectrometry (MS)/MS-based lipidomics approaches were applied to characterize and compare the fasting serum lipidomic profiles of T2DM cynomolgus monkeys and the healthy controls. RESULTS: Multivariate analysis revealed that 196 and 64 lipid molecules differentially expressed in serum samples using untargeted and targeted lipidomics as the comparison between the disease group and healthy group, respectively. Furthermore, the comparative analysis of differential serum lipid metabolites obtained by untargeted and targeted lipidomics approaches, four common serum lipid species (phosphatidylcholine [18:0_22:4], lysophosphatidylcholine [14:0], phosphatidylethanolamine [PE] [16:1_18:2], and PE [18:0_22:4]) were identified as potential biomarkers and all of which were found to be downregulated. By analyzing the metabolic pathway, glycerophospholipid metabolism was associated with the pathogenesis of T2DM cynomolgus monkeys. CONCLUSION: The study found that four downregulated serum lipid species could serve as novel potential biomarkers of T2DM cynomolgus monkeys. Glycerophospholipid metabolism was filtered out as the potential therapeutic target pathway of T2DM progression. Our results showed that the identified biomarkers may offer a novel tool for tracking disease progression and response to therapeutic interventions.

Genetically engineering of Saccharomyces cerevisiae for enhanced oral delivery vaccine vehicle.

As a series of our previous studies reported, recombinant yeast can be the oral vaccines to deliver designed protein and DNA, as well as functional shRNA, into dendritic cells (DCs) in mice for specific immune regulation. Here, we report the further optimization of oral yeast-based vaccine from two aspects (yeast characteristics and recombinant DNA constitution) to improve the effect of immune regulation. After screening four genes in negative regulation of glucan synthesis in yeast (MNN9, GUP1, PBS2 and EXG1), this research combined HDR-based genome editing technology with Cre-loxP technology to acquire 15 gene-knockout strains without drug resistance-gene to exclude biosafety risks; afterward, oral feeding experiments were performed on the mice using 15 oral recombinant yeast-based vaccines constructed by the gene-knockout strains harboring pCMV-MSTN plasmid to screen the target strain with more effective inducing mstn-specific antibody which in turn increasing weight gain effect. And subsequently based on the selected gene-knockout strain, the recombinant DNA in the oral recombinant yeast-based vaccine is optimized via a combination of protein fusion expression (OVA-MSTN) and interfering RNA technology (shRNA-IL21), comparison in terms of both weight gain effect and antibody titer revealed that the selected gene-knockout strain (GUP1ΔEXG1Δ) combined with specific recombinant DNA (pCMV-OVA-MSTN-shIL2) had a better effect of the vaccine. This study provides a useful reference to the subsequent construction of a more efficient oral recombinant yeast-based vaccine in the food and pharmaceutical industry.