Human population history at the crossroads of East and Southeast Asia since 11,000 years ago.

Past human genetic diversity and migration between southern China and Southeast Asia have not been well characterized, in part due to poor preservation of ancient DNA in hot and humid regions. We sequenced 31 ancient genomes from southern China (Guangxi and Fujian), including two ∼12,000- to 10,000-year-old individuals representing the oldest humans sequenced from southern China. We discovered a deeply diverged East Asian ancestry in the Guangxi region that persisted until at least 6,000 years ago. We found that ∼9,000- to 6,000-year-old Guangxi populations were a mixture of local ancestry, southern ancestry previously sampled in Fujian, and deep Asian ancestry related to Southeast Asian Hòabìnhian hunter-gatherers, showing broad admixture in the region predating the appearance of farming. Historical Guangxi populations dating to ∼1,500 to 500 years ago are closely related to Tai-Kadai and Hmong-Mien speakers. Our results show heavy interactions among three distinct ancestries at the crossroads of East and Southeast Asia.

Catalytic Glycosylation for Minimally Protected Donors and Acceptors.

Oligosaccharides have myriad functions throughout biology.1,2 To investigate these functions requires multi-step chemical synthesis of these structurally complex molecules. With a dense concentration of stereocentres and hydroxyl groups, oligosaccharide assembly through O-glycosylation requires simultaneous control of site-, stereo-, and chemoselectivities3,4. Chemists have traditionally relied on protecting group manipulations for this purpose,5-8 adding a lot of synthetic work. Here, we report a glycosylation platform that enables selective coupling between unprotected or minimally protected donor and acceptor sugars, producing 1,2-cis-O-glycosides in a catalyst-controlled, site-selective manner. Radical-based activation9 of allyl glycosyl sulfones forms glycosyl bromides. A designed aminoboronic acid catalysts bring this reactive intermediate close to an acceptor through a network of noncovalent hydrogen bonding and reversible covalent B-O bonding interactions, allowing precise glycosyl transfer. The site of glycosylation can be switched with different aminoboronic acid catalysts by affecting their interaction modes with substrates. The method accommodates a wide range of sugar types, amenable to preparing naturally occurring sugar chains and pentasaccharides containing 11 free hydroxyls. Experimental and computational studies provide insights into the origin of selectivity outcomes.

Ultrathin quantum light source with van der Waals NbOCl2 crystal.

Interlayer electronic coupling in two-dimensional materials enables tunable and emergent properties by stacking engineering. However, it also results in significant evolution of electronic structures and attenuation of excitonic effects in two-dimensional semiconductors as exemplified by quickly degrading excitonic photoluminescence and optical nonlinearities in transition metal dichalcogenides when monolayers are stacked into van der Waals structures. Here we report a van der Waals crystal, niobium oxide dichloride (NbOCl2), featuring vanishing interlayer electronic coupling and monolayer-like excitonic behaviour in the bulk form, along with a scalable second-harmonic generation intensity of up to three orders higher than that in monolayer WS2. Notably, the strong second-order nonlinearity enables correlated parametric photon pair generation, through a spontaneous parametric down-conversion (SPDC) process, in flakes as thin as about 46 nm. To our knowledge, this is the first SPDC source unambiguously demonstrated in two-dimensional layered materials, and the thinnest SPDC source ever reported. Our work opens an avenue towards developing van der Waals material-based ultracompact on-chip SPDC sources as well as high-performance photon modulators in both classical and quantum optical technologies1-4.

Hidden hotspots of amphibian biodiversity in China.

Identifying and protecting hotspots of endemism and species richness is crucial for mitigating the global biodiversity crisis. However, our understanding of spatial diversity patterns is far from complete, which severely limits our ability to conserve biodiversity hotspots. Here, we report a comprehensive analysis of amphibian species diversity in China, one of the most species-rich countries on Earth. Our study combines 20 y of field surveys with new molecular analyses of 521 described species and also identifies 100 potential cryptic species. We identify 10 hotspots of amphibian diversity in China, each with exceptional species richness and endemism and with exceptional phylogenetic diversity and phylogenetic endemism (based on a new time-calibrated, species-level phylogeny for Chinese amphibians). These 10 hotspots encompass 59.6% of China's described amphibian species, 49.0% of cryptic species, and 55.6% of species endemic to China. Only four of these 10 hotspots correspond to previously recognized biodiversity hotspots. The six new hotspots include the Nanling Mountains and other mountain ranges in South China. Among the 186 species in the six new hotspots, only 9.7% are well covered by protected areas and most (88.2%) are exposed to high human impacts. Five of the six new hotspots are under very high human pressure and are in urgent need of protection. We also find that patterns of richness in cryptic species are significantly related to those in described species but are not identical.

The causes of sea-level rise since 1900.

The rate of global-mean sea-level rise since 1900 has varied over time, but the contributing factors are still poorly understood1. Previous assessments found that the summed contributions of ice-mass loss, terrestrial water storage and thermal expansion of the ocean could not be reconciled with observed changes in global-mean sea level, implying that changes in sea level or some contributions to those changes were poorly constrained2,3. Recent improvements to observational data, our understanding of the main contributing processes to sea-level change and methods for estimating the individual contributions, mean another attempt at reconciliation is warranted. Here we present a probabilistic framework to reconstruct sea level since 1900 using independent observations and their inherent uncertainties. The sum of the contributions to sea-level change from thermal expansion of the ocean, ice-mass loss and changes in terrestrial water storage is consistent with the trends and multidecadal variability in observed sea level on both global and basin scales, which we reconstruct from tide-gauge records. Ice-mass loss-predominantly from glaciers-has caused twice as much sea-level rise since 1900 as has thermal expansion. Mass loss from glaciers and the Greenland Ice Sheet explains the high rates of global sea-level rise during the 1940s, while a sharp increase in water impoundment by artificial reservoirs is the main cause of the lower-than-average rates during the 1970s. The acceleration in sea-level rise since the 1970s is caused by the combination of thermal expansion of the ocean and increased ice-mass loss from Greenland. Our results reconcile the magnitude of observed global-mean sea-level rise since 1900 with estimates based on the underlying processes, implying that no additional processes are required to explain the observed changes in sea level since 1900.

Phosphorylation of EZH2 by AMPK Suppresses PRC2 Methyltransferase Activity and Oncogenic Function.

Sustained energy starvation leads to activation of AMP-activated protein kinase (AMPK), which coordinates energy status with numerous cellular processes including metabolism, protein synthesis, and autophagy. Here, we report that AMPK phosphorylates the histone methyltransferase EZH2 at T311 to disrupt the interaction between EZH2 and SUZ12, another core component of the polycomb repressive complex 2 (PRC2), leading to attenuated PRC2-dependent methylation of histone H3 at Lys27. As such, PRC2 target genes, many of which are known tumor suppressors, were upregulated upon T311-EZH2 phosphorylation, which suppressed tumor cell growth both in cell culture and mouse xenografts. Pathologically, immunohistochemical analyses uncovered a positive correlation between AMPK activity and pT311-EZH2, and higher pT311-EZH2 correlates with better survival in both ovarian and breast cancer patients. Our finding suggests that AMPK agonists might be promising sensitizers for EZH2-targeting cancer therapies.

The RNA ligase RNA terminal phosphate cyclase B regulates mRNA alternative splicing and is required for mouse oocyte development and maintenance.

Recent large-scale mRNA sequencing has shown that introns are retained in 5-10% of mRNA, and these events are named intron retention (IR). IR has been recognized as a key mechanism in the regulation of gene expression. However, the role of this mechanism in female reproduction in mammals remains unclear. RNA terminal phosphate cyclase B (RTCB) is a RNA ligase; we found that RTCB conditional knockout mice have premature ovarian failure and that RTCB plays a crucial role in follicular development. RTCB regulated the splicing of transcripts related to DNA methylation and DNA damage repair. In addition, it regulated the resumption of oocyte meiosis by affecting CDK1 activation. Moreover, the loss of RTCB suppressed zygotic genome activation (ZGA) and decreased translation at the global level. In addition, Rtcb deletion resulted in the accumulation of maternal mRNAs containing unspliced introns and in a decline in the overall level of transcripts. As a result, the Rtcb-/- females were sterile. Our study highlights the important role of RTCB-regulated noncanonical alternative splicing in female reproduction.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization.

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and phosphorylated Ser-139 residue of the histone variant H2AX (γH2AX) foci linked to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-µm micropores. To search for pancreatic cancer cell-of-origin, analysis of single-cell data sets revealed that the extracellular matrix shaped an alternate route of acinar-ductal transdifferentiation of acinar cells into topoisomerase II α (TOP2A)-overexpressing cancer cells and derived subclusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in The Cancer Genome Atlas-Pancreatic Adenocarcinoma cohort. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased vav guanine nucleotide exchange factor 1 (VAV1) expression, and prolonged overall survival in the KPC mice. Reduction of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodeled the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induced the onset of cancer cells-of-origin by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes.

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Distinct functional patterns in child and adolescent bipolar and unipolar depression during emotional processing.

Accumulating evidence from functional magnetic resonance imaging studies supported brain dysfunction during emotional processing in bipolar disorder (BD) and major depressive disorder (MDD). However, child and adolescent BD and MDD could display different activation patterns, which have not been fully understood. This study aimed to investigate common and distinct activation patterns of pediatric BD (PBD) and MDD (p-MDD) during emotion processing using meta-analytic approaches. Literature search identified 25 studies, contrasting 252 PBD patients, and 253 healthy controls (HCs) as well as 311 p-MDD patients and 263 HCs. A total of nine meta-analyses were conducted pulling PBD and p-MDD experiments together and separately. The results revealed that PBD and p-MDD showed distinct patterns during negative processing. PBD patients exhibited activity changes in bilateral precuneus, left inferior parietal gyrus, left angular gyrus, and right posterior cingulate cortex while p-MDD patients showed functional disruptions in the left rectus, left triangular part of the inferior frontal gyrus, left orbital frontal cortex, left insula, and left putamen. In conclusion, the activity changes in PBD patients were mainly in regions correlated with emotion perception while the dysfunction among p-MDD patients was in the fronto-limbic circuit and reward-related regions in charge of emotion appraisal and regulation.

Cryo-electron microscopy structure and translocation mechanism of the crenarchaeal ribosome.

Archaeal ribosomes have many domain-specific features; however, our understanding of these structures is limited. We present 10 cryo-electron microscopy (cryo-EM) structures of the archaeal ribosome from crenarchaeota Sulfolobus acidocaldarius (Sac) at 2.7-5.7 Å resolution. We observed unstable conformations of H68 and h44 of ribosomal RNA (rRNA) in the subunit structures, which may interfere with subunit association. These subunit structures provided models for 12 rRNA expansion segments and 3 novel r-proteins. Furthermore, the 50S-aRF1 complex structure showed the unique domain orientation of aRF1, possibly explaining P-site transfer RNA (tRNA) release after translation termination. Sac 70S complexes were captured in seven distinct steps of the tRNA translocation reaction, confirming conserved structural features during archaeal ribosome translocation. In aEF2-engaged 70S ribosome complexes, 3D classification of cryo-EM data based on 30S head domain identified two new translocation intermediates with 30S head domain tilted 5-6° enabling its disengagement from the translocated tRNA and its release post-translocation. Additionally, we observed conformational changes to aEF2 during ribosome binding and switching from three different states. Our structural and biochemical data provide new insights into archaeal translation and ribosome translocation.

Archaeal ribosomes display variations in their ribosomal proteins and ribosomal RNA (rRNA) expansion segments (ESs). Protein translation in archaea combines features in both bacterial and eukaryotic translation. In this study, we present 10 cryo-electron microscopy structures of the archaeal ribosome from crenarchaeota Sulfolobus acidocaldarius (Sac). The 50S and 30S subunit structures present 3 novel ribosomal proteins and 12 rRNA ESs. The 70S Sac ribosome structures were captured in seven distinct functional states, including pre-, intermediate- and post-translocation states. Specifically, we identified two novel translocation intermediates, in which the 30S subunit head domain tilts outward to release the translocated P-site transfer RNA. The structures of archaeal ribosomes provide insights into the archaeal translation and ribosome translocation.

Chemical interference with DSIF complex formation lowers synthesis of mutant huntingtin gene products and curtails mutant phenotypes.

Earlier work has shown that siRNA-mediated reduction of the SUPT4H or SUPT5H proteins, which interact to form the DSIF complex and facilitate transcript elongation by RNA polymerase II (RNAPII), can decrease expression of mutant gene alleles containing nucleotide repeat expansions differentially. Using luminescence and fluorescence assays, we identified chemical compounds that interfere with the SUPT4H-SUPT5H interaction and then investigated their effects on synthesis of mRNA and protein encoded by mutant alleles containing repeat expansions in the huntingtin gene (HTT), which causes the inherited neurodegenerative disorder, Huntington's Disease (HD). Here we report that such chemical interference can differentially affect expression of HTT mutant alleles, and that a prototypical chemical, 6-azauridine (6-AZA), that targets the SUPT4H-SUPT5H interaction can modify the biological response to mutant HTT gene expression. Selective and dose-dependent effects of 6-AZA on expression of HTT alleles containing nucleotide repeat expansions were seen in multiple types of cells cultured in vitro, and in a Drosophila melanogaster animal model for HD. Lowering of mutant HD protein and mitigation of the Drosophila "rough eye" phenotype associated with degeneration of photoreceptor neurons in vivo were observed. Our findings indicate that chemical interference with DSIF complex formation can decrease biochemical and phenotypic effects of nucleotide repeat expansions.

Integrated Manipulation and Addressing of Spin Defect in Diamond.

Scalable and addressable integrated manipulation of qubits is crucial for practical quantum information applications. Different waveguides have been used to transport the optical and electrical driving pulses, which are usually required for qubit manipulation. However, the separated multifields may limit the compactness and efficiency of manipulation and introduce unwanted perturbation. Here, we develop a tapered fiber-nanowire-electrode hybrid structure to realize integrated optical and microwave manipulation of solid-state spins at nanoscale. Visible light and microwave driving pulses are simultaneously transported and concentrated along an Ag nanowire. Studied with spin defects in diamond, the results show that the different driving fields are aligned with high accuracy. The spatially selective spin manipulation is realized. And the frequency-scanning optically detected magnetic resonance (ODMR) of spin qubits is measured, illustrating the potential for portable quantum sensing. Our work provides a new scheme for developing compact, miniaturized quantum sensors and quantum information processing devices.

Cobalt-Catalyzed Enantioselective Reductive Arylation, Heteroarylation, and Alkenylation of Michael Acceptors via an Elementary Mechanism of 1,4-Addition.

Cobalt complexes with chiral quinox ligands effectively promote the enantioselective conjugate addition of enones using aryl, heteroaryl, and alkenyl halides and sulfonates. Additionally, a cobalt complex with a strongly donating diphosphine, BenzP*, successfully catalyzes the asymmetric reductive arylation and alkenylation of α,ß-unsaturated amides. Both catalytic systems show broad scopes and tolerance of sensitive functional groups. Both reactions can be scaled up with low loadings of cobalt catalysts. Experimental results and density functional theory (DFT) calculations suggest a new mechanism of elementary 1,4-addition of aryl cobalt(I) complexes.

Mild Stereoselective Synthesis of Densely Substituted [3]Dendralenes via Ru-Catalyzed Intermolecular Dimerization of 1,1-Disubstituted Allenes.

Described here is a mild and stereoselective protocol for the synthesis of [3]dendralenes via the intermolecular dimerization of allenes. With the proper choice of a ruthenium catalyst, a range of unactivated 1,1-disubstituted allenes, without prefunctionalization in the allylic position, reacted efficiently to provide rapid access to densely substituted [3]dendralenes. An intermolecular C-C bond and three different types of C═C double bonds (di-, tri-, and tetrasubstituted) embedded in an acyclic structure were constructed with good to high E/Z stereocontrol. This is in contrast to the known catalytic protocols that focus on allenes with prefunctionalization at the allylic position and/or monosubstituted allenes, which would proceed by a different mechanism or require less stereocontrol. The silyl-substituted dendralene products are precursors of other useful dendralene molecules. Density functional theory (DFT) studies and control experiments supported a mechanism involving oxidative cyclometalation, ß-H elimination (the rate-determining step), and reductive elimination.

BLMP-1 promotes developmental cell death in C. elegans by timely repression of ced-9 transcription.

Programmed cell death (PCD) is a common cell fate in metazoan development. PCD effectors are extensively studied, but how they are temporally regulated is less understood. Here, we report a mechanism controlling tail-spike cell death onset during Caenorhabditis elegans development. We show that the zinc-finger transcription factor BLMP-1, which controls larval development timing, also regulates embryonic tail-spike cell death initiation. BLMP-1 functions upstream of CED-9 and in parallel to DRE-1, another CED-9 and tail-spike cell death regulator. BLMP-1 expression is detected in the tail-spike cell shortly after the cell is born, and blmp-1 mutations promote ced-9-dependent tail-spike cell survival. BLMP-1 binds ced-9 gene regulatory sequences, and inhibits ced-9 transcription just before cell-death onset. BLMP-1 and DRE-1 function together to regulate developmental timing, and their mammalian homologs regulate B-lymphocyte fate. Our results, therefore, identify roles for developmental timing genes in cell-death initiation, and suggest conservation of these functions.

Discordance of Oncotype DX scores in synchronous bilateral and unilateral multifocal breast cancers.

PURPOSE: Oncotype DX, a 21-gene expression profiling test, has become standard of care in the management of estrogen receptor (ER)-positive breast cancer. In multifocal tumors, it is unclear whether testing of the different foci is necessary. We evaluated the concordance of Oncotype DX recurrence scores (RS) between 2 tumor foci in synchronous bilateral or unilateral multifocal tumors and characterized pathological predictors of discordance. METHODS: We reviewed 713 ER+, HER2- primary invasive breast cancer patients with Oncotype RS and identified 17 bilateral synchronous patients (34 tumors) and 13 unilateral multifocal patients (26 tumors) with available Oncotype RS on all foci. Discordance in Oncotype RS between synchronous tumors was recorded and associations with clinicopathologic features including tumor size, histology, Nottingham histologic grade, progesterone receptor staining, and Ki67 index were analyzed. RESULTS: Bilateral synchronous tumors were present in older patients (median age 59 years) and had larger tumor (median size 17 mm) and more discordant histology (10/17, 59%) as compared to unilateral multifocal tumors (median age 49 years, p < 0.01; median tumor size 12 mm, p = 0.01; discordant histology 2/13, 15%, p = 0.03). Oncotype RS were discordant in 47% (8/17) of bilateral and 54% (7/13) of unilateral multifocal tumors. Concordant Oncotype RS was associated with similar histologic grade and Ki67 index in 78% (7/9) of bilateral and 100% (6/6) of multifocal tumors. In contrast, only 25% (2/8) of bilateral (p = 0.06) and 14% (1/7) of unilateral multifocal (p < 0.01) cases with discordant Oncotype RS had concordant histology grades and Ki67 levels. In synchronous tumors with discordant Oncotype RS and Ki67 index, all (4/4) foci with higher RS had higher Ki67 index. CONCLUSION: Discordance of Oncotype RS is common in both bilateral and unilateral multifocal breast cancer and is likely associated with discordant histologic grade or Ki67.

Immunohistochemical assessment of HER2 low breast cancer: interobserver reproducibility and correlation with digital image analysis.

PURPOSE: The recent findings from the DESTINY-Breast04 trial highlighted the clinical importance of distinguishing between HER2 immunohistochemistry (IHC) scores 0 and 1 + in metastatic breast cancer (BC). However, pathologist interpretation of HER2 IHC scoring is subjective, and standardized methodology is needed. We evaluated the consistency of HER2 IHC scoring among pathologists and the accuracy of digital image analysis (DIA) in interpreting HER2 IHC staining in cases of HER2-low BC. METHODS: Fifty whole-slide biopsies of BC with HER2 IHC staining were evaluated, comprising 25 cases originally reported as IHC score 0 and 25 as 1 +. These slides were digitally scanned. Six pathologists with breast expertise independently reviewed and scored the scanned images, and DIA was applied. Agreement among pathologists and concordance between pathologist scores and DIA results were statistically analyzed using Kendall coefficient of concordance (W) tests. RESULTS: Substantial agreement among at least five of the six pathologists was found for 18 of the score 0 cases (72%) and 15 of the score 1 + cases (60%), indicating excellent interobserver agreement (W = 0.828). DIA scores were highly concordant with pathologist scores in 96% of cases (47/49), indicating excellent concordance (W = 0.959). CONCLUSION: Although breast subspecialty pathologists were relatively consistent in evaluating BC with HER2 IHC scores of 0 and 1 +, DIA may be a reliable supplementary tool to enhance the standardization and quantification of HER2 IHC assessment, especially in challenging cases where results may be ambiguous (i.e., scores 0-1 +). These findings hold promise for improving the accuracy and consistency of HER2 testing.

A functional study reveals CsNAC086 regulated the biosynthesis of flavonols in Camellia sinensis.

MAIN CONCLUSION: CsNAC086 was found to promote the expression of CsFLS, thus promoting the accumulation of flavonols in Camellia sinensis. Flavonols, the main flavonoids in tea plants, play an important role in the taste and quality of tea. In this study, a NAC TF gene CsNAC086 was isolated from tea plants and confirmed its regulatory role in the expression of flavonol synthase which is a key gene involved in the biosynthesis of flavonols in tea plant. Yeast transcription-activity assays showed that CsNAC086 has self-activation activity. The transcriptional activator domain of CsNAC086 is located in the non-conserved C-terminal region (positions 171-550), while the conserved NAC domain (positions 1-170) does not have self-activation activity. Silencing the CsNAC086 gene using antisense oligonucleotides significantly decreased the expression of CsFLS. As a result, the concentration of flavonols decreased significantly. In overexpressing CsNAC086 tobacco leaves, the expression of NtFLS was significantly increased. Compared with wild-type tobacco, the flavonols concentration increased. Yeast one-hybrid assays showed CsNAC086 did not directly regulate the gene expression of CsFLS. These findings indicate that CsNAC086 plays a role in regulating flavonols biosynthesis in tea plants, which has important implications for selecting and breeding of high-flavonols-concentration containing tea-plant cultivars.

Paranasal sinus occupancy assessed from magnetic resonance images-associations with clinical indicators in patients with systemic lupus erythematosus.

OBJECTIVES: Nasal, paranasal sinus and mucosal disorders are common symptoms in autoimmune rheumatic diseases. Soft tissue changes and fluid accumulation in the osteomeatal complexes and paranasal sinuses manifest as opaqueness on radiological images which can be assessed using visual scoring and computational methods on CT scans, but their results do not always correlate. Using MRI, we investigate the applicability of different image analysis methods in SLE. METHODS: We assessed paranasal sinus opaqueness on MRI from 51 SLE patients, using three visual scoring systems and expert-delineated computational volumes, and examined their association with markers of disease activity, inflammation, endothelial dysfunction and common small vessel disease (SVD) indicators, adjusting for age and sex-at-birth. RESULTS: The average paranasal sinus volume occupation was 4.55 (6.47%) [median (interquartile range) = 0.67 (0.25-2.65) ml], mainly in the maxillary and ethmoid sinuses. It was highly correlated with Lund-Mackay (LM) scores modified at 50% opaqueness cut-off (Spearman's ρ: 0.71 maxillary and 0.618 ethmoids, P < 0.001 in all), and with more granular variations of the LM system. The modified LM scores were associated with SVD scores (0: B = 5.078, s.e. = 1.69, P = 0.0026; 2: B = -0.066, s.e. = 0.023, P = 0.0045) and disease activity (anti-dsDNA: B = 4.59, s.e. = 2.22, P = 0.045; SLEDAI 3-7: 2.86 < B < 4.30; 1.38 < s.e. < 1.63; 0.0083 ≤ P ≤ 0.0375). Computationally derived percent opaqueness yielded similar results. CONCLUSION: In patients with SLE, MRI computational assessment of sinuses opaqueness and LM scores modified at a 50% cut-off may be useful tools in understanding the relationships among paranasal sinus occupancy, disease activity and SVD markers.