Human immune organoids to decode B cell response in healthy donors and patients with lymphoma.

Antibodies are produced when naive B cells differentiate into plasma cells within germinal centres (GCs) of lymphoid tissues. Patients with B cell lymphoma on effective immunotherapies exhibit diminished antibody production, leading to higher infection rates and reduced vaccine efficacy, even after B cell recovery. Current ex vivo models fail to sustain long-term GC reactions and effectively test B cell responses. Here we developed synthetic hydrogels mimicking the lymphoid tissue microenvironment, enabling human GCs from tonsils and peripheral blood mononuclear cell-derived B cells. Immune organoids derived from peripheral blood mononuclear cells maintain GC B cells and plasma cells longer than tonsil-derived ones and exhibit unique B cell programming, including GC compartments, somatic hypermutation, immunoglobulin class switching and B cell clones. Chemical inhibition of transcriptional and epigenetic processes enhances plasma cell formation. While integrating polarized CXCL12 protein in a lymphoid organ-on-chip modulates GC responses in healthy donor B cells, it fails with B cells derived from patients with lymphoma. Our system allows rapid, controlled modelling of immune responses and B cell disorders.

Assembly of SARS-CoV-2 ribonucleosomes by truncated N∗ variant of the nucleocapsid protein.

The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) compacts the RNA genome into viral ribonucleoprotein (vRNP) complexes within virions. Assembly of vRNPs is inhibited by phosphorylation of the N protein serine/arginine (SR) region. Several SARS-CoV-2 variants of concern carry N protein mutations that reduce phosphorylation and enhance the efficiency of viral packaging. Variants of the dominant B.1.1 viral lineage also encode a truncated N protein, termed N∗ or Δ(1-209), that mediates genome packaging despite lacking the N-terminal RNA-binding domain and SR region. Here, we use mass photometry and negative stain electron microscopy to show that purified Δ(1-209) and viral RNA assemble into vRNPs that are remarkably similar in size and shape to those formed with full-length N protein. We show that assembly of Δ(1-209) vRNPs requires the leucine-rich helix of the central disordered region and that this helix promotes N protein oligomerization. We also find that fusion of a phosphomimetic SR region to Δ(1-209) inhibits RNA binding and vRNP assembly. Our results provide new insights into the mechanisms by which RNA binding promotes N protein self-association and vRNP assembly, and how this process is modulated by phosphorylation.

Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection.

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

An emerging view of coevolution in the legume-rhizobium mutualism.

Mutualisms are often framed as 'delicately balanced antagonisms' (Bronstein, 1994), with the net fitness benefits to both partners potentially masking underlying conflicts of interest. How commonly symbionts evolve to 'cheat' their hosts and hosts evolve to 'sanction' or 'control' uncooperative symbionts is the subject of debate, especially in legume-rhizobium interactions (Frederickson, 2013; Kiers et al., 2003). This kind of antagonistic coevolution should result in either arms-race dynamics characterized by repeated selective sweeps or fluctuating selection dynamics that leave signatures of balancing selection in host and symbiont genomes (Frederickson, 2013; Kortright et al., 2022; O'Brien et al., 2021). In a From the Cover article in this issue of Molecular Ecology, Epstein et al. (2022) combine GWAS and population genomic analyses to assess the evidence for positive or balancing selection consistent with ongoing, antagonistic coevolution between legumes and rhizobia. They found few genomic signatures of fitness conflicts between mutualistic partners, suggesting that legume and rhizobium fitness interests are largely aligned and symbiotic traits are mostly under stabilizing selection. In combination with other recent work (e.g. Batstone et al., 2020), the results of Epstein et al. (2022) indicate that there is little ongoing fitness conflict between legumes and rhizobia that shapes host and symbiont genomes in this system. It may be time to move beyond symbiont 'cheating' and host 'control' as the dominant paradigm for understanding how partners in mutualism coevolve.

Combinatorial treatment rescues tumour-microenvironment-mediated attenuation of MALT1 inhibitors in B-cell lymphomas.

Activated B-cell-like diffuse large B-cell lymphomas (ABC-DLBCLs) are characterized by constitutive activation of nuclear factor κB driven by the B-cell receptor (BCR) and Toll-like receptor (TLR) pathways. However, BCR-pathway-targeted therapies have limited impact on DLBCLs. Here we used >1,100 DLBCL patient samples to determine immune and extracellular matrix cues in the lymphoid tumour microenvironment (Ly-TME) and built representative synthetic-hydrogel-based B-cell-lymphoma organoids accordingly. We demonstrate that Ly-TME cellular and biophysical factors amplify the BCR-MYD88-TLR9 multiprotein supercomplex and induce cooperative signalling pathways in ABC-DLBCL cells, which reduce the efficacy of compounds targeting the BCR pathway members Bruton tyrosine kinase and mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1). Combinatorial inhibition of multiple aberrant signalling pathways induced higher antitumour efficacy in lymphoid organoids and implanted ABC-DLBCL patient tumours in vivo. Our studies define the complex crosstalk between malignant ABC-DLBCL cells and Ly-TME, and provide rational combinatorial therapies that rescue Ly-TME-mediated attenuation of treatment response to MALT1 inhibitors.

The evolution of partner specificity in mutualisms.

Mutualistic species vary in their level of partner specificity, which has important evolutionary, ecological, and management implications. Yet, the evolutionary mechanisms which underpin partner specificity are not fully understood. Most work on specialization focuses on the trade-off between generalism and specialism, where specialists receive more benefits from preferred partners at the expense of benefits from non-preferred partners, while generalists receive similar benefits from all partners. Because all mutualisms involve some degree of both cooperation and conflict between partners, we highlight that specialization to a mutualistic partner can be cooperative, increasing benefit to a focal species and a partner, or antagonistic, increasing resource extraction by a focal species from a partner. We devise an evolutionary game theoretic model to assess the evolutionary dynamics of cooperative specialization, antagonistic specialization, and generalism. Our model shows that cooperative specialization leads to bistability: stable equilibria with a specialist host and its preferred partner excluding all others. We also show that under cooperative specialization with spatial effects, generalists can thrive at the boundaries between differing specialist patches. Under antagonistic specialization, generalism is evolutionarily stable. We provide predictions for how a cooperation-antagonism continuum may determine the patterns of partner specificity that develop within mutualistic relationships.

Reconstitution of the SARS-CoV-2 ribonucleosome provides insights into genomic RNA packaging and regulation by phosphorylation.

The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 is responsible for compaction of the ∼30-kb RNA genome in the ∼90-nm virion. Previous studies suggest that each virion contains 35 to 40 viral ribonucleoprotein (vRNP) complexes, or ribonucleosomes, arrayed along the genome. There is, however, little mechanistic understanding of the vRNP complex. Here, we show that N protein, when combined in vitro with short fragments of the viral genome, forms 15-nm particles similar to the vRNP structures observed within virions. These vRNPs depend on regions of N protein that promote protein-RNA and protein-protein interactions. Phosphorylation of N protein in its disordered serine/arginine region weakens these interactions to generate less compact vRNPs. We propose that unmodified N protein binds structurally diverse regions in genomic RNA to form compact vRNPs within the nucleocapsid, while phosphorylation alters vRNP structure to support other N protein functions in viral transcription.

Reconstitution of the SARS-CoV-2 ribonucleosome provides insights into genomic RNA packaging and regulation by phosphorylation.

The nucleocapsid (N) protein of coronaviruses is responsible for compaction of the ∼30-kb RNA genome in the ∼100-nm virion. Cryo-electron tomography suggests that each virion contains 35-40 viral ribonucleoprotein (vRNP) complexes, or ribonucleosomes, arrayed along the genome. There is, however, little mechanistic understanding of the vRNP complex. Here, we show that N protein, when combined with viral RNA fragments in vitro, forms cylindrical 15-nm particles similar to the vRNP structures observed within coronavirus virions. These vRNPs form in the presence of stem-loop-containing RNA and depend on regions of N protein that promote protein-RNA and protein-protein interactions. Phosphorylation of N protein in its disordered serine/arginine (SR) region weakens these interactions and disrupts vRNP assembly. We propose that unmodified N binds stem-loop-rich regions in genomic RNA to form compact vRNP complexes within the nucleocapsid, while phosphorylated N maintains uncompacted viral RNA to promote the protein's transcriptional function.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits.

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study.

AIMS/HYPOTHESIS: Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study. METHODS: We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci. RESULTS: Four novel associations were identified (p < 5 × 10-9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis. CONCLUSIONS/INTERPRETATION: Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations. DATA AVAILABILITY: Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog ( https://www.ebi.ac.uk/gwas/downloads/summary-statistics ).

Multisite phosphorylation by Cdk1 initiates delayed negative feedback to control mitotic transcription.

Cell-cycle progression is driven by the phosphorylation of cyclin-dependent kinase (Cdk) substrates.1-3 The order of substrate phosphorylation depends in part on the general rise in Cdk activity during the cell cycle,4-7 together with variations in substrate docking to sites on associated cyclin and Cks subunits.3,6,8-10 Many substrates are modified at multiple sites to provide more complex regulation.10-14 Here, we describe an elegant regulatory circuit based on multisite phosphorylation of Ndd1, a transcriptional co-activator of budding yeast genes required for mitotic progression.11,12 As cells enter mitosis, Ndd1 phosphorylation by Cdk1 is known to promote mitotic cyclin (CLB2) gene transcription, resulting in positive feedback.13-16 Consistent with these findings, we show that low Cdk1 activity promotes CLB2 expression at mitotic entry. We also find, however, that when high Cdk1 activity accumulates in a mitotic arrest, CLB2 expression is inhibited. Inhibition is accompanied by Ndd1 degradation, and we present evidence that degradation is triggered by multisite Ndd1 phosphorylation by high mitotic Cdk1-Clb2 activity. Complete Ndd1 phosphorylation by Clb2-Cdk1-Cks1 requires the phosphothreonine-binding site of Cks1, as well as a recently identified phosphate-binding pocket on the cyclin Clb2.17 We therefore propose that initial phosphorylation by Cdk1 primes Ndd1 for delayed secondary phosphorylation at suboptimal sites that promote degradation. Together, our results suggest that rising levels of mitotic Cdk1 activity act at multiple phosphorylation sites on Ndd1, first triggering rapid positive feedback and then promoting delayed negative feedback, resulting in a pulse of mitotic gene expression.

Role of Endoscopic Retrograde Cholangiopancreatography in the Diagnosis and Management of Cholestatic Liver Diseases.

Cholestatic liver diseases (CLDs) occur as a result of bile duct injury, emanating into duct obstruction and bile stasis. Advances in radiological imaging in the last decade has replaced endoscopic retrograde cholangiopancreatography (ERCP) as the first diagnostic tool, except in certain groups of patients, such as those with ischemic cholangiopathy (IsC) or early stages of primary sclerosing cholangitis (PSC). ERCP provides an opportunity for targeted tissue acquisition for histopathological evaluation and carries a diverse therapeutic profile to restore bile flow. The aim of this review article is to appraise the diagnostic and therapeutic roles of ERCP in CLDs.

Marine heatwaves depress metabolic activity and impair cellular acid-base homeostasis in reef-building corals regardless of bleaching susceptibility.

Ocean warming is causing global coral bleaching events to increase in frequency, resulting in widespread coral mortality and disrupting the function of coral reef ecosystems. However, even during mass bleaching events, many corals resist bleaching despite exposure to abnormally high temperatures. While the physiological effects of bleaching have been well documented, the consequences of heat stress for bleaching-resistant individuals are not well understood. In addition, much remains to be learned about how heat stress affects cellular-level processes that may be overlooked at the organismal level, yet are crucial for coral performance in the short term and ecological success over the long term. Here we compared the physiological and cellular responses of bleaching-resistant and bleaching-susceptible corals throughout the 2019 marine heatwave in Hawai'i, a repeat bleaching event that occurred 4 years after the previous regional event. Relative bleaching susceptibility within species was consistent between the two bleaching events, yet corals of both resistant and susceptible phenotypes exhibited pronounced metabolic depression during the heatwave. At the cellular level, bleaching-susceptible corals had lower intracellular pH than bleaching-resistant corals at the peak of bleaching for both symbiont-hosting and symbiont-free cells, indicating greater disruption of acid-base homeostasis in bleaching-susceptible individuals. Notably, cells from both phenotypes were unable to compensate for experimentally induced cellular acidosis, indicating that acid-base regulation was significantly impaired at the cellular level even in bleaching-resistant corals and in cells containing symbionts. Thermal disturbances may thus have substantial ecological consequences, as even small reallocations in energy budgets to maintain homeostasis during stress can negatively affect fitness. These results suggest concern is warranted for corals coping with ocean acidification alongside ocean warming, as the feedback between temperature stress and acid-base regulation may further exacerbate the physiological effects of climate change.

HP1 proteins compact DNA into mechanically and positionally stable phase separated domains.

In mammals, HP1-mediated heterochromatin forms positionally and mechanically stable genomic domains even though the component HP1 paralogs, HP1α, HP1ß, and HP1γ, display rapid on-off dynamics. Here, we investigate whether phase-separation by HP1 proteins can explain these biological observations. Using bulk and single-molecule methods, we show that, within phase-separated HP1α-DNA condensates, HP1α acts as a dynamic liquid, while compacted DNA molecules are constrained in local territories. These condensates are resistant to large forces yet can be readily dissolved by HP1ß. Finally, we find that differences in each HP1 paralog's DNA compaction and phase-separation properties arise from their respective disordered regions. Our findings suggest a generalizable model for genome organization in which a pool of weakly bound proteins collectively capitalize on the polymer properties of DNA to produce self-organizing domains that are simultaneously resistant to large forces at the mesoscale and susceptible to competition at the molecular scale.

Phosphoregulation of Phase Separation by the SARS-CoV-2 N Protein Suggests a Biophysical Basis for its Dual Functions.

The nucleocapsid (N) protein of coronaviruses serves two major functions: compaction of the RNA genome in the virion and regulation of viral gene transcription. It is not clear how the N protein mediates such distinct functions. The N protein contains two RNA-binding domains surrounded by regions of intrinsic disorder. Phosphorylation of the central disordered region promotes the protein's transcriptional function, but the underlying mechanism is not known. Here, we show that the N protein of SARS-CoV-2, together with viral RNA, forms biomolecular condensates. Unmodified N protein forms partially ordered gel-like condensates and discrete 15-nm particles based on multivalent RNA-protein and protein-protein interactions. Phosphorylation reduces these interactions, generating a more liquid-like droplet. We propose that distinct oligomeric states support the two functions of the N protein: unmodified protein forms a structured oligomer that is suited for nucleocapsid assembly, and phosphorylated protein forms a liquid-like compartment for viral genome processing.

Phosphorylation modulates liquid-liquid phase separation of the SARS-CoV-2 N protein.

The nucleocapsid (N) protein of coronaviruses serves two major functions: compaction of the RNA genome in the virion and regulation of viral gene transcription in the infected cell 1-3 . The N protein contains two globular RNA-binding domains surrounded by regions of intrinsic disorder 4 . Phosphorylation of the central disordered region is required for normal viral genome transcription 5,6 , which occurs in a cytoplasmic structure called the replication transcription complex (RTC) 7-11 . It is not known how phosphorylation controls N protein function. Here we show that the N protein of SARS-CoV-2, together with viral RNA, forms biomolecular condensates 12-15 . Unmodified N protein forms partially ordered gel-like structures that depend on multivalent RNA-protein and protein-protein interactions. Phosphorylation reduces a subset of these interactions, generating a more liquid-like droplet. We speculate that distinct oligomeric states support the two functions of the N protein: unmodified protein forms a structured oligomer that is suited for nucleocapsid assembly, and phosphorylated protein forms a liquid-like compartment for viral genome processing. Inhibitors of N protein phosphorylation could therefore serve as antiviral therapy.

Minority-centric meta-analyses of blood lipid levels identify novel loci in the Population Architecture using Genomics and Epidemiology (PAGE) study.

Lipid levels are important markers for the development of cardio-metabolic diseases. Although hundreds of associated loci have been identified through genetic association studies, the contribution of genetic factors to variation in lipids is not fully understood, particularly in U.S. minority groups. We performed genome-wide association analyses for four lipid traits in over 45,000 ancestrally diverse participants from the Population Architecture using Genomics and Epidemiology (PAGE) Study, followed by a meta-analysis with several European ancestry studies. We identified nine novel lipid loci, five of which showed evidence of replication in independent studies. Furthermore, we discovered one novel gene in a PrediXcan analysis, minority-specific independent signals at eight previously reported loci, and potential functional variants at two known loci through fine-mapping. Systematic examination of known lipid loci revealed smaller effect estimates in African American and Hispanic ancestry populations than those in Europeans, and better performance of polygenic risk scores based on minority-specific effect estimates. Our findings provide new insight into the genetic architecture of lipid traits and highlight the importance of conducting genetic studies in diverse populations in the era of precision medicine.

Genetics of Chronic Kidney Disease Stages Across Ancestries: The PAGE Study.

BACKGROUND: Chronic kidney disease (CKD) is common and disproportionally burdens United States ethnic minorities. Its genetic determinants may differ by disease severity and clinical stages. To uncover genetic factors associated CKD severity among high-risk ethnic groups, we performed genome-wide association studies (GWAS) in diverse populations within the Population Architecture using Genomics and Epidemiology (PAGE) study. METHODS: We assembled multi-ethnic genome-wide imputed data on CKD non-overlapping cases [4,150 mild to moderate CKD, 1,105 end-stage kidney disease (ESKD)] and non-CKD controls for up to 41,041 PAGE participants (African Americans, Hispanics/Latinos, East Asian, Native Hawaiian, and American Indians). We implemented a generalized estimating equation approach for GWAS using ancestry combined data while adjusting for age, sex, principal components, study, and ethnicity. RESULTS: The GWAS identified a novel genome-wide associated locus for mild to moderate CKD nearby NMT2 (rs10906850, p = 3.7 × 10-8) that replicated in the United Kingdom Biobank white British (p = 0.008). Several variants at the APOL1 locus were associated with ESKD including the APOL1 G1 rs73885319 (p = 1.2 × 10-9). There was no overlap among associated loci for CKD and ESKD traits, even at the previously reported APOL1 locus (p = 0.76 for CKD). Several additional loci were associated with CKD or ESKD at p-values below the genome-wide threshold. These loci were often driven by variants more common in non-European ancestry. CONCLUSION: Our genetic study identified a novel association at NMT2 for CKD and showed for the first time strong associations of the APOL1 variants with ESKD across multi-ethnic populations. Our findings suggest differences in genetic effects across CKD severity and provide information for study design of genetic studies of CKD in diverse populations.