ABSTRACT
Long non-coding RNAs (lncRNAs) can act as scaffolds that promote the interaction of proteins, RNA, and DNA. There is increasing evidence of sequence-specific interactions of lncRNAs with DNA via triple-helix (triplex) formation. This process allows lncRNAs to recruit protein complexes to specific genomic regions and regulate gene expression. Here we propose a computational method called Triplex Domain Finder (TDF) to detect triplexes and characterize DNA-binding domains and DNA targets statistically. Case studies showed that this approach can detect the known domains of lncRNAs Fendrr, HOTAIR and MEG3. Moreover, we validated a novel DNA-binding domain in MEG3 by a genome-wide sequencing method. We used TDF to perform a systematic analysis of the triplex-forming potential of lncRNAs relevant to human cardiac differentiation. We demonstrated that the lncRNA with the highest triplex-forming potential, GATA6-AS, forms triple helices in the promoter of genes relevant to cardiac development. Moreover, down-regulation of GATA6-AS impairs GATA6 expression and cardiac development. These data indicate the unique ability of our computational tool to identify novel triplex-forming lncRNAs and their target genes.
Subject(s)
Computational Biology/methods , DNA/metabolism , RNA, Long Noncoding/chemistry , RNA, Long Noncoding/metabolism , Algorithms , Base Sequence , Binding Sites/genetics , DNA/chemistry , Gene Expression , Humans , Nucleic Acid Conformation , Protein Binding , Transcription Factors/metabolismABSTRACT
Human protein-coding genes are often accompanied by divergently transcribed non-coding RNAs whose functions, especially in cell fate decisions, are poorly understood. Using an hESC-based cardiac differentiation model, we define a class of divergent lncRNAs, termed yin yang lncRNAs (yylncRNAs), that mirror the cell-type-specific expression pattern of their protein-coding counterparts. yylncRNAs are preferentially encoded from the genomic loci of key developmental cell fate regulators. Most yylncRNAs are spliced polyadenylated transcripts showing comparable expression patterns in vivo in mouse and in human embryos. Signifying their developmental function, the key mesoderm specifier BRACHYURY (T) is accompanied by yylncT, which localizes to the active T locus during mesoderm commitment. yylncT binds the de novo DNA methyltransferase DNMT3B, and its transcript is required for activation of the T locus, with yylncT depletion specifically abolishing mesodermal commitment. Collectively, we report a lncRNA-mediated regulatory layer safeguarding embryonic cell fate transitions.
