Functional Genetic Variants Revealed by Massively Parallel Precise Genome Editing.

A major challenge in genetics is to identify genetic variants driving natural phenotypic variation. However, current methods of genetic mapping have limited resolution. To address this challenge, we developed a CRISPR-Cas9-based high-throughput genome editing approach that can introduce thousands of specific genetic variants in a single experiment. This enabled us to study the fitness consequences of 16,006 natural genetic variants in yeast. We identified 572 variants with significant fitness differences in glucose media; these are highly enriched in promoters, particularly in transcription factor binding sites, while only 19.2% affect amino acid sequences. Strikingly, nearby variants nearly always favor the same parent's alleles, suggesting that lineage-specific selection is often driven by multiple clustered variants. In sum, our genome editing approach reveals the genetic architecture of fitness variation at single-base resolution and could be adapted to measure the effects of genome-wide genetic variation in any screen for cell survival or cell-sortable markers.

HOTTIP-dependent R-loop formation regulates CTCF boundary activity and TAD integrity in leukemia.

HOTTIP lncRNA is highly expressed in acute myeloid leukemia (AML) driven by MLL rearrangements or NPM1 mutations to mediate HOXA topologically associated domain (TAD) formation and drive aberrant transcription. However, the mechanism through which HOTTIP accesses CCCTC-binding factor (CTCF) chromatin boundaries and regulates CTCF-mediated genome topology remains unknown. Here, we show that HOTTIP directly interacts with and regulates a fraction of CTCF-binding sites (CBSs) in the AML genome by recruiting CTCF/cohesin complex and R-loop-associated regulators to form R-loops. HOTTIP-mediated R-loops reinforce the CTCF boundary and facilitate formation of TADs to drive gene transcription. Either deleting CBS or targeting RNase H to eliminate R-loops in the boundary CBS of ß-catenin TAD impaired CTCF boundary activity, inhibited promoter/enhancer interactions, reduced ß-catenin target expression, and mitigated leukemogenesis in xenograft mouse models with aberrant HOTTIP expression. Thus, HOTTIP-mediated R-loop formation directly reinforces CTCF chromatin boundary activity and TAD integrity to drive oncogene transcription and leukemia development.

Spatially resolved land and grid model of carbon neutrality in China.

China has committed to achieve net carbon neutrality by 2060 to combat global climate change, which will require unprecedented deployment of negative emissions technologies, renewable energies (RE), and complementary infrastructure. At terawatt-scale deployment, land use limitations interact with operational and economic features of power systems. To address this, we developed a spatially resolved resource assessment and power systems planning optimization that models a full year of power system operations, sub-provincial RE siting criteria, and transmission connections. Our modeling results show that wind and solar must be expanded to 2,000 to 3,900 GW each, with one plausible pathway leading to 300 GW/yr combined annual additions in 2046 to 2060, a three-fold increase from today. Over 80% of solar and 55% of wind is constructed within 100 km of major load centers when accounting for current policies regarding land use. Large-scale low-carbon systems must balance key trade-offs in land use, RE resource quality, grid integration, and costs. Under more restrictive RE siting policies, at least 740 GW of distributed solar would become economically feasible in regions with high demand, where utility-scale deployment is limited by competition with agricultural land. Effective planning and policy formulation are necessary to achieve China's climate goals.

Kinetic pathway of HIV-1 TAR cotranscriptional folding.

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

lncRNA BREA2 promotes metastasis by disrupting the WWP2-mediated ubiquitination of Notch1.

Notch has been implicated in human cancers and is a putative therapeutic target. However, the regulation of Notch activation in the nucleus remains largely uncharacterized. Therefore, characterizing the detailed mechanisms governing Notch degradation will identify attractive strategies for treating Notch-activated cancers. Here, we report that the long noncoding RNA (lncRNA) BREA2 drives breast cancer metastasis by stabilizing the Notch1 intracellular domain (NICD1). Moreover, we reveal WW domain containing E3 ubiquitin protein ligase 2 (WWP2) as an E3 ligase for NICD1 at K1821 and a suppressor of breast cancer metastasis. Mechanistically, BREA2 impairs WWP2-NICD1 complex formation and in turn stabilizes NICD1, leading to Notch signaling activation and lung metastasis. BREA2 loss sensitizes breast cancer cells to inhibition of Notch signaling and suppresses the growth of breast cancer patient-derived xenograft tumors, highlighting its therapeutic potential in breast cancer. Taken together, these results reveal the lncRNA BREA2 as a putative regulator of Notch signaling and an oncogenic player driving breast cancer metastasis.

Multiple models for outbreak decision support in the face of uncertainty.

Policymakers must make management decisions despite incomplete knowledge and conflicting model projections. Little guidance exists for the rapid, representative, and unbiased collection of policy-relevant scientific input from independent modeling teams. Integrating approaches from decision analysis, expert judgment, and model aggregation, we convened multiple modeling teams to evaluate COVID-19 reopening strategies for a mid-sized United States county early in the pandemic. Projections from seventeen distinct models were inconsistent in magnitude but highly consistent in ranking interventions. The 6-mo-ahead aggregate projections were well in line with observed outbreaks in mid-sized US counties. The aggregate results showed that up to half the population could be infected with full workplace reopening, while workplace restrictions reduced median cumulative infections by 82%. Rankings of interventions were consistent across public health objectives, but there was a strong trade-off between public health outcomes and duration of workplace closures, and no win-win intermediate reopening strategies were identified. Between-model variation was high; the aggregate results thus provide valuable risk quantification for decision making. This approach can be applied to the evaluation of management interventions in any setting where models are used to inform decision making. This case study demonstrated the utility of our approach and was one of several multimodel efforts that laid the groundwork for the COVID-19 Scenario Modeling Hub, which has provided multiple rounds of real-time scenario projections for situational awareness and decision making to the Centers for Disease Control and Prevention since December 2020.

Decrotonylation of cGAS K254 prompts homologous recombination repair by blocking its DNA binding and releasing PARP1.

Cyclic GMP-AMP synthase (cGAS), a cytosolic DNA sensor, also exhibits nuclear genomic localization and is involved in DNA damage signaling. In this study, we investigated the impact of cGAS crotonylation on the regulation of the DNA damage response, particularly homologous recombination repair, following exposure to ionizing radiation (IR). Lysine 254 of cGAS is constitutively crotonylated by the CREB-binding protein; however, IR-induced DNA damage triggers SIRT3-mediated decrotonylation. Lysine 254 decrotonylation decreased the DNA-binding affinity of cGAS and inhibited its interaction with PARP1, promoting HR repair. Moreover, SIRT3 suppression led to HR repair inhibition and markedly sensitized cancer cells to IR and DNA-damaging chemicals, highlighting SIRT3 as a potential target for cancer therapy. Overall, this study revealed the crucial role of cGAS crotonylation in the DNA damage response. Furthermore, we propose that modulating cGAS and SIRT3 activities could be potential strategies for cancer therapy.

DOCK2 Promotes Atherosclerosis by Mediating the Endothelial Cell Inflammatory Response.

The pathology of atherosclerosis, a leading cause of mortality in patients with cardiovascular disease, involves inflammatory phenotypic changes in vascular endothelial cells. This study explored the role of the dedicator of cytokinesis (DOCK)-2 protein in atherosclerosis. Mice with deficiencies in low-density lipoprotein receptor and Dock2 (Ldlr-/-Dock2-/-) and controls (Ldlr-/-) were fed a high-fat diet (HFD) to induce atherosclerosis. In controls, Dock2 was increased in atherosclerotic lesions, with increased intercellular adhesion molecule (Icam)-1 and vascular cell adhesion molecule (Vcam)-1, after HFD for 4 weeks. Ldlr-/-Dock2-/- mice exhibited significantly decreased oil red O staining in both aortic roots and aortas compared to that in controls after HFD for 12 weeks. In control mice and in humans, Dock2 was highly expressed in the ECs of atherosclerotic lesions. Dock2 deficiency was associated with attenuation of Icam-1, Vcam-1, and monocyte chemoattractant protein (Mcp)-1 in the aortic roots of mice fed HFD. Findings in human vascular ECs in vitro suggested that DOCK2 was required in TNF-α-mediated expression of ICAM-1/VCAM-1/MCP-1. DOCK2 knockdown was associated with attenuated NF-κB phosphorylation with TNF-α, partially accounting for DOCK2-mediated vascular inflammation. With DOCK2 knockdown in human vascular ECs, TNF-α-mediated VCAM-1 promoter activity was inhibited. The findings from this study suggest the novel concept that DOCK2 promotes the pathogenesis of atherosclerosis by modulating inflammation in vascular ECs.

ADAR1 Non-Editing Function in Macrophage Activation and Abdominal Aortic Aneurysm.

BACKGROUND: Macrophage activation plays a critical role in abdominal aortic aneurysm (AAA) development. However, molecular mechanisms controlling macrophage activation and vascular inflammation in AAA remain largely unknown. The objective of the study was to identify novel mechanisms underlying adenosine deaminase acting on RNA (ADAR1) function in macrophage activation and AAA formation. METHODS: Aortic transplantation was conducted to determine the importance of nonvascular ADAR1 in AAA development/dissection. Ang II (Angiotensin II) infusion of ApoE-/- mouse model combined with macrophage-specific knockout of ADAR1 was used to study ADAR1 macrophage-specific role in AAA formation/dissection. The relevance of macrophage ADAR1 to human AAA was examined using human aneurysm specimens. Moreover, a novel humanized AAA model was established to test the role of human macrophages in aneurysm formation in human arteries. RESULTS: Allograft transplantation of wild-type abdominal aortas to ADAR1+/- recipient mice significantly attenuated AAA formation, suggesting that nonvascular ADAR1 is essential for AAA development. ADAR1 deficiency in hematopoietic cells decreased the prevalence and severity of AAA while inhibited macrophage infiltration and aorta wall inflammation. ADAR1 deletion blocked the classic macrophage activation, diminished NF-κB (nuclear factor kappa B) signaling, and enhanced the expression of a number of anti-inflammatory microRNAs. Mechanistically, ADAR1 interacted with Drosha to promote its degradation, which attenuated Drosha-DGCR8 (DiGeorge syndrome critical region 8) interaction, and consequently inhibited pri- to pre-microRNA processing of microRNAs targeting IKKß, resulting in an increased IKKß (inhibitor of nuclear factor kappa-B) expression and enhanced NF-κB signaling. Significantly, ADAR1 was induced in macrophages and interacted with Drosha in human AAA lesions. Reconstitution of ADAR1-deficient, but not the wild type, human monocytes to immunodeficient mice blocked the aneurysm formation in transplanted human arteries. CONCLUSIONS: Macrophage ADAR1 promotes aneurysm formation in both mouse and human arteries through a novel mechanism, that is, Drosha protein degradation, which inhibits the processing of microRNAs targeting NF-kB signaling and thus elicits macrophage-mediated vascular inflammation in AAA.

Loss of BRD4 induces cell senescence in HSC/HPCs by deregulating histone H3 clipping.

Bromodomain-containing protein 4 (BRD4) is overexpressed and functionally implicated in various myeloid malignancies. However, the role of BRD4 in normal hematopoiesis remains largely unknown. Here, utilizing an inducible Brd4 knockout mouse model, we find that deletion of Brd4 (Brd4Δ/Δ ) in the hematopoietic system impairs hematopoietic stem cell (HSC) self-renewal and differentiation, which associates with cell cycle arrest and senescence. ATAC-seq analysis shows increased chromatin accessibility in Brd4Δ/Δ hematopoietic stem/progenitor cells (HSC/HPCs). Genome-wide mapping with cleavage under target and release using nuclease (CUT&RUN) assays demonstrate that increased global enrichment of H3K122ac and H3K4me3 in Brd4Δ/Δ HSC/HPCs is associated with the upregulation of senescence-specific genes. Interestingly, Brd4 deletion increases clipped H3 (cH3) which correlates with the upregulation of senescence-specific genes and results in a higher frequency of senescent HSC/HPCs. Re-expression of BRD4 reduces cH3 levels and rescues the senescence rate in Brd4Δ/Δ HSC/HPCs. This study unveils an important role of BRD4 in HSC/HPC function by preventing H3 clipping and suppressing senescence gene expression.

PAI-1 Regulates the Cytoskeleton and Intrinsic Stiffness of Vascular Smooth Muscle Cells.

BACKGROUND: Plasma concentration of PAI-1 (plasminogen activator inhibitor-1) correlates with arterial stiffness. Vascular smooth muscle cells (SMCs) express PAI-1, and the intrinsic stiffness of SMCs is a major determinant of total arterial stiffness. We hypothesized that PAI-1 promotes SMC stiffness by regulating the cytoskeleton and that pharmacological inhibition of PAI-1 decreases SMC and aortic stiffness. METHODS: PAI-039, a specific inhibitor of PAI-1, and small interfering RNA were used to inhibit PAI-1 expression in cultured human SMCs. Effects of PAI-1 inhibition on SMC stiffness, F-actin (filamentous actin) content, and cytoskeleton-modulating enzymes were assessed. WT (wild-type) and PAI-1-deficient murine SMCs were used to determine PAI-039 specificity. RNA sequencing was performed to determine the effects of PAI-039 on SMC gene expression. In vivo effects of PAI-039 were assessed by aortic pulse wave velocity. RESULTS: PAI-039 significantly reduced intrinsic stiffness of human SMCs, which was accompanied by a significant decrease in cytoplasmic F-actin content. PAI-1 gene knockdown also decreased cytoplasmic F-actin. PAI-1 inhibition significantly increased the activity of cofilin, an F-actin depolymerase, in WT murine SMCs, but not in PAI-1-deficient SMCs. RNA-sequencing analysis suggested that PAI-039 upregulates AMPK (AMP-activated protein kinase) signaling in SMCs, which was confirmed by Western blotting. Inhibition of AMPK prevented activation of cofilin by PAI-039. In mice, PAI-039 significantly decreased aortic stiffness and tunica media F-actin content without altering the elastin or collagen content. CONCLUSIONS: PAI-039 decreases intrinsic SMC stiffness and cytoplasmic stress fiber content. These effects are mediated by AMPK-dependent activation of cofilin. PAI-039 also decreases aortic stiffness in vivo. These findings suggest that PAI-1 is an important regulator of the SMC cytoskeleton and that pharmacological inhibition of PAI-1 has the potential to prevent and treat cardiovascular diseases involving arterial stiffening.

Glioma-associated mesenchymal stem cells.

Recent studies have revealed that glioma-associated mesenchymal stem cells play instrumental roles in tumorigenesis and tumour progression and cannot be ignored as a cellular component of the glioma microenvironment. Nevertheless, the origin of these cells and their roles are poorly understood. The only relevant studies have shown that glioma-associated mesenchymal stem cells play a large role in promoting tumour proliferation, invasion and angiogenesis. This review provides a comprehensive summary of their discovery and definition, origin, differences from other tissue-derived mesenchymal stem cells, spatial distribution, functions and prognostic and therapeutic opportunities to deepen the understanding of these cells and provide new insight into the treatment of glioma.

RNAJP: enhanced RNA 3D structure predictions with non-canonical interactions and global topology sampling.

RNA 3D structures are critical for understanding their functions. However, only a limited number of RNA structures have been experimentally solved, so computational prediction methods are highly desirable. Nevertheless, accurate prediction of RNA 3D structures, especially those containing multiway junctions, remains a significant challenge, mainly due to the complicated non-canonical base pairing and stacking interactions in the junction loops and the possible long-range interactions between loop structures. Here we present RNAJP ('RNA Junction Prediction'), a nucleotide- and helix-level coarse-grained model for the prediction of RNA 3D structures, particularly junction structures, from a given 2D structure. Through global sampling of the 3D arrangements of the helices in junctions using molecular dynamics simulations and in explicit consideration of non-canonical base pairing and base stacking interactions as well as long-range loop-loop interactions, the model can provide significantly improved predictions for multibranched junction structures than existing methods. Moreover, integrated with additional restraints from experiments, such as junction topology and long-range interactions, the model may serve as a useful structure generator for various applications.

An Overlooked and Underrated Endemic Mycosis-Talaromycosis and the Pathogenic Fungus Talaromyces marneffei.

Talaromycosis is an invasive mycosis endemic in tropical and subtropical Asia and is caused by the pathogenic fungus Talaromyces marneffei. Approximately 17,300 cases of T. marneffei infection are diagnosed annually, and the reported mortality rate is extremely high (~1/3). Despite the devastating impact of talaromycosis on immunocompromised individuals, particularly HIV-positive persons, and the increase in reported occurrences in HIV-uninfected persons, diagnostic and therapeutic approaches for talaromycosis have received far too little attention worldwide. In 2021, scientists living in countries where talaromycosis is endemic raised a global demand for it to be recognized as a neglected tropical disease. Therefore, T. marneffei and the infectious disease induced by this fungus must be treated with concern. T. marneffei is a thermally dimorphic saprophytic fungus with a complicated mycological growth process that may produce various cell types in its life cycle, including conidia, hyphae, and yeast, all of which are associated with its pathogenicity. However, understanding of the pathogenic mechanism of T. marneffei has been limited until recently. To achieve a holistic view of T. marneffei and talaromycosis, the current knowledge about talaromycosis and research breakthroughs regarding T. marneffei growth biology are discussed in this review, along with the interaction of the fungus with environmental stimuli and the host immune response to fungal infection. Importantly, the future research directions required for understanding this serious infection and its causative pathogenic fungus are also emphasized to identify solutions that will alleviate the suffering of susceptible individuals worldwide.

Machine learning in RNA structure prediction: Advances and challenges.

RNA molecules play a crucial role in various biological processes, with their functionality closely tied to their structures. The remarkable advancements in machine learning techniques for protein structure prediction have shown promise in the field of RNA structure prediction. In this perspective, we discuss the advances and challenges encountered in constructing machine learning-based models for RNA structure prediction. We explore topics including model building strategies, specific challenges involved in predicting RNA secondary (2D) and tertiary (3D) structures, and approaches to these challenges. In addition, we highlight the advantages and challenges of constructing RNA language models. Given the rapid advances of machine learning techniques, we anticipate that machine learning-based models will serve as important tools for predicting RNA structures, thereby enriching our understanding of RNA structures and their corresponding functions.

Profiling cardiomyocytes at single cell resolution reveals COX7B could be a potential target for attenuating heart failure in cardiac hypertrophy.

Cardiac hypertrophy can develop to end-stage heart failure (HF), which inevitably leading to heart transplantation or death. Preserving cardiac function in cardiomyocytes (CMs) is essential for improving prognosis in hypertrophic cardiomyopathy (HCM) patients. Therefore, understanding transcriptomic heterogeneity of CMs in HCM would be indispensable to aid potential therapeutic targets investigation. We isolated primary CM from HCM patients who had extended septal myectomy, and obtained transcriptomes in 338 human primary CM with single-cell tagged reverse transcription (STRT-seq) approach. Our results revealed that CMs could be categorized into three subsets in nonfailing HCM heart: high energy synthesis cluster, high cellular metabolism cluster and intermediate cluster. The expression of electron transport chain (ETC) was up-regulated in larger-sized CMs from high energy synthesis cluster. Of note, we found the expression of Cytochrome c oxidase subunit 7B (COX7B), a subunit of Complex IV in ETC had trends of positively correlation with CMs size. Further, by assessing COX7B expression in HCM patients, we speculated that COX7B was compensatory up-regulated at early-stage but down-regulated in failing HCM heart. To test the hypothesis that COX7B might participate both in hypertrophy and HF progression, we used adeno associated virus 9 (AAV9) to mediate the expression of Cox7b in pressure overload-induced mice. Mice in vivo data supported that knockdown of Cox7b would accelerate HF and Cox7b overexpression could restore partial cardiac function in hypertrophy. Our result highlights targeting COX7B and preserving energy synthesis in hypertrophic CMs could be a promising translational direction for HF therapeutic strategy.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family.

BACKGROUND: Terpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated. RESULTS: We identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality. CONCLUSIONS: We systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Evaluating genomic inbreeding of two Chinese yak (Bos grunniens) populations.

BACKGROUND: Yaks are a vital livestock in the Qinghai-Tibetan Plateau area for providing food products, maintaining sustainable ecosystems, and promoting cultural heritage. Because of uncontrolled mating, it is impossible to estimate inbreeding level of yak populations using the pedigree-based approaches. With the aims to accurately evaluate inbreeding level of two Chinese yak populations (Maiwa and Jiulong), we obtained genome-wide single nucleotide polymorphisms (SNPs) by DNA sequencing and calculated five SNP-by-SNP estimators ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]), as well as two segment-based estimators of runs of homozygosity (ROH, [Formula: see text]) and homozygous-by-descent (HBD, [Formula: see text]). Functional implications were analyzed for the positional candidate genes located within the related genomic regions. RESULTS: A total of 151,675 and 190,955 high-quality SNPs were obtained from 71 Maiwa and 30 Jiulong yaks, respectively. Jiulong had greater genetic diversity than Maiwa in terms of allele frequency and nucleotide diversity. The two populations could be genetically distinguished by principal component analysis, with the mean differentiation index (Fst) of 0.0054. The greater genomic inbreeding levels of Maiwa yaks were consistently supported by all five SNP-by-SNP estimators. Based on simple proportion of homozygous SNPs ([Formula: see text]), a lower inbreeding level was indicated by three successfully sequenced old leather samples that may represent historical Maiwa yaks about five generations ago. There were 3304 ROH detected among all samples, with mean and median length of 1.97 Mb and 1.0 Mb, respectively. A total of 94 HBD segments were found among all samples, whereas 92 of them belonged to the shortest class with the mean length of 10.9 Kb. Based on the estimates of [Formula: see text] and [Formula: see text], however, there was no difference in inbreeding level between Maiwa and Jiulong yaks. Within the genomic regions with the significant Fst or enriched by ROH, we found several candidate genes and pathways that have been reported to be related to diverse production traits in farm animals. CONCLUSIONS: We successfully evaluated the genomic inbreeding level of two Chinese yak populations. Although different estimators resulted in inconsistent conclusions on their genomic inbreeding levels, our results may be helpful to implement the genetic conservation and utilization programs for the two yak populations.

Potential impact of annual vaccination with reformulated COVID-19 vaccines: Lessons from the US COVID-19 scenario modeling hub.

BACKGROUND: Coronavirus Disease 2019 (COVID-19) continues to cause significant hospitalizations and deaths in the United States. Its continued burden and the impact of annually reformulated vaccines remain unclear. Here, we present projections of COVID-19 hospitalizations and deaths in the United States for the next 2 years under 2 plausible assumptions about immune escape (20% per year and 50% per year) and 3 possible CDC recommendations for the use of annually reformulated vaccines (no recommendation, vaccination for those aged 65 years and over, vaccination for all eligible age groups based on FDA approval). METHODS AND FINDINGS: The COVID-19 Scenario Modeling Hub solicited projections of COVID-19 hospitalization and deaths between April 15, 2023 and April 15, 2025 under 6 scenarios representing the intersection of considered levels of immune escape and vaccination. Annually reformulated vaccines are assumed to be 65% effective against symptomatic infection with strains circulating on June 15 of each year and to become available on September 1. Age- and state-specific coverage in recommended groups was assumed to match that seen for the first (fall 2021) COVID-19 booster. State and national projections from 8 modeling teams were ensembled to produce projections for each scenario and expected reductions in disease outcomes due to vaccination over the projection period. From April 15, 2023 to April 15, 2025, COVID-19 is projected to cause annual epidemics peaking November to January. In the most pessimistic scenario (high immune escape, no vaccination recommendation), we project 2.1 million (90% projection interval (PI) [1,438,000, 4,270,000]) hospitalizations and 209,000 (90% PI [139,000, 461,000]) deaths, exceeding pre-pandemic mortality of influenza and pneumonia. In high immune escape scenarios, vaccination of those aged 65+ results in 230,000 (95% confidence interval (CI) [104,000, 355,000]) fewer hospitalizations and 33,000 (95% CI [12,000, 54,000]) fewer deaths, while vaccination of all eligible individuals results in 431,000 (95% CI: 264,000-598,000) fewer hospitalizations and 49,000 (95% CI [29,000, 69,000]) fewer deaths. CONCLUSIONS: COVID-19 is projected to be a significant public health threat over the coming 2 years. Broad vaccination has the potential to substantially reduce the burden of this disease, saving tens of thousands of lives each year.

Stratifying ICIs-responsive tumor microenvironment in HCC: from parsing out immune-hypoxic crosstalk to clinically applicable MRI-radiomics models.

BACKGROUND: We aimed to redefine Immune checkpoint inhibitors (ICIs)-responsive "hot" TME and develop a corresponding stratification model to maximize ICIs-efficacy in Hepatocellular Carcinoma (HCC). METHODS: Hypoxic scores were designed, and the relevance to immunotherapy responses were validated in pan-cancers through single cell analysis. Multi-omics analysis using the hypoxic scores and immune infiltrate abundance was performed to redefine the ICIs-responsive TME subtype in HCC patients from TCGA (n = 363) and HCCDB database (n = 228). The immune hypoxic stress index (IHSI) was constructed to stratify the ICIs-responsive TME subtype, with exploring biological mechanism in vitro and in vivo. MRI-radiomics models were built for clinical applicability. RESULTS: The hypoxic scores were lower in the dominant cell-subclusters of responders in pan-cancers. The higher immune infiltrate-normoxic (HIN) subtype was redefined as the ICIs-responsive TME. Stratification of the HIN subtype using IHSI effectively identified ICIs-responders in Melanoma (n = 122) and urological cancer (n = 22). TRAF3IP3, the constituent gene of IHSI, was implicated in ICIs-relevant "immune-hypoxic" crosstalk by stimulating MAVS/IFN-I pathway under normoxic condition. MRI-radiomics models assessing TRAF3IP3 with HIF1A expression (AUC > 0.80) screened ICIs-Responders in HCC cohort (n = 75). CONCLUSION: The hypoxic-immune stratification redefined ICIs-responsive TME and provided MRI-Radiomics models for initial ICIs-responders screening, with IHSI facilitating further identification.