Longitudinal cytokine and multi-modal health data of an extremely severe ME/CFS patient with HSD reveals insights into immunopathology, and disease severity.

Introduction: Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) presents substantial challenges in patient care due to its intricate multisystem nature, comorbidities, and global prevalence. The heterogeneity among patient populations, coupled with the absence of FDA-approved diagnostics and therapeutics, further complicates research into disease etiology and patient managment. Integrating longitudinal multi-omics data with clinical, health,textual, pharmaceutical, and nutraceutical data offers a promising avenue to address these complexities, aiding in the identification of underlying causes and providing insights into effective therapeutics and diagnostic strategies. Methods: This study focused on an exceptionally severe ME/CFS patient with hypermobility spectrum disorder (HSD) during a period of marginal symptom improvements. Longitudinal cytokine profiling was conducted alongside the collection of extensive multi-modal health data to explore the dynamic nature of symptoms, severity, triggers, and modifying factors. Additionally, an updated severity assessment platform and two applications, ME-CFSTrackerApp and LexiTime, were introduced to facilitate real-time symptom tracking and enhance patient-physician/researcher communication, and evaluate response to medical intervention. Results: Longitudinal cytokine profiling revealed the significance of Th2-type cytokines and highlighted synergistic activities between mast cells and eosinophils, skewing Th1 toward Th2 immune responses in ME/CFS pathogenesis, particularly in cognitive impairment and sensorial intolerance. This suggests a potentially shared underlying mechanism with major ME/CFS comorbidities such as HSD, Mast cell activation syndrome, postural orthostatic tachycardia syndrome (POTS), and small fiber neuropathy. Additionally, the data identified potential roles of BCL6 and TP53 pathways in ME/CFS etiology and emphasized the importance of investigating adverse reactions to medication and supplements and drug interactions in ME/CFS severity and progression. Discussion: Our study advocates for the integration of longitudinal multi-omics with multi-modal health data and artificial intelligence (AI) techniques to better understand ME/CFS and its major comorbidities. These findings highlight the significance of dysregulated Th2-type cytokines in patient stratification and precision medicine strategies. Additionally, our results suggest exploring the use of low-dose drugs with partial agonist activity as a potential avenue for ME/CFS treatment. This comprehensive approach emphasizes the importance of adopting a patient-centered care approach to improve ME/CFS healthcare management, disease severity assessment, and personalized medicine. Overall, these findings contribute to our understanding of ME/CFS and offer avenues for future research and clinical practice.

Phenotypic characteristics of peripheral immune cells of Myalgic encephalomyelitis/chronic fatigue syndrome via transmission electron microscopy: A pilot study.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a complex chronic multi-systemic disease characterized by extreme fatigue that is not improved by rest, and worsens after exertion, whether physical or mental. Previous studies have shown ME/CFS-associated alterations in the immune system and mitochondria. We used transmission electron microscopy (TEM) to investigate the morphology and ultrastructure of unstimulated and stimulated ME/CFS immune cells and their intracellular organelles, including mitochondria. PBMCs from four participants were studied: a pair of identical twins discordant for moderate ME/CFS, as well as two age- and gender- matched unrelated subjects-one with an extremely severe form of ME/CFS and the other healthy. TEM analysis of CD3/CD28-stimulated T cells suggested a significant increase in the levels of apoptotic and necrotic cell death in T cells from ME/CFS patients (over 2-fold). Stimulated Tcells of ME/CFS patients also had higher numbers of swollen mitochondria. We also found a large increase in intracellular giant lipid droplet-like organelles in the stimulated PBMCs from the extremely severe ME/CFS patient potentially indicative of a lipid storage disorder. Lastly, we observed a slight increase in platelet aggregation in stimulated cells, suggestive of a possible role of platelet activity in ME/CFS pathophysiology and disease severity. These results indicate extensive morphological alterations in the cellular and mitochondrial phenotypes of ME/CFS patients' immune cells and suggest new insights into ME/CFS biology.

Candidate variants in TUB are associated with familial tremor.

Essential tremor (ET) is the most common adult-onset movement disorder. In the present study, we performed whole exome sequencing of a large ET-affected family (10 affected and 6 un-affected family members) and identified a TUB p.V431I variant (rs75594955) segregating in a manner consistent with autosomal-dominant inheritance. Subsequent targeted re-sequencing of TUB in 820 unrelated individuals with sporadic ET and 630 controls revealed significant enrichment of rare nonsynonymous TUB variants (e.g. rs75594955: p.V431I, rs1241709665: p.Ile20Phe, rs55648406: p.Arg49Gln) in the ET cohort (SKAT-O test p-value = 6.20e-08). TUB encodes a transcription factor predominantly expressed in neuronal cells and has been previously implicated in obesity. ChIP-seq analyses of the TUB transcription factor across different regions of the mouse brain revealed that TUB regulates the pathways responsible for neurotransmitter production as well thyroid hormone signaling. Together, these results support the association of rare variants in TUB with ET.

Heterogeneity in old fibroblasts is linked to variability in reprogramming and wound healing.

Age-associated chronic inflammation (inflammageing) is a central hallmark of ageing1, but its influence on specific cells remains largely unknown. Fibroblasts are present in most tissues and contribute to wound healing2,3. They are also the most widely used cell type for reprogramming to induced pluripotent stem (iPS) cells, a process that has implications for regenerative medicine and rejuvenation strategies4. Here we show that fibroblast cultures from old mice secrete inflammatory cytokines and exhibit increased variability in the efficiency of iPS cell reprogramming between mice. Variability between individuals is emerging as a feature of old age5-8, but the underlying mechanisms remain unknown. To identify drivers of this variability, we performed multi-omics profiling of fibroblast cultures from young and old mice that have different reprogramming efficiencies. This approach revealed that fibroblast cultures from old mice contain 'activated fibroblasts' that secrete inflammatory cytokines, and that the proportion of activated fibroblasts in a culture correlates with the reprogramming efficiency of that culture. Experiments in which conditioned medium was swapped between cultures showed that extrinsic factors secreted by activated fibroblasts underlie part of the variability between mice in reprogramming efficiency, and we have identified inflammatory cytokines, including TNF, as key contributors. Notably, old mice also exhibited variability in wound healing rate in vivo. Single-cell RNA-sequencing analysis identified distinct subpopulations of fibroblasts with different cytokine expression and signalling in the wounds of old mice with slow versus fast healing rates. Hence, a shift in fibroblast composition, and the ratio of inflammatory cytokines that they secrete, may drive the variability between mice in reprogramming in vitro and influence wound healing rate in vivo. This variability may reflect distinct stochastic ageing trajectories between individuals, and could help in developing personalized strategies to improve iPS cell generation and wound healing in elderly individuals.

Microfluidic isoform sequencing shows widespread splicing coordination in the human transcriptome.

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k-200k partitions of 10-200 molecules at a time, enabling analysis of 10-100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT Coordination of internal exons is found to be extensive among protein coding genes: 23.5%-59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein-protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.

Cell Type-Specific Chromatin Signatures Underline Regulatory DNA Elements in Human Induced Pluripotent Stem Cells and Somatic Cells.

RATIONALE: Regulatory DNA elements in the human genome play important roles in determining the transcriptional abundance and spatiotemporal gene expression during embryonic heart development and somatic cell reprogramming. It is not well known how chromatin marks in regulatory DNA elements are modulated to establish cell type-specific gene expression in the human heart. OBJECTIVE: We aimed to decipher the cell type-specific epigenetic signatures in regulatory DNA elements and how they modulate heart-specific gene expression. METHODS AND RESULTS: We profiled genome-wide transcriptional activity and a variety of epigenetic marks in the regulatory DNA elements using massive RNA-seq (n=12) and ChIP-seq (chromatin immunoprecipitation combined with high-throughput sequencing; n=84) in human endothelial cells (CD31+CD144+), cardiac progenitor cells (Sca-1+), fibroblasts (DDR2+), and their respective induced pluripotent stem cells. We uncovered 2 classes of regulatory DNA elements: class I was identified with ubiquitous enhancer (H3K4me1) and promoter (H3K4me3) marks in all cell types, whereas class II was enriched with H3K4me1 and H3K4me3 in a cell type-specific manner. Both class I and class II regulatory elements exhibited stimulatory roles in nearby gene expression in a given cell type. However, class I promoters displayed more dominant regulatory effects on transcriptional abundance regardless of distal enhancers. Transcription factor network analysis indicated that human induced pluripotent stem cells and somatic cells from the heart selected their preferential regulatory elements to maintain cell type-specific gene expression. In addition, we validated the function of these enhancer elements in transgenic mouse embryos and human cells and identified a few enhancers that could possibly regulate the cardiac-specific gene expression. CONCLUSIONS: Given that a large number of genetic variants associated with human diseases are located in regulatory DNA elements, our study provides valuable resources for deciphering the epigenetic modulation of regulatory DNA elements that fine-tune spatiotemporal gene expression in human cardiac development and diseases.

The lncRNA BORG Drives Breast Cancer Metastasis and Disease Recurrence.

Long noncoding RNAs (lncRNAs) have emerged as potent regulators of breast cancer development and progression, including the metastatic spread of disease. Through in silico and biological analyses, we identified a novel lncRNA, BMP/OP-Responsive Gene (BORG), whose expression directly correlates with aggressive breast cancer phenotypes, as well as with metastatic competence and disease recurrence in multiple clinical cohorts. Mechanistically, BORG elicits the metastatic outgrowth of latent breast cancer cells by promoting the localization and transcriptional repressive activity of TRIM28, which binds BORG and induces substantial alterations in carcinoma proliferation and survival. Moreover, inhibiting BORG expression in metastatic breast cancer cells impedes their metastatic colonization of the lungs of mice, implying that BORG acts as a novel driver of the genetic and epigenetic alterations that underlie the acquisition of metastatic and recurrent phenotypes by breast cancer cells.

Association of AHSG with alopecia and mental retardation (APMR) syndrome.

Alopecia with mental retardation syndrome (APMR) is a very rare autosomal recessive condition that is associated with total or partial absence of hair from the scalp and other parts of the body as well as variable intellectual disability. Here we present whole-exome sequencing results of a large consanguineous family segregating APMR syndrome with seven affected family members. Our study revealed a novel predicted pathogenic, homozygous missense mutation in the AHSG (OMIM 138680) gene (AHSG: NM_001622:exon7:c.950G>A:p.Arg317His). The variant is predicted to affect a region of the protein required for protein processing and disrupts a phosphorylation motif. In addition, the altered protein migrates with an aberrant size relative to healthy individuals. Consistent with the phenotype, AHSG maps within APMR linkage region 1 (APMR 1) as reported before, and falls within runs of homozygosity (ROH). Previous families with APMR syndrome have been studied through linkage analyses and the linkage resolution did not allow pointing out to a single gene candidate. Our study is the first report to identify a homozygous missense mutation for APMR syndrome through whole-exome sequencing.

iPSC-derived cardiomyocytes reveal abnormal TGF-ß signalling in left ventricular non-compaction cardiomyopathy.

Left ventricular non-compaction (LVNC) is the third most prevalent cardiomyopathy in children and its pathogenesis has been associated with the developmental defect of the embryonic myocardium. We show that patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) generated from LVNC patients carrying a mutation in the cardiac transcription factor TBX20 recapitulate a key aspect of the pathological phenotype at the single-cell level and this was associated with perturbed transforming growth factor beta (TGF-ß) signalling. LVNC iPSC-CMs have decreased proliferative capacity due to abnormal activation of TGF-ß signalling. TBX20 regulates the expression of TGF-ß signalling modifiers including one known to be a genetic cause of LVNC, PRDM16, and genome editing of PRDM16 caused proliferation defects in iPSC-CMs. Inhibition of TGF-ß signalling and genome correction of the TBX20 mutation were sufficient to reverse the disease phenotype. Our study demonstrates that iPSC-CMs are a useful tool for the exploration of pathological mechanisms underlying poorly understood cardiomyopathies including LVNC.

Effects of cellular origin on differentiation of human induced pluripotent stem cell-derived endothelial cells.

Human induced pluripotent stem cells (iPSCs) can be derived from various types of somatic cells by transient overexpression of 4 Yamanaka factors (OCT4, SOX2, C-MYC, and KLF4). Patient-specific iPSC derivatives (e.g., neuronal, cardiac, hepatic, muscular, and endothelial cells [ECs]) hold great promise in drug discovery and regenerative medicine. In this study, we aimed to evaluate whether the cellular origin can affect the differentiation, in vivo behavior, and single-cell gene expression signatures of human iPSC-derived ECs. We derived human iPSCs from 3 types of somatic cells of the same individuals: fibroblasts (FB-iPSCs), ECs (EC-iPSCs), and cardiac progenitor cells (CPC-iPSCs). We then differentiated them into ECs by sequential administration of Activin, BMP4, bFGF, and VEGF. EC-iPSCs at early passage (10 < P < 20) showed higher EC differentiation propensity and gene expression of EC-specific markers (PECAM1 and NOS3) than FB-iPSCs and CPC-iPSCs. In vivo transplanted EC-iPSC-ECs were recovered with a higher percentage of CD31+ population and expressed higher EC-specific gene expression markers (PECAM1, KDR, and ICAM) as revealed by microfluidic single-cell quantitative PCR (qPCR). In vitro EC-iPSC-ECs maintained a higher CD31+ population than FB-iPSC-ECs and CPC-iPSC-ECs with long-term culturing and passaging. These results indicate that cellular origin may influence lineage differentiation propensity of human iPSCs; hence, the somatic memory carried by early passage iPSCs should be carefully considered before clinical translation.

IGF2BP3 Modulates the Interaction of Invasion-Associated Transcripts with RISC.

Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with malignancy, but its role(s) in pathogenesis remains enigmatic. We interrogated the IGF2BP3-RNA interaction network in pancreatic ductal adenocarcinoma (PDAC) cells. Using a combination of genome-wide approaches, we have identified 164 direct mRNA targets of IGF2BP3. These transcripts encode proteins enriched for functions such as cell migration, proliferation, and adhesion. Loss of IGF2BP3 reduced PDAC cell invasiveness and remodeled focal adhesion junctions. Individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) revealed significant overlap of IGF2BP3 and microRNA (miRNA) binding sites. IGF2BP3 promotes association of the RNA-induced silencing complex (RISC) with specific transcripts. Our results show that IGF2BP3 influences a malignancy-associated RNA regulon by modulating miRNA-mRNA interactions.

Synthetic long-read sequencing reveals intraspecies diversity in the human microbiome.

Identifying bacterial strains in metagenome and microbiome samples using computational analyses of short-read sequences remains a difficult problem. Here, we present an analysis of a human gut microbiome using TruSeq synthetic long reads combined with computational tools for metagenomic long-read assembly, variant calling and haplotyping (Nanoscope and Lens). Our analysis identifies 178 bacterial species, of which 51 were not found using shotgun reads alone. We recover bacterial contigs that comprise multiple operons, including 22 contigs of >1 Mbp. Furthermore, we observe extensive intraspecies variation within microbial strains in the form of haplotypes that span up to hundreds of Kbp. Incorporation of synthetic long-read sequencing technology with standard short-read approaches enables more precise and comprehensive analyses of metagenomic samples.

Comprehensive transcriptome analysis using synthetic long-read sequencing reveals molecular co-association of distant splicing events.

Alternative splicing shapes mammalian transcriptomes, with many RNA molecules undergoing multiple distant alternative splicing events. Comprehensive transcriptome analysis, including analysis of exon co-association in the same molecule, requires deep, long-read sequencing. Here we introduce an RNA sequencing method, synthetic long-read RNA sequencing (SLR-RNA-seq), in which small pools (≤1,000 molecules/pool, ≤1 molecule/gene for most genes) of full-length cDNAs are amplified, fragmented and short-read-sequenced. We demonstrate that these RNA sequences reconstructed from the short reads from each of the pools are mostly close to full length and contain few insertion and deletion errors. We report many previously undescribed isoforms (human brain: ∼13,800 affected genes, 14.5% of molecules; mouse brain ∼8,600 genes, 18% of molecules) and up to 165 human distant molecularly associated exon pairs (dMAPs) and distant molecularly and mutually exclusive pairs (dMEPs). Of 16 associated pairs detected in the mouse brain, 9 are conserved in human. Our results indicate conserved mechanisms that can produce distant but phased features on transcript and proteome isoforms.

Exome sequencing and genome-wide copy number variant mapping reveal novel associations with sensorineural hereditary hearing loss.

BACKGROUND: The genetic diversity of loci and mutations underlying hereditary hearing loss is an active area of investigation. To identify loci associated with predominantly non-syndromic sensorineural hearing loss, we performed exome sequencing of families and of single probands, as well as copy number variation (CNV) mapping in a case-control cohort. RESULTS: Analysis of three distinct families revealed several candidate loci in two families and a single strong candidate gene, MYH7B, for hearing loss in one family. MYH7B encodes a Type II myosin, consistent with a role for cytoskeletal proteins in hearing. High-resolution genome-wide CNV analysis of 150 cases and 157 controls revealed deletions in genes known to be involved in hearing (e.g. GJB6, OTOA, and STRC, encoding connexin 30, otoancorin, and stereocilin, respectively), supporting CNV contributions to hearing loss phenotypes. Additionally, a novel region on chromosome 16 containing part of the PDXDC1 gene was found to be frequently deleted in hearing loss patients (OR=3.91, 95% CI: 1.62-9.40, p=1.45×10(-7)). CONCLUSIONS: We conclude that many known as well as novel loci and distinct types of mutations not typically tested in clinical settings can contribute to the etiology of hearing loss. Our study also demonstrates the challenges of exome sequencing and genome-wide CNV mapping for direct clinical application, and illustrates the need for functional and clinical follow-up as well as curated open-access databases.

Genome-wide map of regulatory interactions in the human genome.

Increasing evidence suggests that interactions between regulatory genomic elements play an important role in regulating gene expression. We generated a genome-wide interaction map of regulatory elements in human cells (ENCODE tier 1 cells, K562, GM12878) using Chromatin Interaction Analysis by Paired-End Tag sequencing (ChIA-PET) experiments targeting six broadly distributed factors. Bound regions covered 80% of DNase I hypersensitive sites including 99.7% of TSS and 98% of enhancers. Correlating this map with ChIP-seq and RNA-seq data sets revealed cohesin, CTCF, and ZNF143 as key components of three-dimensional chromatin structure and revealed how the distal chromatin state affects gene transcription. Comparison of interactions between cell types revealed that enhancer-promoter interactions were highly cell-type-specific. Construction and comparison of distal and proximal regulatory networks revealed stark differences in structure and biological function. Proximal binding events are enriched at genes with housekeeping functions, while distal binding events interact with genes involved in dynamic biological processes including response to stimulus. This study reveals new mechanistic and functional insights into regulatory region organization in the nucleus.

A novel RNA motif mediates the strict nuclear localization of a long noncoding RNA.

The ubiquitous presence of long noncoding RNAs (lncRNAs) in eukaryotes points to the importance of understanding how their sequences impact function. As many lncRNAs regulate nuclear events and thus must localize to nuclei, we analyzed the sequence requirements for nuclear localization in an intergenic lncRNA named BORG (BMP2-OP1-responsive gene), which is both spliced and polyadenylated but is strictly localized in nuclei. Subcellular localization of BORG was not dependent on the context or level of its expression or decay but rather depended on the sequence of the mature, spliced transcript. Mutational analyses indicated that nuclear localization of BORG was mediated through a novel RNA motif consisting of the pentamer sequence AGCCC with sequence restrictions at positions -8 (T or A) and -3 (G or C) relative to the first nucleotide of the pentamer. Mutation of the motif to a scrambled sequence resulted in complete loss of nuclear localization, while addition of even a single copy of the motif to a cytoplasmically localized RNA was sufficient to impart nuclear localization. Further, the presence of this motif in other cellular RNAs showed a direct correlation with nuclear localization, suggesting that the motif may act as a general nuclear localization signal for cellular RNAs.