Brain-Enriched Coding and Long Non-coding RNA Genes Are Overrepresented in Recurrent Neurodevelopmental Disorder CNVs

A Critical Review of the Impact of Candidate Copy Number Variants on Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder (NDD) influenced by genetic, epigenetic, and environmental factors. Recent advancements in genomic analysis have shed light on numerous genes associated with ASD, highlighting the significant role of both common and rare genetic mutations, as well as copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and unique de novo variants. These genetic variations disrupt neurodevelopmental pathways, contributing to the disorder's complexity. Notably, CNVs are present in 10%-20% of individuals with autism, with 3%-7% detectable through cytogenetic methods. While the role of submicroscopic CNVs in ASD has been recently studied, their association with genomic loci and genes has not been thoroughly explored. In this review, we focus on 47 CNV regions linked to ASD, encompassing 1,632 genes, including protein-coding genes and long non-coding RNAs (lncRNAs), of which 659 show significant brain expression. Using a list of ASD-associated genes from SFARI, we detect 17 regions harboring at least one known ASD-related protein-coding gene. Of the remaining 30 regions, we identify 24 regions containing at least one protein-coding gene with brain-enriched expression and a nervous system phenotype in mouse mutants, and one lncRNA with both brain-enriched expression and upregulation in iPSC to neuron differentiation. This review not only expands our understanding of the genetic diversity associated with ASD but also underscores the potential of lncRNAs in contributing to its etiology. Additionally, the discovered CNVs will be a valuable resource for future diagnostic, therapeutic, and research endeavors aimed at prioritizing genetic variations in ASD.

Long non-coding RNAs in schizophrenia: Genetic variations, treatment markers and potential targeted signaling pathways

Schizophrenia (SZ), a complex and debilitating spectrum of psychiatric disorders, is now mainly attributed to multifactorial etiology that includes genetic and environmental factors. Long non-coding RNAs (lncRNAs) are gaining popularity as a way to better understand the comprehensive mechanisms beneath the clinical manifestation of SZ. Only in recent years has it been elucidated that mammalian genomes encode thousands of lncRNAs. Strikingly, roughly 30-40% of these lncRNAs are extensively expressed in different regions across the brain, which may be closely associated with SZ. The therapeutic and adverse effects of atypical antipsychotic drugs (AAPDs) are partially reflected by their role in the regulation of lncRNAs. This begs the question directly, do any lncRNAs exist as biomarkers for AAPDs treatment? Furthermore, we comprehend a range of mechanistic investigations that have revealed the regulatory roles for lncRNAs both involved in the brain and the periphery of SZ. More crucially, we also combine insights from a variety of signaling pathways to argue that lncRNAs probably play critical roles in SZ via their interactive downstream factors. This review provides a thorough understanding regarding dysregulation of lncRNAs, corresponding genetic alternations, as well as their potential regulatory roles in the pathology of SZ, which might help reveal useful therapeutic targets in SZ.

Pan-glioma analyses reveal species- and tumor-specific regulation of neuron-glioma synapse genes by lncRNAs

Gliomas are highly heterogeneous and aggressive. Malignant cells in gliomas can contact normal neurons through a synapse-like structure (called neuron-to-glioma synapse, NGS) to promote their proliferation, but it is unclear whether NGS gene expression and regulation show species- and tumor-specificity. This question is important in that many anti-cancer drugs are developed upon mouse models. To address this question, we conducted a pan-glioma analysis using nine scRNA-seq datasets from humans and mice. We also experimentally validated the key element of our methods and verified a key result using TCGA datasets of the same glioma types. Our analyses revealed that NGS gene expression and regulation by lncRNAs are highly species- and tumor-specific. Importantly, simian-specific lncRNAs are more involved in NGS gene regulation than lncRNAs conserved in mammals, and transgenic mouse gliomas have little in common with PDX mouse models and human gliomas in terms of NGS gene regulation. The analyses suggest that simian-specific lncRNAs are a new and rich class of potential targets for tumor-specific glioma treatment, and provide pertinent data for further experimentally and clinically exmining the targets.

Integrated analysis of copy number variation-associated lncRNAs identifies candidates contributing to the etiologies of congenital kidney anomalies

Congenital anomalies of the kidney and urinary tract (CAKUT) are disorders resulting from defects in the development of the kidneys and their outflow tract. Copy number variations (CNVs) have been identified as important genetic variations leading to CAKUT, whereas most CAKUT-associated CNVs cannot be attributed to a specific pathogenic gene. Here we construct coexpression networks involving long noncoding RNAs (lncRNAs) within these CNVs (CNV-lncRNAs) using human kidney developmental transcriptomic data. The results show that CNV-lncRNAs encompassed in recurrent CAKUT associated CNVs have highly correlated expression with CAKUT genes in the developing kidneys. The regulatory effects of two hub CNV-lncRNAs (HSALNG0134318 in 22q11.2 and HSALNG0115943 in 17q12) in the module most significantly enriched in known CAKUT genes (CAKUT_sig1, P = 1.150 × 10^-6) are validated experimentally. Our results indicate that the reduction of CNV-lncRNAs can downregulate CAKUT genes as predicted by our computational analyses. Furthermore, knockdown of HSALNG0134318 would downregulate HSALNG0115943 and affect kidney development related pathways. The results also indicate that the CAKUT_sig1 module has function significance involving multi-organ development. Overall, our findings suggest that CNV-lncRNAs play roles in regulating CAKUT genes, and the etiologies of CAKUT-associated CNVs should take account of effects on the noncoding genome.

DeepGenePrior: A deep learning model for prioritizing genes affected by copy number variants

The genetic etiology of brain disorders is highly heterogeneous, characterized by abnormalities in the development of the central nervous system that lead to diminished physical or intellectual capabilities. The process of determining which gene drives disease, known as "gene prioritization," is not entirely understood. Genome-wide searches for gene-disease associations are still underdeveloped due to reliance on previous discoveries and evidence sources with false positive or negative relations. This paper introduces DeepGenePrior, a model based on deep neural networks that prioritizes candidate genes in genetic diseases. Using the well-studied Variational AutoEncoder (VAE), we developed a score to measure the impact of genes on target diseases. Unlike other methods that use prior data to select candidate genes, based on the "guilt by association" principle and auxiliary data sources like protein networks, our study exclusively employs copy number variants (CNVs) for gene prioritization. By analyzing CNVs from 74,811 individuals with autism, schizophrenia, and developmental delay, we identified genes that best distinguish cases from controls. Our findings indicate a 12% increase in fold enrichment in brain-expressed genes compared to previous studies and a 15% increase in genes associated with mouse nervous system phenotypes. Furthermore, we identified common deletions in ZDHHC8, DGCR5, and CATG00000022283 among the top genes related to all three disorders, suggesting a common etiology among these clinically distinct conditions. DeepGenePrior is publicly available online at http://git.dml.ir/z_rahaie/DGP to address obstacles in existing gene prioritization studies identifying candidate genes.

The contribution of copy number variants to psychiatric symptoms and cognitive ability

Copy number variants (CNVs) are deletions and duplications of DNA sequence. The most frequently studied CNVs, which are described in this review, are recurrent CNVs that occur in the same locations on the genome. These CNVs have been strongly implicated in neurodevelopmental disorders, namely autism spectrum disorder (ASD), intellectual disability (ID), and developmental delay (DD), but also in schizophrenia. More recent work has also shown that CNVs increase risk for other psychiatric disorders, namely, depression, bipolar disorder, and post-traumatic stress disorder. Many of the same CNVs are implicated across all of these disorders, and these neuropsychiatric CNVs are also associated with cognitive ability in the general population, as well as with structural and functional brain alterations. Neuropsychiatric CNVs also show incomplete penetrance, such that carriers do not always develop any psychiatric disorder, and may show only mild symptoms, if any. Variable expressivity, whereby the same CNVs are associated with many different phenotypes of varied severity, also points to highly complex mechanisms underlying disease risk in CNV carriers. Comprehensive and longitudinal phenotyping studies of individual CNVs have provided initial insights into these mechanisms. However, more work is needed to estimate and predict the effect of non-recurrent, ultra-rare CNVs, which also contribute to psychiatric and cognitive outcomes. Moreover, delineating the broader phenotypic landscape of neuropsychiatric CNVs in both clinical and general population cohorts may also offer important mechanistic insights.

Copy number variation-associated lncRNAs may contribute to the etiologies of congenital heart disease

Copy number variations (CNVs) have long been recognized as pathogenic factors for congenital heart disease (CHD). Few CHD associated CNVs could be interpreted as dosage effect due to disruption of coding sequences. Emerging evidences have highlighted the regulatory roles of long noncoding RNAs (lncRNAs) in cardiac development. Whereas it remains unexplored whether lncRNAs within CNVs (CNV-lncRNAs) could contribute to the etiology of CHD associated CNVs. Here we constructed coexpression networks involving CNV-lncRNAs within CHD associated CNVs and protein coding genes using the human organ developmental transcriptomic data, and showed that CNV-lncRNAs within 10 of the non-syndromic CHD associated CNVs clustered in the most significant heart correlated module, and had highly correlated coexpression with multiple key CHD genes. HSALNG0104472 within 15q11.2 region was identified as a hub CNV-lncRNA with heart-biased expression and validated experimentally. Our results indicated that HSALNG0104472 should be a main effector responsible for cardiac defects of 15q11.2 deletion through regulating cardiomyocytes differentiation. Our findings suggested that CNV-lncRNAs could potentially contribute to the pathologies of a maximum proportion of 68.4% (13/19) of non-syndromic CHD associated CNVs. These results indicated that explaining the pathogenesis of CHD associated CNVs should take account of the noncoding regions.

A Comprehensive Investigation of Genomic Variants in Prostate Cancer Reveals 30 Putative Regulatory Variants

Prostate cancer (PC) is the most frequently diagnosed non-skin cancer in the world. Previous studies have shown that genomic alterations represent the most common mechanism for molecular alterations responsible for the development and progression of PC. This highlights the importance of identifying functional genomic variants for early detection in high-risk PC individuals. Great efforts have been made to identify common protein-coding genetic variations; however, the impact of non-coding variations, including regulatory genetic variants, is not well understood. Identification of these variants and the underlying target genes will be a key step in improving the detection and treatment of PC. To gain an understanding of the functional impact of genetic variants, and in particular, regulatory variants in PC, we developed an integrative pipeline (AGV) that uses whole genome/exome sequences, GWAS SNPs, chromosome conformation capture data, and ChIP-Seq signals to investigate the potential impact of genomic variants on the underlying target genes in PC. We identified 646 putative regulatory variants, of which 30 significantly altered the expression of at least one protein-coding gene. Our analysis of chromatin interactions data (Hi-C) revealed that the 30 putative regulatory variants could affect 131 coding and non-coding genes. Interestingly, our study identified the 131 protein-coding genes that are involved in disease-related pathways, including Reactome and MSigDB, for most of which targeted treatment options are currently available. Notably, our analysis revealed several non-coding RNAs, including RP11-136K7.2 and RAMP2-AS1, as potential enhancer elements of the protein-coding genes CDH12 and EZH1, respectively. Our results provide a comprehensive map of genomic variants in PC and reveal their potential contribution to prostate cancer progression and development.

PeakCNV: A multi-feature ranking algorithm-based tool for genome-wide copy number variation-association study

Copy Number Variation (CNV) refers to a type of structural genomic alteration in which a segment of chromosome is duplicated or deleted. To date, many CNVs have been identified as causative genetic elements for several diseases and phenotypes. However, performing a CNV-based genome-wide association study is challenging due to inconsistency in length and occurrence of CNVs across different individuals under investigation. One of the most efficient strategies to address this issue is building CNV regions (genomic regions in which CNVs are overlapping - CNVRs). However, this approach is susceptible to a high false positive rate due to overlapping and co-occurring of confounding CNVRs with true positive CNVRs. Here, we develop PeakCNV that differentiates false-positive CNVRs from true positives by calculating a new metric, independence ranking score, (IR-score) via a feature ranking approach. We compared the performance of PeakCNV with other current existing tools by carrying out two case studies one using the CNV genotype data for individuals with prostate cancer (194 cases and 2,392 healthy individuals) and the second one for individuals with neurodevelopmental disorders (19,642 cases and 6,451 healthy individuals). Crucially, our benchmarking analyses on prostate cancer cohort indicated that PeakCNV identifies a fewer risk candidate CNVRs with shorter lengths compared to other tools. Importantly, these CNVRs cover a greater proportion of case over healthy individuals compared to other tools. The accuracy of PeakCNV in identifying relevant candidate CNVRs was reproducible in the case study on neurodevelopmental disorders. Using data from the FANTOM5 expression atlas and the Clinical Genomic Database, we show that the candidate CNVRs identified by PeakCNV for neurodevelopmental disorders overlap with a greater number of genes with the brain-enriched expression, and a greater number of genes that are associated with neurological conditions compared to candidate CNVRs identified by other tools. Taken together, PeakCNV outperformed current existing CNV association study tools by identifying more biologically meaningful CNVRs relevant to the phenotype of interest. PeakCNV is publicly available for the analysis of CNV-associated diseases and is accessible from https://rdrr.io/github/mahdieh1/PeakCNV.

Somatic point mutations are enriched in non-coding RNAs with possible regulatory function in breast cancer

Non-coding RNAs (ncRNAs) form a large portion of the mammalian genome. However, their biological functions are poorly characterized in cancers. In this study, using a newly developed tool, SomaGene, we analyze de novo somatic point mutations from the International Cancer Genome Consortium (ICGC) whole-genome sequencing data of 1,855 breast cancer samples. We identify 1030 candidates of ncRNAs that are significantly and explicitly mutated in breast cancer samples. By integrating data from the ENCODE regulatory features and FANTOM5 expression atlas, we show that the candidate ncRNAs significantly enrich active chromatin histone marks (1.9 times), CTCF binding sites (2.45 times), DNase accessibility (1.76 times), HMM predicted enhancers (2.26 times) and eQTL polymorphisms (1.77 times). Importantly, we show that the 1030 ncRNAs contain a much higher level (3.64 times) of breast cancer-associated genome-wide association (GWAS) single nucleotide polymorphisms (SNPs) than genome-wide expectation. Such enrichment has not been seen with GWAS SNPs from other cancers. Using breast cell line related Hi-C data, we then show that 82% of our candidate ncRNAs (1.9 times) significantly interact with the promoter of protein-coding genes, including previously known cancer-associated genes, suggesting the critical role of candidate ncRNA genes in the activation of essential regulators of development and differentiation in breast cancer. We provide an extensive web-based resource (https://www.ihealthe.unsw.edu.au/research) to communicate our results with the research community. Our list of breast cancer-specific ncRNA genes has the potential to provide a better understanding of the underlying genetic causes of breast cancer. Lastly, the tool developed in this study can be used to analyze somatic mutations in all cancers.

The SomaGene tool is developed to identify non-coding RNAs (ncRNAs) mutated in breast cancer but can be used for other cancers. Candidate ncRNAs are shown to be enriched for regulatory features and to contain specific trait loci polymorphisms.

The Emerging Roles of Long Non-Coding RNAs in Intellectual Disability and Related Neurodevelopmental Disorders

In the human brain, long non-coding RNAs (lncRNAs) are widely expressed in an exquisitely temporally and spatially regulated manner, thus suggesting their contribution to normal brain development and their probable involvement in the molecular pathology of neurodevelopmental disorders (NDD). Bypassing the classic protein-centric conception of disease mechanisms, some studies have been conducted to identify and characterize the putative roles of non-coding sequences in the genetic pathogenesis and diagnosis of complex diseases. However, their involvement in NDD, and more specifically in intellectual disability (ID), is still poorly documented and only a few genomic alterations affecting the lncRNAs function and/or expression have been causally linked to the disease endophenotype. Considering that a significant fraction of patients still lacks a genetic or molecular explanation, we expect that a deeper investigation of the non-coding genome will unravel novel pathogenic mechanisms, opening new translational opportunities. Here, we present evidence of the possible involvement of many lncRNAs in the etiology of different forms of ID and NDD, grouping the candidate disease-genes in the most frequently affected cellular processes in which ID-risk genes were previously collected. We also illustrate new approaches for the identification and prioritization of NDD-risk lncRNAs, together with the current strategies to exploit them in diagnosis.

Enhancing DLG2 Implications in Neuropsychiatric Disorders: Analysis of a Cohort of Eight Patients with 11q14.1 Imbalances

Neurodevelopmental disorders (NDDs) are considered synaptopathies, as they are due to anomalies in neuronal connectivity during development. DLG2 is a gene involved insynaptic function; the phenotypic effect of itsalterations in NDDs has been underestimated since few cases have been thoroughly described.We report on eight patients with 11q14.1 imbalances involving DLG2, underlining its potential effects on clinical presentation and its contribution to NDD comorbidity by accurate neuropsychiatric data collection. DLG2 is a very large gene in 11q14.1, extending over 2.172 Mb, with alternative splicing that gives rise to numerous isoforms differentially expressed in brain tissues. A thorough bioinformatic analysis of the altered transcripts was conducted for each patient. The different expression profiles of the isoforms of this gene and their influence on the excitatory–inhibitory balance in crucial brain structures could contribute to the phenotypic variability related to DLG2 alterations. Further studies on patients would be helpful to enrich clinical and neurodevelopmental findings and elucidate the molecular mechanisms subtended to NDDs.

Integrative analysis of mutated genes and mutational processes reveals novel mutational biomarkers in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the leading causes of cancer-related deaths worldwide. Recent studies have observed causative mutations in susceptible genes related to colorectal cancer in 10 to 15% of the patients. This highlights the importance of identifying mutations for early detection of this cancer for more effective treatments among high risk individuals. Mutation is considered as the key point in cancer research. Many studies have performed cancer subtyping based on the type of frequently mutated genes, or the proportion of mutational processes. However, to the best of our knowledge, combination of these features has never been used together for this task. This highlights the potential to introduce better and more inclusive subtype classification approaches using wider range of related features to enable biomarker discovery and thus inform drug development for CRC.

RESULTS

In this study, we develop a new pipeline based on a novel concept called 'gene-motif', which merges mutated gene information with tri-nucleotide motif of mutated sites, for colorectal cancer subtype identification. We apply our pipeline to the International Cancer Genome Consortium (ICGC) CRC samples and identify, for the first time, 3131 gene-motif combinations that are significantly mutated in 536 ICGC colorectal cancer samples. Using these features, we identify seven CRC subtypes with distinguishable phenotypes and biomarkers, including unique cancer related signaling pathways, in which for most of them targeted treatment options are currently available. Interestingly, we also identify several genes that are mutated in multiple subtypes but with unique sequence contexts.

CONCLUSION

Our results highlight the importance of considering both the mutation type and mutated genes in identification of cancer subtypes and cancer biomarkers. The new CRC subtypes presented in this study demonstrates distinguished phenotypic properties which can be effectively used to develop new treatments. By knowing the genes and phenotypes associated with the subtypes, a personalized treatment plan can be developed that considers the specific phenotypes associated with their genomic lesion.

Current and Future Perspectives of Noncoding RNAs in Brain Function and Neuropsychiatric Disease

Noncoding RNAs (ncRNAs) represent the majority of the transcriptome and play important roles in regulating neuronal functions. ncRNAs are exceptionally diverse in both structure and function and include enhancer RNAs, long ncRNAs, and microRNAs, all of which demonstrate specific temporal and regional expression in the brain. Here, we review recent studies demonstrating that ncRNAs modulate chromatin structure, act as chaperone molecules, and contribute to synaptic remodeling and behavior. In addition, we discuss ncRNA function within the context of neuropsychiatric diseases, particularly focusing on addiction and schizophrenia, and the recent methodological developments that allow for better understanding of ncRNA function in the brain. Overall, ncRNAs represent an underrecognized molecular contributor to complex neuronal processes underlying neuropsychiatric disorders.

17q12 Recurrent Deletions and Duplications: Description of a Case Series with Neuropsychiatric Phenotype

Syndromic neurodevelopmental disorders are usually investigated through genetics technologies, within which array comparative genomic hybridization (Array-CGH) is still considered the first-tier clinical diagnostic test. Among recurrent syndromic imbalances, 17q12 deletions and duplications are characterized by neurodevelopmental disorders associated with visceral developmental disorders, although expressive variability is common. Here we describe a case series of 12 patients with 17q12 chromosomal imbalances, in order to expand the phenotypic characterization of these recurrent syndromes whose diagnosis is often underestimated, especially if only mild traits are present. Gene content and genotype-phenotype correlations have been discussed, with special regard to neuropsychiatric features, whose impact often requires etiologic analysis.

Neurodevelopmental Disorders in Patients With Complex Phenotypes and Potential Complex Genetic Basis Involving Non-Coding Genes, and Double CNVs

Neurodevelopmental disorders (NDDs) are a heterogeneous class of brain diseases, with a complex genetic basis estimated to account for up to 50% of cases. Nevertheless, genetic diagnostic yield is about 20%. Array-comparative genomic hybridization (array-CGH) is an established first-level diagnostic test able to detect pathogenic copy number variants (CNVs), however, most identified variants remain of uncertain significance (VUS). Failure of interpretation of VUSs may depend on various factors, including complexity of clinical phenotypes and inconsistency of genotype-phenotype correlations. Indeed, although most NDD-associated CNVs are de novo, transmission from unaffected parents to affected children of CNVs with high risk for NDDs has been observed. Moreover, variability of genetic components overlapped by CNVs, such as long non-coding genes, genomic regions with long-range effects, and additive effects of multiple CNVs can make CNV interpretation challenging. We report on 12 patients with complex phenotypes possibly explained by complex genetic mechanisms, including involvement of antisense genes and boundaries of topologically associating domains. Eight among the 12 patients carried two CNVs, either de novo or inherited, respectively, by each of their healthy parents, that could additively contribute to the patients’ phenotype. CNVs overlapped either known NDD-associated or novel candidate genes (PTPRD, BUD13, GLRA3, MIR4465, ABHD4, and WSCD2). Bioinformatic enrichment analyses showed that genes overlapped by the co-occurring CNVs have synergistic roles in biological processes fundamental in neurodevelopment. Double CNVs could concur in producing deleterious effects, according to a two-hit model, thus explaining the patients’ phenotypes and the incomplete penetrance, and variable expressivity, associated with the single variants. Overall, our findings could contribute to the knowledge on clinical and genetic diagnosis of complex forms of NDD.

Quantitative neurogenetics: applications in understanding disease

Neurodevelopmental and neurodegenerative disorders (NNDs) are a group of conditions with a broad range of core and co-morbidities, associated with dysfunction of the central nervous system. Improvements in high throughput sequencing have led to the detection of putative risk genetic loci for NNDs, however, quantitative neurogenetic approaches need to be further developed in order to establish causality and underlying molecular genetic mechanisms of pathogenesis. Here, we discuss an approach for prioritizing the contribution of genetic risk loci to complex-NND pathogenesis by estimating the possible impacts of these loci on gene regulation. Furthermore, we highlight the use of a tissue-specificity gene expression index and the application of artificial intelligence (AI) to improve the interpretation of the role of genetic risk elements in NND pathogenesis. Given that NND symptoms are associated with brain dysfunction, risk loci with direct, causative actions would comprise genes with essential functions in neural cells that are highly expressed in the brain. Indeed, NND risk genes implicated in brain dysfunction are disproportionately enriched in the brain compared with other tissues, which we refer to as brain-specific expressed genes. In addition, the tissue-specificity gene expression index can be used as a handle to identify non-brain contexts that are involved in NND pathogenesis. Lastly, we discuss how using an AI approach provides the opportunity to integrate the biological impacts of risk loci to identify those putative combinations of causative relationships through which genetic factors contribute to NND pathogenesis.

Non-Syndromic Intellectual Disability and Its Pathways: A Long Noncoding RNA Perspective

Non-syndromic intellectual disability (NS-ID or idiopathic) is a complex neurodevelopmental disorder that represents a global health issue. Although many efforts have been made to characterize it and distinguish it from syndromic intellectual disability (S-ID), the highly heterogeneous aspect of this disorder makes it difficult to understand its etiology. Long noncoding RNAs (lncRNAs) comprise a large group of transcripts that can act through various mechanisms and be involved in important neurodevelopmental processes. In this sense, comprehending the roles they play in this intricate context is a valuable way of getting new insights about how NS-ID can arise and develop. In this review, we attempt to bring together knowledge available in the literature about lncRNAs involved with molecular and cellular pathways already described in intellectual disability and neural function, to better understand their relevance in NS-ID and the regulatory complexity of this disorder.

A systematic review of long non-coding RNAs with a potential role in breast cancer

The human transcriptome contains many non-coding RNAs (ncRNAs), which play important roles in gene regulation. Long noncoding RNAs (lncRNAs) are an important class of ncRNAs with lengths between 200 and 200,000 bases. Unlike mRNA, lncRNA lacks protein-coding features, specifically, open-reading frames, and start and stop codons. LncRNAs have been reported to play a role in the pathogenesis and progression of many cancers, including breast cancer (BC), acting as tumor suppressors or oncogenes. In this review, we systematically mined the literature to identify 65 BC-related lncRNAs. We then perform an integrative bioinformatics analysis to identify 14 lncRNAs with a potential regulatory role in BC. The biological function of these 14 lncRNAs, their regulatory mechanisms, and roles in the initiation and progression of BC are discussed in this review. Additionally, we elaborate on the current and future applications of lncRNAs as diagnostic and/or therapeutic biomarkers in BC.