HOXs and lincRNAs: Two sides of the same coin

The Potential Links between lncRNAs and Drug Tolerance in Lung Adenocarcinoma

Lung cancer patients treated with targeted therapies frequently respond well but invariably relapse due to the development of drug resistance. Drug resistance is in part mediated by a subset of cancer cells termed "drug-tolerant persisters" (DTPs), which enter a dormant, slow-cycling state that enables them to survive drug exposure. DTPs also exhibit stem cell-like characteristics, broad epigenetic reprogramming, altered metabolism, and a mutagenic phenotype mediated by adaptive mutability. While several studies have characterised the transcriptional changes that lead to the altered phenotypes exhibited in DTPs, these studies have focused predominantly on protein coding changes. As long non-coding RNAs (lncRNAs) are also implicated in the phenotypes altered in DTPs, it is likely that they play a role in the biology of drug tolerance. In this review, we outline how lncRNAs may contribute to the key characteristics of DTPs, their potential roles in tolerance to targeted therapies, and the emergence of genetic resistance in lung adenocarcinoma.

Deciphering the emerging landscape of HOX genes in cardiovascular biology, atherosclerosis and beyond (Review)

Atherosclerosis, a dominant driving force underlying multiple cardiovascular events, is an intertwined and chronic inflammatory disease characterized by lipid deposition in the arterial wall, which leads to diverse cardiovascular problems. Despite unprecedented advances in understanding the pathogenesis of atherosclerosis and the substantial decline in cardiovascular mortality, atherosclerotic cardiovascular disease remains a global public health issue. Understanding the molecular landscape of atherosclerosis is imperative in the field of molecular cardiology. Recently, compelling evidence has shown that an important family of homeobox (HOX) genes endows causality in orchestrating the interplay between various cardiovascular biological processes and atherosclerosis. Despite seemingly scratching the surface, such insight into the realization of biology promises to yield extraordinary breakthroughs in ameliorating atherosclerosis. Primarily recapitulated herein are the contributions of HOX in atherosclerosis, including diverse cardiovascular biology, knowledge gaps, remaining challenges and future directions. A snapshot of other cardiovascular biological processes was also provided, including cardiac/vascular development, cardiomyocyte pyroptosis/apoptosis, cardiac fibroblast proliferation and cardiac hypertrophy, which are responsible for cardiovascular disorders. Further in‑depth investigation of HOX promises to provide a potential yet challenging landscape, albeit largely undetermined to date, for partially pinpointing the molecular mechanisms of atherosclerosis. A plethora of new targeted therapies may ultimately emerge against atherosclerosis, which is rapidly underway. However, translational undertakings are crucially important but increasingly challenging and remain an ongoing and monumental conundrum in the field.

Ablation of Long Noncoding RNA Hoxb3os Exacerbates Cystogenesis in Mouse Polycystic Kidney Disease

SIGNIFICANCE STATEMENT

Long noncoding RNAs (lncRNAs) are a class of nonprotein coding RNAs with pivotal functions in development and disease. They have emerged as an exciting new drug target category for many common conditions. However, the role of lncRNAs in autosomal dominant polycystic kidney disease (ADPKD) has been understudied. This study provides evidence implicating a lncRNA in the pathogenesis of ADPKD. We report that Hoxb3os is downregulated in ADPKD and regulates mammalian target of rapamycin (mTOR)/Akt pathway in the in vivo mouse kidney. Ablating the expression of Hoxb3os in mouse polycystic kidney disease (PKD) activated mTOR complex 2 (mTORC2) signaling and exacerbated the cystic phenotype. The results from our study provide genetic proof of concept for future studies that focus on targeting lncRNAs as a treatment option in PKD.

BACKGROUND

ADPKD is a monogenic disorder characterized by the formation of kidney cysts and is primarily caused by mutations in two genes, PKD1 and PKD2 .

METHODS

In this study, we investigated the role of lncRNA Hoxb3os in ADPKD by ablating its expression in the mouse.

RESULTS

Hoxb3os -null mice were viable and had grossly normal kidney morphology but displayed activation of mTOR/Akt signaling and subsequent increase in kidney cell proliferation. To determine the role of Hoxb3os in cystogenesis, we crossed the Hoxb3os -null mouse to two orthologous Pkd1 mouse models: Pkhd1/Cre; Pkd1F/F (rapid cyst progression) and Pkd1RC/RC (slow cyst progression). Ablation of Hoxb3os exacerbated cyst growth in both models. To gain insight into the mechanism whereby Hoxb3os inhibition promotes cystogenesis, we performed western blot analysis of mTOR/Akt pathway between Pkd1 single-knockout and Pkd1 - Hoxb3os double-knockout (DKO) mice. Compared with single-knockout, DKO mice presented with enhanced levels of total and phosphorylated Rictor. This was accompanied by increased phosphorylation of Akt at Ser 473 , a known mTORC2 effector site. Physiologically, kidneys from DKO mice displayed between 50% and 60% increase in cell proliferation and cyst number.

CONCLUSIONS

The results from this study indicate that ablation of Hoxb3os in mouse PKD exacerbates cystogenesis and dysregulates mTORC2.

Spotlight on HOX cluster‑embedded antisense lncRNAs in cardiovascular diseases (Review)

Atherosclerosis is a complex and chronic inflammatory disease driven by multiple pathophysiological processes that are responsible for diverse cardiovascular events. Atherosclerotic cardiovascular disease, despite substantial triumphs in primary and secondary prevention, remains a dominant epidemic that impairs human health. Therefore, deciphering the pathogenesis of atherosclerosis will provide a real‑world translational understanding. Homeobox cluster‑embedded antisense long non‑coding RNAs (HOX‑lncRNAs), a nascent class of lncRNA molecules with versatile roles in cancer, can also orchestrate various cell functions in cardiovascular disorders and have thus captured the attention of many researchers. Subsequently, numerous studies have demonstrated the role of HOX‑lncRNAs as potential modulators of atherosclerosis. Nevertheless, given that the understanding of HOX‑lncRNAs in atherosclerosis is only just emerging, ongoing research must be initiated to thoroughly pinpoint such causal roles. The present review aimed to highlight the important contributions of HOX‑lncRNAs to atherosclerosis and other pivotal biological processes related to cardiovascular disease. The review concludes with a discussion of the limitations, outlook, challenges and possible solutions associated with HOX‑lncRNAs in atherosclerosis. Looking forward, this may lead to extraordinary breakthroughs in revealing the molecular underpinnings of HOX‑lncRNAs and may offer a promising yet challenging landscape for robust therapeutic strategies for atherosclerosis and/or associated cardiovascular disorders.

The role of the 5' HoxA genes in the development of the hindgut, vent, and a novel sphincter in a derived teleost (bluebanded goby, Lythrypnus dalli)

Unique expression patterns of the 5' HoxA genes are associated with the evolution and development of novel features including claspers in cartilaginous fishes, modified pectoral fins in batoids, and the yolk sac extension in Cypriniformes. Here, we demonstrate a role for HoxA11a and HoxA13a in demarcating the hindgut in fishes of the family Gobiidae, including a novel sphincter called the intestinal rectal sphincter (IRS). Disruption of 5' HoxA expression, via manipulation of retinoic acid signaling, results in failure of the IRS and/or vent to develop. Furthermore, exposure to HoxA disruptors alters 5' HoxA expression, in association with developmental phenotypes, demonstrating a functional link between 5' HoxA expression and development of a novel feature in the bluebanded goby, Lythrypnus dalli.

DNA methylation differences in noncoding regions in ER negative breast tumors between Black and White women

Introduction

Incidence of estrogen receptor (ER)-negative breast cancer, an aggressive tumor subtype associated with worse prognosis, is higher among African American/Black women than other US racial and ethnic groups. The reasons for this disparity remain poorly understood but may be partially explained by differences in the epigenetic landscape.

Methods

We previously conducted genome-wide DNA methylation profiling of ER- breast tumors from Black and White women and identified a large number of differentially methylated loci (DML) by race. Our initial analysis focused on DML mapping to protein-coding genes. In this study, motivated by increasing appreciation for the biological importance of the non-protein coding genome, we focused on 96 DMLs mapping to intergenic and noncoding RNA regions, using paired Illumina Infinium Human Methylation 450K array and RNA-seq data to assess the relationship between CpG methylation and RNA expression of genes located up to 1Mb away from the CpG site.

Results

Twenty-three (23) DMLs were significantly correlated with the expression of 36 genes (FDR<0.05), with some DMLs associated with the expression of single gene and others associated with more than one gene. One DML (cg20401567), hypermethylated in ER- tumors from Black versus White women, mapped to a putative enhancer/super-enhancer element located 1.3 Kb downstream of HOXB2. Increased methylation at this CpG correlated with decreased expression of HOXB2 (Rho=-0.74, FDR<0.001) and other HOXB/HOXB-AS genes. Analysis of an independent set of 207 ER- breast cancers from TCGA similarly confirmed hypermethylation at cg20401567 and reduced HOXB2 expression in tumors from Black versus White women (Rho=-0.75, FDR<0.001).

Discussion

Our findings indicate that epigenetic differences in ER- tumors between Black and White women are linked to altered gene expression and may hold functional significance in breast cancer pathogenesis.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

A single-cell transcriptome atlas profiles early organogenesis in human embryos

The early window of human embryogenesis is largely a black box for developmental biologists. Here we probed the cellular diversity of 4-6 week human embryos when essentially all organs are just laid out. On the basis of over 180,000 single-cell transcriptomes, we generated a comprehensive atlas of 313 clusters in 18 developmental systems, which were annotated with a collection of ontology and markers from 157 publications. Together with spatial transcriptome on embryonic sections, we characterized the molecule and spatial architecture of previously unappreciated cell types. Combined with data from other vertebrates, the rich information shed light on spatial patterning of axes, systemic temporal regulation of developmental progression and potential human-specific regulation. Our study provides a compendium of early progenitor cells of human organs, which can serve as the root of lineage analysis in organogenesis.

HOX cluster-embedded lncRNAs and epithelial-mesenchymal transition in cancer: Molecular mechanisms and therapeutic opportunities

Although there has been substantial improvement in the treatment modalities, cancer remains the major cause of fatality worldwide. Metastasis, recurrence, and resistance to oncological therapies are the leading causes of cancer mortality. Epithelial-mesenchymal transition (EMT) is a complex biological process that allows cancer cells to undergo morphological transformation into a mesenchymal phenotype to acquire invasive potential. It encompasses reversible and dynamic ontogenesis by neoplastic cells during metastatic dissemination. Hence, understanding the molecular landscape of EMT is imperative to identify a reliable clinical biomarker to combat metastatic spread. Accumulating evidence reveals the role of HOX (homeobox) cluster-embedded long non-coding RNAs (lncRNAs) in EMT and cancer metastasis. They play a crucial role in the induction of EMT, modulating diverse biological targets. The present review emphasizes the involvement of HOX cluster-embedded lncRNAs in EMT as a molecular sponge, chromatin remodeler, signaling regulator, and immune system modulator. Furthermore, the molecular mechanisms behind therapy resistance and the potential use of novel drugs targeting HOX cluster-embedded lncRNAs in the clinical management of distant metastasis will be discussed.

Small hairpin inhibitory RNA delivery in the metanephric organ culture identifies long noncoding RNA Pvt1 as a modulator of cyst growth

Autosomal dominant polycystic kidney disease (ADPKD) is a monogenic disorder characterized by the formation of kidney cysts that originate from the epithelial tubules of the nephron and primarily results from mutations in polycystin-1 (PKD1) and polycystin-2 (PKD2). The metanephric organ culture (MOC) is an ex vivo system in which explanted embryonic kidneys undergo tubular differentiation and kidney development. MOC has been previously used to study polycystic kidney disease as treatment with 8-bromo-cAMP induces the formation of kidney cysts. However, the inefficiency of manipulating gene expression in MOC has limited its utility for identifying genes and pathways that are involved in cystogenesis. Here, we used a lentivirus and three serotypes of self-complementary adeno-associated viral (scAAV) plasmids that express green fluorescent protein and found that scAAV serotype D/J transduces the epithelial compartment of MOC at an efficiency of 68%. We used scAAV/DJ to deliver shRNA to knockdown Pvt1, a long noncoding RNA, which was upregulated in kidneys from Pkd1 and Pkd2 mutant mice and humans with ADPKD. shRNA delivery by scAAV/DJ downregulated expression of Pvt1 by 45% and reduced the cyst index by 53% in wild-type MOCs and 32% in Pkd1-null MOCs. Knockdown of Pvt1 decreased the level of c-MYC protein by 60% without affecting Myc mRNA, indicating that Pvt1 regulation of c-MYC was posttranscriptional. These results identify Pvt1 as a long noncoding RNA that modulates cyst progression in MOC.NEW & NOTEWORTHY This study identified scAAV/DJ as effective in transducing epithelial cells of the metanephric organ culture (MOC). We used scAAV/DJ shRNA to knockdown Pvt1 in cystic MOCs derived from Pkd1-null embryos. Downregulation of Pvt1 reduced cyst growth and decreased levels of c-MYC protein. These data suggest that suppression of Pvt1 activity in autosomal dominant polycystic kidney disease might reduce cyst growth.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

A coordinated function of lncRNA HOTTIP and miRNA-196b underpinning leukemogenesis by targeting FAS signaling

MicroRNAs (miRNAs) may modulate more than 60% of human coding genes and act as negative regulators, whereas long noncoding RNAs (lncRNAs) regulate gene expression on multiple levels by interacting with chromatin, functional proteins, and RNAs such as mRNAs and microRNAs. However, the crosstalk between HOTTIP lncRNA and miRNAs in leukemogenesis remains elusive. Using combined integrated analyses of global miRNA expression profiling and state-of-the-art genomic analyses of chromatin such as ChIRP-seq (HOTTIP binding in genomewide), ChIP-seq, and ATAC-seq, we found that some miRNA genes are directly controlled by HOTTIP. Specifically, the HOX cluster miRNAs (miR-196a, miR-196b, miR-10a, and miR-10b), located cis and trans, were most dramatically regulated and significantly decreased in HOTTIP-/- AML cells. HOTTIP bound to the miR-196b promoter and HOTTIP deletion reduced chromatin accessibility and enrichment of active histone modifications at HOX cluster-associated miRNAs in AML cells, whereas reactivation of HOTTIP restored miR gene expression and chromatin accessibility in the CTCF-boundary-attenuated AML cells. Inactivation of HOTTIP or miR-196b promotes apoptosis by altering the chromatin signature at the FAS promoter and increasing FAS expression. Transplantation of miR-196b knockdown MOLM13 cells in NSG mice increased overall survival of mice compared to wild-type cells transplanted into mice. Thus, HOTTIP remodels the chromatin architecture around miRNAs to promote their transcription and consequently represses tumor suppressors and promotes leukemogenesis.

The expression of Drosophila melanogaster Hox gene Ultrabithorax is not overtly regulated by the intronic long non-coding RNA lncRNA: PS4 in a wild type genetic background

Long noncoding RNAs (lncRNAs) have been implicated in a variety of processes in development, differentiation, and disease. In Drosophila melanogaster, the bithorax Hox cluster contains three Hox genes [Ultrabithorax (Ubx), abdominal-A, and Abdominal-B], along with a number of lncRNAs, most with unknown functions. Here, we investigated the function of a lncRNA, lncRNA:PS4 that originates in the second intron of Ubx and is transcribed in the antisense orientation to Ubx. The expression pattern of lncRNA:PS4 is complementary to Ubx in the thoracic primordia, and the lncRNA:PS4 coding region overlaps the location of the large insertion that causes the dominant homeotic mutation Contrabithorax-1 (Ubx^Cbx-1), which partially transforms Drosophila wings into halteres via ectopic activation of Ubx. This led us to investigate the potential role of this lncRNA in regulation of Ubx expression. The Ubx^Cbx-1 mutation dramatically changes the pattern of lncRNA:PS4, eliminating the expression of most lncRNA:PS4 sequences from parasegment 4 (where Ubx protein is normally absent) and ectopically activating lncRNA:PS4 at high levels in the abdomen (where Ubx is normally expressed). These changes, however, did not lead to changes in the Ubx embryonic transcription pattern. Targeted deletion of the two promoters of lncRNA:PS4 did not result in the change of Ubx expression in the embryos. In the genetic background of a Ubx^Cbx-1 mutation, the lncRNA:PS4 mutation does slightly enhance the ectopic activation of Ubx protein expression in wing discs and also slightly enhances the wing phenotype seen in Ubx^Cbx-1 heterozygotes.

Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

The long non-coding RNAs (lncRNAs) have been known to play important roles in several biological processes as well as in reproduction. This study aimed to identify lncRNA in ovary female banana shrimp, Fenneropenaeus merguiensis, and investigate the potential role of lncPV13 in the vitellogenesis. After the in silico identification of the ovarian transcriptome, a total of 24,733 putative lncRNAs were obtained, and only 147 putative lncRNAs were significantly differentially expressed among the ovarian development stages. To validate the in silico identification of lncRNAs, the 16 lncRNAs with the highest differential expression in the transcriptome analysis were evaluated by RT-qPCR. The 6 lncRNAs showed higher expression levels in the mature stage than in the previtellogenic stage and were found in several tissues such as in eyestalks, brains, thoracic ganglia, gills, and muscle. Furthermore, most candidate lncRNAs were amplifiable in Litopenaeus vannamei's and Penaeus monodon's DNA but not in Macrobrachium rosenbergii's DNA, suggesting some lncRNAs are expressed in a species-specific manner among penaeid shrimp. In this study, the lncPV13 was investigated for its vitellogenin regulating function by RNA interference. The result indicates that the lncPV13 expression was suppressed in the ovary on day 7 after the injection of double-stranded RNA specific to lncPV13 (dslncPV13), while vitellogenin (Vg) expression was significantly decreased. In contrast, the gonad inhibiting hormone (GIH) expression was significantly increased in the lncPV13 knockdown shrimp. However, the oocyte proliferation was not significantly different between control and lncPV13 knockdown shrimp. This suggests that lncPV13 regulate Vg synthesis through GIH inhibition. Finally, our findings provide lncRNA information and potential lncRNAs involved in the vitellogenesis of female banana shrimp.

HOXB7 Overexpression Leads Triple-Negative Breast Cancer Cells to a Less Aggressive Phenotype

HOX genes appear to play a role in breast cancer progression in a molecular subtype-dependent way. The altered expression of HOXB7, for example, was reported to promote breast cancer progression in specific subtypes. Here we induced HOXB7 overexpression in MDA-MB-231 cells, a cellular model of the Triple-Negative breast cancer molecular subtype, and evaluated the phenotypic changes in cell viability, morphogenesis, migration, invasion, and colony formation. During the phenotypic characterization of the HOXB7-overexpressing cells, we consistently found less aggressive behavior represented by lower cell viability, inhibition of cell migration, invasion, and attachment-independent colony formation capacities added to the more compact and organized spheroid growth in 3D cultures. We then evaluated the expression of putative downstream targets and their direct binding to HOXB7 comparing ChIP-qPCR data generated from HOXB7-overexpressing cells and controls. In the manipulated cells, we found enriched biding of HOXB7 to CTNNB1, EGFR, FGF2, CDH1, DNMT3B, TGFB2, and COMMD7. Taken together, these results highlight the plasticity of the HOXB7 function in breast cancer, according to the cellular genetic background and expression levels, and provide evidence that in Triple-Negative breast cancer cells, HOXB7 overexpression has the potential to promote less aggressive phenotypes.

Temporal mechanisms of myogenic specification in human induced pluripotent stem cells

Understanding the mechanisms of myogenesis in human induced pluripotent stem cells (hiPSCs) is a prerequisite to achieving patient-specific therapy for diseases of skeletal muscle. hiPSCs of different origin show distinctive kinetics and ability to differentiate into myocytes. To address the unique cellular and temporal context of hiPSC differentiation, we perform a longitudinal comparison of the transcriptomic profiles of three hiPSC lines that display differential myogenic specification, one robust and two blunted. We detail temporal differences in mechanisms that lead to robust myogenic specification. We show gene expression signatures of putative cell subpopulations and extracellular matrix components that may support myogenesis. Furthermore, we show that targeted knockdown of ZIC3 at the outset of differentiation leads to improved myogenic specification in blunted hiPSC lines. Our study suggests that β-catenin transcriptional cofactors mediate cross-talk between multiple cellular processes and exogenous cues to facilitate specification of hiPSCs to mesoderm lineage, leading to robust myogenesis.

Comprehensive analysis of the homeobox family genes in breast cancer demonstrates their similar roles in cancer and development

BACKGROUND

The homeobox (HOX) family consists of 39 genes whose expressions are tightly controlled and coordinated within the family, during development. We performed a comprehensive analysis of this gene family in cancer settings.

METHODS

Gene correlation analysis was performed using breast cancer data available in The Cancer Genome Atlas (TCGA) and data from the patients admitted to our hospital. We also analyzed the data of normal breast tissue (GSE20437). We next collected gene expression and prognosis data of breast cancer patients (GSE11121, GSE7390, GSE3494, and GSE2990) and performed unsupervised hierarchal clustering by the HOX gene expression pattern and compared prognosis. We additionally performed this analysis to leukemia (available in TCGA) and sarcoma (GSE20196) data.

RESULTS

Gene correlation analysis showed that the proximal HOX genes exhibit strong interactions and are expressed together in breast cancer, similar to the expression observed during development. However, in normal breast tissue, less interactions were observed. Breast cancer microarray meta-data classified by the HOX gene expression pattern predicted the prognosis of luminal B breast cancer patients (p = 0.016). Leukemia (p = 0.00016) and sarcoma (p = 0.018) presented similar results. The Wnt signaling pathway, one of the major upstream signals of HOX genes in development, was activated in the poor prognostic group. Interestingly, poor prognostic cancer presented stronger correlation in the gene family compared to favorable prognostic cancer.

CONCLUSION

Comprehensive analysis of the HOX family demonstrated their similar roles in cancer and development, and indicated that the strong interaction of HOX genes might be specific to malignancies, especially in the case of poor prognostic cancer.

Functional Long Non-coding RNAs Evolve from Junk Transcripts

Transcriptome studies reveal pervasive transcription of complex genomes, such as those of mammals. Despite popular arguments for functionality of most, if not all, of these transcripts, genome-wide analysis of selective constraints indicates that most of the produced RNA are junk. However, junk is not garbage. On the contrary, junk transcripts provide the raw material for the evolution of diverse long non-coding (lnc) RNAs by non-adaptive mechanisms, such as constructive neutral evolution. The generation of many novel functional entities, such as lncRNAs, that fuels organismal complexity does not seem to be driven by strong positive selection. Rather, the weak selection regime that dominates the evolution of most multicellular eukaryotes provides ample material for functional innovation with relatively little adaptation involved.

Polycomb Repressive Complexes in Hox Gene Regulation: Silencing and Beyond: The Functional Dynamics of Polycomb Repressive Complexes in Hox Gene Regulation

The coordinated expression of the Hox gene family encoding transcription factors is critical for proper embryonic development and patterning. Major efforts have thus been dedicated to understanding mechanisms controlling Hox expression. In addition to the temporal and spatial sequential activation of Hox genes, proper embryonic development requires that Hox genes get differentially silenced in a cell-type specific manner as development proceeds. Factors contributing to Hox silencing include the polycomb repressive complexes (PRCs), which control gene expression through epigenetic modifications. This review focuses on PRC-dependent regulation of the Hox genes and is aimed at integrating the growing complexity of PRC functional properties in the context of Hox regulation. In particular, mechanisms underlying PRC binding dynamics as well as a series of studies that have revealed the impact of PRC on the 3D organization of the genome is discussed, which has a significant role on Hox regulation during development.

Long Noncoding RNAs and Repetitive Elements: Junk or Intimate Evolutionary Partners?

Our recent ability to sequence entire genomes, along with all of their transcribed RNAs, has led to the surprising finding that only ∼1% of the human genome is used to encode proteins. This finding has led to vigorous debate over the functional importance of the transcribed but untranslated portions of the genome. Currently, scientists tend to assume coding genes are functional until proven not to be, while the opposite is true for noncoding genes. This review takes a new look at the evidence for and against widespread noncoding gene functionality. We focus in particular on long noncoding RNA (noncoding RNAs longer than 200 nucleotides) genes and their 'junk' associates, transposable elements, and satellite repeats. Taken together, the suggestion put forward is that more of this junk DNA may be functional than nonfunctional and that noncoding RNAs and transposable elements act symbiotically to drive evolution.

The various and shared roles of lncRNAs during development

lncRNAs, genes transcribed but not translated, longer than 200 nucleotides, are classified as a separate class of nonprotein coding genes. Since their discovery, largely from RNAseq data, a number of pioneer studies have begun to unravel its myriad functions, centered on gene expression regulation, suggesting developmental and evolutionary conservation. Since they do not code for proteins and have no open reading frames, their functional constraints likely differ from that of protein coding genes, or of genes where the majority of the nucleotide sequence is required for function, like tRNAs. This has complicated assessment of both developmental and evolutionary conservation, and the identification of homologs in different species. Here we argue that other characteristics: general synteny and particular chromosomal placement regardless of sequence, sequence micro-motifs, and secondary structure allow for "homologs" to be identified and compared, confirming developmental and evolutionary conservation of lncRNAs. We conclude exemplifying a case in point: that of the evolutionarily conserved lncRNA acal, characterized and required for embryogenesis in Drosophila.

HOXA cluster gene expression during osteoblast differentiation involves epigenetic control

The HOXA gene cluster is generally recognized as a pivotal mediator of positional identity in the skeletal system, expression of different orthologues conferring alternative locational phenotype of the vertebrate bone. Strikingly, however, the molecular mechanisms that regulate orthologue-specific expression of different HOXA cluster members in gestating osteoblasts remain largely obscure, but in analogy to the processes observed in acute lymphatic leukemia it is assumed that alternative methylation of HOXA promoter regions drives position specific expression patterns. In an effort to understand HOXA cluster gene expression in osteogenesis we characterize both expression and the epigenetic landscape of the HOXA gene cluster during in vitro osteoblast formation from mesenchymal precursors. We observe that osteoblast formation per se provokes strong upregulation of HOXA gene cluster expression, in particular of midcluster genes, and paradoxal downregulation of HOXA7 and HOXA10. These differences in expression appear related to promoter methylation. LnRNAs HOTAIR and HOTTIP, known to modulate HOXA expression, are also regulated by their promoter methylation processing, but do not correlate with HOXA cluster expression profile. We thus conclude that HOXA expression is profoundly regulated during osteoblast differentiation through canonical methylation-dependent mechanisms but not through the flanking lnRNAs.

Identification of an acute myeloid leukaemia associated noncoding somatic mutation at 3' end of HOXA cluster

Noncoding somatic mutations have been demonstrated to play important role in tumourigenesis. Here we show that there exists an acute myeloid leukaemia associated noncoding somatic mutation at 3' terminal of conserved HOXA cluster. The mutation was identified in the bone marrow blasts but not peripheral blood mononuclear cells or buccal cells of two M3 (acute promyelocytic leukaemia, APL) type patients from 45 acute myeloid leukaemia patients. The mutation also existed in a pair of twins one of them developed acute myeloid leukaemia M4 (acute myelomonocytic leukaemia) type. The mutation resides in about 2-kb downstream of HOXA1 gene where a functional retinoic acid response element is located and also bound by histone demethylase KDM3B. Reporter assay showed that the mutation results in the upregulation of transcriptional activity and unresponsiveness to retinoic acid receptor. To sum up, we identified a new acute myeloid leukaemia associated noncoding somatic mutation.

Hox genes: Downstream "effectors" of retinoic acid signaling in vertebrate embryogenesis

One of the major regulatory challenges of animal development is to precisely coordinate in space and time the formation, specification, and patterning of cells that underlie elaboration of the basic body plan. How does the vertebrate plan for the nervous and hematopoietic systems, heart, limbs, digestive, and reproductive organs derive from seemingly similar population of cells? These systems are initially established and patterned along the anteroposterior axis (AP) by opposing signaling gradients that lead to the activation of gene regulatory networks involved in axial specification, including the Hox genes. The retinoid signaling pathway is one of the key signaling gradients coupled to the establishment of axial patterning. The nested domains of Hox gene expression, which provide a combinatorial code for axial patterning, arise in part through a differential response to retinoic acid (RA) diffusing from anabolic centers established within the embryo during development. Hence, Hox genes are important direct effectors of retinoid signaling in embryogenesis. This review focuses on describing current knowledge on the complex mechanisms and regulatory processes, which govern the response of Hox genes to RA in several tissue contexts including the nervous system during vertebrate development.

Overexpression of BP1, an isoform of Homeobox Gene DLX4, promotes cell proliferation, migration and predicts poor prognosis in endometrial cancer

The expression of homeobox gene DLX4 has been verified in some tumors, but not in endometrial cancer. We found that expression of DLX7, a splicing isoform of DLX4, did not show any significant difference in expression between endometrial cancer and endometrium. However, BP1, another splicing isoform of DLX4, was highly expressed in endometrial cancer, and its expression was positively correlated with patient prognosis, cancer pathological grade, tumor invasion and metastasis. Lentiviral-mediated expression of BP1 in HEC-1-B cells accelerated the cell cycle progression from G0/G1 into S phase, and promoted cell proliferation and migration both in vitro and in vivo. Real-time PCR and western blotting showed that the expression levels of p15, p21 and E-cadherin significantly decreased, and levels of cyclinD1 and MMP-2 increased in endometrial cancer cells. In conclusion, our results demonstrate that high expression of BP1 is associated with poor prognosis in patients with endometrial cancer and promotes cell proliferation and migration.

Exploring Differential Transcriptome between Jejunal and Cecal Tissue of Broiler Chickens

The study proposed an exploratory functional analysis on differential gene expression of the jejunum and of cecum in chickens. For this study, 150 Ross 308 male chickens were randomly allotted in six pens (25 birds/pen) and fed the same commercial diet. From 19 birds of 42 days of age, jejunum and cecum mucosae were collected for RNA extraction for transcriptome microarray analysis. Differentially expressed genes (DEGs) submitted to DAVID (Database for Annotation, Visualization, and Integrated Discovery) and Gene Set Enrichment Analysis (GSEA) software evidenced enriched gene clusters for biological functions differentiated in the tissues. DAVID analysis in the jejunum showed enriched annotations for cell membrane integral components, PPAR (peroxisome proliferator-activated receptor) signaling pathway, and peroxisome and lipid metabolism, and showed DEGs for gluconeogenesis, not previously reported in chicken jejunum. The cecum showed enriched annotations for disulfide bond category, cysteine and methionine metabolism, glycoprotein category, cell cycle, and extracellular matrix (ECM). GSEA analysis in the jejunum showed peroxisome and PPAR signaling pathway-related gene sets, as found with DAVID, and gene sets for immune regulation, tryptophan and histidine metabolism, and renin-angiotensin system, like in mammals. The cecum showed cell cycle and regulation processes, as well as ECM receptor interaction and focal adhesion-related gene sets. Typical intestinal functions specific for the gut site and interesting functional genes groups emerged, revealing tissue-related key aspects which future studies might take advantage of.

Emerging role of HOX genes and their related long noncoding RNAs in lung cancer

The transcription factor homeobox (Hox) proteins are the master regulator for the embryonic development. Studies have identified new functions for HOX in the regulation of metabolism and other primary cellular processes in humans. Their dysregulation has been observed in a variety of cancers and accumulating evidence has revealed the crucial role of HOX in cancer progression, metastasis, and resistance to therapy. HOX-related long non-coding RNAs (lncRNAs) became the most attracting lncRNAs recently that play critical role in gene regulation and chromatin dynamics in cancers. In this review, we explore the roles of HOX and their related lncRNAs in lung cancer, indicating HOX genes as potential therapeutic targets in lung cancer.

Identification of long non-coding RNAs expressed in knee and hip osteoarthritic cartilage

OBJECTIVE

Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in gene expression; however, little is known about the lincRNA expression changes that occur in osteoarthritis (OA). Here we aimed to define a transcriptome of lncRNAs in OA cartilage, specifically comparing the lincRNA transcriptome of knee and hip cartilage.

METHOD

RNA-seq was performed on nucleic acid extracted from hip cartilage from patients undergoing joint replacement surgery because of either OA (n = 10) or because of a neck of femur fracture (NOF; n = 6). After transcript alignment, counts were performed using Salmon and differential expression for ENSEMBL lincRNAs determined using DESeq2. Hip RNA-seq lincRNA expression was compared to a knee dataset (ArrayExpress; E-MTAB-4304). ChIP-seq data from ENCODE was used to determine whether lincRNAs were associated with promoters (plncRNA) or unidirectional enhancer-like regulatory elements (elncRNAs).

RESULTS

Our analysis of the hip transcriptome identified 1692 expressed Transcripts Per Million (TPM ≥1) Ensembl lincRNAs, of which 198 were significantly (FDR ≤0.05) differentially expressed in OA vs normal (NOF) cartilage. Similar analysis of knee cartilage transcriptome identified 648 Emsembl lincRNAs with 93 significantly (FDR ≤0.05) differentially expressed in intact vs damaged cartilage. In total, 1834 lincRNAs were expressed in both hip and knee cartilage, with a highly significant correlation in expression between the two cartilages.

CONCLUSION

This is the first study to use RNA-seq to map and compare the lincRNA transcriptomes of hip and knee cartilage. We propose that lincRNAs expressed selectively in cartilage, or showing differential expression in OA, will play a role in cartilage homoeostasis.

Hox gene expression determines cell fate of adult periosteal stem/progenitor cells

Hox genes are evolutionarily conserved transcription factors that during embryonic development function as master regulators of positional identity. In postnatal life, the function of Hox proteins is less clear: Hox genes are expressed during tissue repair, but in this context their function(s) are largely unknown. Here we show that Hox genes are expressed in periosteal stem/progenitor cells in a distribution similar to that during embryonic development. Using unbiased sequencing, we established that periosteal stem/progenitor cells from distinct anatomic sites within the skeleton significantly differ in their transcriptome, and that Hox expression status best defines these differences. Lastly, we provide evidence that Hox gene expression is one potential mechanism that maintains periosteal stem/progenitor cells in a more primitive, tripotent state, while suppression of Hox genes leads to fate changes with loss of tripotency. Together, our data describe an adult role of Hox genes other than positional identity, and the modulatory role of Hox genes in fate decisions may offer potential druggable targets for the treatment of fractures, non-unions and bone defects.

Obesity impacts the regulation of miR-10b and its targets in primary breast tumors

BACKGROUND

Obesity increases breast cancer (BC) risk in post-menopausal women by mostly unknown molecular mechanisms which may partly be regulated by microRNAs (miRNAs).

METHODS

We isolated RNA from paired benign and malignant biopsies from 83 BC patients and determined miRNA profiles in samples from 12 women at the extremes of the BMI distribution by RNA-seq. Candidates were validated in all samples. Associations between miR-10b expression and validated target transcript levels, and effects of targeted manipulation of miR-10b levels in a primary BC cell line on proliferation and invasion potential, were explored.

RESULTS

Of the 148 miRNAs robustly expressed in breast tissues, the levels of miR-21, miR-10b, miR-451a, miR-30c, and miR-378d were significantly associated with presence of cancer. Of these, miR-10b showed a stronger down-regulation in the tumors of the obese subjects, as opposed to the lean. In ductal but not lobular tumors, significant inverse correlations were observed between the tumor levels of miR-10b and miR-30c and the mRNA levels of cancer-relevant target genes SRSF1, PIEZO1, MAPRE1, CDKN2A, TP-53 and TRA2B, as well as tumor grade. Suppression of miR-10b levels in BT-549 primary BC-derived cells increased cell proliferation and invasive capacity, while exogenous miR-10b mimic decreased invasion. Manipulation of miR-10b levels also inversely affected the mRNA levels of miR-10b targets BCL2L11, PIEZO1 and NCOR2.

CONCLUSIONS

Our findings suggest that miR-10b may be a mediator between obesity and cancer in post-menopausal women, regulating several known cancer-relevant genes. MiR-10b expression may have diagnostic and therapeutic implications for the incidence and prognosis of BC in obese women.

Long noncoding RNA Hoxb3os is dysregulated in autosomal dominant polycystic kidney disease and regulates mTOR signaling

Autosomal dominant polycystic kidney disease (ADPKD) is a debilitating disease that is characterized by the accumulation of numerous fluid-filled cysts in the kidney. ADPKD is primarily caused by mutations in two genes, PKD1 and PKD2 Long noncoding RNAs (lncRNA), defined by a length >200 nucleotides and absence of a long ORF, have recently emerged as epigenetic regulators of development and disease; however, their involvement in PKD has not been explored previously. Here, we performed deep RNA-Seq to identify lncRNAs that are dysregulated in two orthologous mouse models of ADPKD (kidney-specific Pkd1 and Pkd2 mutant mice). We identified a kidney-specific, evolutionarily conserved lncRNA called Hoxb3os that was down-regulated in cystic kidneys from Pkd1 and Pkd2 mutant mice. The human ortholog HOXB3-AS1 was down-regulated in cystic kidneys from ADPKD patients. Hoxb3os was highly expressed in renal tubules in adult WT mice, whereas its expression was lost in the cyst epithelium of mutant mice. To investigate the function of Hoxb3os, we utilized CRISPR/Cas9 to knock out its expression in mIMCD3 cells. Deletion of Hoxb3os resulted in increased phosphorylation of mTOR and its downstream targets, including p70 S6 kinase, ribosomal protein S6, and the translation repressor 4E-BP1. Consistent with activation of mTORC1 signaling, Hoxb3os mutant cells displayed increased mitochondrial respiration. The Hoxb3os mutant phenotype was partially rescued upon re-expression of Hoxb3os in knockout cells. These findings identify Hoxb3os as a novel lncRNA that is down-regulated in ADPKD and regulates mTOR signaling and mitochondrial respiration.

From Heterochromatin to Long Noncoding RNAs in Drosophila: Expanding the Arena of Gene Function and Regulation

Recent years have witnessed a remarkable interest in exploring the significance of pervasive noncoding transcripts in diverse eukaryotes. Classical cytogenetic studies using the Drosophila model system unraveled the perplexing attributes and "functions" of the "gene"-poor heterochromatin. Recent molecular studies in the fly model are likewise revealing the very diverse and significant roles played by long noncoding RNAs (lncRNAs) in development, gene regulation, chromatin organization, cell and nuclear architecture, etc. There has been a rapid increase in the number of identified lncRNAs, although a much larger number still remains unknown. The diversity of modes of actions and functions of the limited number of Drosophila lncRNAs, which have been examined, already reflects the profound roles of such RNAs in generating and sustaining the biological complexities of eukaryotes. Several of the known Drosophila lncRNAs originate as independent sense or antisense transcripts from promoter or intergenic, intronic, or 5'/3'-UTR regions, while many of them are independent genes that produce only lncRNAs or coding as well as noncoding RNAs. The different lncRNAs affect chromatin organization (local or large-scale pan-chromosomal), transcription, RNA processing/stability, or translation either directly through interaction with their target DNA sequences or indirectly by acting as intermediary molecules for specific regulatory proteins or may act as decoys/sinks, or storage sites for specific proteins or groups of proteins, or may provide a structural framework for the assembly of substructures in nucleus/cytoplasm. It is interesting that many of the "functions" alluded to heterochromatin in earlier cytogenetic studies appear to find correlates with the known subtle as well as far-reaching actions of the different small and long noncoding RNAs. Further studies exploiting the very rich and powerful genetic and molecular resources available for the Drosophila model are expected to unravel the mystery underlying the long reach of ncRNAs.

MicroRNA filters Hox temporal transcription noise to confer boundary formation in the spinal cord

The initial rostrocaudal patterning of the neural tube leads to differential expression of Hox genes that contribute to the specification of motor neuron (MN) subtype identity. Although several 3' Hox mRNAs are expressed in progenitors in a noisy manner, these Hox proteins are not expressed in the progenitors and only become detectable in postmitotic MNs. MicroRNA biogenesis impairment leads to precocious expression and propagates the noise of Hoxa5 at the protein level, resulting in an imprecise Hoxa5-Hoxc8 boundary. Here we uncover, using in silico simulation, two feed-forward Hox-miRNA loops accounting for the precocious and noisy Hoxa5 expression, as well as an ill-defined boundary phenotype in Dicer mutants. Finally, we identify mir-27 as a major regulator coordinating the temporal delay and spatial boundary of Hox protein expression. Our results provide a novel trans Hox-miRNA circuit filtering transcription noise and controlling the timing of protein expression to confer robust individual MN identity.

Chromosome folding and its regulation in health and disease

There are many ways in which cells may not adequately behave or respond to their environment, and the molecular mechanisms leading to these defects are as diverse as they are many. In this review, we report on how spatial chromatin organization contributes to the proper expression of genes, relating how CTCF-one of its main architects-contributes to gene regulation. We also touch on the emerging role of long noncoding RNAs in shaping chromatin organization and activity. The HOX gene clusters have been used as paradigm in the study of various biological pathways, and the overview we provide gives emphasis to what research on these loci has revealed about chromatin architecture and its regulation in the control of gene expression.

Origin and evolution of the metazoan non-coding regulatory genome

Animals rely on genomic regulatory systems to direct the dynamic spatiotemporal and cell-type specific gene expression that is essential for the development and maintenance of a multicellular lifestyle. Although it is widely appreciated that these systems ultimately evolved from genomic regulatory mechanisms present in single-celled stem metazoans, it remains unclear how this occurred. Here, we focus on the contribution of the non-coding portion of the genome to the evolution of animal gene regulation, specifically on recent insights from non-bilaterian metazoan lineages, and unicellular and colonial holozoan sister taxa. High-throughput next-generation sequencing, largely in bilaterian model species, has led to the discovery of tens of thousands of non-coding RNA genes (ncRNAs), including short, long and circular forms, and uncovered the central roles they play in development. Based on the analysis of non-bilaterian metazoan, unicellular holozoan and fungal genomes, the evolution of some ncRNAs, such as Piwi-interacting RNAs, correlates with the emergence of metazoan multicellularity, while others, including microRNAs, long non-coding RNAs and circular RNAs, appear to be more ancient. Analysis of non-coding regulatory DNA and histone post-translational modifications have revealed that some cis-regulatory mechanisms, such as those associated with proximal promoters, are present in non-animal holozoans, while others appear to be metazoan innovations, most notably distal enhancers. In contrast, the cohesin-CTCF system for regulating higher-order chromatin structure and enhancer-promoter long-range interactions appears to be restricted to bilaterians. Taken together, most bilaterian non-coding regulatory mechanisms appear to have originated before the divergence of crown metazoans. However, differential expansion of non-coding RNA and cis-regulatory DNA repertoires in bilaterians may account for their increased regulatory and morphological complexity relative to non-bilaterians.

HOXA5 localization in postnatal and adult mouse brain is suggestive of regulatory roles in postmitotic neurons

Hoxa5 is a member of the Hox gene family, which plays critical roles in successive steps of the central nervous system formation during embryonic and fetal development. Hoxa5 expression in the adult mouse brain has been reported, suggesting that this gene may be functionally required in the brain after birth. To provide further insight into the Hoxa5 expression pattern and potential functions in the brain, we have characterized its neuroanatomical profile from embryonic stages to adulthood. While most Hox mapping studies have been based solely on transcript analysis, we extended our analysis to HOXA5 protein localization in adulthood using specific antibodies. Our results show that Hoxa5 expression appears in the most caudal part of the hindbrain at fetal stages, where it is maintained until adulthood. In the medulla oblongata and pons, we detected Hoxa5 expression in many precerebellar neurons and in several nuclei implicated in the control of autonomic functions. In these territories, the HOXA5 protein is present solely in neurons, specifically in γ-aminobutyric acid (GABA)ergic, glutamatergic, and catecholaminergic neurons. Finally, we also detected Hoxa5 transcripts, but not the HOXA5 protein, in the thalamus and the cortex, from postnatal stages to adult stages, and in the cerebellum at adulthood. We provide evidence that some larger variants of Hoxa5 transcripts are present in these territories. Our mapping analysis allowed us to build hypotheses regarding HOXA5 functions in the nervous system after birth, such as a potential role in the establishment and refinement/plasticity of precerebellar circuits during postnatal and adult life. J. Comp. Neurol. 525:1155-1175, 2017. © 2016 Wiley Periodicals, Inc.