Loperamide positive deaths in Sweden 2012-2022 and Finland 2017-2022: Fatal loperamide intoxication exclusively for Sweden.

Loperamide, a potent µ-opioid receptor agonist used as an antidiarrheal drug, exhibits increased bioavailability at supratherapeutic doses, causing potential central nervous system effects. Its misuse for opioid withdrawal relief and euphoria can lead to dangerously elevated blood levels, causing severe cardiac dysrhythmias and death. This study aimed to compare loperamide positive autopsy cases in Sweden and Finland after the introduction of postmortem toxicological analysis of loperamide, focusing on loperamide's role in fatalities and identifying common characteristics among those affected. All cases with detected loperamide in femoral blood at forensic autopsies in Sweden (2012-2022) and Finland (2017-2022) were included. In Sweden, loperamide was detected in 126 individuals, and in Finland, in 111 individuals. The incidence of individuals positive for loperamide in postmortem femoral blood increased steadily over the study duration in both Sweden and Finland. Loperamide related fatalities were observed exclusively in Sweden (n=80), predominantly involving younger males with histories of substance abuse, typically classified as accidental deaths. The group of loperamide nonrelated deaths in Sweden mirrored the entirety of cases in Finland. The concentration of loperamide in postmortem femoral blood was significantly higher in cases where loperamide was considered the cause of death (median 0.140 µg/g) compared to cases where loperamide contributed (median 0.080 µg/g), as well as in deaths unrelated to loperamide in both countries (Sweden: median 0.029 µg/g; Finland: median 0.010 µg/ml). The high limit of quantification for loperamide in Sweden may underestimate therapeutic users in epidemiological assessments. This study underscores the absence of loperamide misuse in Finland and indicates a rising trend of loperamide abuse in Sweden.

Clonal associations between lymphocyte subsets and functional states in rheumatoid arthritis synovium.

Rheumatoid arthritis (RA) is an autoimmune disease involving antigen-specific T and B cells. Here, we perform single-cell RNA and repertoire sequencing on paired synovial tissue and blood samples from 12 seropositive RA patients. We identify clonally expanded CD4 + T cells, including CCL5+ cells and T peripheral helper (Tph) cells, which show a prominent transcriptomic signature of recent activation and effector function. CD8 + T cells show higher oligoclonality than CD4 + T cells, with the largest synovial clones enriched in GZMK+ cells. CD8 + T cells with possibly virus-reactive TCRs are distributed across transcriptomic clusters. In the B cell compartment, NR4A1+ activated B cells, and plasma cells are enriched in the synovium and demonstrate substantial clonal expansion. We identify synovial plasma cells that share BCRs with synovial ABC, memory, and activated B cells. Receptor-ligand analysis predicted IFNG and TNFRSF members as mediators of synovial Tph-B cell interactions. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate the synovium of patients with RA.

Genome-wide association study identifies host genetic variants influencing oral microbiota diversity and metabolic health.

The microbial communities of the oral cavity are important elements of oral and systemic health. With emerging evidence highlighting the heritability of oral bacterial microbiota, this study aimed to identify host genome variants that influence oral microbial traits. Using data from 16S rRNA gene amplicon sequencing, we performed genome-wide association studies with univariate and multivariate traits of the salivary microbiota from 610 unrelated adults from the Danish ADDITION-PRO cohort. We identified six single nucleotide polymorphisms (SNPs) in human genomes that showed associations with abundance of bacterial taxa at different taxonomical tiers (P < 5 × 10-8). Notably, SNP rs17793860 surpassed our study-wide significance threshold (P < 1.19 × 10-9). Additionally, rs4530093 was linked to bacterial beta diversity (P < 5 × 10-8). Out of these seven SNPs identified, six exerted effects on metabolic traits, including glycated hemoglobin A1c, triglyceride and high-density lipoprotein cholesterol levels, the risk of type 2 diabetes and stroke. Our findings highlight the impact of specific host SNPs on the composition and diversity of the oral bacterial community. Importantly, our results indicate an intricate interplay between host genetics, the oral microbiota, and metabolic health. We emphasize the need for integrative approaches considering genetic, microbial, and metabolic factors.

The chromatin landscape of pathogenic transcriptional cell states in rheumatoid arthritis.

Synovial tissue inflammation is a hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. Here, we examine genome-wide open chromatin at single-cell resolution in 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identify 24 chromatin classes and predict their associated transcription factors, including a CD8 + GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating with an RA tissue transcriptional atlas, we propose that these chromatin classes represent 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrate the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance, as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

A harmonized database of European forest simulations under climate change.

Process-based forest models combine biological, physical, and chemical process understanding to simulate forest dynamics as an emergent property of the system. As such, they are valuable tools to investigate the effects of climate change on forest ecosystems. Specifically, they allow testing of hypotheses regarding long-term ecosystem dynamics and provide means to assess the impacts of climate scenarios on future forest development. As a consequence, numerous local-scale simulation studies have been conducted over the past decades to assess the impacts of climate change on forests. These studies apply the best available models tailored to local conditions, parameterized and evaluated by local experts. However, this treasure trove of knowledge on climate change responses remains underexplored to date, as a consistent and harmonized dataset of local model simulations is missing. Here, our objectives were (i) to compile existing local simulations on forest development under climate change in Europe in a common database, (ii) to harmonize them to a common suite of output variables, and (iii) to provide a standardized vector of auxiliary environmental variables for each simulated location to aid subsequent investigations. Our dataset of European stand- and landscape-level forest simulations contains over 1.1 million simulation runs representing 135 million simulation years for more than 13,000 unique locations spread across Europe. The data were harmonized to consistently describe forest development in terms of stand structure (dominant height), composition (dominant species, admixed species), and functioning (leaf area index). Auxiliary variables provided include consistent daily climate information (temperature, precipitation, radiation, vapor pressure deficit) as well as information on local site conditions (soil depth, soil physical properties, soil water holding capacity, plant-available nitrogen). The present dataset facilitates analyses across models and locations, with the aim to better harness the valuable information contained in local simulations for large-scale policy support, and for fostering a deeper understanding of the effects of climate change on forest ecosystems in Europe.

Fatal Intoxications with Zopiclone-A Cause for Concern?

INTRODUCTION: Zopiclone, a controlled substance prescribed for insomnia, has become a common toxicological finding in forensic autopsy cases. This study investigated the role and extent of zopiclone use in fatal intoxications in Sweden. METHODS: All forensic autopsy cases positive for zopiclone in femoral blood during 2012-2020 were selected. Among these cases, fatalities caused by intoxication according to the cause of death certificates issued by the forensic pathologist were identified. Intoxications where zopiclone contributed to the cause of death were included in the study. The Swedish Prescribed Drug Register was utilized to examine whether the included cases were prescribed zopiclone or not. RESULTS: In total 7320 fatal intoxications underwent a forensic autopsy during the study period, 573 of them were caused by zopiclone. Among the zopiclone fatalities, 87% (n = 494) had a prescription for zopiclone, and 8% (n = 43) were monointoxications. Most fatalities, 62% (n = 354) were suicides, and zopiclone was involved in about 17% (n = 354) of all intoxication suicides in Sweden. Women were significantly (p < 0.01) overrepresented in suicides with zopiclone, comprising 56% (n = 291) of fatalities. The median age was 55 years among zopiclone intoxications compared with 44 years amongst all fatal intoxications. CONCLUSION: This study demonstrates that the toxicity of zopiclone can be lethal both in combination with other substances and on its own. Most individuals dying in fatal zopiclone intoxications were prescribed zopiclone, which potentially indicates that a more restrictive prescribing rate could prevent future intoxication deaths, especially when caring for patients with an increased suicide risk.

Granzyme K+ CD8 T cells in autoimmunity.

CD8 T cells expressing granzyme K are enriched in synovial tissue from patients with rheumatoid arthritis and in tissues affected by several other autoimmune diseases. The roles these cells play in autoimmune disease is under active investigation, and several recent studies have begun to shed light on this question. Putting this cell type into functional perspective is especially important given their enrichment at the sites of disease. This review summarizes available evidence for the presence of CD8 T cells and other granzyme K-expressing cells in tissues in autoimmune diseases and discusses the effects these cells may have on the pathogenesis of autoimmune conditions.

Associations Between Rheumatoid Arthritis Clinical Factors and Synovial Cell Types and States.

OBJECTIVE: Recent studies have uncovered diverse cell types and states in the rheumatoid arthritis (RA) synovium; however, limited data exist correlating these findings with patient-level clinical information. Using the largest cohort to date with clinical and multicell data, we determined associations between RA clinical factors with cell types and states in the RA synovium. METHODS: The Accelerated Medicines Partnership Rheumatoid Arthritis study recruited patients with active RA who were not receiving disease-modifying antirheumatic drugs (DMARDs) or who had an inadequate response to methotrexate (MTX) or tumor necrosis factor inhibitors. RA clinical factors were systematically collected. Biopsies were performed on an inflamed joint, and tissue were disaggregated and processed with a cellular indexing of transcriptomes and epitopes sequencing pipeline from which the following cell type percentages and cell type abundance phenotypes (CTAPs) were derived: endothelial, fibroblast, and myeloid (EFM); fibroblasts; myeloid; T and B cells; T cells and fibroblasts (TF); and T and myeloid cells. Correlations were measured between RA clinical factors, cell type percentage, and CTAPs. RESULTS: We studied 72 patients (mean age 57 years, 75% women, 83% seropositive, mean RA duration 6.6 years, mean Disease Activity Score-28 C-reactive Protein 3 [DAS28-CRP3] score 4.8). Higher DAS28-CRP3 correlated with a higher T cell percentage (P < 0.01). Those receiving MTX and not a biologic DMARD (bDMARD) had a higher percentage of B cells versus those receiving no DMARDs (P < 0.01). Most of those receiving bDMARDs were categorized as EFM (57%), whereas none were TF. No significant difference was observed across CTAPs for age, sex, RA disease duration, or DAS28-CRP3. CONCLUSION: In this comprehensive screen of clinical factors, we observed differential associations between DMARDs and cell phenotypes, suggesting that RA therapies, more than other clinical factors, may impact cell type/state in the synovium and ultimately influence response to subsequent therapies.

Synovial Tissue Insights into Heterogeneity of Rheumatoid Arthritis.

PURPOSE OF REVIEW: Rheumatoid arthritis is one of the most common rheumatic and autoimmune diseases. While it can affect many different organ systems, RA primarily involves inflammation in the synovium, the tissue that lines joints. Patients with RA exhibit significant clinical heterogeneity in terms of presence or absence of autoantibodies, degree of permanent deformities, and most importantly, treatment response. These clinical characteristics point to heterogeneity in the cellular and molecular pathogenesis of RA, an area that several recent studies have begun to address. RECENT FINDINGS: Single-cell RNA-sequencing initiatives and deeper focused studies have revealed several RA-associated cell populations in synovial tissues, including peripheral helper T cells, autoimmunity-associated B cells (ABCs), and NOTCH3+ sublining fibroblasts. Recent large transcriptional studies and translational clinical trials present frameworks to capture cellular and molecular heterogeneity in RA synovium. Technological developments, such as spatial transcriptomics and machine learning, promise to further elucidate the different types of RA synovitis and the biological mechanisms that characterize them, key elements of precision medicine to optimize patient care and outcomes in RA. This review recaps the findings of those recent studies and puts our current knowledge and future challenges into scientific and clinical perspective.

Altered glucagon and GLP-1 responses to oral glucose in children and adolescents with obesity and insulin resistance.

CONTEXT: Pediatric obesity is characterized by insulin resistance, yet it remains unclear whether insulin resistance contributes to abnormalities in glucagon and incretin secretion. OBJECTIVE: To examine whether fasting and stimulated glucagon, GLP-1, and GIP concentrations differ between children and adolescents with obesity and insulin resistance (OIR), obesity and normal insulin sensitivity (OIS), and controls with normal weight (NW). METHODS: 80 (34 boys) children and adolescents, aged 7-17 years with OIR (n=22), OIS (n=22), and NW (n=36) underwent an oral glucose tolerance test with measurements of serum insulin, plasma glucose, glucagon, total GLP-1, and total GIP. Homeostatic model assessment of insulin resistance (HOMA-IR), single point insulin sensitivity estimator (SPISE), Matsuda index, insulinogenic index (IGI), and oral disposition index (ODI) were calculated. RESULTS: Fasting concentrations of glucagon and GLP-1 were higher in the OIR-group, with no significant differences for GIP. The OIR-group had higher glucagon total area under the curve (tAUC0-120) and lower GLP-1 incremental AUC (iAUC0-120), with no significant differences for GIP iAUC0-120. Higher fasting glucagon was associated with higher HOMA-IR, lower Matsuda index, lower SPISE, higher IGI, and higher plasma alanine transaminase, whereas higher fasting GLP-1 was associated with higher HOMA-IR, lower Matsuda index, and lower ODI. Higher glucagon tAUC0-120 was associated lower SPISE and lower Matsuda index, whereas lower GLP-1 iAUC0-120 was associated with a higher HOMA-IR, lower Matsuda index, and lower ODI. CONCLUSIONS: The OIR-group had elevated fasting concentrations of glucagon and GLP-1, and higher glucagon, but lower GLP-1 responses during an OGTT compared to the OIS- and NW-groups. In contrast, the OIS-group had similar hormone responses to the NW-group.

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes.

Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction1. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments varies across patients, suggesting an undefined pathogenic diversity1,2. Here, to deconstruct the cell states and pathways that characterize this pathogenic heterogeneity, we profiled the full spectrum of cells in inflamed synovium from patients with rheumatoid arthritis. We used multi-modal single-cell RNA-sequencing and surface protein data coupled with histology of synovial tissue from 79 donors to build single-cell atlas of rheumatoid arthritis synovial tissue that includes more than 314,000 cells. We stratified tissues into six groups, referred to as cell-type abundance phenotypes (CTAPs), each characterized by selectively enriched cell states. These CTAPs demonstrate the diversity of synovial inflammation in rheumatoid arthritis, ranging from samples enriched for T and B cells to those largely lacking lymphocytes. Disease-relevant cell states, cytokines, risk genes, histology and serology metrics are associated with particular CTAPs. CTAPs are dynamic and can predict treatment response, highlighting the clinical utility of classifying rheumatoid arthritis synovial phenotypes. This comprehensive atlas and molecular, tissue-based stratification of rheumatoid arthritis synovial tissue reveal new insights into rheumatoid arthritis pathology and heterogeneity that could inform novel targeted treatments.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Genetics of Plasma Bilirubin and Associations between Bilirubin and Cardiometabolic Risk Profiles in Danish Children and Adolescents.

Bilirubin is the end product of heme catabolism, mainly produced by the breakdown of mature red blood cells. Due to its anti-inflammatory, antioxidant, antidiabetic, and antilipemic properties, circulating bilirubin concentrations are inversely associated with the risk of cardiovascular disease, type 2 diabetes, and all-cause mortality in adults. Some genetic loci associated with circulating bilirubin concentrations have been identified by genome-wide association studies in adults. We aimed to examine the relationship between circulating bilirubin, cardiometabolic risk factors, and inflammation in children and adolescents and the genetic architecture of plasma bilirubin concentrations. We measured fasting plasma bilirubin, cardiometabolic risk factors, and inflammatory markers in a sample of Danish children and adolescents with overweight or obesity (n = 1530) and in a population-based sample (n = 1820) of Danish children and adolescents. Linear and logistic regression analyses were performed to analyze the associations between bilirubin, cardiometabolic risk factors, and inflammatory markers. A genome-wide association study (GWAS) of fasting plasma concentrations of bilirubin was performed in children and adolescents with overweight or obesity and in a population-based sample. Bilirubin is associated inversely and significantly with a number of cardiometabolic risk factors, including body mass index (BMI) standard deviation scores (SDS), waist circumference, high-sensitivity C-reactive protein (hs-CRP), homeostatic model assessment for insulin resistance (HOMA-IR), hemoglobin A1c (HbA1c), low-density lipoprotein cholesterol (LDL-C), triglycerides, and the majority of measured inflammatory markers. In contrast, bilirubin was positively associated with fasting plasma concentrations of alanine transaminase (ALT), high-density lipoprotein cholesterol (HDL-C), systolic blood pressure (SDS), and the inflammatory markers GH, PTX3, THBS2, TNFRSF9, PGF, PAPPA, GT, CCL23, CX3CL1, SCF, and TRANCE. The GWAS showed that two loci were positively associated with plasma bilirubin concentrations at a p-value threshold of <5 × 10-8 (rs76999922: ß = -0.65 SD; p = 4.3 × 10-8, and rs887829: ß = 0.78 SD; p = 2.9 × 10-247). Approximately 25% of the variance in plasma bilirubin concentration was explained by rs887829. The rs887829 was not significantly associated with any of the mentioned cardiometabolic risk factors except for hs-CRP. Our findings suggest that plasma concentrations of bilirubin non-causally associates with cardiometabolic risk factors in children and adolescents.

Deep immunophenotyping reveals circulating activated lymphocytes in individuals at risk for rheumatoid arthritis.

Rheumatoid arthritis (RA) is a systemic autoimmune disease with currently no universally highly effective prevention strategies. Identifying pathogenic immune phenotypes in 'At-Risk' populations prior to clinical disease onset is crucial to establishing effective prevention strategies. Here, we applied mass cytometry to deeply characterize the immunophenotypes in blood from At-Risk individuals identified through the presence of serum antibodies to citrullinated protein antigens (ACPA) and/or first-degree relative (FDR) status (n=52), as compared to established RA (n=67), and healthy controls (n=48). We identified significant cell expansions in At-Risk individuals compared with controls, including CCR2+CD4+ T cells, T peripheral helper (Tph) cells, type 1 T helper cells, and CXCR5+CD8+ T cells. We also found that CD15+ classical monocytes were specifically expanded in ACPA-negative FDRs, and an activated PAX5 low naïve B cell population was expanded in ACPA-positive FDRs. Further, we developed an "RA immunophenotype score" classification method based on the degree of enrichment of cell states relevant to established RA patients. This score significantly distinguished At-Risk individuals from controls. In all, we systematically identified activated lymphocyte phenotypes in At-Risk individuals, along with immunophenotypic differences among both ACPA+ and ACPA-FDR At-Risk subpopulations. Our classification model provides a promising approach for understanding RA pathogenesis with the goal to further improve prevention strategies and identify novel therapeutic targets.

Genome-wide association study and functional characterization identifies candidate genes for insulin-stimulated glucose uptake.

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10-8) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits.

Adipocytes regulate fibroblast function, and their loss contributes to fibroblast dysfunction in inflammatory diseases.

Fibroblasts play critical roles in tissue homeostasis, but in pathologic states can drive fibrosis, inflammation, and tissue destruction. In the joint synovium, fibroblasts provide homeostatic maintenance and lubrication. Little is known about what regulates the homeostatic functions of fibroblasts in healthy conditions. We performed RNA sequencing of healthy human synovial tissue and identified a fibroblast gene expression program characterized by enhanced fatty acid metabolism and lipid transport. We found that fat-conditioned media reproduces key aspects of the lipid-related gene signature in cultured fibroblasts. Fractionation and mass spectrometry identified cortisol in driving the healthy fibroblast phenotype, confirmed using glucocorticoid receptor gene ( NR3C1 ) deleted cells. Depletion of synovial adipocytes in mice resulted in loss of the healthy fibroblast phenotype and revealed adipocytes as a major contributor to active cortisol generation via Hsd11 ß 1 expression. Cortisol signaling in fibroblasts mitigated matrix remodeling induced by TNFα- and TGFß, while stimulation with these cytokines repressed cortisol signaling and adipogenesis. Together, these findings demonstrate the importance of adipocytes and cortisol signaling in driving the healthy synovial fibroblast state that is lost in disease.

The Chromatin Landscape of Pathogenic Transcriptional Cell States in Rheumatoid Arthritis.

Synovial tissue inflammation is the hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. We measured genome-wide open chromatin at single cell resolution from 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identified 24 chromatin classes and predicted their associated transcription factors, including a CD8+ GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating an RA tissue transcriptional atlas, we found that the chromatin classes represented 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrated the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

Type 2 Diabetes risk alleles in Peptidyl-glycine Alpha-amidating Monooxygenase influence GLP-1 levels and response to GLP-1 Receptor Agonists.

Patients with type 2 diabetes vary in their response to currently available therapeutic agents (including GLP-1 receptor agonists) leading to suboptimal glycemic control and increased risk of complications. We show that human carriers of hypomorphic T2D-risk alleles in the gene encoding peptidyl-glycine alpha-amidating monooxygenase (PAM), as well as Pam-knockout mice, display increased resistance to GLP-1 in vivo. Pam inactivation in mice leads to reduced gastric GLP-1R expression and faster gastric emptying: this persists during GLP-1R agonist treatment and is rescued when GLP-1R activity is antagonized, indicating resistance to GLP-1's gastric slowing properties. Meta-analysis of human data from studies examining GLP-1R agonist response (including RCTs) reveals a relative loss of 44% and 20% of glucose lowering (measured by glycated hemoglobin) in individuals with hypomorphic PAM alleles p.S539W and p.D536G treated with GLP-1R agonist. Genetic variation in PAM has effects on incretin signaling that alters response to medication used commonly for treatment of T2D.

Convergence and divergence emerging in climatic controls of polynomial trends for northern ecosystem productivity over 2000-2018.

Recent rapid warming has caused uneven impacts on the composition, structure, and functioning of northern ecosystems. It remains unknown how climatic drivers control linear and non-linear trends in ecosystem productivity. Based on a plant phenology index (PPI) product at a spatial resolution of 0.05° over 2000-2018, we used an automated polynomial fitting scheme to detect and characterize trend types (i.e., polynomial trends and no-trends) in the yearly-integrated PPI (PPIINT) for northern (> 30°N) ecosystems and their dependence on climatic drivers and ecosystem types. The averaged slope for the linear trends (p < 0.05) of PPIINT was positive across all the ecosystems, among which deciduous broadleaved forests and evergreen needle-leaved forests (ENF) showed the highest and lowest mean slopes, respectively. More than 50% of the pixels in ENF, arctic and boreal shrublands, and permanent wetlands (PW) had linear trends. A large fraction of PW also showed quadratic and cubic trends. These trend patterns agreed well with estimates of global vegetation productivity based on solar-induced chlorophyll fluorescence. Across all the biomes, PPIINT in pixels with linear trends showed lower mean values and higher partial correlation coefficients with temperature or precipitation than in pixels without linear trends. Overall, our study revealed the emergence of latitudinal convergence and divergence in climatic controls on the linear and non-linear trends of PPIINT, implying that northern shifts of vegetation and climate change may potentially increase the non-linear nature of climatic controls on ecosystem productivity. These results can improve our understanding and prediction of climate-induced changes in plant phenology and productivity and facilitate sustainable management of ecosystems by accounting for their resilience and vulnerability to future climate change.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.