The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called k-mers. We found that lncRNAs of related function often had similar k-mer profiles despite lacking linear homology, and that k-mer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. K-mer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.

Patients with late-onset Alzheimer's disease (LOAD) frequently manifest comorbid neuropsychiatric symptoms with depression and anxiety being most frequent, and individuals with major depressive disorder (MDD) have an increased prevalence of LOAD. This suggests shared etiologies and intersecting pathways between LOAD and MDD. We performed pleiotropy analyses using LOAD and MDD GWAS data sets from the International Genomics of Alzheimer's Project (IGAP) and the Psychiatric Genomics Consortium (PGC), respectively. We found a moderate enrichment for SNPs associated with LOAD across increasingly stringent levels of significance with the MDD GWAS association (LOAD| MDD), of maximum four and eightfolds, including and excluding the APOE-region, respectively. Association analysis excluding the APOE-region identified numerous SNPs corresponding to 40 genes, 9 of which are known LOAD-risk loci primarily in chromosome 11 regions that contain the SPI1 gene and MS4A genes cluster, and others were novel pleiotropic risk-loci for LOAD conditional with MDD. The most significant associated SNPs on chromosome 11 overlapped with eQTLs found in whole-blood and monocytes, suggesting functional roles in gene regulation. The reverse conditional association analysis (MDD|LOAD) showed a moderate level,~sevenfold, of polygenic overlap, however, no SNP showed significant association. Pathway analyses replicated previously reported LOAD biological pathways related to immune response and regulation of endocytosis. In conclusion, we provide insights into the overlapping genetic signatures underpinning the common phenotypic manifestations and inter-relationship between LOAD and MDD. This knowledge is crucial to the development of actionable targets for novel therapies to treat depression preceding dementia, in an effort to delay or ultimately prevent the onset of LOAD.

The marsupial inactive X chromosome expresses a long noncoding RNA (lncRNA) called Rsx that has been proposed to be the functional analog of eutherian Xist. Despite the possibility that Xist and Rsx encode related functions, the two lncRNAs harbor no linear sequence similarity. However, both lncRNAs harbor domains of tandemly repeated sequence. In Xist, these repeat domains are known to be critical for function. Using k-mer based comparison, we show that the repeat domains of Xist and Rsx unexpectedly partition into two major clusters that each harbor substantial levels of nonlinear sequence similarity. Xist Repeats B, C, and D were most similar to each other and to Rsx Repeat 1, whereas Xist Repeats A and E were most similar to each other and to Rsx Repeats 2, 3, and 4. Similarities at the level of k-mers corresponded to domain-specific enrichment of protein-binding motifs. Within individual domains, protein-binding motifs were often enriched to extreme levels. Our data support the hypothesis that Xist and Rsx encode similar functions through different spatial arrangements of functionally analogous protein-binding domains. We propose that the two clusters of repeat domains in Xist and Rsx function in part to cooperatively recruit PRC1 and PRC2 to chromatin. The physical manner in which these domains engage with protein cofactors may be just as critical to the function of the domains as the protein cofactors themselves. The general approaches we outline in this report should prove useful in the study of any set of RNAs.