Interplay Between GH-regulated, Sex-biased Liver Transcriptome and Hepatic Zonation Revealed by Single-Nucleus RNA Sequencing

The zonation of liver metabolic processes is well-characterized; however, little is known about the cell type-specificity and zonation of sexually dimorphic gene expression or its growth hormone (GH)-dependent transcriptional regulators. We address these issues using single-nucleus RNA-sequencing of 32 000 nuclei representing 9 major liver cell types. Nuclei were extracted from livers from adult male and female mice; from males infused with GH continuously, mimicking the female plasma GH pattern; and from mice exposed to TCPOBOP, a xenobiotic agonist ligand of the nuclear receptor CAR that perturbs sex-biased gene expression. Analysis of these rich transcriptomic datasets revealed the following: 1) expression of sex-biased genes and their GH-dependent transcriptional regulators is primarily restricted to hepatocytes and is not a feature of liver nonparenchymal cells; 2) many sex-biased transcripts show sex-dependent zonation within the liver lobule; 3) gene expression is substantially feminized both in periportal and pericentral hepatocytes when male mice are infused with GH continuously; 4) sequencing nuclei increases the sensitivity for detecting thousands of nuclear-enriched long-noncoding RNAs (lncRNAs) and enables determination of their liver cell type-specificity, sex-bias and hepatocyte zonation profiles; 5) the periportal to pericentral hepatocyte cell ratio is significantly higher in male than female liver; and 6) TCPOBOP exposure disrupts both sex-specific gene expression and hepatocyte zonation within the liver lobule. These findings highlight the complex interconnections between hepatic sexual dimorphism and zonation at the single-cell level and reveal how endogenous hormones and foreign chemical exposure can alter these interactions across the liver lobule with large effects both on protein-coding genes and lncRNAs.

Global analysis of expression, maturation and subcellular localization of mouse liver transcriptome identifies novel sex-biased and TCPOBOP-responsive long non-coding RNAs

BACKGROUND

While nuclear transcription and RNA processing and localization are well established for protein coding genes (PCGs), these processes are poorly understood for long non-coding (lnc)RNAs. Here, we characterize global patterns of transcript expression, maturation and localization for mouse liver RNA, including more than 15,000 lncRNAs. PolyA-selected liver RNA was isolated and sequenced from four subcellular fractions (chromatin, nucleoplasm, total nucleus, and cytoplasm), and from the chromatin-bound fraction without polyA selection.

RESULTS

Transcript processing, determined from normalized intronic to exonic sequence read density ratios, progressively increased for PCG transcripts in going from the chromatin-bound fraction to the nucleoplasm and then on to the cytoplasm. Transcript maturation was similar for lncRNAs in the chromatin fraction, but was significantly lower in the nucleoplasm and cytoplasm. LncRNA transcripts were 11-fold more likely to be significantly enriched in the nucleus than cytoplasm, and 100-fold more likely to be significantly chromatin-bound than nucleoplasmic. Sequencing chromatin-bound RNA greatly increased the sensitivity for detecting lowly expressed lncRNAs and enabled us to discover and localize hundreds of novel regulated liver lncRNAs, including lncRNAs showing sex-biased expression or responsiveness to TCPOBOP a xenobiotic agonist ligand of constitutive androstane receptor (Nr1i3).

CONCLUSIONS

Integration of our findings with prior studies and lncRNA annotations identified candidate regulatory lncRNAs for a variety of hepatic functions based on gene co-localization within topologically associating domains or transcription divergent or antisense to PCGs associated with pathways linked to hepatic physiology and disease.