Identification of the cis-regulatory elements controlling cell-type specific gene expression patterns is essential for understanding the origin of cellular diversity. Conventional assays to map regulatory elements via open chromatin analysis of primary tissues is hindered by sample heterogeneity. Single cell analysis of accessible chromatin (scATAC-seq) can overcome this limitation. However, the high-level noise of each single cell profile and the large volume of data pose unique computational challenges. Here, we introduce SnapATAC, a software package for analyzing scATAC-seq datasets. SnapATAC dissects cellular heterogeneity in an unbiased manner and map the trajectories of cellular states. Using the Nyström method, SnapATAC can process data from up to a million cells. Furthermore, SnapATAC incorporates existing tools into a comprehensive package for analyzing single cell ATAC-seq dataset. As demonstration of its utility, SnapATAC is applied to 55,592 single-nucleus ATAC-seq profiles from the mouse secondary motor cortex. The analysis reveals ~370,000 candidate regulatory elements in 31 distinct cell populations in this brain region and inferred candidate cell-type specific transcriptional regulators.

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

We still know very little about how the human immune system responds to SARS-CoV-2. Here we construct a SARS-CoV-2 proteome microarray containing 18 out of the 28 predicted proteins and apply it to the characterization of the IgG and IgM antibodies responses in the sera from 29 convalescent patients. We find that all these patients had IgG and IgM antibodies that specifically bind SARS-CoV-2 proteins, particularly the N protein and S1 protein. Besides these proteins, significant antibody responses to ORF9b and NSP5 are also identified. We show that the S1 specific IgG signal positively correlates with age and the level of lactate dehydrogenase (LDH) and negatively correlates with lymphocyte percentage. Overall, this study presents a systemic view of the SARS-CoV-2 specific IgG and IgM responses and provides insights to aid the development of effective diagnostic, therapeutic and vaccination strategies.

Tandem repeats (TRs) represent one of the largest sources of genetic variation in humans and are implicated in a range of phenotypes. Here we present a deep characterization of TR variation based on high coverage whole genome sequencing from 3550 diverse individuals from the 1000 Genomes Project and H3Africa cohorts. We develop a method, EnsembleTR, to integrate genotypes from four separate methods resulting in high-quality genotypes at more than 1.7 million TR loci. Our catalog reveals novel sequence features influencing TR heterozygosity, identifies population-specific trinucleotide expansions, and finds hundreds of novel eQTL signals. Finally, we generate a phased haplotype panel which can be used to impute most TRs from nearby single nucleotide polymorphisms (SNPs) with high accuracy. Overall, the TR genotypes and reference haplotype panel generated here will serve as valuable resources for future genome-wide and population-wide studies of TRs and their role in human phenotypes.

Carbapenem-resistant Escherichia coli (CREC) ST410 has recently emerged as a major global health problem. Here, we report a shift in CREC prevalence in Chinese hospitals between 2017 and 2021 with ST410 becoming the most commonly isolated sequence type. Genomic analysis identifies a hypervirulent CREC ST410 clone, B5/H24RxC, which caused two separate outbreaks in a children's hospital. It may have emerged from the previously characterised B4/H24RxC in 2006 and has been isolated in ten other countries from 2015 to 2021. Compared with B4/H24RxC, B5/H24RxC lacks the blaOXA-181-bearing X3 plasmid, but carries a F-type plasmid containing blaNDM-5. Most of B5/H24RxC also carry a high pathogenicity island and a novel O-antigen gene cluster. We find that B5/H24RxC grew faster in vitro and is more virulent in vivo. The identification of this newly emerged but already globally disseminated hypervirulent CREC clone, highlights the ongoing evolution of ST410 towards increased resistance and virulence.

X-chromosome inactivation (XCI), i.e., the inactivation of one of the female X chromosomes, restores equal expression of X-chromosomal genes between females and males. However, ~10% of genes show variable degrees of escape from XCI between females, although little is known about the causes of variable XCI. Using a discovery data-set of 1867 females and 1398 males and a replication sample of 3351 females, we show that genetic variation at three autosomal loci is associated with female-specific changes in X-chromosome methylation. Through cis-eQTL expression analysis, we map these loci to the genes SMCHD1/METTL4, TRIM6/HBG2, and ZSCAN9. Low-expression alleles of the loci are predominantly associated with mild hypomethylation of CpG islands near genes known to variably escape XCI, implicating the autosomal genes in variable XCI. Together, these results suggest a genetic basis for variable escape from XCI and highlight the potential of a population genomics approach to identify genes involved in XCI.

Human cytomegalovirus (HCMV) is a widespread pathogen that in immunocompromised hosts can cause life-threatening disease. Studying HCMV-exposed monocyte-derived dendritic cells by single-cell RNA sequencing, we observe that most cells are entered by the virus, whereas less than 30% of them initiate viral gene expression. Increased viral gene expression is associated with activation of the stimulator of interferon genes (STING) that usually induces anti-viral interferon responses, and with the induction of several pro- (RHOB, HSP1A1, DNAJB1) and anti-viral (RNF213, TNFSF10, IFI16) genes. Upon progression of infection, interferon-beta but not interferon-lambda transcription is inhibited. Similarly, interferon-stimulated gene expression is initially induced and then shut off, thus further promoting productive infection. Monocyte-derived dendritic cells are composed of 3 subsets, with one being especially susceptible to HCMV. In conclusion, HCMV permissiveness of monocyte-derived dendritic cells depends on complex interactions between virus sensing, regulation of the interferon response, and viral gene expression.

Hepatocellular carcinoma (HCC) cells often invade the portal venous system and subsequently develop into portal vein tumour thrombosis (PVTT). Long noncoding RNAs (lncRNAs) have been associated with HCC, but a comprehensive analysis of their specific association with HCC metastasis has not been conducted. Here, by analysing 60 clinical samples' RNA-seq data from 20 HCC patients, we have identified and characterized 8,603 candidate lncRNAs. The expression patterns of 917 recurrently deregulated lncRNAs are correlated with clinical data in a TCGA cohort and published liver cancer data. Matched array data from the 60 samples show that copy number variations (CNVs) and alterations in DNA methylation contribute to the observed recurrent deregulation of 235 lncRNAs. Many recurrently deregulated lncRNAs are enriched in co-expressed clusters of genes related to cell adhesion, immune response and metabolic processes. Candidate lncRNAs related to metastasis, such as HAND2-AS1, were further validated using RNAi-based loss-of-function assays. Thus, we provide a valuable resource of functional lncRNAs and biomarkers associated with HCC tumorigenesis and metastasis.

Cancer stem cells (CSCs) play a critical role in invasive growth and metastasis of human head and neck squamous cell carcinoma (HNSCC). Although significant progress has been made in understanding the self-renewal and pro-tumorigenic potentials of CSCs, a key challenge remains on how to eliminate CSCs and halt metastasis effectively. Here we show that super-enhancers (SEs) play a critical role in the transcription of cancer stemness genes as well as pro-metastatic genes, thereby controlling their tumorigenic potential and metastasis. Mechanistically, we find that bromodomain-containing protein 4 (BRD4) recruits Mediators and NF-κB p65 to form SEs at cancer stemness genes such as TP63, MET and FOSL1, in addition to oncogenic transcripts. In vivo lineage tracing reveals that disrupting SEs by BET inhibitors potently inhibited CSC self-renewal and eliminated CSCs in addition to elimination of proliferating non-stem tumor cells in a mouse model of HNSCC. Moreover, disrupting SEs also inhibits the invasive growth and lymph node metastasis of human CSCs isolated from human HNSCC. Taken together, our results suggest that targeting SEs may serve as an effective therapy for HNSCC by eliminating CSCs.

Chronic lymphocytic leukaemia (CLL) is a frequent B-cell malignancy, characterized by recurrent somatic chromosome alterations and a low level of point mutations. Here we present single-nucleotide polymorphism microarray analyses of a single CLL patient over 29 years of observation and treatment, and transcriptome and whole-genome sequencing at selected time points. We identify chromosome alterations 13q14-, 6q- and 12q+ in early cell clones, elimination of clonal populations following therapy, and subsequent appearance of a clone containing trisomy 12 and chromosome 10 copy-neutral loss of heterogeneity that marks a major population dominant at death. Serial single-cell RNA sequencing reveals an expression pattern with high FOS, JUN and KLF4 at disease acceleration, which resolves following therapy, but reoccurs following relapse and death. Transcriptome evolution indicates complex changes in expression occur over time. In conclusion, CLL can evolve gradually during indolent phases, and undergo rapid changes following therapy.