Search

BACKGROUND: The incidence of early-onset colorectal cancer (EOCRC) has been rising at an alarming rate in the USA, and EOCRC disproportionately affects racial/ethnic minorities. Here, we construct comprehensive profiles of EOCRC DNA methylomes at base-pair resolution for a cohort of Hispanic and African American patients. RESULTS: We show the epigenetic landscape of these EOCRC patients differs from that of late-onset colorectal cancer patients, and methylation canyons in EOCRC tumor tissue preferentially overlapped genes in cancer-related pathways. Furthermore, we identify epigenetic alterations in metabolic genes that are specific to our racial/ethnic minority EOCRC cohort but not Caucasian patients from TCGA. Top genes differentially methylated between these cohorts included the obesity-protective MFAP2 gene as well as cancer risk susceptibility genes APOL3 and RNASEL. CONCLUSIONS: In this study, we provide to the scientific community high-resolution DNA methylomes for a cohort of EOCRC patients from underrepresented populations. Our exploratory findings in this cohort highlight epigenetic mechanisms underlying the pathogenesis of EOCRC and nominate novel biomarkers for EOCRC in underrepresented populations.

Although rare non-coding variants (RVs) play crucial roles in complex traits and diseases, understanding their mechanisms and identifying disease-associated RVs continue to be major challenges. Here we constructed a comprehensive atlas of alternative polyadenylation (APA) outliers (aOutliers), including 1334 3 UTR and 200 intronic aOutliers, from 15,201 samples across 49 human tissues. These aOutliers exhibit unique characteristics from transcription or splicing outliers, with a pronounced RV enrichment. Mechanistically, aOutlier-RVs alter poly(A) signals and splicing sites, and perturbation indeed triggers APA events. Furthermore, we developed a Bayesian-based APA RV prediction model, which successfully pinpointed a specific set of 1799 RVs impacting 278 genes with significantly large disease effect sizes. Notably, we observed a convergence effect between rare and common cancer variants, exemplified by regulation in the DDX18 gene. Together, this study introduced an APA-enhanced framework for genome annotation, underscoring APAs role in uncovering functional RVs linked to complex traits and diseases.

BACKGROUND: Plasma cell-free DNA (cfDNA) is derived from cellular death in various tissues. Investigating the tissue origin of cfDNA through cell type deconvolution, we can detect changes in tissue homeostasis that occur during disease progression or in response to treatment. Consequently, cfDNA has emerged as a valuable noninvasive biomarker for disease detection and treatment monitoring. Although there are many methylation-based methods for cfDNA cell type deconvolution, a comprehensive and systematic evaluation of these methods has yet to be conducted. RESULTS: In this study, we benchmark five methods: MethAtlas, cfNOMe toolkit, CelFiE, CelFEER, and UXM. Utilizing deep whole-genome bisulfite sequencing data from 35 human cell types, we generate in silico cfDNA samples with ground truth cell type proportions to assess the deconvolution performance of the five methods under multiple scenarios. Our findings indicate that multiple factors, including reference marker selection, sequencing depth, and reference atlas completeness, jointly influence the deconvolution performance. Notably, an incomplete reference with missing markers or cell types leads to suboptimal results. We observe performance differences among methods under varying conditions, underscoring the importance of tailoring cfDNA deconvolution analyses. To increase the clinical relevance of our findings, we further evaluate each methods performance in potential clinical applications using real-world datasets. CONCLUSIONS: Based on the benchmark results, we propose general guidelines to choose the suitable methods based on sequencing depth of the cfDNA data and completeness of the reference atlas to maximize the performance of methylation-based cfDNA cell type deconvolution.

Two correspondences raised concerns or comments about our analyses regarding exaggerated false positives found by differential expression (DE) methods. Here, we discuss the points they raise and explain why we agree or disagree with these points. We add new analysis to confirm that the Wilcoxon rank-sum test remains the most robust method compared to the other five DE methods (DESeq2, edgeR, limma-voom, dearseq, and NOISeq) in two-condition DE analyses after considering normalization and winsorization, the data preprocessing steps discussed in the two correspondences.

Transfer-RNA-derived fragments (tRFs) are a novel class of small non-coding RNAs that have been implicated in oncogenesis. tRFs may act as post-transcriptional regulators by recruiting AGO proteins and binding to highly complementary regions of mRNA at seed regions, resulting in the knockdown of the transcript. Therefore, tRFs may be critical to tumorigenesis and warrant investigation as potential biomarkers. Meanwhile, the incidence of papillary thyroid carcinoma (PTC) has increased in recent decades and current diagnostic technology stands to benefit from new detection methods. Although small non-coding RNAs have been studied for their role in oncogenesis, there is currently no standard for their use as PTC biomarkers, and tRFs are especially underexplored. Accordingly, we aim to identify dysregulated tRFs in PTC that may serve as biomarker candidates. We identified dysregulated tRFs and driver genes between PTC primary tumor samples (n = 511) and adjacent normal tissue samples (n = 59). Expression data were obtained from MINTbase v2.0 and The Cancer Genome Atlas. Dysregulated tRFs and genes were analyzed in tandem to find pairs with anticorrelated expression. Significantly anticorrelated tRF-gene pairs were then tested for potential binding affinity using RNA22-if a heteroduplex can form via complementary binding, this would support the hypothesized RNA silencing mechanism. Four tRFs were significantly dysregulated in PTC tissue (p < 0.05), with only AsnGTT 3-tRF being upregulated. Binding affinity analysis revealed that tRF-30-RY73W0K5KKOV (AsnGTT 3-tRF) exhibits sufficient complementarity to potentially bind to and regulate transcripts of SLC26A4, SLC5A8, DIO2, and TPO, which were all found to be downregulated in PTC tissue. In the present study, we identified dysregulated tRFs in PTC and found that AsnGTT 3-tRF is a potential post-transcriptional regulator and biomarker.

We present a TALEN-based workflow to generate and maintain dual-edited (IL-15+/+/TGFβR2-/-) iPSCs that produce enhanced iPSC-derived natural killer (iNK) cells for cancer immunotherapy. It involves using a cell lineage promoter for knocking in (KI) gene(s) to minimize the potential effects of expression of any exogenous genes on iPSCs. As a proof-of-principle, we KI IL-15 under the endogenous B2M promoter and show that it results in high expression of the sIL-15 in iNK cells but minimal expression in iPSCs. Furthermore, given that it is known that knockout (KO) of TGFβR2 in immune cells can enhance resistance to the suppressive TGF-β signaling in the tumor microenvironment, we develop a customized medium containing Nodal that can maintain the pluripotency of iPSCs with TGFβR2 KO, enabling banking of these iPSC clones. Ultimately, we show that the dual-edited IL-15+/+/TGFβR2-/- iPSCs can be efficiently differentiated into NK cells that show enhanced autonomous growth and are resistant to the suppressive TGF-β signaling.

BACKGROUND: RNA sequencing (RNA-Seq) offers profound insights into the complex transcriptomes of diverse biological systems. However, standard differential expression analysis pipelines based on DESeq2 and edgeR encounter challenges when applied to the immediate early transcriptomes of Chlamydia spp., obligate intracellular bacteria. These challenges arise from their reliance on assumptions that do not hold in scenarios characterized by extensive transcriptomic activation and limited repression. RESULTS: Standard analyses using unique chlamydial RNA-Seq reads alone identify nearly 300 upregulated and about 300 downregulated genes, significantly deviating from actual RNA-Seq read trends. By incorporating both chlamydial and host reads or adjusting for total sequencing depth, the revised normalization methods each detected over 700 upregulated genes and 30 or fewer downregulated genes, closely aligned with observed RNA-Seq data. Further validation through qRT-PCR analysis confirmed the effectiveness of these adjusted approaches in capturing the true extent of transcriptomic activation during the immediate early phase of chlamydial infection. CONCLUSIONS: This study highlights the limitations of standard RNA-Seq analysis tools in scenarios with extensive transcriptomic activation, such as in Chlamydia spp. during early infection. Our revised normalization methods, incorporating host reads or total sequencing depth, provide a more accurate representation of gene expression dynamics. These approaches may inform similar adjustments in other systems with unbalanced gene expression dynamics, enhancing the accuracy of transcriptomic analysis.

We identify a population of Protogenin-positive (PRTG+ve) MYChigh NESTINlow stem cells in the four-week-old human embryonic hindbrain that subsequently localizes to the ventricular zone of the rhombic lip (RLVZ). Oncogenic transformation of early Prtg+ve rhombic lip stem cells initiates group 3 medulloblastoma (Gr3-MB)-like tumors. PRTG+ve stem cells grow adjacent to a human-specific interposed vascular plexus in the RLVZ, a phenotype that is recapitulated in Gr3-MB but not in other types of medulloblastoma. Co-culture of Gr3-MB with endothelial cells promotes tumor stem cell growth, with the endothelial cells adopting an immature phenotype. Targeting the PRTGhigh compartment of Gr3-MB in vivo using either the diphtheria toxin system or chimeric antigen receptor T cells constitutes effective therapy. Human Gr3-MBs likely arise from early embryonic RLVZ PRTG+ve stem cells inhabiting a specific perivascular niche. Targeting the PRTGhigh compartment and/or the perivascular niche represents an approach to treat children with Gr3-MB.

The Genome Aggregation Database (gnomAD), widely recognized as the gold-standard reference map of human genetic variation, has largely overlooked tandem repeat (TR) expansions, despite the fact that TRs constitute ∼6% of our genome and are linked to over 50 human diseases. Here, we introduce the TR-gnomAD (https://wlcb.oit.uci.edu/TRgnomAD), a biobank-scale reference of 0.86 million TRs derived from 338,963 whole-genome sequencing (WGS) samples of diverse ancestries (39.5% non-European samples). TR-gnomAD offers critical insights into ancestry-specific disease prevalence using disparities in TR unit number frequencies among ancestries. Moreover, TR-gnomAD is able to differentiate between common, presumably benign TR expansions, which are prevalent in TR-gnomAD, from those potentially pathogenic TR expansions, which are found more frequently in disease groups than within TR-gnomAD. Together, TR-gnomAD is an invaluable resource for researchers and physicians to interpret TR expansions in individuals with genetic diseases.

BACKGROUND: Although many molecules have been investigated as biomarkers for spinal cord injury (SCI) or ischemic stroke, none of them are specifically induced in central nervous system (CNS) neurons following injuries with low baseline expression. However, neuronal injury constitutes a major pathology associated with SCI or stroke and strongly correlates with neurological outcomes. Biomarkers characterized by low baseline expression and specific induction in neurons post-injury are likely to better correlate with injury severity and recovery, demonstrating higher sensitivity and specificity for CNS injuries compared to non-neuronal markers or pan-neuronal markers with constitutive expressions. METHODS: In animal studies, young adult wildtype and global Atf3 knockout mice underwent unilateral cervical 5 (C5) SCI or permanent distal middle cerebral artery occlusion (pMCAO). Gene expression was assessed using RNA-sequencing and qRT-PCR, while protein expression was detected through immunostaining. Serum ATF3 levels in animal models and clinical human samples were measured using commercially available enzyme-linked immune-sorbent assay (ELISA) kits. RESULTS: Activating transcription factor 3 (ATF3), a molecular marker for injured dorsal root ganglion sensory neurons in the peripheral nervous system, was not expressed in spinal cord or cortex of naïve mice but was induced specifically in neurons of the spinal cord or cortex within 1 day after SCI or ischemic stroke, respectively. Additionally, ATF3 protein levels in mouse blood significantly increased 1 day after SCI or ischemic stroke. Importantly, ATF3 protein levels in human serum were elevated in clinical patients within 24 hours after SCI or ischemic stroke. Moreover, Atf3 knockout mice, compared to the wildtype mice, exhibited worse neurological outcomes and larger damage regions after SCI or ischemic stroke, indicating that ATF3 has a neuroprotective function. CONCLUSIONS: ATF3 is an easily measurable, neuron-specific biomarker for clinical SCI and ischemic stroke, with neuroprotective properties. HIGHLIGHTS: ATF3 was induced specifically in neurons of the spinal cord or cortex within 1 day after SCI or ischemic stroke, respectively. Serum ATF3 protein levels are elevated in clinical patients within 24 hours after SCI or ischemic stroke. ATF3 exhibits neuroprotective properties, as evidenced by the worse neurological outcomes and larger damage regions observed in Atf3 knockout mice compared to wildtype mice following SCI or ischemic stroke.

The obligate intracellular bacterium Chlamydia has a unique developmental cycle that alternates between two contrasting cell types. With a hardy envelope and highly condensed genome, the small elementary body (EB) maintains limited metabolic activities yet survives in extracellular environments and is infectious. After entering host cells, EBs differentiate into larger and proliferating reticulate bodies (RBs). Progeny EBs are derived from RBs in late developmental stages and eventually exit host cells. How expression of the chlamydial genome consisting of nearly 1,000 genes governs the chlamydial developmental cycle is unclear. A previous microarray study identified only 29 Chlamydia trachomatis immediate early genes, defined as genes with increased expression during the first hour postinoculation in cultured cells. In this study, we performed more sensitive RNA sequencing (RNA-Seq) analysis for C. trachomatis cultures with high multiplicities of infection. Remarkably, we observed well over 700 C. trachomatis genes that underwent 2- to 900-fold activation within 1 hour postinoculation. Quantitative reverse transcription real-time PCR analysis was further used to validate the activated expression of a large subset of the genes identified by RNA-Seq. Importantly, our results demonstrate that the immediate early transcriptome is over 20 times more extensive than previously realized. Gene ontology analysis indicates that the activated expression spans all functional categories. We conclude that over 70% of C. trachomatis genes are activated in EBs almost immediately upon entry into host cells, thus implicating their importance in initiating rapid differentiation into RBs and establishing an intracellular niche conducive with chlamydial development and growth.

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3 untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3 untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

The MYC oncogene is often dysregulated in human cancer, including hepatocellular carcinoma (HCC). MYC is considered undruggable to date. Here, we comprehensively identify genes essential for survival of MYChigh but not MYClow cells by a CRISPR/Cas9 genome-wide screen in a MYC-conditional HCC model. Our screen uncovers novel MYC synthetic lethal (MYC-SL) interactions and identifies most MYC-SL genes described previously. In particular, the screen reveals nucleocytoplasmic transport to be a MYC-SL interaction. We show that the majority of MYC-SL nucleocytoplasmic transport genes are upregulated in MYChigh murine HCC and are associated with poor survival in HCC patients. Inhibiting Exportin-1 (XPO1) in vivo induces marked tumor regression in an autochthonous MYC-transgenic HCC model and inhibits tumor growth in HCC patient-derived xenografts. XPO1 expression is associated with poor prognosis only in HCC patients with high MYC activity. We infer that MYC may generally regulate and require altered expression of nucleocytoplasmic transport genes for tumorigenesis.

Hallmarks of the developmental cycle of the obligate intracellular pathogenic bacterium Chlamydia are the primary differentiation of the infectious elementary body (EB) into the proliferative reticulate body (RB) and the secondary differentiation of RBs back into EBs. The mechanisms regulating these transitions remain unclear. In this report, we developed an effective novel strategy termed dependence on plasmid-mediated expression (DOPE) that allows for the knockdown of essential genes in Chlamydia. We demonstrate that GrgA, a Chlamydia-specific transcription factor, is essential for the secondary differentiation and optimal growth of RBs. We also show that GrgA, a chromosome-encoded regulatory protein, controls the maintenance of the chlamydial virulence plasmid. Transcriptomic analysis further indicates that GrgA functions as a critical regulator of all three sigma factors that recognize different promoter sets at developmental stages. The DOPE strategy outlined here should provide a valuable tool for future studies examining chlamydial growth, development, and pathogenicity.

BACKGROUND & AIMS: Lacticaseibacillus rhamnosus GG (LGG) is the worlds most consumed probiotic but its mechanism of action on intestinal permeability and differentiation along with its interactions with an essential source of signaling metabolites, dietary tryptophan (trp), are unclear. METHODS: Untargeted metabolomic and transcriptomic analyses were performed in LGG monocolonized germ-free mice fed trp-free or -sufficient diets. LGG-derived metabolites were profiled in vitro under anaerobic and aerobic conditions. Multiomic correlations using a newly developed algorithm discovered novel metabolites tightly linked to tight junction and cell differentiation genes whose abundances were regulated by LGG and dietary trp. Barrier-modulation by these metabolites were functionally tested in Caco2 cells, mouse enteroids, and dextran sulfate sodium experimental colitis. The contribution of these metabolites to barrier protection is delineated at specific tight junction proteins and enterocyte-promoting factors with gain and loss of function approaches. RESULTS: LGG, strictly with dietary trp, promotes the enterocyte program and expression of tight junction genes, particularly Ocln. Functional evaluations of fecal and serum metabolites synergistically stimulated by LGG and trp revealed a novel vitamin B3 metabolism pathway, with methylnicotinamide (MNA) unexpectedly being the most robust barrier-protective metabolite in vitro and in vivo. Reduced serum MNA is significantly associated with increased disease activity in patients with inflammatory bowel disease. Exogenous MNA enhances gut barrier in homeostasis and robustly promotes colonic healing in dextran sulfate sodium colitis. MNA is sufficient to promote intestinal epithelial Ocln and RNF43, a master inhibitor of Wnt. Blocking trp or vitamin B3 absorption abolishes barrier recovery in vivo. CONCLUSIONS: Our study uncovers a novel LGG-regulated dietary trp-dependent production of MNA that protects the gut barrier against colitis.

Huntingtons disease (HD) is a progressive neurodegenerative brain disorder associated with uncontrolled body movements, cognitive decline, and reduced circulating melatonin levels. Melatonin is a potent antioxidant and exogenous melatonin treatment is neuroprotective in experimental HD models. In neurons, melatonin is exclusively synthesized in the mitochondrial matrix. Thus, we investigated the integrity of melatonin biosynthesis pathways in pineal and extrapineal brain areas in human HD brain samples, in the R6/2 mouse model of HD and in full-length mutant huntingtin knock-in cells. Aralkylamine N-acetyltransferase (AANAT) is the rate-limiting step enzyme in the melatonin biosynthetic pathway. We found that AANAT expression is significantly decreased in the pineal gland and the striatum of HD patients compared to normal controls. In the R6/2 mouse forebrain, AANAT protein expression was decreased in synaptosomal, but not nonsynaptosomal, mitochondria and was associated with decreased synaptosomal melatonin levels compared to wild type mice. We also demonstrate sequestration of AANAT in mutant-huntingtin protein aggregates likely resulting in decreased AANAT bioavailability. Paradoxically, AANAT mRNA expression is increased in tissues where AANAT protein expression is decreased, suggesting a potential feedback loop that is, ultimately unsuccessful. In conclusion, we demonstrate that pineal, extrapineal, and synaptosomal melatonin levels are compromised in the brains of HD patients and R6/2 mice due, at least in part, to protein aggregation.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

Paneth cells (PCs), a specialized secretory cell type in the small intestine, are increasingly recognized as having an essential role in host responses to microbiome and environmental stresses. Whether and how commensal and pathogenic microbes modify PC composition to modulate inflammation remain unclear. Using newly developed PC-reporter mice under conventional and gnotobiotic conditions, we determined PC transcriptomic heterogeneity in response to commensal and invasive microbes at single cell level. Infection expands the pool of CD74+ PCs, whose number correlates with auto or allogeneic inflammatory disease progressions in mice. Similar correlation was found in human inflammatory disease tissues. Infection-stimulated cytokines increase production of reactive oxygen species (ROS) and expression of a PC-specific mucosal pentraxin (Mptx2) in activated PCs. A PC-specific ablation of MyD88 reduced CD74+ PC population, thus ameliorating pathogen-induced systemic disease. A similar phenotype was also observed in mice lacking Mptx2. Thus, infection stimulates expansion of a PC subset that influences disease progression.

The majority of oncogenic drivers are intracellular proteins, constraining their immunotherapeutic targeting to mutated peptides (neoantigens) presented by individual human leukocyte antigen (HLA) allotypes1. However, most cancers have a modest mutational burden that is insufficient for generating responses using neoantigen-based therapies2,3. Neuroblastoma is a paediatric cancer that harbours few mutations and is instead driven by epigenetically deregulated transcriptional networks4. Here we show that the neuroblastoma immunopeptidome is enriched with peptides derived from proteins essential for tumorigenesis. We focused on targeting the unmutated peptide QYNPIRTTF discovered on HLA-A*24:02, which is derived from the neuroblastoma-dependency gene and master transcriptional regulator PHOX2B. To target QYNPIRTTF, we developed peptide-centric chimeric antigen receptors (PC-CARs) through a counter panning strategy using predicted potentially cross-reactive peptides. We further proposed that PC-CARs can recognize peptides on additional HLA allotypes when presenting a similar overall molecular surface. Informed by our computational modelling results, we show that PHOX2B PC-CARs also recognize QYNPIRTTF presented by HLA-A*23:01, the most common non-A2 allele in people with African ancestry. Finally, we demonstrate potent and specific killing of neuroblastoma cells expressing these HLAs in vitro and complete tumour regression in mice. These data suggest that PC-CARs have the potential to expand the pool of immunotherapeutic targets to include non-immunogenic intracellular oncoproteins and allow targeting through additional HLA allotypes in a clinical setting.

Chronic infection of hepatitis B virus (HBV) is the major cause of hepatocellular carcinoma (HCC). Notably, 90% of HBV-positive HCC cases exhibit detectable HBV integrations, hinting at the potential early entanglement of these viral integrations in tumorigenesis and their subsequent oncogenic implications. Nevertheless, the precise chronology of integration events during HCC tumorigenesis, alongside their sequential structural patterns, has remained elusive thus far. In this study, we applied whole-genome sequencing to multiple biopsies extracted from six HBV-positive HCC cases. Through this approach, we identified point mutations and viral integrations, offering a blueprint for the intricate tumor phylogeny of these samples. The emergent narrative paints a rich tapestry of diverse evolutionary trajectories characterizing the analyzed tumors. We uncovered oncogenic integration events in some samples that appear to happen before and during the initiation stage of tumor development based on their locations in reconstituted trajectories. Furthermore, we conducted additional long-read sequencing of selected samples and unveiled integration-bridged chromosome rearrangements and tandem repeats of the HBV sequence within integrations. In summary, this study revealed premalignant oncogenic and sequential complex integrations and highlighted the contributions of HBV integrations to HCC development and genome instability.

The advent of single-cell RNA sequencing (scRNA-seq) technologies has enabled gene expression profiling at the single-cell resolution, thereby enabling the quantification and comparison of transcriptional variability among individual cells. Although alterations in transcriptional variability have been observed in various biological states, statistical methods for quantifying and testing differential variability between groups of cells are still lacking. To identify the best practices in differential variability analysis of single-cell gene expression data, we propose and compare 12 statistical pipelines using different combinations of methods for normalization, feature selection, dimensionality reduction and variability calculation. Using high-quality synthetic scRNA-seq datasets, we benchmarked the proposed pipelines and found that the most powerful and accurate pipeline performs simple library size normalization, retains all genes in analysis and uses denSNE-based distances to cluster medoids as the variability measure. By applying this pipeline to scRNA-seq datasets of COVID-19 and autism patients, we have identified cellular variability changes between patients with different severity status or between patients and healthy controls.

IntroductionAdverse drug events (ADEs) are associated with poor outcomes and increased costs but may be prevented with prediction tools. With the National Institute of Health All of Us (AoU) database, we employed machine learning (ML) to predict selective serotonin reuptake inhibitor (SSRI)-associated bleeding.MethodsThe AoU program, beginning in 05/2018, continues to recruit ≥ 18 years old individuals across the United States. Participants completed surveys and consented to contribute electronic health record (EHR) for research. Using the EHR, we determined participants who were exposed to SSRIs (citalopram, escitalopram, fluoxetine, fluvoxamine, paroxetine, sertraline, vortioxetine). Features (n = 88) were selected with clinicians' input and comprised sociodemographic, lifestyle, comorbidities, and medication use information. We identified bleeding events with validated EHR algorithms and applied logistic regression, decision tree, random forest, and extreme gradient boost to predict bleeding during SSRI exposure. We assessed model performance with area under the receiver operating characteristic curve statistic (AUC) and defined clinically significant features as resulting in > 0.01 decline in AUC after removal from the model, in three of four ML models.ResultsThere were 10,362 participants exposed to SSRIs, with 9.6% experiencing a bleeding event during SSRI exposure. For each SSRI, performance across all four ML models was relatively consistent. AUCs from the best models ranged 0.632-0.698. Clinically significant features included health literacy for escitalopram, and bleeding history and socioeconomic status for all SSRIs.ConclusionsWe demonstrated feasibility of predicting ADEs using ML. Incorporating genomic features and drug interactions with deep learning models may improve ADE prediction.

BACKGROUND: Methylation of the p16 promoter resulting in epigenetic gene silencing-known as p16 epimutation-is frequently found in human colorectal cancer and is also common in normal-appearing colonic mucosa of aging individuals. Thus, to improve clinical care of colorectal cancer (CRC) patients, we explored the role of age-related p16 epimutation in intestinal tumorigenesis. METHODS: We established a mouse model that replicates two common genetic and epigenetic events observed in human CRCs: Apc mutation and p16 epimutation. We conducted long-term survival and histological analysis of tumor development and progression. Colonic epithelial cells and tumors were collected from mice and analyzed by RNA sequencing (RNA-seq), quantitative PCR, and flow cytometry. We performed single-cell RNA sequencing (scRNA-seq) to characterize tumor-infiltrating immune cells throughout tumor progression. We tested whether anti-PD-L1 immunotherapy affects overall survival of tumor-bearing mice and whether inhibition of both epigenetic regulation and immune checkpoint is more efficacious. RESULTS: Mice carrying combined Apc mutation and p16 epimutation had significantly shortened survival and increased tumor growth compared to those with Apc mutation only. Intriguingly, colon tumors with p16 epimutation exhibited an activated interferon pathway, increased expression of programmed death-ligand 1 (Pdl1), and enhanced infiltration of immune cells. scRNA-seq further revealed the presence of Foxp3+ Tregs and γδT17 cells, which contribute to an immunosuppressive tumor microenvironment (TME). Furthermore, we showed that a combined therapy using an inhibitor of DNA methylation and a PD-L1 immune checkpoint inhibitor is more effective for improving survival in tumor-bearing mice than blockade of either pathway alone. CONCLUSIONS: Our study demonstrated that age-dependent p16 epimutation creates a permissive microenvironment for malignant transformation of polyps to colon cancer. Our findings provide a mechanistic rationale for future targeted therapy in patients with p16 epimutation.

Multiple randomized, controlled clinical trials have yielded discordant results regarding the efficacy of convalescent plasma in outpatients, with some showing an approximately 2-fold reduction in risk and others showing no effect. We quantified binding and neutralizing antibody levels in 492 of the 511 participants from the Clinical Trial of COVID-19 Convalescent Plasma in Outpatients (C3PO) of a single unit of COVID-19 convalescent plasma (CCP) versus saline infusion. In a subset of 70 participants, peripheral blood mononuclear cells were obtained to define the evolution of B and T cell responses through day 30. Binding and neutralizing antibody responses were approximately 2-fold higher 1 hour after infusion in recipients of CCP compared with saline plus multivitamin, but levels achieved by the native immune system by day 15 were almost 10-fold higher than those seen immediately after CCP administration. Infusion of CCP did not block generation of the host antibody response or skew B or T cell phenotype or maturation. Activated CD4+ and CD8+ T cells were associated with more severe disease outcome. These data show that CCP leads to a measurable boost in anti-SARS-CoV-2 antibodies but that the boost is modest and may not be sufficient to alter disease course.

Sex or gender differences in the risk of Alzheimer's disease and related dementias (ADRD) differ by world region, suggesting that there are potentially modifiable risk factors for intervention. However, few epidemiological or clinical ADRD studies examine sex differences; even fewer evaluate gender in the context of ADRD risk. The goals of this perspective are to: (1) provide definitions of gender, biologic sex, and sexual orientation. and the limitations of examining these as binary variables; (2) provide an overview of what is known with regard to sex and gender differences in the risk, prevention, and diagnosis of ADRD; and (3) discuss these sex and gender differences from a global, worldwide perspective. Identifying drivers of sex and gender differences in ADRD throughout the world is a first step in developing interventions unique to each geographical and sociocultural area to reduce these inequities and to ultimately reduce global ADRD risk. HIGHLIGHTS: The burden of dementia is unevenly distributed geographically and by sex and gender. Scientific advances in genetics and biomarkers challenge beliefs that sex is binary. Discrimination against women and sex and gender minority (SGM) populations contributes to cognitive decline. Sociocultural factors lead to gender inequities in Alzheimer's disease and related dementias (ADRD) worldwide.

When identifying differentially expressed genes between two conditions using human population RNA-seq samples, we found a phenomenon by permutation analysis: two popular bioinformatics methods, DESeq2 and edgeR, have unexpectedly high false discovery rates. Expanding the analysis to limma-voom, NOISeq, dearseq, and Wilcoxon rank-sum test, we found that FDR control is often failed except for the Wilcoxon rank-sum test. Particularly, the actual FDRs of DESeq2 and edgeR sometimes exceed 20% when the target FDR is 5%. Based on these results, for population-level RNA-seq studies with large sample sizes, we recommend the Wilcoxon rank-sum test.

Head and neck squamous cell carcinoma (HNSCC) tumor phenotypes and clinical outcomes are significantly influenced by etiological agents, such as HPV infection, smoking, and alcohol consumption. Accordingly, the intratumor microbiome has been increasingly implicated in cancer progression and metastasis. However, few studies characterize the intratumor microbial landscape of HNSCC with respect to these etiological agents. In this study, we aimed to investigate the bacterial and fungal landscape of HNSCC in association with HPV infection, smoking, and alcohol consumption. RNA-sequencing data were extracted from The Cancer Genome Atlas (TCGA) regarding 449 tissue samples and 44 normal samples. Pathoscope 2.0 was used to extract the microbial reads. Microbe abundance was compared to clinical variables, oncogenic signatures, and immune-associated pathways. Our results demonstrated that a similar number of dysregulated microbes was overabundant in smokers and nonsmokers, while heavy drinkers were characterized by an underabundance of dysregulated microbes. Conversely, the majority of dysregulated microbes were overabundant in HPV+ tumor samples when compared to HPV- tumor samples. Moreover, we observed that many dysregulated microbes were associated with oncogenic and metastatic pathways, suggesting their roles in influencing carcinogenesis. These microbes provide insights regarding potential mechanisms for tumor pathogenesis and progression with respect to the three etiological agents.

The limited performance of guideline-recommended abdominal ultrasound and serum alpha-fetoprotein (AFP) highlights the urgent, unmet need for new biomarkers for more accurate detection of early hepatocellular carcinoma (HCC). To this end, we have conducted a prospective clinical validation study to evaluate the performance of the HelioLiver Test, a multi-analyte blood test combining cell-free DNA methylation patterns, clinical variables, and protein tumor markers. A blinded, multicenter validation study was performed with 247 subjects, including 122 subjects with HCC and 125 control subjects with chronic liver disease. The performance of the HelioLiver Test was compared with AFP and the GALAD score as established HCC surveillance blood tests. The performance of the HelioLiver Test (area under the receiver operating characteristic curve [AUROC] = 0.944) was superior to both AFP (AUROC = 0.851; p < 0.0001) and GALAD (AUROC = 0.899; p < 0.0001). Using a prespecified diagnostic algorithm, the HelioLiver Test showed sensitivities of 85% (95% confidence interval [CI], 78%-90%) for HCC of any stage and 76% (95% CI, 60%-87%) for early stage (American Joint Committee on Cancer [AJCC] I and II) HCC. In contrast, AFP (≥20 ng/mL) alone and the GALAD score (≥-0.63) showed lower sensitivities of 62% (95% CI, 54%-70%) and 75% (95% CI, 67%-82%) for HCC overall, and 57% (95% CI, 40%-71%) and 65% (95% CI, 49%-79%) for early stage (AJCC I and II) HCC, respectively. The specificities of the HelioLiver Test (91%; 95% CI, 85%-95%), AFP (97%; 95% CI, 92%-99%), and the GALAD score (94%; 95% CI, 88%-97%) were similar for control subjects. The HelioLiver Test showed superior performance for HCC detection compared to with both AFP and the GALAD score and warrants further evaluation in HCC surveillance settings.

Pontine autosomal dominant microangiopathy and leukoencephalopathy (PADMAL) is a rare hereditary cerebral small vessel disease. We report a novel collagen type IV alpha 1 (COL4A1) gene mutation in a Chinese family with PADMAL. The index case was followed up for 6 years. Neuroimaging, whole-exome sequencing, skin biopsy, and pedigree analysis were performed. She initially presented with minor head injury at age 38. MRI brain showed chronic lacunar infarcts in the pons, left thalamus, and right centrum semiovale. Extensive workup was unremarkable except for a patent foramen ovale (PFO). Despite anticoagulation, PFO closure, and antiplatelet therapy, the patient had recurrent lacunar infarcts in the pons and deep white matter, as well as subcortical microhemorrhages. Whole-exome sequencing demonstrated a novel c.*34G > T mutation in the 3' untranslated region of COL4A1 gene. Skin biopsy subsequently demonstrated thickening of vascular basement membrane, proliferation of endothelial cells, and stenosis of vascular lumen. Three additional family members had gene testing and 2 of them were found to have the same heterozygous mutation. Of the 18 individuals in the pedigree of 3 generations, 12 had clinical and MRI evidence of PADMAL. The mechanisms of both ischemic and hemorrhagic stroke are likely the overexpression of COLT4A1 in the basement membrane and frugality of the vessel walls. Our findings suggest that the novel c.*34G > T mutation appears to have the same functional consequences as the previously reported COL4A1 gene mutations in patients with PADMAL and multi-infarct dementia of Swedish type.

A comprehensive evaluation of the clear cell renal cell carcinoma (ccRCC) immune landscape was found using 584 RNA-sequencing datasets from The Cancer Genome Atlas (TCGA), we identified 17 key dysregulated immune-associated genes in ccRCC based on association with clinical variables and important immune pathways. Of the numerous findings from our analyses, we found that several of the 17 key dysregulated genes are heavily involved in interleukin and NF-kB signaling and that somatic copy number alteration (SCNA) hotspots may be causally associated with gene dysregulation. More importantly, we also found that key immune-associated genes and pathways are strongly upregulated in ccRCC. Our study may lend novel insights into the clinical implications of immune dysregulation in ccRCC and suggests potential immunotherapeutic targets for further evaluation.

Immune checkpoints (ICPs) consist of paired receptor-ligand molecules that exert inhibitory or stimulatory effects on immune defense, surveillance, regulation, and self-tolerance. ICPs exist in both membrane and soluble forms in vivo and in vitro. Imbalances between inhibitory and stimulatory membrane-bound ICPs (mICPs) in malignant cells and immune cells in the tumor immune microenvironment (TIME) have been well documented. Blockades of inhibitory mICPs have emerged as an immense breakthrough in cancer therapeutics. However, the origin, structure, production regulation, and biological significance of soluble ICPs (sICPs) in health and disease largely remains elusive. Soluble ICPs can be generated through either alternative mRNA splicing and secretion or protease-mediated shedding from mICPs. Since sICPs are found in the bloodstream, they likely form a circulating immune regulatory system. In fact, there is increasing evidence that sICPs exhibit biological functions including (1) regulation of antibacterial immunity, (2) interaction with their mICP compartments to positively or negatively regulate immune responses, and (3) competition with their mICP compartments for binding to the ICP blocking antibodies, thereby reducing the efficacy of ICP blockade therapies. Here, we summarize current data of sICPs in cancer and infectious diseases. We particularly focus on sICPs in COVID-19 and HIV infection as they are the two ongoing global pandemics and have created the world's most serious public health challenges. A "storm" of sICPs occurs in the peripheral circulation of COVID-19 patients and is associated with the severity of COVID-19. Similarly, sICPs are highly dysregulated in people living with HIV (PLHIV) and some sICPs remain dysregulated in PLHIV on antiretroviral therapy (ART), indicating these sICPs may serve as biomarkers of incomplete immune reconstitution in PLHIV on ART. We reveal that HIV infection in the setting of alcohol misuse exacerbates sICP dysregulation as PLHIV with heavy alcohol consumption have significantly elevated plasma levels of many sICPs. Thus, both stimulatory and inhibitory sICPs are present in the bloodstream of healthy people and their balance can be disrupted under pathophysiological conditions such as cancer, COVID-19, HIV infection, and alcohol misuse. There is an urgent need to study the role of sICPs in immune regulation in health and disease.

Phagocytic clearance of dying cells, termed efferocytosis, is essential for maintaining tissue homeostasis, yet our understanding of efferocytosis regulation remains incomplete. Here we perform a FACS-based, genome-wide CRISPR knockout screen in primary mouse macrophages to search for novel regulators of efferocytosis. The results show that Wdfy3 knockout in macrophages specifically impairs uptake, but not binding, of apoptotic cells due to defective actin disassembly. Additionally, WDFY3 interacts with GABARAP, thus facilitating LC3 lipidation and subsequent lysosomal acidification to permit the degradation of apoptotic cell components. Mechanistically, while the C-terminus of WDFY3 is sufficient to rescue the impaired degradation induced by Wdfy3 knockout, full-length WDFY3 is required to reconstitute the uptake of apoptotic cells. Finally, WDFY3 is also required for efficient efferocytosis in vivo in mice and in vitro in primary human macrophages. This work thus expands our knowledge of the mechanisms of macrophage efferocytosis, as well as supports genome-wide CRISPR screen as a platform for interrogating complex functional phenotypes in primary macrophages.

Mood disorders are an important public health issue and recent advances in genomic studies have indicated that molecules involved in neurodevelopment are causally related to mood disorders. BLM-s (BCL-2-like molecule, small transcript isoform), a BH3-only proapoptotic BCL-2 family member, mediates apoptosis of postmitotic immature neurons during embryonic cortical development, but its role in the adult brain is unknown. To better understand the physiological role of Blm-s gene in vivo, we generated a Blm-s-knockout (Blm-s-/-) mouse. The Blm-s-/- mice breed normally and exhibit grossly normal development. However, global depletion of Blm-s is highly associated with depression- and anxiety-related behaviors in adult mutant mice with intact learning and memory capacity. Functional magnetic resonance imaging of adult Blm-s-/- mice reveals reduced connectivity mainly in the ventral dentate gyrus (vDG) of the hippocampus with no alteration in the dorsal DG connectivity and in total hippocampal volume. At the cellular level, BLM-s is expressed in DG granule cells (GCs), and Blm-s-/- mice show reduced dendritic complexity and decreased spine density in mature GCs. Electrophysiology study uncovers that mature vGCs in adult Blm-s-/- DG are intrinsically more excitable. Interestingly, certain genetic variants of the human Blm homologue gene (VPS50) are significantly associated with depression traits from publicly resourced UK Biobank data. Taken together, BLM-s is required for the hippocampal mood control function. Loss of BLM-s causes abnormality in the electrophysiology and morphology of GCs and a disrupted vDG neural network, which could underlie Blm-s-null-associated anxiety and depression.

Lung development is precisely controlled by underlying gene regulatory networks (GRN). Disruption of genes in the network can interrupt normal development and cause diseases such as bronchopulmonary dysplasia (BPD) - a chronic lung disease in preterm infants with morbid and sometimes lethal consequences characterized by lung immaturity and reduced alveolarization. Here, we generated a transgenic mouse exhibiting a moderate severity BPD phenotype by blocking IGF1 signaling in secondary crest myofibroblasts (SCMF) at the onset of alveologenesis. Using approaches mirroring the construction of the model GRN in sea urchin's development, we constructed the IGF1 signaling network underlying alveologenesis using this mouse model that phenocopies BPD. The constructed GRN, consisting of 43 genes, provides a bird's eye view of how the genes downstream of IGF1 are regulatorily connected. The GRN also reveals a mechanistic interpretation of how the effects of IGF1 signaling are transduced within SCMF from its specification genes to its effector genes and then from SCMF to its neighboring alveolar epithelial cells with WNT5A and FGF10 signaling as the bridge. Consistently, blocking WNT5A signaling in mice phenocopies BPD as inferred by the network. A comparative study on human samples suggests that a GRN of similar components and wiring underlies human BPD. Our network view of alveologenesis is transforming our perspective to understand and treat BPD. This new perspective calls for the construction of the full signaling GRN underlying alveologenesis, upon which targeted therapies for this neonatal chronic lung disease can be viably developed.

Asthma currently affects more than 339 million people worldwide. In the present preliminary study, we examined the efficacy of a new, inhalable soluble epoxide hydrolase inhibitor (sEHI), 1-trifluoromethoxyphenyl-3-(1-propionylpiperidin-4-yl) urea (TPPU), to attenuate airway inflammation, mucin secretion, and hyper-responsiveness (AHR) in an ovalbumin (OVA)-sensitized murine model. Male BALB/c mice were divided into phosphate-buffered saline (PBS), OVA, and OVA+TPPU (2- or 6-h) exposure groups. On days 0 and 14, the mice were administered PBS or sensitized to OVA in PBS. From days 26-38, seven challenge exposures were performed with 30 min inhalation of filtered air or OVA alone. In the OVA+TPPU groups, a 2- or 6-h TPPU inhalation preceded each 30-min OVA exposure. On day 39, pulmonary function tests (PFTs) were performed, and biological samples were collected. Lung tissues were used to semi-quantitatively evaluate the severity of inflammation and airway constriction and the volume of stored intracellular mucosubstances. Bronchoalveolar lavage (BAL) and blood samples were used to analyze regulatory lipid mediator profiles. Significantly (p < 0.05) attenuated alveolar, bronchiolar, and pleural inflammation; airway resistance and constriction; mucosubstance volume; and inflammatory lipid mediator levels were observed with OVA+TPPU relative to OVA alone. Cumulative findings indicated TPPU inhalation effectively inhibited inflammation, suppressed AHR, and prevented mucosubstance accumulation in the murine asthmatic model. Future studies should determine the pharmacokinetics (i.e., absorption, distribution, metabolism, and excretion) and pharmacodynamics (i.e., concentration/dose responses) of inhaled TPPU to explore its potential as an asthma-preventative or -rescue treatment.

Epidemiological studies show strong associations between fine particulate matter (PM2.5) air pollution and adverse pulmonary effects. In the present study, wintertime PM2.5 samples were collected from three geographically similar regions-Sacramento, California, USA; Jinan, Shandong, China; and Taiyuan, Shanxi, China-and extracted to form PMCA, PMSD, and PMSX, respectively, for comparison in a BALB/c mouse model. Each of four groups was oropharyngeally administered Milli-Q water vehicle control (50 μL) or one type of PM extract (20 μg/50 μL) five times over two weeks. Mice were necropsied on post-exposure days 1, 2, and 4 and examined using bronchoalveolar lavage (BAL), histopathology, and assessments of cytokine/chemokine mRNA and protein expression. Chemical analysis demonstrated all three extracts contained black carbon, but PMSX contained more sulfates and polycyclic aromatic hydrocarbons (PAHs) associated with significantly greater neutrophil numbers and greater alveolar/bronchiolar inflammation on post-exposure days 1 and 4. On day 4, PMSX-exposed mice also exhibited significant increases in interleukin-1 beta, tumor necrosis factor-alpha, and chemokine C-X-C motif ligands-3 and -5 mRNA, and monocyte chemoattractant protein-1 protein. These combined findings suggest greater sulfate and PAH content contributed to a more intense and progressive inflammatory response with repeated PMSX compared to PMCA or PMSD exposure.

Most eukaryotic genes express alternative polyadenylation (APA) isoforms. A growing number of RNA sequencing methods, especially those used for single-cell transcriptome analysis, generate reads close to the polyadenylation site (PAS), termed nearSite reads, hence inherently containing information about APA isoform abundance. Here, we present a probabilistic model-based method named MAAPER to utilize nearSite reads for APA analysis. MAAPER predicts PASs with high accuracy and sensitivity and examines different types of APA events with robust statistics. We show MAAPER's performance with both bulk and single-cell data and its applicability in unpaired or paired experimental designs.

Abstract Background Allopurinol is a xanthine oxidase inhibitor that lowers serum uric acid and is used to prevent acute gout flares in patients with gout. Observational and small interventional studies have suggested beneficial cardiovascular effects of allopurinol. Objective To determine whether allopurinol improves major cardiovascular outcomes in patients with ischaemic heart disease. Design Prospective, randomised, open-label, blinded endpoint multicentre clinical trial. Setting Four hundred and twenty-four UK primary care practices. Participants Aged 60 years and over with ischaemic heart disease but no gout. Interventions Participants were randomised (1 : 1) using a central web-based randomisation system to receive allopurinol up to 600 mg daily that was added to usual care or to continue usual care. Main outcome measures The primary outcome was the composite of non-fatal myocardial infarction, non-fatal stroke or cardiovascular death. Secondary outcomes were non-fatal myocardial infarction, non-fatal stroke, cardiovascular death, all-cause mortality, hospitalisation for heart failure, hospitalisation for acute coronary syndrome, coronary revascularisation, hospitalisation for acute coronary syndrome or coronary revascularisation, all cardiovascular hospitalisations, quality of life and cost-effectiveness. The hazard ratio (allopurinol vs. usual care) in a Cox proportional hazards model was assessed for superiority in a modified intention-to-treat analysis. Results From 7 February 2014 to 2 October 2017, 5937 participants were enrolled and randomised to the allopurinol arm (n = 2979) or the usual care arm (n = 2958). A total of 5721 randomised participants (2853 allopurinol; 2868 usual care) were included in the modified intention-to-treat analysis population (mean age 72.0 years; 75.5% male). There was no difference between the allopurinol and usual care arms in the primary endpoint, 314 (11.0%) participants in the allopurinol arm (2.47 events per 100 patient-years) and 325 (11.3%) in the usual care arm (2.37 events per 100 patient-years), hazard ratio 1.04 (95% confidence interval 0.89 to 1.21); p = 0.65. Two hundred and eighty-eight (10.1%) participants in the allopurinol arm and 303 (10.6%) participants in the usual care arm died, hazard ratio 1.02 (95% confidence interval 0.87 to 1.20); p = 0.77. The pre-specified health economic analysis plan was to perform a ‘within trial’ cost-utility analysis if there was no statistically significant difference in the primary endpoint, so NHS costs and quality-adjusted life-years were estimated over a 5-year period. The difference in costs between treatment arms was +£115 higher for allopurinol (95% confidence interval £17 to £210) with no difference in quality-adjusted life-years (95% confidence interval −0.061 to +0.060). We conclude that there is no evidence that allopurinol used in line with the study protocol is cost-effective. Limitations The results may not be generalisable to younger populations, other ethnic groups or patients with more acute ischaemic heart disease. One thousand six hundred and thirty-seven participants (57.4%) in the allopurinol arm withdrew from randomised treatment, but an on-treatment analysis gave similar results to the main analysis. Conclusions The ALL-HEART study showed that treatment with allopurinol 600 mg daily did not improve cardiovascular outcomes compared to usual care in patients with ischaemic heart disease. We conclude that allopurinol should not be recommended for the secondary prevention of cardiovascular events in patients with ischaemic heart disease but no gout. Future work The effects of allopurinol on cardiovascular outcomes in patients with ischaemic heart disease and co-existing hyperuricaemia or clinical gout could be explored in future studies. Trial registration This trial is registered as EU Clinical Trials Register (EudraCT 2013-003559-39) and ISRCTN (ISRCTN 32017426). Funding This award was funded by the National Institute for Health and Care Research (NIHR) Health Technology Assessment programme (NIHR award ref: 11/36/41) and is published in full in Health Technology Assessment; Vol. 28, No. 18. See the NIHR Funding and Awards website for further award information. Plain language summary What was the question? The purpose of the ALL-HEART study was to determine whether giving allopurinol to people with ischaemic heart disease (also commonly known as coronary heart disease) would reduce their risk of having a heart attack, stroke or of dying from cardiovascular disease. Allopurinol is a medication usually given to patients with gout to prevent acute gout flares. It is not currently used to treat ischaemic heart disease. What did we do? We randomly allocated people aged over 60 years with ischaemic heart disease to take up to 600 mg of allopurinol daily (in addition to their usual care) or to continue with their usual care. We then monitored participants for several years and recorded any major health events such as heart attacks, strokes and deaths. We obtained most of the follow-up data from centrally held electronic hospital admissions and death records, making the study easier for participants and more cost-efficient. We asked participants in both groups to complete questionnaires to assess their quality of life during the study. We also collected data to determine whether there was any economic benefit to the NHS of using allopurinol in patients with ischaemic heart disease. What did we find? There was no difference in the risk of heart attacks, strokes or death from cardiovascular disease between the participants given allopurinol and those in the group continuing their usual care. We also found no difference in the risks of other cardiovascular events, deaths from any cause or quality-of-life measurements between the allopurinol and usual care groups. What does this mean? The results of the ALL-HEART study suggest that we should not recommend that allopurinol be given to people with ischaemic heart disease to prevent further cardiovascular events or deaths. Scientific summary Background Allopurinol is a xanthine oxidase inhibitor that lowers serum uric acid (SUA) and is widely used in patients with gout to prevent acute gout flares. Xanthine oxidase promotes inflammation and atherosclerosis via the production of reactive oxygen species and xanthine oxidase levels are raised in several conditions including coronary artery disease. The role of SUA in cardiovascular (CV) disease is controversial, with some studies associating higher SUA levels with worse CV outcomes, but more recent genome-wide association studies suggest no major role of uric acid levels in determining CV outcomes. Some observational studies have suggested that allopurinol therapy may improve CV outcomes, while others have not found an association. Small interventional studies have shown that allopurinol therapy improves some CV parameters, including endothelial function, left ventricular hypertrophy, blood pressure, carotid intimal media thickness and arterial stiffness. Allopurinol therapy was also found to improve outcomes after acute coronary syndrome (ACS) in one study and to improve chest pain in patients with chronic stable angina with documented coronary artery disease in another. However, results have not been consistent across different studies. Before the ALL-HEART study, no large, randomised trial of the effects of allopurinol therapy on CV outcomes in patients with ischaemic heart disease (IHD) had been performed. Objectives Primary Does allopurinol therapy added to usual care improve major CV outcomes in patients aged over 60 years with IHD but no gout? Secondary Does allopurinol therapy added to usual care improve all-cause mortality or other CV outcomes in patients with IHD? What is the cost-effectiveness of adding allopurinol up to 600 mg daily to usual care in patients with IHD? Does allopurinol therapy added to usual care improve quality of life assessed by general health survey [EuroQol-5 Dimensions, five-level version (EQ-5D-5L)] or coronary heart disease-specific questionnaire (Seattle angina questionnaire)? Methods Design and participants The ALL-HEART study was a prospective, randomised, open-label, blinded endpoint (PROBE) multicentre trial undertaken in patients with IHD. Participants were primarily recruited from 424 primary care practices via 18 regional centres in the UK, with a small number also referred into the study from secondary care centres. Eligible patients were aged 60 years or over, with a history of IHD [myocardial infarction (MI), angina or other evidence of IHD]. Exclusion criteria were: history of gout; known severe renal impairment [estimated glomerular filtration rate (eGFR) < 30 ml/minute/1.73 m2]; moderate-to-severe heart failure (HF) [New York Heart Association (NYHA) III–IV]; significant hepatic disease [e.g. alanine transaminase (ALT) > 3 × upper limit of normal, cirrhosis, ascites] (investigator opinion); currently taking part in another interventional clinical trial of an investigational medicinal product or medical device (or taken part in one within the last 3 months); previous allergy to allopurinol; previous serious adverse cutaneous (skin) reaction to any drug (e.g. Stevens–Johnson syndrome, toxic epidermal necrolysis, hospitalisation due to skin reaction to drug) (investigator opinion); already taking urate-lowering therapy (including allopurinol, febuxostat, sulfinpyrazone, benzbromarone, probenecid, rasburicase); taking azathioprine, mercaptopurine, ciclosporin or theophylline; malignancy (except non-metastatic, non-melanoma skin cancers, cervical in situ carcinoma, breast ductal carcinoma in situ, or stage 1 prostate carcinoma) within the last 5 years (investigator opinion). The exclusion criterion relating to renal impairment was originally ‘known renal impairment eGFR < 60 ml/minute/1.73 m2’ for patients recruited from the start of the trial (7 February 2014) until 4 April 2016 when an updated version of the protocol was implemented at all study sites to allow the inclusion of patients with moderate renal impairment in the study with the purpose of making the study results more generalisable. Fifty-two per cent of the target number of patients had been randomised by this date. Randomisation At a single screening and randomisation visit usually held at the patient’s primary care practice by a research nurse, patients were consented, screened and randomised to receive allopurinol therapy or to continue their usual care. Baseline demographics, medical history, CV risk factors and concomitant medications were recorded. Blood pressure, height and weight were measured. Baseline blood tests were taken for urea, creatinine and electrolytes, full blood count and SUA. Participants were randomised before screening blood results were available. When the screening results were available, a nurse telephoned the patient to advise them to start taking randomised therapy. If the screening visit eGFR result was below the exclusion limit, the participant did not receive any allopurinol and was excluded from the modified intention-to-treat (mITT) analysis population, whichever arm of the study they had been randomised to. Randomisation was via a web-based randomisation facility located at the Robertson Centre for Biostatistics, University of Glasgow, accessed using a web-based application or an interactive voice response system. Randomisation was based on randomised permuted blocks of variable size, stratified by history of MI, history of stroke and primary care practice. Randomised intervention Allopurinol 100 or 300 mg tablets were prescribed to participants by their primary care physicians via the usual NHS prescription system. Participants with a screening eGFR result ≥ 60 ml/minute/1.73 m2 were prescribed allopurinol at 100 mg daily for 2 weeks, then 300 mg daily for 2 weeks then 600 mg daily (given as 300 mg twice daily) thereafter if tolerated for the duration of the study. After 4 April 2016, participants with a screening eGFR result 30–59 ml/minute/1.73 m2 were prescribed allopurinol at 100 mg daily for 2 weeks, then 300 mg daily thereafter if tolerated for the duration of the study. If there were tolerability issues, the dose could be decreased at the discretion of a physician. Participants who stopped randomised therapy were encouraged to continue with study follow-up. The comparison arm of the study was ‘usual care’; no placebo was given as this was a pragmatic open-label study. Randomised therapy was blinded to the endpoint adjudication committee, but not to participants, study staff or treating healthcare professionals. If allopurinol was started in participants randomised to the usual care group at any point in the study (e.g. for clinical reasons such as developing gout), this was recorded. Other study procedures At the screening visit, participants completed the EQ-5D-5L questionnaire to assess general health outcomes and the Seattle angina questionnaire to assess coronary artery disease-specific quality of life. Participants who had taken any allopurinol therapy attended a 6-week study visit. Blood samples were taken for urea, creatinine and electrolytes, full blood count and SUA. Participants were then followed up remotely by electronic record-linkage with centralised databases of hospitalisations, cancers and deaths held by Public Health Scotland (PHS) and NHS Digital and by annual questionnaires (online, by post or by telephone). Data on adverse events, skin rashes, gout flares and self-reported adherence to randomised therapy were collected at the 6-week visit, then by annual questionnaire (but could also be reported at any time by participants or health professionals). Participants completed the EQ-5D-5L and Seattle angina questionnaires again after 1 year and at the end of the trial. The trial recruitment period was extended once, and the trial follow-up period was extended twice. Recruitment exceeded the target of 5215 randomised participants to a final total of 5937 randomised participants due to an increase in recruitment rate in the last weeks of recruitment. The follow-up period of the trial ended on 30 September 2021. After this date, participants stopped randomised therapy and continued to receive their usual care. Study outcomes The primary and secondary outcomes were as follows. Primary outcome: Composite of non-fatal MI, non-fatal stroke or CV death. Secondary outcomes: Non-fatal MI. Non-fatal stroke. CV death. All-cause mortality. Hospitalisation for ACS. Coronary revascularisation. Hospitalisation for ACS or coronary revascularisation. Hospitalisation for HF. All CV hospitalisations. Quality of life. Cost-effectiveness. Adverse events Serious adverse events (SAEs) occurring during the study (until 30 September 2021) were recorded and any events ongoing at that time were followed up for a further 30 days, unless consent had been withdrawn. Treatment-related adverse events (adverse reactions), gout flares and rashes were also recorded. Allopurinol therapy was stopped if participants developed a rash that may have been due to allopurinol. For any SAEs that were potential study endpoints, detailed information was obtained from medical records and death certificates to support the production of an anonymised endpoint package that was adjudicated by an independent clinical endpoint adjudication committee, blinded to randomised therapy allocation. Statistical analysis The power calculation suggested that 5215 participants would need to be randomised 1 : 1 to give 80% power to detect a 20% reduction in the primary outcome for the intervention (allowing for a 4% dropout for withdrawal of consent and non-CV deaths). A primary event rate of 14% over an average of 4 years follow-up was estimated from other trials in similar patient groups. The study ended when 631 adjudicated first primary events had occurred. Baseline characteristics are presented by treatment group as means [standard deviation (SD)] and medians [interquartile range (IQR)] for continuous variables and as numbers (%) for categorical variables. Clinical outcomes were analysed on a time-to-first event basis using Cox proportional hazards models. Treatment effects (allopurinol vs. usual care) were estimated as hazard ratios (HR) and 95% confidence intervals (CIs) for the Cox models. Analyses were adjusted for the stratification variables history of MI and history of stroke and p-values were calculated from Wald statistics. The primary analysis was a mITT analysis. Health economic analysis A health economic analysis was planned to determine whether allopurinol was a cost-effective intervention in the context of the NHS setting in the UK in patients with IHD. The analysis plan was based around NHS costs and estimation of quality-adjusted life-years (QALYs) and was to include a ‘within trial’ cost–utility analysis if the primary outcome was not statistically different between randomised study arms. Trial approvals and committees The study was approved by an ethics committee, Medicines and Healthcare Products Regulatory Agency (MHRA) and Health Research Authority (HRA). An Independent Data Monitoring Committee oversaw trial safety. A Trial Steering Committee including independent and patient/lay members oversaw trial progress. Results From 7 February 2014 to 29 September 2017, 6134 patients consented to enrol in the trial and were assessed for eligibility. Also, 167 of the 6134 consented participants were not eligible and 30 others were not randomised. The final randomisation was completed on 2 October 2017. Of the 5937 randomised participants, 216 were excluded post randomisation from the mITT analysis population (184 did not meet the eGFR entry criteria once their screening blood results were available; 32 were later found to have not met all of the inclusion/exclusion criteria). Five thousand seven hundred and twenty-one participants (2853 in the allopurinol arm and 2868 in the usual care arm) were included in the mITT analysis population for the efficacy analysis. The population for safety analyses consisted of 2805 participants in the allopurinol arm (excluding the 48 participants in the allopurinol arm who never received any randomised therapy) and 2868 participants in the usual care arm. The mean duration of follow-up in the trial was 4.8 years. Two hundred and fifty-eight participants (9.0%) in the allopurinol group and 76 participants (2.6%) in the usual care group withdrew consent for all follow-up. One thousand six hundred and thirty-seven participants (57.4%) in the allopurinol arm withdrew from randomised treatment. Baseline characteristics were well balanced in the two groups. The mean age at study entry was 72.0 years (SD 6.8), 4321 participants (75.5%) were male, 5676 (99.2%) were white and 1241 (21.7%) had a history of diabetes mellitus. The median duration of IHD at study entry was 10.1 years (IQR 5.1–16.1). Three thousand four hundred and sixty-four (60.5%) participants were recruited in England and 2257 (39.5%) participants in Scotland. The most commonly taken dose of allopurinol was 600 mg daily. In the 2447 participants in the allopurinol arm with both a baseline and 6-week SUA result, SUA fell from a mean of 0.34 (SD 0.08) mmol/l at baseline to a mean of 0.18 (SD 0.09) mmol/l 6 weeks after randomisation. Forty-five participants in the usual care arm started allopurinol therapy during follow-up (for clinical reasons, mainly gout). Primary and secondary outcomes There was no significant difference between the randomised treatment groups in the rates of the primary outcome or any of the secondary time-to-event outcomes. Three hundred and fourteen (11.0%) participants in the allopurinol arm (2.47 events per 100 patient-years) and 325 (11.3%) in the usual care arm (2.37 events per 100 patient-years) experienced a primary outcome, HR 1.04 (95% CI 0.89 to 1.21); p = 0.65. Two hundred and eighty-eight (10.1%) participants in the allopurinol arm and 303 (10.6%) participants in the usual care arm died, HR 1.02 (95% CI 0.87 to 1.20); p = 0.77. Results for the primary outcome were consistent across all pre-specified subgroups. In a supporting on-treatment analysis, results for the time-to-event clinical outcomes were broadly similar to those in the mITT analysis. There was limited evidence of any effect of allopurinol on quality-of-life outcomes, with no differences in EQ-5D-5L outcomes or Seattle angina questionnaire outcomes at the end of the first year or at the final visit, except for a nominally significant but only slightly greater fall in the physical domain score of the Seattle angina questionnaire at the end of the first year [treatment difference = 1.219 (95% CI 0.027 to 2.410); p = 0.045] in the allopurinol arm. Health economic analysis There was strong evidence that allopurinol treatment was associated with incremental costs relative to usual care [incremental cost per patient £115.4, 95% CI (£17.0 to £210.2)], with little evidence of improvement in relation to incremental QALYs [incremental QALYs −0.000, 95% CI (−0.061 to 0.060)]. The cost-effectiveness acceptability curve asymptoted at a probability of 0.456, meaning that even with a willingness to pay of an infinite amount, the probability of cost-effectiveness was only 0.465. At a willingness to pay of £20,000 the probability of cost-effectiveness was 0.41. Adverse events There was no difference in SAE rates between treatment arms, except for the grouping of endocrine disorders where no events occurred in the allopurinol treatment group but 14 in the usual care group. However, these endocrine events included events with a spread of different types. Fifteen participants had SAEs that were considered potentially treatment related and none of these were fatal. There was no difference between treatment arms in the rates of incident cancers. Adjudicated causes of death were well balanced between the treatment groups, including deaths from COVID-19 and COVID-19 pneumonia. Conclusions The ALL-HEART study showed that treatment with allopurinol 600 mg daily did not improve CV outcomes compared to usual care in patients with IHD. There were also no benefits on quality of life. There was no evidence that allopurinol used in line with the study protocol is cost-effective within the NHS system. We conclude that allopurinol should not be recommended for the secondary prevention of CV events in patients with IHD but no gout. Recommendations for research Future research should explore other therapeutic options for the improvement of CV outcomes and quality of life in patients with IHD. Further exploration of the effects of allopurinol on CV outcomes in patients with co-existing IHD and clinical gout or hyperuricaemia could be considered (patients with gout were excluded from this study). Trial registration The ALL-HEART trial is registered with the EU Clinical Trials Register (EudraCT 2013-003559-39) and ISRCTN (ISRCTN 32017426). Funding This award was funded by the National Institute for Health and Care Research (NIHR) Health Technology Assessment programme (NIHR award ref: 11/36/41) and is published in full in Health Technology Assessment; Vol. 28, No. 18. See the NIHR Funding and Awards website for further award information.

Smoking and HPV infection are known causes for the vast majority of head and neck squamous cell carcinomas (HNSCC) due to their likelihood of causing gene dysregulation and genomic alterations. Enhancer RNAs (eRNAs) are non-coding RNAs that are known to increase nearby and target gene expression, and activity that has been suggested to be affected by genetic and epigenetic alterations. Here we sought to identify the effects of smoking and HPV status on eRNA expression in HNSCC tumors. We focused on four patient cohorts including smoking/HPV+, smoking/HPV-, non-smoking/HPV+, and non-smoking/HPV- patients. We used TCGA RNA-seq data from cancer tumors and adjacent normal tissue, extracted eRNA read counts, and correlated these to survival, clinical variables, immune infiltration, cancer pathways, and genomic alterations. We found a large number of differentially expressed eRNA in each patient cohort. We also found several dysregulated eRNA correlated to patient survival, clinical variables, immune pathways, and genomic alterations. Additionally, we were able to find dysregulated eRNA nearby seven key HNSCC-related oncogenes. For example, we found eRNA chr14:103272042-103272430 (eRNA-24036), which is located close to the TRAF3 gene to be differentially expressed and correlated with the pathologic N stage and immune cell populations. Using a separate validation dataset, we performed differential expression and immune infiltration analysis to validate our results from the TCGA data. Our findings may explain the association between eRNA expression, enhancer activity, and nearby gene dysregulation.

Single-cell RNA sequencing (scRNA-seq) technologies facilitate the characterization of transcriptomic landscapes in diverse species, tissues, and cell types with unprecedented molecular resolution. In order to evaluate various biological hypotheses using high-dimensional single-cell gene expression data, most computational and statistical methods depend on a gene feature selection step to identify genes with high biological variability and reduce computational complexity. Even though many gene selection methods have been developed for scRNA-seq analysis, there lacks a systematic comparison of the assumptions, statistical models, and selection criteria used by these methods. In this article, we summarize and discuss 17 computational methods for selecting gene features in unsupervised analysis of single-cell gene expression data, with unified notations and statistical frameworks. Our discussion provides a useful summary to help practitioners select appropriate methods based on their assumptions and applicability, and to assist method developers in designing new computational tools for unsupervised learning of scRNA-seq data.

UTX (also known as KDM6A) encodes a histone H3K27 demethylase and is an important tumour suppressor that is frequently mutated in human cancers1. However, as the demethylase activity of UTX is often dispensable for mediating tumour suppression and developmental regulation2-8, the underlying molecular activity of UTX remains unknown. Here we show that phase separation of UTX underlies its chromatin-regulatory activity in tumour suppression. A core intrinsically disordered region (cIDR) of UTX forms phase-separated liquid condensates, and cIDR loss caused by the most frequent cancer mutation of UTX is mainly responsible for abolishing tumour suppression. Deletion, mutagenesis and replacement assays of the intrinsically disordered region demonstrate a critical role of UTX condensation in tumour suppression and embryonic stem cell differentiation. As shown by reconstitution in vitro and engineered systems in cells, UTX recruits the histone methyltransferase MLL4 (also known as KMT2D) to the same condensates and enriches the H3K4 methylation activity of MLL4. Moreover, UTX regulates genome-wide histone modifications and high-order chromatin interactions in a condensation-dependent manner. We also found that UTY, the Y chromosome homologue of UTX with weaker tumour-suppressive activity, forms condensates with reduced molecular dynamics. These studies demonstrate a crucial biological function of liquid condensates with proper material states in enabling the tumour-suppressive activity of a chromatin regulator.

Although parental HIV disclosure has benefits for parents and children, the disclosure rate among parents remains low. This study aims to qualitatively examine parental concerns regarding disclosure of their HIV status to their children. Eighty parents were enrolled in a randomized controlled trial of a three-session disclosure-support intervention, with forty receiving the intervention and forty receiving treatment as usual. Intervention sessions were audio recorded, and transcriptions were qualitatively coded for content related to concerns of disclosure. Four themes emerged: Intention to disclose, disclosure approach, indicators for disclosure, and fears about disclosure. These themes reveal struggles that parents experience when considering HIV disclosure suggesting that an effective disclosure intervention must help parents assess pros and cons, discuss the emotions of the children after the disclosure, and monitor the impact on children's lives after disclosure over time. Future research is needed to implement interventions supporting HIV-positive parents' disclosure decision-making and actions.

Current understanding holds that Klinefelter syndrome (KS) is not inherited, but arises randomly during meiosis. Whether there is any genetic basis for the origin of KS is unknown. Here, guided by our identification of some USP26 variations apparently associated with KS, we found that knockout of Usp26 in male mice resulted in the production of 41, XXY offspring. USP26 protein is localized at the XY body, and the disruption of Usp26 causes incomplete sex chromosome pairing by destabilizing TEX11. The unpaired sex chromosomes then result in XY aneuploid spermatozoa. Consistent with our mouse results, a clinical study shows that some USP26 variations increase the proportion of XY aneuploid spermatozoa in fertile men, and we identified two families with KS offspring wherein the father of the KS patient harbored a USP26-mutated haplotype, further supporting that paternal USP26 mutation can cause KS offspring production. Thus, some KS should originate from XY spermatozoa, and paternal USP26 mutations increase the risk of producing KS offspring.

OBJECTIVE: To measure the global impact of COVID-19 pandemic on volumes of IV thrombolysis (IVT), IVT transfers, and stroke hospitalizations over 4 months at the height of the pandemic (March 1 to June 30, 2020) compared with 2 control 4-month periods. METHODS: We conducted a cross-sectional, observational, retrospective study across 6 continents, 70 countries, and 457 stroke centers. Diagnoses were identified by their ICD-10 codes or classifications in stroke databases. RESULTS: There were 91,373 stroke admissions in the 4 months immediately before compared to 80,894 admissions during the pandemic months, representing an 11.5% (95% confidence interval [CI] -11.7 to -11.3, p < 0.0001) decline. There were 13,334 IVT therapies in the 4 months preceding compared to 11,570 procedures during the pandemic, representing a 13.2% (95% CI -13.8 to -12.7, p < 0.0001) drop. Interfacility IVT transfers decreased from 1,337 to 1,178, or an 11.9% decrease (95% CI -13.7 to -10.3, p = 0.001). Recovery of stroke hospitalization volume (9.5%, 95% CI 9.2-9.8, p < 0.0001) was noted over the 2 later (May, June) vs the 2 earlier (March, April) pandemic months. There was a 1.48% stroke rate across 119,967 COVID-19 hospitalizations. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection was noted in 3.3% (1,722/52,026) of all stroke admissions. CONCLUSIONS: The COVID-19 pandemic was associated with a global decline in the volume of stroke hospitalizations, IVT, and interfacility IVT transfers. Primary stroke centers and centers with higher COVID-19 inpatient volumes experienced steeper declines. Recovery of stroke hospitalization was noted in the later pandemic months.

Patients with underlying cardiovascular conditions are particularly vulnerable to severe COVID-19. In this project, we aimed to characterize similarities in dysregulated immune pathways between COVID-19 patients and patients with cardiomyopathy, venous thromboembolism (VTE), or coronary artery disease (CAD). We hypothesized that these similarly dysregulated pathways may be critical to how cardiovascular diseases (CVDs) exacerbate COVID-19. To evaluate immune dysregulation in different diseases, we used four separate datasets, including RNA-sequencing data from human left ventricular cardiac muscle samples of patients with dilated or ischemic cardiomyopathy and healthy controls; RNA-sequencing data of whole blood samples from patients with single or recurrent event VTE and healthy controls; RNA-sequencing data of human peripheral blood mononuclear cells (PBMCs) from patients with and without obstructive CAD; and RNA-sequencing data of platelets from COVID-19 subjects and healthy controls. We found similar immune dysregulation profiles between patients with CVDs and COVID-19 patients. Interestingly, cardiomyopathy patients display the most similar immune landscape to COVID-19 patients. Additionally, COVID-19 patients experience greater upregulation of cytokine- and inflammasome-related genes than patients with CVDs. In all, patients with CVDs have a significant overlap of cytokine- and inflammasome-related gene expression profiles with that of COVID-19 patients, possibly explaining their greater vulnerability to severe COVID-19.

A pressing challenge in single-cell transcriptomics is to benchmark experimental protocols and computational methods. A solution is to use computational simulators, but existing simulators cannot simultaneously achieve three goals: preserving genes, capturing gene correlations, and generating any number of cells with varying sequencing depths. To fill this gap, we propose scDesign2, a transparent simulator that achieves all three goals and generates high-fidelity synthetic data for multiple single-cell gene expression count-based technologies. In particular, scDesign2 is advantageous in its transparent use of probabilistic models and its ability to capture gene correlations via copulas.

Promoter DNA methylation is a well-established mechanism of transcription repression, though its global correlation with gene expression is weak. This weak correlation can be attributed to the failure of current methylation quantification methods to consider the heterogeneity among sequenced bulk cells. Here, we introduce Cell Heterogeneity-Adjusted cLonal Methylation (CHALM) as a methylation quantification method. CHALM improves understanding of the functional consequences of DNA methylation, including its correlations with gene expression and H3K4me3. When applied to different methylation datasets, the CHALM method enables detection of differentially methylated genes that exhibit distinct biological functions supporting underlying mechanisms.

The mammalian DNA methylome is formed by two antagonizing processes, methylation by DNA methyltransferases (DNMT) and demethylation by ten-eleven translocation (TET) dioxygenases. Although the dynamics of either methylation or demethylation have been intensively studied in the past decade, the direct effects of their interaction on gene expression remain elusive. Here, we quantify the concurrence of DNA methylation and demethylation by the percentage of unmethylated CpGs within a partially methylated read from bisulfite sequencing. After verifying 'methylation concurrence' by its strong association with the co-localization of DNMT and TET enzymes, we observe that methylation concurrence is strongly correlated with gene expression. Notably, elevated methylation concurrence in tumors is associated with the repression of 40~60% of tumor suppressor genes, which cannot be explained by promoter hypermethylation alone. Furthermore, methylation concurrence can be used to stratify large undermethylated regions with negligible differences in average methylation into two subgroups with distinct chromatin accessibility and gene regulation patterns. Together, methylation concurrence represents a unique methylation metric important for transcription regulation and is distinct from conventional metrics, such as average methylation and methylation variation.

High-throughput biological data analysis commonly involves identifying features such as genes, genomic regions, and proteins, whose values differ between two conditions, from numerous features measured simultaneously. The most widely used criterion to ensure the analysis reliability is the false discovery rate (FDR), which is primarily controlled based on p-values. However, obtaining valid p-values relies on either reasonable assumptions of data distribution or large numbers of replicates under both conditions. Clipper is a general statistical framework for FDR control without relying on p-values or specific data distributions. Clipper outperforms existing methods for a broad range of applications in high-throughput data analysis.

Perspiration evaporation plays an indispensable role in human body heat dissipation. However, conventional textiles tend to focus on sweat removal and pay little attention to the basic thermoregulation function of sweat, showing limited evaporation ability and cooling efficiency in moderate/profuse perspiration scenarios. Here, we propose an integrated cooling (i-Cool) textile with unique functional structure design for personal perspiration management. By integrating heat conductive pathways and water transport channels decently, i-Cool exhibits enhanced evaporation ability and high sweat evaporative cooling efficiency, not merely liquid sweat wicking function. In the steady-state evaporation test, compared to cotton, up to over 100% reduction in water mass gain ratio, and 3 times higher skin power density increment for every unit of sweat evaporation are demonstrated. Besides, i-Cool shows about 3 °C cooling effect with greatly reduced sweat consumption than cotton in the artificial sweating skin test. The practical application feasibility of i-Cool design principles is well validated based on commercial fabrics. Owing to its exceptional personal perspiration management performance, we expect the i-Cool concept can provide promising design guidelines for next-generation perspiration management textiles.