The RNAs bound to the genome and their binding sites are detected with GRID-seq. Higher eukaryotic genomes are bound by a large number of coding and non-coding RNAs, but approaches to comprehensively map the identity and binding sites of these RNAs are lacking. Here we report a method to capture in situ global RNA interactions with DNA by deep sequencing (GRID-seq), which enables the comprehensive identification of the entire repertoire of chromatin-interacting RNAs and their respective binding sites. In human, mouse, and Drosophila cells, we detected a large set of tissue-specific coding and non-coding RNAs that are bound to active promoters and enhancers, especially super-enhancers. Assuming that most mRNA–chromatin interactions indicate the physical proximity of a promoter and an enhancer, we constructed a three-dimensional global connectivity map of promoters and enhancers, revealing transcription-activity-linked genomic interactions in the nucleus.

Development of liver fibrosis is associated with activation of quiescent hepatic stellate cells (HSCs) into collagen type I-producing myofibroblasts (activated HSCs). Cessation of liver injury often results in fibrosis resolution and inactivation of activated HSCs/myofibroblasts into a quiescent-like state (inactivated HSCs). We aimed to identify molecular features of phenotypes of HSCs from mice and humans. We performed studies with Lrat^Cre, Ets1-floxed, Nf1-floxed, Pparγ-floxed, Gata6-floxed, Rag2^-/-γc^-/-, and C57/Bl6 (control) mice. Some mice were given carbon tetrachloride (CCl₄) to induce liver fibrosis, with or without a peroxisome proliferator-activated receptor-γ (PPARγ) agonist. Livers from mice were analyzed by immunohistochemistry. Quiescent, activated, and inactivated HSCs were isolated from livers of Col1α1^YFP mice and analyzed by chromatin immunoprecipitation and sequencing. Human HSCs were isolated from livers denied for transplantation. We compared changes in gene expression patterns and epigenetic modifications (histone H3 lysine 4 dimethylation and histone H3 lysine 27 acetylation) in primary mouse and human HSCs. Transcription factors were knocked down with small hairpin RNAs in mouse HSCs. Motif enrichment identified E26 transcription-specific transcription factors (ETS) 1, ETS2, GATA4, GATA6, interferon regulatory factor (IRF) 1, and IRF2 transcription factors as regulators of the mouse and human HSC lineage. Small hairpin RNA-knockdown of these transcription factors resulted in increased expression of genes that promote fibrogenesis and inflammation, and loss of HSC phenotype. Disruption of Gata6 or Ets1, or Nf1 or Pparγ (which are regulated by ETS1), increased the severity of CCl₄-induced liver fibrosis in mice compared to control mice. Only mice with disruption of Gata6 or Pparγ had defects in fibrosis resolution after CCl₄ administration was stopped, associated with persistent activation of HSCs. Administration of a PPARγ agonist accelerated regression of liver fibrosis after CCl₄ administration in control mice but not in mice with disruption of Pparγ. Phenotypes of HSCs from humans and mice are regulated by transcription factors, including ETS1, ETS2, GATA4, GATA6, IRF1, and IRF2. Activated mouse and human HSCs can revert to a quiescent-like, inactivated phenotype. We found GATA6 and PPARγ to be required for inactivation of human HSCs and regression of liver fibrosis in mice.

Advanced hepatic fibrosis, driven by the activation of hepatic stellate cells (HSCs), affects millions worldwide and is the strongest predictor of mortality in nonalcoholic steatohepatitis (NASH); however, there are no approved antifibrotic therapies. To identify antifibrotic drug targets, we integrated progressive transcriptomic and morphological responses that accompany HSC activation in advanced disease using single-nucleus RNA sequencing and tissue clearing in a robust murine NASH model. In advanced fibrosis, we found that an autocrine HSC signaling circuit emerged that was composed of 68 receptor-ligand interactions conserved between murine and human NASH. These predicted interactions were supported by the parallel appearance of markedly increased direct stellate cell-cell contacts in murine NASH. As proof of principle, pharmacological inhibition of one such autocrine interaction, neurotrophic receptor tyrosine kinase 3–neurotrophin 3, inhibited human HSC activation in culture and reversed advanced murine NASH fibrosis. In summary, we uncovered a repertoire of antifibrotic drug targets underlying advanced fibrosis in vivo. The findings suggest a therapeutic paradigm in which stage-specific therapies could yield enhanced antifibrotic efficacy in patients with advanced hepatic fibrosis.

T and B cell receptors (TCRs and BCRs) play a pivotal role in the adaptive immune system by recognizing an enormous variety of external and internal antigens. Understanding these receptors is critical for exploring the process of immunoreaction and exploiting potential applications in immunotherapy and antibody drug design. Although a large number of samples have had their TCR and BCR repertoires sequenced using high-throughput sequencing in recent years, very few databases have been constructed to store these kinds of data. To resolve this issue, we developed a database. We developed a database, the Pan Immune Repertoire Database (PIRD), located in China National GeneBank (CNGBdb), to collect and store annotated TCR and BCR sequencing data, including from Homo sapiens and other species. In addition to data storage, PIRD also provides functions of data visualization and interactive online analysis. Additionally, a manually curated database of TCRs and BCRs targeting known antigens (TBAdb) was also deposited in PIRD. PIRD can be freely accessed at https://db.cngb.org/pird.

Background: Structural variants (SVs) are less common than single nucleotide polymorphisms and indels in the population, but collectively account for a significant fraction of genetic polymorphism and diseases. Base pair differences arising from SVs are on a much higher order (>100 fold) than point mutations; however, none of the current detection methods are comprehensive, and currently available methodologies are incapable of providing sufficient resolution and unambiguous information across complex regions in the human genome. To address these challenges, we applied a high-throughput, cost-effective genome mapping technology to comprehensively discover genome-wide SVs and characterize complex regions of the YH genome using long single molecules (>150 kb) in a global fashion. Results: Utilizing nanochannel-based genome mapping technology, we obtained 708 insertions/deletions and 17 inversions larger than 1 kb. Excluding the 59 SVs (54 insertions/deletions, 5 inversions) that overlap with N-base gaps in the reference assembly hg19, 666 non-gap SVs remained, and 396 of them (60%) were verified by paired-end data from whole-genome sequencing-based re-sequencing or de novo assembly sequence from fosmid data. Of the remaining 270 SVs, 260 are insertions and 213 overlap known SVs in the Database of Genomic Variants. Overall, 609 out of 666 (90%) variants were supported by experimental orthogonal methods or historical evidence in public databases. At the same time, genome mapping also provides valuable information for complex regions with haplotypes in a straightforward fashion. In addition, with long single-molecule labeling patterns, exogenous viral sequences were mapped on a whole-genome scale, and sample heterogeneity was analyzed at a new level. Conclusion: Our study highlights genome mapping technology as a comprehensive and cost-effective method for detecting structural variation and studying complex regions in the human genome, as well as deciphering viral integration into the host genome.