Temporal transcriptomics identifies transregulation by multiple lncRNAs in human iPSCs

Some long noncoding RNAs (lncRNAs) are known to regulate gene expression. However, the underlying temporal dynamics of lncRNAs influencing gene and epigenetic regulation, and mechanisms of lncRNA regulation in trans, are less understood. To investigate this, we genetically engineer 17 doxycycline-inducible lncRNA transgenes for ectopic expression at the H11 safe harbor locus in human pluripotent stem cells (hiPSCs), and we generate high-density temporal RNA-seq and ATAC-seq profiles. Most lncRNA transgenes are induced at 2 hours and maintain expression through the 96-hour time course. Surprisingly, when we seek to identify gene expression changes due to the lncRNAs, we find that the global transcriptional landscape is dominated by a strong systemic response triggered by doxycycline exposure. We rigorously define this cohort of genes as a Doxycycline-Responsive Gene Signature (DRGS). The DRGS is also present in at least 28 public datasets from dox-inducible transgene studies involving diverse cell types. Next, we determine which lncRNAs exhibit trans-regulatory events. We identify DANCR, FENDRR, LINC00667, LINC00847, LNCPRESS1, and LNCPNKY as lncRNAs that regulate specific transcript expression in trans The downstream target genes encode 53 mRNAs and 10 lncRNAs. None of the target lncRNAs alter gene expression proximal to their own loci (i.e., triggering secondary cis effects). Surprisingly, the target genes of LINC00847 (transcribed from Chromosome 22) are substantially enriched on Chromosome 19, with a preponderance of target genes encoding RNA metabolism and RNA splicing factors. Collectively, our study provides a resource to discern artifacts in the doxycycline-inducible system and identifies temporally regulated targets of 6 lncRNAs for future mechanistic studies.