Pairs Of Sequences Research Articles

Resolving the chromatin impact of mosaic variants with targeted Fiber-seq.

Accurately quantifying the functional consequences of noncoding mosaic variants requires the pairing of DNA sequences with both accessible and closed chromatin architectures along individual DNA molecules-a pairing that cannot be achieved using traditional fragmentation-based chromatin assays. We demonstrate that targeted single-molecule chromatin fiber sequencing (Fiber-seq) achieves this, permitting single-molecule, long-read genomic, and epigenomic profiling across targeted >100 kb loci with ∼10-fold enrichment over untargeted sequencing. Targeted Fiber-seq reveals that pathogenic expansions of the DMPK CTG repeat that underlie Myotonic Dystrophy 1 are characterized by somatic instability and disruption of multiple nearby regulatory elements, both of which are repeat length-dependent. Furthermore, we reveal that therapeutic adenine base editing of the segmentally duplicated γ-globin (HBG1/HBG2) promoters in primary human hematopoietic cells induced toward an erythroblast lineage increases the accessibility of the HBG1 promoter as well as neighboring regulatory elements. Overall, we find that these non-protein coding mosaic variants can have complex impacts on chromatin architectures, including extending beyond the regulatory element harboring the variant.

Large scale paired antibody language models.

Antibodies are proteins produced by the immune system that can identify and neutralise a wide variety of antigens with high specificity and affinity, and constitute the most successful class of biotherapeutics. With the advent of next-generation sequencing, billions of antibody sequences have been collected in recent years, though their application in the design of better therapeutics has been constrained by the sheer volume and complexity of the data. To address this challenge, we present IgBert and IgT5, the best performing antibody-specific language models developed to date which can consistently handle both paired and unpaired variable region sequences as input. These models are trained comprehensively using the more than two billion unpaired sequences and two million paired sequences of light and heavy chains present in the Observed Antibody Space dataset. We show that our models outperform existing antibody and protein language models on a diverse range of design and regression tasks relevant to antibody engineering. This advancement marks a significant leap forward in leveraging machine learning, large scale data sets and high-performance computing for enhancing antibody design for therapeutic development.

A new species of Myxobolus (Cnidaria: Myxosporea: Myxobolidae) from the gills of the green sunfish, Lepomis cyanellus (Perciformes: Centrarchidae), from the Ouachita River Drainage of Western Arkansas.

Forty-two species of myxozoans, including 31 species of Myxobolus Bütschli, 1882, have been described from centrarchid fishes. One species, the green sunfish (Lepomis cyanellus), has been reported to host at least three species of Myxobolus. Between March 2023 and June 2024, 42 L. cyanellus were collected from watersheds in Montgomery and Polk counties, Arkansas, and their gills, gallbladders, urinary bladders, fins, integument, other major organs, and musculature were examined for myxozoans. Two (5%) L. cyanellus from Polk County were found to harbor a new species of Myxobolus infecting the gill lamellae. A qualitative and quantitative morphological description was based on formalin-fixed preserved plasmodia and myxospores. Plasmodia of Myxobolus polkensis n. sp. are 133 µm long × 123 µm wide, and myxospores are 19.8 µm long × 6.5 µm wide, with two narrowly pyriform unequal polar capsules, one usually longer (9.2 µm) than the other (8.5 µm). Molecular data consisted of a 2,025 base pair sequence of the partial small subunit rRNA gene (SSU). Variably sized, polysporic plasmodia were randomly distributed throughout the gill lamellae in intralamellar locations. Large plasmodia displaced adjacent lamellae inducing mild epithelial proliferation but minimal inflammatory changes. Phylogenetic analysis revealed that M. polkensis n. sp. is a member of a clade of myxozoan species that predominately infect centrarchids from North America. This is the fourth report of a Myxobolus from L. cyanellus but the first report of a species infecting the gill lamellae. This article was registered in the Official Register of Zoological Nomenclature (ZooBank) as urn:lsid:zoobank.org:pub:E9BA08E0-0A7E-487B-80AE-C5C723DBB661.

Gene and Its Promoter Cloning, and Functional Validation of JmSOC1 Revealed Its Role in Promoting Early Flowering and the Interaction with the JmSVP Protein

Juglans mandshurica, a notable woody oil tree species, possesses both fruit and timber value. However, the complete heterodichogamous flowering mechanism in this species remains elusive. SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) is a crucial regulator of flower bud development in Arabidopsis thaliana. In this study, we cloned the coding DNA sequence (CDS) of the JmSOC1 gene, revealing a 705 base pair (bp) sequence that encodes a protein of 234 amino acids. The JmSOC1 protein contains a highly conserved MADS-box domain, indicating its role as a transcription factor, and is predominantly localized in the nucleus. The JmSOC1 gene expressed the highest in flower buds. The peak expression level of JmSOC1 during the physiological differentiation phase occurred earlier in male flower buds of protandry (MPD) on April 10th compared to female flower buds of protandry (FPD) on April 14th; similarly, the peak expression in female flower buds of protogyny (FPG) on April 2nd preceded that in male flower buds of protogyny (MPG) on April 14th. This may be the primary reason for the earlier differentiation of the male flowers in protandry individuals and the female flowers in protogyny individuals. Overexpression of JmSOC1 in wild-type A. thaliana resulted in earlier flowering, accompanied by an upregulation of key flowering-related genes such as LEAFY (LFY), APETALA1 (AP1), and FLOWERING LOCUS T (FT). To further explore the function of JmSOC1, a 782 bp promoter sequence of JmSOC1 gene was cloned, which has been verified to have promoter activity by GUS staining. Furthermore, the interaction between the JmSOC1 gene promoter and its upstream regulatory protein JmSVP was verified using a yeast one-hybrid. These results offer valuable insights into the molecular mechanisms underpinning the promotion of early flowering in J. mandshurica and hold promise for laying a theoretical foundation for the flowering regulation network of this species.

Doppler resilient MIMO radar waveform design from zero correlation zone complementary sequence pairs

The aperiodic auto-correlation sum of Golay complementary pair (GCP) is zero at all non-zero delay, which makes Golay complementary pair have good performance in radar systems. In the literature, GCPs have been widely used in single user (or single antenna) radar system in which waveforms with only good auto-correlation are considered. While in the multi-user (or MIMO) scenario, radar waveforms with both good auto-correlation and cross-correlation are desired to eliminate mutual interference in target detection. To address this issue, constructions of waveforms with both good auto-ambiguity property and cross-ambiguity property are proposed based on zero correlation zone complementary pair (ZCP) sets. The simulation results show that the proposed waveform set perform as well as Golay complementary waveforms in target detection.

Why Symptoms Linger in Quiescent Crohn's Disease: Investigating the Impact of Sulfidogenic Microbes and Sulfur Metabolic Pathways.

Even in the absence of inflammation, persistent symptoms in patients with Crohn's disease (CD) are prevalent and worsen quality of life. We previously demonstrated enrichment in sulfidogenic microbes in quiescent Crohn's disease patients with (qCD + S) vs without persistent GI symptoms (qCD-S). Thus, we hypothesized that sulfur metabolic pathways would be enriched in stool while differentially abundant microbes would be associated with important sulfur metabolic pathways in qCD + S. We performed a multicenter observational study nested within SPARC IBD. Quiescent inflammation was defined by fecal calprotectin level < 150 mcg/g. Persistent symptoms were defined by CD-PRO2. Active CD (aCD) and non-IBD diarrhea-predominant irritable bowel syndrome (IBS-D) were included as controls. Thirty-nine patients with qCD + S, 274 qCD-S, 21 aCD, and 40 IBS-D underwent paired shotgun metagenomic sequencing and untargeted metabolomic profiling. The fecal metabolome in qCD + S was significantly different relative to qCD-S and IBS-D but not aCD. Patients with qCD + S were enriched in sulfur-containing amino acid pathways, including cysteine and methionine, as well as serine, glycine, and threonine. Glutathione and nicotinate/nicotinamide pathways were also enriched in qCD + S relative to qCD-S, suggestive of mitochondrial dysfunction, a downstream target of H2S signaling. Multi-omic integration demonstrated that enriched microbes in qCD + S were associated with important sulfur metabolic pathways. Bacterial sulfur metabolic genes, including CTH, isfD, sarD, and asrC, were dysregulated in qCD + S. Finally, sulfur metabolites with and without sulfidogenic microbes showed good accuracy in predicting the presence of qCD + S. Microbial-derived sulfur pathways and downstream mitochondrial function are perturbed in qCD + S, which implicate H2S signaling in the pathogenesis of this condition. Future studies will determine whether targeting H2S pathways results in improved quality of life in qCD + S.

Abstract IA005: Tumor cell-based liquid biopsy to characterize prostate cancer

Abstract Circulating Tumor Cells (CTCs) hold the promise of comprehensive yet noninvasive molecular characterization of cancer cells during the course of tumor evolution and progression. To dissect the epigenetic profile of prostate CTCs, we undertook paired DNA methylation and RNA sequencing of single CTCs isolated from blood specimens of patients with prostate cancer. Across the genome, our study identified 40 core Partial Methylation Domains (PMDs) that arise early in tumorigenesis and are shared by all individual CTCs across multiple patients. Key among these is a single genomic locus harboring multiple interferon-inducible genes along with all the CD1 gene family members, implicated in innate immunity. Based on mouse models, we propose that early biallelic silencing of multiple genes residing within a single PMD may contribute to the escape from immune surveillance in prostate and potentially other cancers. Studies of CTC heterogeneity within individual patients are limited by the very small number of cells that can be isolated from standard 10ml blood volumes. To enable such analyses, we applied a new high throughput microfluidic platform, capable of processing an entire leukapheresis product, interrogating 2-3 L of blood volume equivalent. Applying this technology to patients with metastatic prostate cancer, we successfully enriched thousands of CTCs from individual cases, enabling the characterization of intra-patient heterogeneity at the level of cell morphology and marker expression, chromosome copy number variation and exome-wide genotyping, and transcriptionally-defined tumor cell subpopulations. Such high-volume microfluidic enrichment of CTCs constitutes a new dimension in liquid biopsies. Citation Format: Daniel A Haber, Shyamala Maheswaran, Hongshan Guo, Joanna Vuille, Shih-Bo Huang, Ben S Wittner, Lecia V Sequist, Richard J Lee, Patricia AR Brunker, David T Miyamoto, Avanish Mishra, Mehmet Toner. Tumor cell-based liquid biopsy to characterize prostate cancer [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr IA005.

EPCO-12. MENINGIOMA COPY NUMBER ANALYSIS IDENTIFIES OPTIMIZED SIZE THRESHOLDS AND CO-OCCURRENCE MODELS FOR INDIVIDUALIZED RISK-STRATIFICATION

Abstract Arm-level copy number alterations (CNAs) are common in meningiomas and are useful for risk-stratification, but there is no consensus regarding the optimal size threshold to define CNAs for prognostic models. Moreover, the prognostic importance of co-occurring CNA pairs and overall CNA burden in meningiomas remains incompletely understood. To address this, CNAs were analyzed in 691 retrospective meningiomas with clinical data and DNA methylation profiles from 3 international institutions. Arm-level CNAs were defined from DNA methylation profiles using iterative size thresholds ranging from 1-99% of each chromosome arm. The area under the curve (AUC) for 5-year local freedom from recurrence (LFFR) or overall survival (OS) was calculated using CNAs across size thresholds, and AUC standard deviation across thresholds was used as a measure of size-dependence. LASSO and Elastic Net Cox regressions were used to generate multi-CNA models for LFFR or OS and to identify prognostic co-occurring CNA pairs. These analyses identified 13 size-dependent prognostic CNAs in meningiomas, including loss of 1p, 3q, 4p, 4q, 6p, 6q, 9p, 10q, 12q, 14q,18q, and 22q, and gain of 1q. Regression models restricted to size-dependent CNAs improved risk stratification within strata of previously published models (Integrated score, Integrated grade), and identified prognostic CNAs not included in previous models, such as 1q gain, 7p loss, and 12q loss. Ontology analysis of genes in focal prognostic CNAs identified dysregulated STAT signaling and decreased expression of interferon-related genes from matched paired RNA sequencing. Models of size-dependent, prognostic, co-occurring CNA pairs identified 1p/22q codeletion, 1q gain/22q loss, and 9p/18q codeletion as significant predictors of LFFR. Increased CNA burden correlated with higher tumor grade and was significantly associated with worse LFFR and OS. In summary, chromosome-specific size thresholds identify prognostic arm-level CNAs and co-occurring CNA pairs in meningiomas. Thus, granular copy number analysis may improve risk stratification models for meningioma.

CTNI-15. A PERIOPERATIVE STUDY OF SAFUSIDENIB IN PATIENTS WITH IDH1 MUTATED GLIOMA

Abstract BACKGROUND IDH mutations cause aberrant accumulation of (R)-2-hydroxyglutarate (2-HG), an oncometabolite that promotes tumorigenesis through abnormal hypermethylation of histones and DNA as well as suppression of normal cellular differentiation. IDH inhibition significantly improves progression free survival in patients with IDH mutant WHO grade 2 glioma, however most patients progress, and there is a lack of understanding on mechanistic reasons for failure. METHODS We tested Safusidenib, an oral, blood brain barrier penetrant IDH1 R132X enzyme inhibitor in a single arm, open label perioperative study. Patients with IDH1 mutated glioma, treatment naïve to radiation or chemotherapy, received Safusidenib (250 mg BID) for four weeks following surgical biopsy and prior to definitive resection (Part A). After resection patients continued Safusidenib until toxicity or progression (Part B). Translational endpoints include paired whole genome transcriptome sequencing (WGTS), single nuclei RNA sequencing (snRNAseq), spatial transcriptomics and spatial metabolomics as well as longitudinal sampling of blood and CSF for circulating tumour DNA (ctDNA). Here we are reporting Part A of the study. RESULTS As of May/27/2024 10 patients (mean age 37, 40% female) have been enrolled: seven with astrocytoma, three with oligodendroglioma; six had prior surgical resection. Pharmacokinetic sampling showed a mean tumor concentration 2457 ng/mL (range 957.0 – 4810) with tumor:plasma ratio 0.33. Consistent with mechanism, tumor 2-HG reduced 86% from biopsy (mean 75.0μM, range 7.8-203) to surgery (10.3μM, range 2.3 – 30.8), as well as identification of novel biomarkers of response. Paired tumour snRNAseq reveals transcriptional shift with changing cell-to-cell communication, and increased immune infiltration post-treatment, validated by spatial transcriptomics. Additionally, we find reduced heterogeneity consistent with a more differentiated tumour state. Treatment was well tolerated with ongoing follow-up. CONCLUSIONS This innovative trial facilitates intra-patient comparisons and confirmation of the on-target effects of Safusidenib, enhancing understanding of the mechanisms of response and resistance. The approach serves as a proof of concept for advancing drug development in glioma through perioperative trials.

TMIC-33. MICROENVIRONMENT-DRIVEN CHANGES IN GLIOBLASTOMA CELL STATE SIGNATURES EXAMINED USING REGION-SPECIFIC HUMAN BRAIN ORGANOIDS

Abstract Cell state plasticity is retained by all stem-like and differentiated cells within glioblastoma (GBM) tumors. GBM cells can exploit these state-shifting capacities to evade and resist therapeutics. Thus, there is a need to identify and target the regulators of GBM cell state transitions. To address this question, we use a fully human platform where hiSPC-derived organoids serve as a recipient environment to host cells engrafted directly from surgically resected GBMs. This allows us to precisely examine how microenvironmental variables impact the molecular signatures of GBM cells. We acutely isolate GBM cells from surgical resections, label cells with a GFP-lentivirus for 4 hours, and then engraft these tumor cells into parallelly patterned region-specific organoids that mimic dorsal forebrain, ventral forebrain, midbrain, and hindbrain identities. These regions were chosen as they capture derivatives of the three primary brain vesicles and because these regions contain unique cell types and signaling molecules. We performed paired single-cell RNA sequencing on tumor cells pre-engraftment and 14 days post-engraftment. Universal to all regional organoids, we find that GBM cells shift towards a neural progenitor-like and neuron-like signature within organoids, regardless of the cell state of the parent tumor prior to engraftment. To support our findings, we also observed that patient-derived glioma cell lines converge on a neural progenitor-like state within organoids. Thus, the dominance of this shared GBM cell state within organoids implies the existence of shared extrinsic factors in these immature (neurogenic) organoids that favor a progenitor-like or neuron-like identity across all four regions. We are now asking whether this bias towards neuronal and/or neural progenitor identity is a result of engrafting into young neurogenic organoids, or if the same phenomenon could be true in late-stage (post-gliogenic) mature organoids. We are also using this valuable platform and data to identify the extrinsic signals that GBM cells rely on to fulfill their plastic potential, adapt to the microenvironment, and resist therapeutics.

The symmetric 2-adic complexity of Tang-Gong interleaved sequences from Legendre sequence pair

The symmetric 2-adic complexity of Tang-Gong interleaved sequences from Legendre sequence pair

Fine-scale spatial and social patterns of SARS-CoV-2 transmission from identical pathogen sequences.

Pathogen genomics can provide insights into underlying infectious disease transmission patterns, but new methods are needed to handle modern large-scale pathogen genome datasets and realize this full potential. In particular, genetically proximal viruses should be highly informative about transmission events as genetic proximity indicates epidemiological linkage. Here, we leverage pairs of identical sequences to characterise fine-scale transmission patterns using 114,298 SARS-CoV-2 genomes collected through Washington State (USA) genomic sentinel surveillance with associated age and residence location information between March 2021 and December 2022. This corresponds to 59,660 sequences with another identical sequence in the dataset. We find that the location of pairs of identical sequences is highly consistent with expectations from mobility and social contact data. Outliers in the relationship between genetic and mobility data can be explained by SARS-CoV-2 transmission between postal codes with male prisons, consistent with transmission between prison facilities. We find that transmission patterns between age groups vary across spatial scales. Finally, we use the timing of sequence collection to understand the age groups driving transmission. Overall, this work improves our ability to leverage large pathogen genome datasets to understand the determinants of infectious disease spread.

Comparative Analyses of the Complete Mitogenomes of Two Oxyria Species (Polygonaceae) Provide Insights into Understanding the Mitogenome Evolution Within the Family.

Oxyria (Polygonaceae) is a small genus only comprising two species, Oxyria digyna and O. sinensis. Both species have well-documented usage in Chinese herbal medicine. We sequenced and assembled the complete mitogenomes of these two species and conducted a comparative analysis of the mitogenomes within Polygonaceae. Both O. digyna and O. sinensis displayed distinctive multi-branched conformations, consisting of one linear and one circular molecule. These two species shared similar gene compositions and exhibited distinct codon preferences, with mononucleotides as the most abundant type of simple sequence repeats. In the mitogenome of O. sinensis, a pair of long forward repeat sequences can mediate the division of molecule 1 into two sub-genomic circular molecules. Homologous sequence analysis revealed the occurrence of gene transfer between the chloroplast and mitochondrial genomes within Oxyria species. Additionally, a substantial number of homologous collinear blocks with varied arrangements were observed across different Polygonaceae species. Phylogenetic analysis suggested that mitogenome genes can serve as reliable markers for constructing phylogenetic relationships within Polygonaceae. Comparative analysis of eight species revealed Polygonaceae mitogenomes exhibited variability in gene presence, and most protein-coding genes (PCGs) have undergone negative selection. Overall, our study provided a comprehensive overview of the structural, functional, and evolutionary characteristics of the Polygonaceae mitogenomes.

Enzyme-Assisted Fluorescence Biosensor Based on Circular Single-Stranded DNA Without Group Modification for MicroRNA Detection.

Fluorescent biosensor, which has the characteristics of high sensitivity, specificity, and low cost, can be directly detected in physiological fluids such as blood and serum. Therefore, the development of fluorescence sensor platforms for miRNA detection has a positive effect on the prevention and treatment of various diseases. In this paper, miR-34a was selected as a biological indicator of Alzheimer's disease (AD). We designed a circular single-stranded DNA (CSSD) biosensor, which uses two unmodified single-stranded DNA (ssDNA) with complementary ends, DNAa and DNAb, to form CSSD by DNA sequence pairing to improve thermal stability and achieve signal amplification. At the same time, CSSD can react with miR-34a, and then the DNA of the DNA-RNA chain is hydrolyzed by duplex-specific nuclease (DSN enzyme). Finally, miR-34a is released to partake in the subsequent step, thus realizing cycle amplification. By evaluating the change in fluorescence signal under the optimized conditions, we discovered that this approach exhibits impressive sensitivity, with a detection threshold reaching as low as 0.36 nM. This surpasses the performance of numerous preceding miRNA detection biosensors. Furthermore, the system displays excellent detection capabilities even in intricate settings like serum, showcasing a strong ability to differentiate and choose effectively. In summary, this is a signal-off fluorescent biosensor, which realizes the purpose of double amplification of biosensor signal by using CSSD and enzyme assistance so that it can be used as a valuable tool for early diagnosis of diseases.

TIRTL-seq: Deep, quantitative, and affordable paired TCR repertoire sequencing.

ɑ/β T cells are key players in adaptive immunity. The specificity of T cells is determined by the sequences of the hypervariable T cell receptor (TCR) ɑ and β chains. Although bulk TCR sequencing offers a cost-effective approach for in-depth TCR repertoire profiling, it does not provide chain pairings, which are essential for determining T cell specificity. In contrast, single-cell TCR sequencing technologies produce paired chain data, but are limited in throughput to thousands of cells and are cost-prohibitive for cohort-scale studies. Here, we present TIRTL-seq (Throughput-Intensive Rapid TCR Library sequencing), a novel approach that generates ready-to-sequence TCR libraries from live cells in less than 7 hours. The protocol is optimized for use with non-contact liquid handlers in an automation-friendly 384-well plate format. Reaction volume miniaturization reduces library preparation costs to <$0.50 per well. The core principle of TIRTL-seq is the parallel generation of hundreds of libraries providing multiple biological replicates from a single sample that allows precise inference of both frequencies of individual clones and TCR chain pairings from well-occurrence patterns. We demonstrate scalability of our approach up to 1 million unique paired αβTCR clonotypes corresponding to over 30 million T cells per sample at a cost of less than $2000. For a sample of 10 million cells the cost is ~$200. We benchmarked TIRTL-seq against state-of-the-art 5'RACE bulk TCR-seq and 10x Genomics Chromium technologies on longitudinal samples. We show that TIRTL-seq is able to quantitatively identify expanding and contracting clonotypes between timepoints while providing accurate TCR chain pairings, including distinct temporal dynamics of SARS-CoV-2-specific and EBV-specific CD8+ T cell responses after infection. While clonal expansion was followed by sharp contraction for SARS-CoV-2 specific TCRs, EBV-specific TCRs remained stable once established. The sequences of both ɑ and β TCR chains are essential for determining T cell specificity. As the field moves towards greater applications in diagnostics and immunotherapy that rely on TCR specificity, we anticipate that our scalable paired TCR sequencing methodology will be instrumental for collecting large paired-chain datasets and ultimately extracting therapeutically relevant information from the TCR repertoire.

MiR-644a is a tumor cell-intrinsic mediator of sex bias in glioblastoma.

Biological sex is an important risk factor for glioblastoma (GBM), with males having a higher incidence and poorer prognosis. The mechanisms for this sex bias are thought to be both tumor intrinsic and tumor extrinsic. MicroRNAs (miRNAs), key posttranscriptional regulators of gene expression, have been previously linked to sex differences in various cell types and diseases, but their role in the sex bias of GBM remains unknown. We leveraged previously published paired miRNA and mRNA sequencing of 39 GBM patients (22 male, 17 female) to identify sex-biased miRNAs. We further interrogated a separate single-cell RNA-sequencing dataset of 110 GBM patients to examine whether differences in miRNA target gene expression were tumor cell-intrinsic or tumor cell extrinsic. Results were validated in a panel of patient-derived cell models. We identified 10 sex-biased miRNAs (p adjusted < .1), of which 3 were more highly expressed in males and 7 more highly expressed in females. Of these, miR-644a was higher in females, and increased expression of miR-644a target genes was significantly associated with decreased overall survival (HR 1.3, P = .02). Furthermore, analysis of an independent single-cell RNA-sequencing dataset confirmed sex-specific expression of miR-644a target genes in tumor cells (P < 10-15). Among patient-derived models, miR-644a was expressed a median of 4.8-fold higher in females compared to males. Our findings implicate miR-644a as a candidate tumor cell-intrinsic regulator of sex-biased gene expression in GBM.

Paired single-cell RNA sequencing and T cell receptor sequencing reveals both pro- and anti-inflammatory roles of T cells in patients with calcific aortic valve disease

Abstract Background Calcified aortic valve disease (CAVD) is a chronic unresolvable inflammatory disorder that affects 25% of adults over 65 years old, leading to aortic stenosis, heart failure, and death [1]. Since the molecular mechanisms of its pathogenesis are not well understood, therapies to treat root causes of the disease remain to be developed, with valve replacement as the major intervention strategy currently. Increasing evidence indicates that substantial T cells infiltrate in the aortic valve during calcification, and clonal expanded CD8+ T cells have been identified in the calcificed aortic valve by bulk T cell receptor (TCR) repertoire sequencing [2]. Purpose To decipher T cell heterogenity and understand the precise mechanisms of T cell subsets involving CAVD progression as a basis for novel therapeutic targets. Methods In Cohort 1, we collected aortic valves from CAVD (n=3) and healthy (n=2) patients for single-cell RNA sequencing (scRNA-seq). In Cohort 2, we performed paired scRNA-seq and TCR sequencing (scTCR-seq) on aortic valves obtained from CAVD (n=3) and healthy (n=3) patients to evaluate both gene expression and clonality. Pathway enrichment analysis and cell-cell interaction analysis elucidated the function of T cells in the calcificied valves. Immunohistochemistry and fluorescence-activated cell sorting were recruited to verify the key findings. Results scRNA-seq demonstrated a conversion from anti-inflammatory regulatory T (Treg) cells towards pro-inflammatory T helper 17 (Th17) cells in the calcificied aortic valve. Paired scRNA-seq and scTCR-seq revealed CAVD-related signaling pathways (osteoclast differentiation, NOD-like receptor signaling pathways, IL-17 signaling pathways) significantly upregulated in highly expanded CD8+ T cells in the calcified valve, with upregulation of pro-inflammatory mediators such as IFNG, CCL4 and CCL5, suggesting pro-inflammation functions in clonally expanded T cells. Moreover, we characterized tissue-resident memory T cells, which resembled similar transcriptome and TCR clonality as exhausted T cells in calcified valves. These cells up-regulated anti-inflammation-related transcriptomes but their cellular fraction decreased in calcified valves, implying anti-inflammatory functions of these tissue-resident memory and exhausted T cells. Conclusion This study elucidates transcriptomic and immune profiling of T cells in the calcified aortic valve, therefore shedding light on the pathophysiological mechanisms underlying aortic valve calcification from the novel perspective of tissue-specific T lymphocytes, providing promising targets for immune checkpoint-based therapeutic interventions.

Overloading And unpacKing (OAK) - droplet-based combinatorial indexing for ultra-high throughput single-cell multiomic profiling

Multiomic profiling of single cells by sequencing is a powerful technique for investigating cellular diversity. Existing droplet-based microfluidic methods produce many cell-free droplets, underutilizing bead barcodes and reagents. Combinatorial indexing on microplates is more efficient for barcoding but labor-intensive. Here we present Overloading And unpacKing (OAK), which uses a droplet-based barcoding system for initial compartmentalization followed by a second aliquoting round to achieve combinatorial indexing. We demonstrate OAK’s versatility with single-cell RNA sequencing as well as paired single-nucleus RNA sequencing and accessible chromatin profiling. We further showcase OAK’s performance on complex samples, including differentiated bronchial epithelial cells and primary retinal tissue. Finally, we examine transcriptomic responses of over 400,000 melanoma cells to a RAF inhibitor, belvarafenib, discovering a rare resistant cell population (0.12%). OAK’s ultra-high throughput, broad compatibility, high sensitivity, and simplified procedures make it a powerful tool for large-scale molecular analysis, even for rare cells.

The brain hierarchically represents the past and future during multistep anticipation

Memory for temporal structure enables both planning of future events and retrospection of past events. We investigated how the brain flexibly represents extended temporal sequences into the past and future during anticipation. Participants learned sequences of environments in immersive virtual reality. Pairs of sequences had the same environments in a different order, enabling context-specific learning. During fMRI, participants anticipated upcoming environments multiple steps into the future in a given sequence. Temporal structure was represented in the hippocampus and across higher-order visual regions (1) bidirectionally, with graded representations into the past and future and (2) hierarchically, with further events into the past and future represented in successively more anterior brain regions. In hippocampus, these bidirectional representations were context-specific, and suppression of far-away environments predicted response time costs in anticipation. Together, this work sheds light on how we flexibly represent sequential structure to enable planning over multiple timescales.

Assessment of Genetic Diversity and Population Structure of Exotic Sugar Beet (Beta vulgaris L.) Varieties Using Three Molecular Markers.

Sugar beet (Beta vulgaris L.) is a biennial herb belonging to the Amaranthaceae family. It contributes to approximately 30% of the world's total sucrose production and is an economically important crop. In this study, we analyzed the genetic diversity and population structure of 132 exotic sugar beet varieties using three molecular makers: four pairs of simple sequence repeat (SSR) primers, three pairs of insertion-deletion sequence (InDel) primers, and 20 cis-element amplification polymorphism (CEAP) primers. The results indicated that the number of alleles (Na) was 298, among which the number of effective alleles (Ne) was 182.426 (accounting for approximately 61.2%). The mean value of the genetic diversity index was 0.836. The polymorphic information content (PIC) was 0.639-0.907 (mean = 0.819), indicating a high level of polymorphism. These sugar beet varieties were classified into six clusters using the UPGMA method of cluster analysis. Population structure analysis revealed that the most ideal K value was 6. This indicated that the test materials could be divided into six categories, consistent with the clustering results. The clustering results indicated that most sugar beet varieties from the same breeding company clustered together, and the genetic distance between them was small, indicating that they may share the same male and/or female parent. Some varieties from different companies clustered together, indicating a narrow genetic base and potential exchange of germplasm resources between breeding companies. This study revealed the genetic differences among exotic sugar beet varieties and characteristics of the population structure. It provided a scientific basis for the identification of sugar beet varieties and markers-assisted breeding in China in the future.