A conserved complex lipid signature marks human muscle aging and responds to short-term exercise.

Studies in preclinical models suggest that complex lipids, such as phospholipids, play a role in the regulation of longevity. However, identification of universally conserved complex lipid changes that occur during aging, and how these respond to interventions, is lacking. Here, to comprehensively map how complex lipids change during aging, we profiled ten tissues in young versus aged mice using a lipidomics platform. Strikingly, from >1,200 unique lipids, we found a tissue-wide accumulation of bis(monoacylglycero)phosphate (BMP) during mouse aging. To investigate translational value, we assessed muscle tissue of young and older people, and found a similar marked BMP accumulation in the human aging lipidome. Furthermore, we found that a healthy-aging intervention consisting of moderate-to-vigorous exercise was able to lower BMP levels in postmenopausal female research participants. Our work implicates complex lipid biology as central to aging, identifying a conserved aging lipid signature of BMP accumulation that is modifiable upon a short-term healthy-aging intervention.

Systematic evaluation of B-cell clonal family inference approaches.

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

Understanding repertoire sequencing data through a multiscale computational model of the germinal center.

Sequencing of B-cell and T-cell immune receptor repertoires helps us to understand the adaptive immune response, although it only provides information about the clonotypes (lineages) and their frequencies and not about, for example, their affinity or antigen (Ag) specificity. To further characterize the identified clones, usually with special attention to the particularly abundant ones (dominant), additional time-consuming or expensive experiments are generally required. Here, we present an extension of a multiscale model of the germinal center (GC) that we previously developed to gain more insight in B-cell repertoires. We compare the extent that these simulated repertoires deviate from experimental repertoires established from single GCs, blood, or tissue. Our simulations show that there is a limited correlation between clonal abundance and affinity and that there is large affinity variability among same-ancestor (same-clone) subclones. Our simulations suggest that low-abundance clones and subclones, might also be of interest since they may have high affinity for the Ag. We show that the fraction of plasma cells (PCs) with high B-cell receptor (BcR) mRNA content in the GC does not significantly affect the number of dominant clones derived from single GCs by sequencing BcR mRNAs. Results from these simulations guide data interpretation and the design of follow-up experiments.

B-cell receptor profiling before and after IVIG monotherapy in newly diagnosed idiopathic inflammatory myopathies.

OBJECTIVE: To unravel B-cell receptor (BcR) characteristics in muscle tissues and peripheral blood and gain more insight into BcR repertoire changes in peripheral blood in idiopathic inflammatory myopathies (IIMs), and study how this correlates to the clinical response to IVIG. METHODS: Nineteen treatment-naive patients with newly diagnosed IIM were prospectively treated with IVIG monotherapy. RNA-based BcR repertoire sequencing was performed in muscle biopsies collected before, and in peripheral blood (PB) collected before and nine weeks after IVIG treatment. Results were correlated to patients' clinical improvement based on the total improvement score (TIS). RESULTS: Prior to IVIG treatment, BcR clones found in muscle tissue could be retrieved in peripheral blood. Nine weeks after IVIG treatment, new patient-specific dominant BcR clones appeared in peripheral blood while pre-treatment dominant BcR clones disappeared. The cumulative frequency of all dominant BcR clones before treatment was significantly higher in individuals who responded to IVIG compared with those who did not respond to IVIG, and correlated with a higher CK. During follow-up, a decrease in the cumulative frequency of all dominant clones correlated with a higher TIS. CONCLUSION: In treatment-naive patients with newly diagnosed IIM, muscle tissue and peripheral blood share expanded BcR clones. In our study a higher cumulative frequency of dominant BcR clones in blood before treatment was associated with a higher CK and better treatment response, suggesting that response to IVIG may depend on the composition of the pre-treatment BcR repertoire.

TCRß clones in muscle tissue share structural features in patients with idiopathic inflammatory myopathy and are associated with disease activity.

Objectives: To characterize the T cell receptor (TCRß) repertoire in peripheral blood and muscle tissues of treatment naïve patients with newly diagnosed idiopathic inflammatory myopathies (IIMs). Methods: High throughput RNA sequencing of the TCRß chain was performed in peripheral blood and muscle tissue in twenty newly-diagnosed treatment-naïve IIM patients (9 DM, 5 NM/OM, 5 IMNM and 1 ASyS) and healthy controls. Results thereof were correlated with markers of disease activity. Results: Muscle tissue of IIM patients shows more expansion of TCRß clones and decreased diversity when compared to peripheral blood of IIM as well as healthy controls (both p=0.0001). Several expanded TCRß clones in muscle are tissue restricted and cannot be retrieved in peripheral blood. These clones have significantly longer CDR3 regions when compared to clones (also) found in circulation (p=0.0002), while their CDR3 region is more hydrophobic (p<0.01). Network analysis shows that clonal TCRß signatures are shared between patients. Increased clonal expansion in muscle tissue is significantly correlated with increased CK levels (p=0.03), while it tends to correlate with decreased muscle strength (p=0.08). Conclusion: Network analysis of clones in muscle of IIM patients shows shared clusters of sequences across patients. Muscle-restricted CDR3 TCRß clones show specific structural features in their T cell receptor. Our results indicate that clonal TCRß expansion in muscle tissue might be associated with disease activity. Collectively, these findings support a role for specific clonal T cell responses in muscle tissue in the pathogenesis of the IIM subtypes studied.

In rheumatoid arthritis inflamed joints share dominant patient-specific B-cell clones.

Background: In patients with rheumatoid arthritis (RA) different joints were shown to share the same dominant T-cell clones, suggesting shared characteristics of the inflammatory process and indicating that strategies to selectively target the antigen receptor might be feasible. Since T- and B-lymphocytes closely interact in adaptive responses, we analysed to what extent different joints also share dominant B-cell clones. Methods: In 11 RA patients, quantitative B-cell receptor (BCR) repertoire analysis was performed in simultaneously obtained samples from inflamed synovial tissue (ST) from distinct locations within one joint, from multiple joints, from synovial fluid (SF) and peripheral blood (PB). Results: ST biopsies from different locations in the same joint showed clear overlap in the top-25 dominant BCR clones (16.7%, SD 12.5), in the same range as the overlap between ST and SF in the same joint (8.0%, SD 8.8) and the overlap between ST-ST between different joints (9.1%, SD 8.2), but clearly higher than the overlap between ST and PB (1.7%, SD 2.4; p<0.05) and SF and PB (2.7%, SD 4.1; p<0.05). Interestingly, these figures were substantially lower than the overlap observed in previous T-cell clonality studies. Conclusions: We conclude that in RA BCR clonal responses may be more localized than TCR clonal responses, pointing to antigen-selective influx, proliferation and/or maturation of B-cells. B lineage cells in the SF may adequately represent the dominant BCR clones of the ST, which is in contrast to T-cells. Collectively, the presence of shared B- and especially T-cells in different joints from the same patient suggests that approaches might be feasible that aim to develop antigen-receptor specific targeting of lymphocyte clones in RA as an alternative to more generalized immunosuppressive strategies.

Identification of IOMA-class neutralizing antibodies targeting the CD4-binding site on the HIV-1 envelope glycoprotein.

A major goal of current HIV-1 vaccine design efforts is to induce broadly neutralizing antibodies (bNAbs). The VH1-2-derived bNAb IOMA directed to the CD4-binding site of the HIV-1 envelope glycoprotein is of interest because, unlike the better-known VH1-2-derived VRC01-class bNAbs, it does not require a rare short light chain complementarity-determining region 3 (CDRL3). Here, we describe three IOMA-class NAbs, ACS101-103, with up to 37% breadth, that share many characteristics with IOMA, including an average-length CDRL3. Cryo-electron microscopy revealed that ACS101 shares interactions with those observed with other VH1-2 and VH1-46-class bNAbs, but exhibits a unique binding mode to residues in loop D. Analysis of longitudinal sequences from the patient suggests that a transmitter/founder-virus lacking the N276 glycan might have initiated the development of these NAbs. Together these data strengthen the rationale for germline-targeting vaccination strategies to induce IOMA-class bNAbs and provide a wealth of sequence and structural information to support such strategies.

A public antibody class recognizes an S2 epitope exposed on open conformations of SARS-CoV-2 spike.

Delineating the origins and properties of antibodies elicited by SARS-CoV-2 infection and vaccination is critical for understanding their benefits and potential shortcomings. Therefore, we investigate the SARS-CoV-2 spike (S)-reactive B cell repertoire in unexposed individuals by flow cytometry and single-cell sequencing. We show that ∼82% of SARS-CoV-2 S-reactive B cells harbor a naive phenotype, which represents an unusually high fraction of total human naive B cells (∼0.1%). Approximately 10% of these naive S-reactive B cells share an IGHV1-69/IGKV3-11 B cell receptor pairing, an enrichment of 18-fold compared to the complete naive repertoire. Following SARS-CoV-2 infection, we report an average 37-fold enrichment of IGHV1-69/IGKV3-11 B cell receptor pairing in the S-reactive memory B cells compared to the unselected memory repertoire. This class of B cells targets a previously undefined non-neutralizing epitope on the S2 subunit that becomes exposed on S proteins used in approved vaccines when they transition away from the native pre-fusion state because of instability. These findings can help guide the improvement of SARS-CoV-2 vaccines.

Longitudinal analysis of anti-drug antibody development in multiple sclerosis patients treated with interferon beta-1a (Rebif™) using B cell receptor repertoire analysis.

A significant proportion of multiple sclerosis (MS) patients treated with interferon beta-1a (Rebif™) develop anti-drug antibodies (ADA) with a negative impact on treatment efficacy. We hypothesized that high-throughput B-cell receptor (BCR) repertoire analysis could be used to predict and monitor ADA development. To study this we analyzed 228 peripheral blood samples from 68 longitudinally followed patients starting on interferon beta-1a. Our results show that whole blood BCR analysis does not reflect, and does not predict ADA development in MS patients treated with interferon beta-1a. We propose that BCR analysis of phenotypically selected cell subsets or tissues might be more informative.

From affinity selection to kinetic selection in Germinal Centre modelling.

Affinity maturation is an evolutionary process by which the affinity of antibodies (Abs) against specific antigens (Ags) increases through rounds of B-cell proliferation, somatic hypermutation, and positive selection in germinal centres (GC). The positive selection of B cells depends on affinity, but the underlying mechanisms of affinity discrimination and affinity-based selection are not well understood. It has been suggested that selection in GC depends on both rapid binding of B-cell receptors (BcRs) to Ags which is kinetically favourable and tight binding of BcRs to Ags, which is thermodynamically favourable; however, it has not been shown whether a selection bias for kinetic properties is present in the GC. To investigate the GC selection bias towards rapid and tight binding, we developed an agent-based model of GC and compared the evolution of founder B cells with initially identical low affinities but with different association/dissociation rates for Ag presented by follicular dendritic cells in three Ag collection mechanisms. We compared an Ag collection mechanism based on association/dissociation rates of B-cell interaction with presented Ag, which includes a probabilistic rupture of bonds between the B-cell and Ag (Scenario-1) with a reference scenario based on an affinity-based Ag collection mechanism (Scenario-0). Simulations showed that the mechanism of Ag collection affects the GC dynamics and the GC outputs concerning fast/slow (un)binding of B cells to FDC-presented Ags. In particular, clones with lower dissociation rates outcompete clones with higher association rates in Scenario-1, while remaining B cells from clones with higher association rates reach higher affinities. Accordingly, plasma cell and memory B cell populations were biased towards B-cell clones with lower dissociation rates. Without such probabilistic ruptures during the Ag extraction process (Scenario-2), the selective advantage for clones with very low dissociation rates diminished, and the affinity maturation level of all clones decreased to the reference level.

Cell Type-Selective Loss of Peroxisomal ß-Oxidation Impairs Bipolar Cell but Not Photoreceptor Survival in the Retina.

Retinal degeneration is a common feature in peroxisomal disorders leading to blindness. Peroxisomes are present in the different cell types of the retina; however, their precise contribution to retinal integrity is still unclear. We previously showed that mice lacking the central peroxisomal ß-oxidation enzyme, multifunctional protein 2 (MFP2), develop an early onset retinal decay including photoreceptor cell death. To decipher the function of peroxisomal ß-oxidation in photoreceptors, we generated cell type selective Mfp2 knockout mice, using the Crx promotor targeting photoreceptors and bipolar cells. Surprisingly, Crx-Mfp2-/- mice maintained photoreceptor length and number until the age of 1 year. A negative electroretinogram was indicative of preserved photoreceptor phototransduction, but impaired downstream bipolar cell signaling from the age of 6 months. The photoreceptor ribbon synapse was affected, containing free-floating ribbons and vesicles with altered size and density. The bipolar cell interneurons sprouted into the ONL and died. Whereas docosahexaenoic acid levels were normal in the neural retina, levels of lipids containing very long chain polyunsaturated fatty acids were highly increased. Crx-Pex5-/- mice, in which all peroxisomal functions are inactivated in photoreceptors and bipolar cells, developed the same phenotype as Crx-Mfp2-/- mice. In conclusion, the early photoreceptor death in global Mfp2-/- mice is not driven cell autonomously. However, peroxisomal ß-oxidation is essential for the integrity of photoreceptor ribbon synapses and of bipolar cells.

Coupled Antigen and BLIMP1 Asymmetric Division With a Large Segregation Between Daughter Cells Recapitulates the Temporal Transition From Memory B Cells to Plasma Cells and a DZ-to-LZ Ratio in the Germinal Center.

Memory B cells and antibody-secreting plasma cells are generated within germinal centers during affinity maturation in which B-cell proliferation, selection, differentiation, and self-renewal play important roles. The mechanisms behind memory B cell and plasma cell differentiation in germinal centers are not well understood. However, it has been suggested that cell fate is (partially) determined by asymmetric cell division, which involves the unequal distribution of cellular components to both daughter cells. To investigate what level and/or probability of asymmetric segregation of several fate determinant molecules, such as the antigen and transcription factors (BCL6, IRF4, and BLIMP1) recapitulates the temporal switch and DZ-to-LZ ratio in the germinal center, we implemented a multiscale model that combines a core gene regulatory network for plasma cell differentiation with a model describing the cellular interactions and dynamics in the germinal center. Our simulations show that BLIMP1 driven plasma cell differentiation together with coupled asymmetric division of antigen and BLIMP1 with a large segregation between the daughter cells results in a germinal center DZ-to-LZ ratio and a temporal switch from memory B cells to plasma cells that have been observed in experiments.

Circadian misalignment disturbs the skeletal muscle lipidome in healthy young men.

Circadian misalignment, as seen in shift work, is associated with an increased risk to develop type 2 diabetes. In an experimental setting, we recently showed that a rapid day-night shift for 3 consecutive nights leads to misalignment of the core molecular clock, induction of the PPAR pathway, and insulin resistance in skeletal muscle of young, healthy men. Here, we investigated if circadian misalignment affects the skeletal muscle lipidome and intramyocellular lipid droplet characteristics, explaining the misalignment-induced insulin resistance. Fourteen healthy men underwent one aligned and one circadian misalignment period, both consisting of ~3.5 days. In the misaligned condition, day and night were rapidly shifted by 12 hours leading to opposite eating, sleep, and activity times compared with the aligned condition. For each condition, two muscle biopsies were taken from the m. vastus lateralis in the morning and evening and subjected to semi-targeted lipidomics and confocal microscopy analysis. We found that only 2% of detected lipids were different between morning and evening in the aligned condition, whereas 12% displayed a morning-evening difference upon misalignment. Triacylglycerols, in particular species of a carbon length ≥55, were the most abundant lipid species changed upon misalignment. Cardiolipins were decreased upon misalignment, whereas phosphatidylcholines consistently followed the same morning-evening pattern, suggesting regulation by the circadian clock. Cholesteryl esters adjusted to the shifted behavior. Lipid droplet characteristics remained unaltered upon misalignment. Together, these findings show that simulated shift work disturbs the skeletal muscle lipidome, which may contribute to misalignment-induced insulin resistance.

Metabolomics and lipidomics in Caenorhabditis elegans using a single-sample preparation.

Comprehensive metabolomic and lipidomic mass spectrometry methods are in increasing demand; for instance, in research related to nutrition and aging. The nematode Caenorhabditis elegans is a key model organism in these fields, owing to the large repository of available C. elegans mutants and their convenient natural lifespan. Here, we describe a robust and sensitive analytical method for the semi-quantitative analysis of >100 polar (metabolomics) and >1000 apolar (lipidomics) metabolites in C. elegans, using a single-sample preparation. Our method is capable of reliably detecting a wide variety of biologically relevant metabolic aberrations in, for example, glycolysis and the tricarboxylic acid cycle, pyrimidine metabolism and complex lipid biosynthesis. In conclusion, we provide a powerful analytical tool that maximizes metabolic data yield from a single sample. This article has an associated First Person interview with the joint first authors of the paper.

Peroxisomal Multifunctional Protein 2 Deficiency Perturbs Lipid Homeostasis in the Retina and Causes Visual Dysfunction in Mice.

Patients lacking multifunctional protein 2 (MFP2), the central enzyme of the peroxisomal ß-oxidation pathway, develop retinopathy. This pathway is involved in the metabolism of very long chain (VLCFAs) and polyunsaturated (PUFAs) fatty acids, which are enriched in the photoreceptor outer segments (POS). The molecular mechanisms underlying the retinopathy remain, however, elusive. Here, we report that mice with MFP2 inactivation display decreased retinal function already at the age of 3 weeks, which is accompanied by a profound shortening of the photoreceptor outer and inner segments, but with preserved photoreceptor ultrastructure. Furthermore, MFP2 deficient retinas exhibit severe changes in gene expression with downregulation of genes involved in the phototransduction pathway and upregulation of inflammation related genes. Lipid profiling of the mutant retinas revealed a profound reduction of DHA-containing phospholipids. This was likely due to a hampered systemic supply and retinal traffic of this PUFA, although we cannot exclude that the local defect of peroxisomal ß-oxidation contributes to this DHA decrease. Moreover, very long chain PUFAs were also reduced, with the exception of those containing ≥ 34 carbons that accumulated. The latter suggests that there is an uncontrollable elongation of retinal PUFAs. In conclusion, our data reveal that intact peroxisomal ß-oxidation is indispensable for retinal integrity, most likely by maintaining PUFA homeostasis.

An autosomal dominant neurological disorder caused by de novo variants in FAR1 resulting in uncontrolled synthesis of ether lipids.

PURPOSE: In this study we investigate the disease etiology in 12 patients with de novo variants in FAR1 all resulting in an amino acid change at position 480 (p.Arg480Cys/His/Leu). METHODS: Following next-generation sequencing and clinical phenotyping, functional characterization was performed in patients' fibroblasts using FAR1 enzyme analysis, FAR1 immunoblotting/immunofluorescence, and lipidomics. RESULTS: All patients had spastic paraparesis and bilateral congenital/juvenile cataracts, in most combined with speech and gross motor developmental delay and truncal hypotonia. FAR1 deficiency caused by biallelic variants results in defective ether lipid synthesis and plasmalogen deficiency. In contrast, patients' fibroblasts with the de novo FAR1 variants showed elevated plasmalogen levels. Further functional studies in fibroblasts showed that these variants cause a disruption of the plasmalogen-dependent feedback regulation of FAR1 protein levels leading to uncontrolled ether lipid production. CONCLUSION: Heterozygous de novo variants affecting the Arg480 residue of FAR1 lead to an autosomal dominant disorder with a different disease mechanism than that of recessive FAR1 deficiency and a diametrically opposed biochemical phenotype. Our findings show that for patients with spastic paraparesis and bilateral cataracts, FAR1 should be considered as a candidate gene and added to gene panels for hereditary spastic paraplegia, cerebral palsy, and juvenile cataracts.

Cerebellar and hepatic alterations in ACBD5-deficient mice are associated with unexpected, distinct alterations in cellular lipid homeostasis.

ACBD5 deficiency is a novel peroxisome disorder with a largely uncharacterized pathology. ACBD5 was recently identified in a tethering complex mediating membrane contacts between peroxisomes and the endoplasmic reticulum (ER). An ACBD5-deficient mouse was analyzed to correlate ACBD5 tethering functions with the disease phenotype. ACBD5-deficient mice exhibit elevated very long-chain fatty acid levels and a progressive cerebellar pathology. Liver did not exhibit pathologic changes but increased peroxisome abundance and drastically reduced peroxisome-ER contacts. Lipidomics of liver and cerebellum revealed tissue-specific alterations in distinct lipid classes and subspecies. In line with the neurological pathology, unusual ultra-long chain fatty acids (C > 32) were elevated in phosphocholines from cerebelli but not liver indicating an organ-specific imbalance in fatty acid degradation and elongation pathways. By contrast, ether lipid formation was perturbed in liver towards an accumulation of alkyldiacylglycerols. The alterations in several lipid classes suggest that ACBD5, in addition to its acyl-CoA binding function, might maintain peroxisome-ER contacts in order to contribute to the regulation of anabolic and catabolic cellular lipid pathways.

Time-Resolved Profiling Reveals ATF3 as a Novel Mediator of Endocrine Resistance in Breast Cancer.

Breast cancer is one of the leading causes of death for women worldwide. Patients whose tumors express Estrogen Receptor α account for around 70% of cases and are mostly treated with targeted endocrine therapy. However, depending on the degree of severity of the disease at diagnosis, 10 to 40% of these tumors eventually relapse due to resistance development. Even though recent novel approaches as the combination with CDK4/6 inhibitors increased the overall survival of relapsing patients, this remains relatively short and there is a urgent need to find alternative targetable pathways. In this study we profiled the early phases of the resistance development process to uncover drivers of this phenomenon. Time-resolved analysis revealed that ATF3, a member of the ATF/CREB family of transcription factors, acts as a novel regulator of the response to therapy via rewiring of central signaling processes towards the adaptation to endocrine treatment. ATF3 was found to be essential in controlling crucial processes such as proliferation, cell cycle, and apoptosis during the early response to treatment through the regulation of MAPK/AKT signaling pathways. Its essential role was confirmed in vivo in a mouse model, and elevated expression of ATF3 was verified in patient datasets, adding clinical relevance to our findings. This study proposes ATF3 as a novel mediator of endocrine resistance development in breast cancer and elucidates its role in the regulation of downstream pathways activities.