Systematic analysis of head-to-head gene organization: evolutionary conservation and potential biological relevance

RHNO1: at the crossroads of DNA replication stress, DNA repair, and cancer

The DNA replication stress (DRS) response is a crucial homeostatic mechanism for maintaining genome integrity in the face of intrinsic and extrinsic barriers to DNA replication. Importantly, DRS is often significantly increased in tumor cells, making tumors dependent on the cellular DRS response for growth and survival. Rad9-Hus1-Rad1 Interacting Nuclear Orphan 1 (RHNO1), a protein involved in the DRS response, has recently emerged as a potential therapeutic target in cancer. RHNO1 interacts with the 9-1-1 checkpoint clamp and TopBP1 to activate the ATR/Chk1 signaling pathway, the crucial mediator of the DRS response. Moreover, RHNO1 was also recently identified as a key facilitator of theta-mediated end joining (TMEJ), a DNA repair mechanism implicated in cancer progression and chemoresistance. In this literature review, we provide an overview of our current understanding of RHNO1, including its structure, function in the DRS response, and role in DNA repair, and discuss its potential as a cancer therapeutic target. Therapeutic targeting of RHNO1 holds promise for tumors with elevated DRS as well as tumors with DNA repair deficiencies, including homologous recombination DNA repair deficient (HRD) tumors. Further investigation into RHNO1 function in cancer, and development of approaches to target RHNO1, are expected to yield novel strategies for cancer treatment.

Novel insights on the positive correlation between sense and antisense pairs on gene expression

A considerable proportion of the eukaryotic genome undergoes transcription, leading to the generation of noncoding RNA molecules that lack protein-coding information and are not subjected to translation. These noncoding RNAs (ncRNAs) are well recognized to have essential roles in several biological processes. Long noncoding RNAs (lncRNAs) represent the most extensive category of ncRNAs found in the human genome. Much research has focused on investigating the roles of cis-acting lncRNAs in the regulation of specific target gene expression. In the majority of instances, the regulation of sense gene expression by its corresponding antisense pair occurs in a negative (discordant) manner, resulting in the suppression of the target genes. The notion that a negative correlation exists between sense and antisense pairings is, however, not universally valid. In fact, several recent studies have reported a positive relationship between corresponding cis antisense pairs within plants, budding yeast, and mammalian cancer cells. The positive (concordant) correlation between anti-sense and sense transcripts leads to an increase in the level of the sense transcript within the same genomic loci. In addition, mechanisms such as altering chromatin structure, the formation of R loops, and the recruitment of transcription factors can either enhance transcription or stabilize sense transcripts through their antisense pairs. The primary objective of this work is to provide a comprehensive understanding of both aspects of antisense regulation, specifically focusing on the positive correlation between sense and antisense transcripts in the context of eukaryotic gene expression, including its implications towards cancer progression. This article is categorized under: RNA Processing > 3' End Processing Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

A plant cell death-inducing protein from litchi interacts with Peronophythora litchii pectate lyase and enhances plant resistance

Cell wall degrading enzymes, including pectate lyases (PeLs), released by plant pathogens, break down protective barriers and/or activate host immunity. The direct interactions between PeLs and plant immune-related proteins remain unclear. We identify two PeLs, PlPeL1 and PlPeL1-like, critical for full virulence of Peronophythora litchii on litchi (Litchi chinensis). These proteins enhance plant susceptibility to oomycete pathogens in a PeL enzymatic activity-dependent manner. However, LcPIP1, a plant immune regulator secreted by litchi, binds to PlPeL1/PlPeL1-like, and attenuates PlPeL1/PlPeL1-like induced plant susceptibility to Phytophthora capsici. LcPIP1 also induces cell death and various immune responses in Nicotiana benthamiana. Conserved in plants, LcPIP1 homologs bear a conserved "VDMASG" motif and exhibit immunity-inducing activity. Furthermore, SERK3 interacts with LcPIP1 and is required for LcPIP1-induced cell death. NbPIP1 participates in immune responses triggered by the PAMP protein INF1. In summary, our study reveals the dual roles of PlPeL1/PlPeL1-like in plant-pathogen interactions: enhancing pathogen virulence through PeL enzymatic activity while also being targeted by LcPIP1, thus enhancing plant immunity.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Cysteine Enrichment Mediates Co-Option of Uricase in Reptilian Skin and Transition to Uricotelism

Uric acid is the main means of nitrogen excretion in uricotelic vertebrates (birds and reptiles) and the end product of purine catabolism in humans and a few other mammals. While uricase is inactivated in mammals unable to degrade urate, the presence of orthologous genes without inactivating mutations in avian and reptilian genomes is unexplained. Here we show that the Gallus gallus gene we name cysteine-rich urate oxidase (CRUOX) encodes a functional protein representing a unique case of cysteine enrichment in the evolution of vertebrate orthologous genes. CRUOX retains the ability to catalyze urate oxidation to hydrogen peroxide and 5-hydroxyisourate (HIU), albeit with a 100-fold reduced efficiency. However, differently from all uricases hitherto characterized, it can also facilitate urate regeneration from HIU, a catalytic property that we propose depends on its enrichment in cysteine residues. X-ray structural analysis highlights differences in the active site compared to known orthologs and suggests a mechanism for cysteine-mediated self-aggregation under H2O2-oxidative conditions. Cysteine enrichment was concurrent with the transition to uricotelism and a shift in gene expression from the liver to the skin where CRUOX is co-expressed with β-keratins. Therefore, the loss of urate degradation in amniotes has followed opposite evolutionary trajectories: while uricase has been eliminated by pseudogenization in some mammals, it has been repurposed as a redox-sensitive enzyme in the reptilian skin.

Pan-Cancer Analysis Reveals the Prognostic Potential of the THAP9/THAP9-AS1 Sense-Antisense Gene Pair in Human Cancers

Human THAP9, which encodes a domesticated transposase of unknown function, and lncRNA THAP9-AS1 (THAP9-antisense1) are arranged head-to-head on opposite DNA strands, forming a sense and antisense gene pair. We predict that there is a bidirectional promoter that potentially regulates the expression of THAP9 and THAP9-AS1. Although both THAP9 and THAP9-AS1 are reported to be involved in various cancers, their correlative roles on each other's expression has not been explored. We analyzed the expression levels, prognosis, and predicted biological functions of the two genes across different cancer datasets (TCGA, GTEx). We observed that although the expression levels of the two genes, THAP9 and THAP9-AS1, varied in different tumors, the expression of the gene pair was strongly correlated with patient prognosis; higher expression of the gene pair was usually linked to poor overall and disease-free survival. Thus, THAP9 and THAP9-AS1 may serve as potential clinical biomarkers of tumor prognosis. Further, we performed a gene co-expression analysis (using WGCNA) followed by a differential gene correlation analysis (DGCA) across 22 cancers to identify genes that share the expression pattern of THAP9 and THAP9-AS1. Interestingly, in both normal and cancer samples, THAP9 and THAP9-AS1 often co-express; moreover, their expression is positively correlated in each cancer type, suggesting the coordinated regulation of this H2H gene pair.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development

KEY MESSAGE

Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

Bidirectional promoters: an enigmatic genome architecture and their roles in cancers

Bidirectional promoters are the transcription regulatory regions of genes positioned head-to-head on opposite strands. Specific sequence signals, chromatin modifications and three-dimensional structures of the transcription site facilitate the unconventional yet tightly regulated transcription proceeding in both directions from these promoters. Mutations or aberrant epigenetic changes can lead to abnormal enhanced or reduced expression from either of the bidirectionally transcribed genes resulting in tumorigenesis. Moreover, bidirectionally transcribed genes might also contribute towards the immune regulation in tumor microenvironment. In this review, we aimed to expound the characteristic features of bidirectional promoters alongside their transcriptional regulations, and ultimately, the association of these enigmatic genomic elements in different cancers.

COL4A3 expression in asthmatic epithelium depends on intronic methylation and ZNF263 binding

Background

Reduction of COL4A3, one of the six isoforms of collagen 4, in asthmatic airways results in increased inflammation and angiogenesis, implicating it as a central part of asthma pathogenesis. However, to date, the path underlying these diminished COL4A3 levels has been elusive. This study investigated a possible mechanism underlying the reduction of COL4A3 expression.

Methods

Bronchial biopsies of 76 patients with asthma and 83 controls were subjected to RNA-sequencing and DNA methylation bead arrays to identify expression and methylation changes. The binding of ZNF263 was analysed by chromatin-immunoprecipitation sequencing coupled with quantitative (q)PCR. Effects of ZNF263 silencing, using small interfering RNA, on the COL4A3 expression were studied using qPCR.

Results

COL4A3 expression was significantly reduced in bronchial biopsies compared to healthy controls, whereas DNA methylation levels at cg11797365 were increased. COL4A3 expression levels were significantly low in asthmatics without inhaled corticosteroid (ICS) use, whereas the expression was not statistically different between asthmatics using ICS and controls. Methylation levels at cg11797365 in vitro were increased upon consecutive rhinovirus infections.

Conclusion

Our data indicate an epigenetic modification as a contributing factor for the loss of COL4A3 expression in asthmatic airway epithelium.

An epigenetic modification interrupts ZNF263 binding, which may contribute to the loss of COL4A3 expression in asthmatic airway epitheliumhttps://bit.ly/39cZbyn

Tissue-specific activation of gene expression by the Synergistic Activation Mediator (SAM) CRISPRa system in mice

CRISPR-based transcriptional activation is a powerful tool for functional gene interrogation; however, delivery difficulties have limited its applications in vivo. Here, we created a mouse model expressing all components of the CRISPR-Cas9 guide RNA-directed Synergistic Activation Mediator (SAM) from a single transcript that is capable of activating target genes in a tissue-specific manner. We optimized Lipid Nanoparticles and Adeno-Associated Virus guide RNA delivery approaches to achieve expression modulation of one or more genes in vivo. We utilized the SAM mouse model to generate a hypercholesteremia disease state that we could bidirectionally modulate with various guide RNAs. Additionally, we applied SAM to optimize gene expression in a humanized Transthyretin mouse model to recapitulate human expression levels. These results demonstrate that the SAM gene activation platform can facilitate in vivo research and drug discovery.

Co-regulation and function of FOXM1/RHNO1 bidirectional genes in cancer

The FOXM1 transcription factor is an oncoprotein and a top biomarker of poor prognosis in human cancer. Overexpression and activation of FOXM1 is frequent in high-grade serous carcinoma (HGSC), the most common and lethal form of human ovarian cancer, and is linked to copy number gains at chromosome 12p13.33. We show that FOXM1 is co-amplified and co-expressed with RHNO1, a gene involved in the ATR-Chk1 signaling pathway that functions in the DNA replication stress response. We demonstrate that FOXM1 and RHNO1 are head-to-head (i.e., bidirectional) genes (BDG) regulated by a bidirectional promoter (BDP) (named F/R-BDP). FOXM1 and RHNO1 each promote oncogenic phenotypes in HGSC cells, including clonogenic growth, DNA homologous recombination repair, and poly-ADP ribosylase inhibitor resistance. FOXM1 and RHNO1 are one of the first examples of oncogenic BDG, and therapeutic targeting of FOXM1/RHNO1 BDG is a potential therapeutic approach for ovarian and other cancers.

Molecular evolution and regulation of DHN melanin-related gene clusters are closely related to adaptation of different melanin-producing fungi

Many genes responsible for melanin biosynthesis in fungi were physically linked together. The PKS gene clusters in most of the melanin-producing fungi were regulated by the Cmr1. It was found that a close rearrangement of the PKS gene clusters had evolved in most of the melanin-producing fungi and various functions of melanin in them were beneficial to their adaptation to the changing environments. The melanin-producing fungi had undergone at least five large-scale differentiations, making their PKS gene clusters be quickly evolved and the fungi be adapted to different harsh environments. The recent gene losses and expansion were remarkably frequent in the PKS gene clusters, leading to their rapid evolution and adaptation of their hosts to different environments. The PKS gene and the CMR1 gene in them were subject to a strong co-evolution, but the horizontal gene transfer events might have occurred in the genome-duplicated species, Aspergillus and Penicillium.

Functional Expression, Purification and Identification of Interaction Partners of PACRG

PACRG (Parkin co-regulated gene) shares a bi-directional promoter with the Parkinson’s disease-associated gene Parkin, but the physiological roles of PACRG have not yet been fully elucidated. Recombinant expression methods are indispensable for protein structural and functional studies. In this study, the coding region of PACRG was cloned to a conventional vector pQE80L, as well as two cold-shock vectors pCold II and pCold-GST, respectively. The constructs were transformed into Escherichia coli (DE3), and the target proteins were overexpressed. The results showed that the cold-shock vectors are more suitable for PACRG expression. The soluble recombinant proteins were purified with Ni²⁺ chelating column, glutathione S-transferase (GST) affinity chromatography and gel filtration. His₆ pull down assay and LC-MS/MS were carried out for identification of PACRG-binding proteins in HEK293T cell lysates, and a total number of 74 proteins were identified as potential interaction partners of PACRG. GO (Gene ontology) enrichment analysis (FunRich) of the 74 proteins revealed multiple molecular functions and biological processes. The highest proportion of the 74 proteins functioned as transcription regulator and transcription factor activity, suggesting that PACRG may play important roles in regulation of gene transcription.

Pan-Cancer Analysis of Head-to-Head Gene Pairs in Terms of Transcriptional Activity, Co-expression and Regulation

Background

Head-to-Head (H2H) gene pairs are regulated by bidirectional promoters and divergently transcribed from opposite DNA strands with transcription start sites (TSSs) separated within 1 kb. H2H organization is ancient and conserved, and H2H pairs tend to exhibit similar expression patterns. Although some H2H genes have been reported to be associated with disease and cancer, there is a lack of systematic studies on H2H organization in the scenario of cancer development.

Methods

Human H2H gene pairs were identified based on GENCODE hg19 and the functional relevance of H2H pairs was explored through function enrichment and semantic similarity analysis. To investigate the association between H2H organization and carcinogenesis, pan-cancer differential analysis of H2H genes about transcriptional activity, co-expression and transcriptional regulation by transcription factors and enhancers were performed based on data from The Cancer Genome Atlas. Cox proportional hazards regression model and log-rank test were used to determine the prognostic powers of H2H pairs.

Results

In the present study, we first updated H2H genes from 1,447 to 3,150 pairs, from which the peak group with TSS distance of 1–100 was observed as expected in our previous work. It was found that housekeeping genes, mitochondrial-functional associated genes and cancer genes tend to be organized in H2H arrangement. Pan-cancer analysis indicates that H2H genes are transcriptionally active than random genes in both normal and cancer tissues, but H2H pairs display higher correlation in cancer than in normal. Particularly, housekeeping H2H pairs are differentially correlated much more significantly than non-housekeeping H2H pairs are. Some of differentially correlated H2H pairs were found to be associated with prognosis. The alteration of TF similarity seems to contribute to differential co-expression of H2H pairs during carcinogenesis; meanwhile remote enhancers also at least partly explain the differential co-expression and co-regulation of H2H pairs.

Conclusion

H2H pairs tend to show much stronger positive expression correlation in cancer than in normal due to differential regulation of bidirectional promoters. The study provides insights into the significance of H2H organization in carcinogenesis and the underlying dysfunctional regulation mechanisms. Those differentially correlated H2H pairs associated with survival have the potential to be prognostic biomarkers and therapeutic targets for cancer.

Unexplored regulatory sequences of divergently paired GLA and HNRNPH2 loci pertinent to Fabry disease in human kidney and skin cells: Presence of an active bidirectional promoter

Fabry disease (FD) is a rare hereditary disorder characterized by a wide range of symptoms caused by a variety of mutations in the galactosidase α (GLA) gene. The heterogeneous nuclear ribonucleoprotein (HNRNPH2) gene is divergently paired with GLA on chromosome X and is thought to be implicated in FD. However, insufficient information is available on the regulatory mechanisms associated with the expression of HNRNPH2 and the GLA loci. Therefore, the current study performed bioinformatics analyses to assess the GLA and HNRNPH2 loci and investigate the regulatory mechanisms involved in the expression of each gene. The regulatory mechanisms underlying GLA and HNRNPH2 were revealed. The expression of each gene was associated with a bidirectional promoter (BDP) characterized by the absence of TATA box motifs and the presence of specific transcription factor binding sites (TFBSs) and a CpG Island (CGI). The nuclear run-on transcription assay confirmed the activity of BDP GLA and HNRNPH2 transcription in 293T. Methylation-specific PCR analysis demonstrated a statistically significant variation in the DNA methylation pattern of BDP in several cell lines, including human adult epidermal keratinocytes (AEKs), human renal glomerular endothelial cells, human renal epithelial cells and 293T cells. The highest observed significance was demonstrated in AEKs (P<0.05). The results of the chromatin-immunoprecipitation assay using 293T cells identified specific TFBS motifs for Yin Yang 1 and nuclear respiratory factor 1 transcription factors in BDPs. The National Center for Biotechnology Information-single nucleotide polymorphism database revealed pathogenic variants in the BDP sequence. Additionally, a previously reported variant associated with a severe heterozygous female case of GLA FD was mapped in BDP. The results of the present study suggested that the expression of the divergent paired loci, GLA and HNRNPH2, were controlled by BDP. Mutations in BDP may also serve a role in FD and may explain clinical disease diversity.

The parkin-coregulated gene product PACRG promotes TNF signaling by stabilizing LUBAC

The Parkin-coregulated gene (PACRG), which encodes a protein of unknown function, shares a bidirectional promoter with Parkin (PRKN), which encodes an E3 ubiquitin ligase. Because PRKN is important in mitochondrial quality control and protection against stress, we tested whether PACRG also affected these pathways in various cultured human cell lines and in mouse embryonic fibroblasts. PACRG did not play a role in mitophagy but did play a role in tumor necrosis factor (TNF) signaling. Similarly to Parkin, PACRG promoted nuclear factor κB (NF-κB) activation in response to TNF. TNF-induced nuclear translocation of the NF-κB subunit p65 and NF-κB-dependent transcription were decreased in PACRG-deficient cells. Defective canonical NF-κB activation in the absence of PACRG was accompanied by a decrease in linear ubiquitylation mediated by the linear ubiquitin chain assembly complex (LUBAC), which is composed of the two E3 ubiquitin ligases HOIP and HOIL-1L and the adaptor protein SHARPIN. Upon TNF stimulation, PACRG was recruited to the activated TNF receptor complex and interacted with LUBAC components. PACRG functionally replaced SHARPIN in this context. In SHARPIN-deficient cells, PACRG prevented LUBAC destabilization, restored HOIP-dependent linear ubiquitylation, and protected cells from TNF-induced apoptosis. This function of PACRG in positively regulating TNF signaling may help to explain the association of PACRG and PRKN polymorphisms with an increased susceptibility to intracellular pathogens.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

Gene Repositioning Is Under Constraints After Evolutionary Conserved Gene Neighborhood Separate

Genes are not randomly distributed on eukaryotic chromosomes. Some neighboring genes show order conservation among species, while some neighboring genes separate during evolution even though their neighborhoods are conserved in some species. Here, I investigated whether after-separation gene repositioning is under natural selection for evolutionary conserved gene neighborhoods compared with nonconserved neighborhoods. After separation, genes with conserved neighborhoods show low-expression divergence between the after-separation species and the before-separation species. After genes separate from their conserved gene neighbors, their after-separation gene neighbors tend to show coexpression and coprotein complex with their before-separation gene neighbors. These results indicate evolutionary constraints on the selection of neighboring genes after evolutionary conserved gene neighborhoods separate.

Large-scale mammalian genome rearrangements coincide with chromatin interactions

Motivation

Genome rearrangements drastically change gene order along great stretches of a chromosome. There has been initial evidence that these apparently non-local events in the 1D sense may have breakpoints that are close in the 3D sense. We harness the power of the Double Cut and Join model of genome rearrangement, along with Hi-C chromosome conformation capture data to test this hypothesis between human and mouse.

Results

We devise novel statistical tests that show that indeed, rearrangement scenarios that transform the human into the mouse gene order are enriched for pairs of breakpoints that have frequent chromosome interactions. This is observed for both intra-chromosomal breakpoint pairs, as well as for inter-chromosomal pairs. For intra-chromosomal rearrangements, the enrichment exists from close (<20 Mb) to very distant (100 Mb) pairs. Further, the pattern exists across multiple cell lines in Hi-C data produced by different laboratories and at different stages of the cell cycle. We show that similarities in the contact frequencies between these many experiments contribute to the enrichment. We conclude that either (i) rearrangements usually involve breakpoints that are spatially close or (ii) there is selection against rearrangements that act on spatially distant breakpoints.

Availability and implementation

Our pipeline is freely available at https://bitbucket.org/thekswenson/locality.

Supplementary information

Supplementary data are available at Bioinformatics online.

Neighboring genes are closely related to whole genome duplications after their separation

BACKGROUND

The gene order in a eukaryotic genome is not random. Some neighboring genes show specific similarities, while others become separated during evolution. Whole genome duplication events (WGDs) have been recognized as an important evolutionary force. The potential relationship between the separation of neighboring genes and WGDs needs to be investigated. In this study, we investigated whether there is a potential relationship between separated neighboring gene pairs and WGDs, and the mechanism by which neighboring genes are separated. Additionally, we studied whether neighboring genes tend to show intrachromosomal colocalization after their neighborhood was disrupted and the factors facilitating the intrachromosomal colocalization of separated neighboring genes.

RESULTS

The separation of neighboring gene pairs is closely related to whole genome duplication events. Furthermore, we found that there is a double linear relationship between separated neighboring genes, total genes, and WGDs. The process of separation of neighboring genes caused by WGDs is also not random but abides by the double linear model. Separated neighboring gene pairs tend to show intrachromosomal colocalization. The conservativism of separated neighboring genes and histone modification facilitate the intrachromosomal colocalization of neighboring genes after their separation.

CONCLUSIONS

These results provide new insight into the understanding of evolutionary roles of locations and the relationship of neighboring gene pairs with whole genome duplications. Furthermore, understanding the proposed mechanism for intrachromosomal colocalization of separated genes benefits our knowledge of chromosomal interactions in the nucleus.

Kynurenine 3-Monooxygenase Gene Associated With Nicotine Initiation and Addiction: Analysis of Novel Regulatory Features at 5' and 3'-Regions

Tobacco smoking is widespread behavior in Qatar and worldwide and is considered one of the major preventable causes of ill health and death. Nicotine is part of tobacco smoke that causes numerous health risks and is incredibly addictive; it binds to the α7 nicotinic acetylcholine receptor (α7nAChR) in the brain. Recent studies showed α7nAChR involvement in the initiation and addiction of smoking. Kynurenic acid (KA), a significant tryptophan metabolite, is an antagonist of α7nAChR. Inhibition of kynurenine 3-monooxygenase enzyme encoded by KMO enhances the KA levels. Modulating KMO gene expression could be a useful tactic for the treatment of tobacco initiation and dependence. Since KMO regulation is still poorly understood, we aimed to investigate the 5' and 3'-regulatory factors of KMO gene to advance our knowledge to modulate KMO gene expression. In this study, bioinformatics methods were used to identify the regulatory sequences associated with expression of KMO. The displayed differential expression of KMO mRNA in the same tissue and different tissues suggested the specific usage of the KMO multiple alternative promoters. Eleven KMO alternative promoters identified at 5'-regulatory region contain TATA-Box, lack CpG Island (CGI) and showed dinucleotide base-stacking energy values specific to transcription factor binding sites (TFBSs). The structural features of regulatory sequences can influence the transcription process and cell type-specific expression. The uncharacterized LOC105373233 locus coding for non-coding RNA (ncRNA) located on the reverse strand in a convergent manner at the 3'-side of KMO locus. The two genes likely expressed by a promoter that lacks TATA-Box harbor CGI and two TFBSs linked to the bidirectional transcription, the NRF1, and ZNF14 motifs. We identified two types of microRNA (miR) in the uncharacterized LOC105373233 ncRNA, which are like hsa-miR-5096 and hsa-miR-1285-3p and can target the miR recognition element (MRE) in the KMO mRNA. Pairwise sequence alignment identified 52 nucleotides sequence hosting MRE in the KMO 3' UTR untranslated region complementary to the ncRNA LOC105373233 sequence. We speculate that the identified miRs can modulate the KMO expression and together with alternative promoters at the 5'-regulatory region of KMO might contribute to the development of novel diagnostic and therapeutic algorithm for tobacco smoking.

Tissue-specific impact of FADS cluster variants on FADS1 and FADS2 gene expression

Omega-6 (n-6) and omega-3 (n-3) long (≥ 20 carbon) chain polyunsaturated fatty acids (LC-PUFAs) play a critical role in human health and disease. Biosynthesis of LC-PUFAs from dietary 18 carbon PUFAs in tissues such as the liver is highly associated with genetic variation within the fatty acid desaturase (FADS) gene cluster, containing FADS1 and FADS2 that encode the rate-limiting desaturation enzymes in the LC-PUFA biosynthesis pathway. However, the molecular mechanisms by which FADS genetic variants affect LC-PUFA biosynthesis, and in which tissues, are unclear. The current study examined associations between common single nucleotide polymorphisms (SNPs) within the FADS gene cluster and FADS1 and FADS2 gene expression in 44 different human tissues (sample sizes ranging 70-361) from the Genotype-Tissue Expression (GTEx) Project. FADS1 and FADS2 expression were detected in all 44 tissues. Significant cis-eQTLs (within 1 megabase of each gene, False Discovery Rate, FDR<0.05, as defined by GTEx) were identified in 12 tissues for FADS1 gene expression and 23 tissues for FADS2 gene expression. Six tissues had significant (FDR< 0.05) eQTLs associated with both FADS1 and FADS2 (including artery, esophagus, heart, muscle, nerve, and thyroid). Interestingly, the identified eQTLs were consistently found to be associated in opposite directions for FADS1 and FADS2 expression. Taken together, findings from this study suggest common SNPs within the FADS gene cluster impact the transcription of FADS1 and FADS2 in numerous tissues and raise important questions about how the inverse expression of these two genes impact intermediate molecular (such a LC-PUFA and LC-PUFA-containing glycerolipid levels) and ultimately clinical phenotypes associated with inflammatory diseases and brain health.

Chicken CCDC152 shares an NFYB-regulated bidirectional promoter with a growth hormone receptor antisense transcript and inhibits cells proliferation and migration

The chicken coiled-coil domain-containing protein 152 (CCDC152) recently has been identified as a novel one implicated in cell cycle regulation, cellular proliferation and migration by us. Here we demonstrate that CCDC152 is oriented in a head-to-head configuration with the antisense transcript of growth hormone receptor (GHR) gene. Through serial luciferase reporter assays, we firstly identified a minimal 102 bp intergenic region as a core bidirectional promoter to drive basal transcription in divergent orientations. And site mutation and transient transfected assays showed that nuclear transcription factor Y subunit beta (NFYB) could bind to the CCAAT box and directly transactivate this bidirectional promoter. SiRNA-mediated NFYB depletion could significantly down-regulate the expression of both GHR-AS-I6 and CCDC152. Additionally, the expression of GHR-AS-I6 was significantly up-regulated after CCDC152 overexpression. Overexpression of CCDC152 remarkably reduced cell proliferation and migration through JAK2/STAT signaling pathway. Thus, the GHR-AS-I6-CCDC152 bidirectional transcription unit, as a novel direct target of NFYB, is possibly essential for the accelerated proliferation and motility of different cells.

Pervasive coexpression of spatially proximal genes is buffered at the protein level

Genes are not randomly distributed in the genome. In humans, 10% of protein-coding genes are transcribed from bidirectional promoters and many more are organised in larger clusters. Intriguingly, neighbouring genes are frequently coexpressed but rarely functionally related. Here we show that coexpression of bidirectional gene pairs, and closeby genes in general, is buffered at the protein level. Taking into account the 3D architecture of the genome, we find that co-regulation of spatially close, functionally unrelated genes is pervasive at the transcriptome level, but does not extend to the proteome. We present evidence that non-functional mRNA coexpression in human cells arises from stochastic chromatin fluctuations and direct regulatory interference between spatially close genes. Protein-level buffering likely reflects a lack of coordination of post-transcriptional regulation of functionally unrelated genes. Grouping human genes together along the genome sequence, or through long-range chromosome folding, is associated with reduced expression noise. Our results support the hypothesis that the selection for noise reduction is a major driver of the evolution of genome organisation.

Genetic variability of the activity of bidirectional promoters: a pilot study in bovine muscle

Bidirectional promoters are regulatory regions co-regulating the expression of two neighbouring genes organized in a head-to-head orientation. In recent years, these regulatory regions have been studied in many organisms; however, no investigation to date has been done to analyse the genetic variation of the activity of this type of promoter regions. In our study, we conducted an investigation to first identify bidirectional promoters sharing genes expressed in bovine Longissimus thoracis and then to find genetic variants affecting the activity of some of these bidirectional promoters. Combining bovine gene information and expression data obtained using RNA-Seq, we identified 120 putative bidirectional promoters active in bovine muscle. We experimentally validated in vitro 16 of these bidirectional promoters. Finally, using gene expression and whole-genome genotyping data, we explored the variability of the activity in muscle of the identified bidirectional promoters and discovered genetic variants affecting their activity. We found that the expression level of 77 genes is correlated with the activity of 12 bidirectional promoters. We also identified 57 single nucleotide polymorphisms associated with the activity of 5 bidirectional promoters. To our knowledge, our study is the first analysis in any species of the genetic variability of the activity of bidirectional promoters.

Jmjd6, a JmjC Dioxygenase with Many Interaction Partners and Pleiotropic Functions

Lysyl hydroxylation and arginyl demethylation are post-translational events that are important for many cellular processes. The jumonji domain containing protein 6 (JMJD6) has been reported to catalyze both lysyl hydroxylation and arginyl demethylation on diverse protein substrates. It also interacts directly with RNA. This review summarizes knowledge of JMJD6 functions that have emerged in the last 15 years and considers how a single Jumonji C (JmjC) domain-containing enzyme can target so many different substrates. New links and synergies between the three main proposed functions of Jmjd6 in histone demethylation, promoter proximal pause release of polymerase II and RNA splicing are discussed. The physiological context of the described molecular functions is considered and recently described novel roles for JMJD6 in cancer and immune biology are reviewed. The increased knowledge of JMJD6 functions has wider implications for our general understanding of the JmjC protein family of which JMJD6 is a member.

Mining the Cicer arietinum genome for the mildew locus O (Mlo) gene family and comparative evolutionary analysis of the Mlo genes from Medicago truncatula and some other plant species

The mildew locus O (Mlo) gene family is ubiquitous in land plants. Some members of this gene family are involved in negative regulation of powdery mildew resistance, while others are involved in several other biological functions. Mlo proteins have characteristic seven transmembrane domains and a calmodulin-binding domain at their C-termini, and are associated with plasma membrane. The Mlo gene family has been studied in several economically important cereals, but little information is available on this gene family in the important legumes, Medicago truncatula Gaertn. and Cicer arietinum L. We carried out a comprehensive and comparative investigation of the Mlo gene family in these two species using the genome sequences available at the M. truncatula genome database (Mt v4.0) and NCBI (C. arietinum). A genome-wide homology-based search using Arabidopsis Mlo proteins as query identified 16 MtMlo (M. truncatula Mlo) and 14 CarMlo (C. arietinum Mlo) genes. The MtMlo and CarMlo genes had comparable gene structure, protein sequence and topology. Their chromosomal locations indicated the occurrence of extensive reorganization in the genomes of the two species after their divergence from the common ancestor. A multiple sequence alignment of 53 Mlo proteins from these two and several other species showed a highly conserved sequence block of seven amino acids, viz., L-ETPTW, towards their N-termini. The evolutionary phylogenetic analysis grouped the MtMlo and CarMlo members into four clusters, and most of the MtMlo and CarMlo members formed one-to-one ortholog pairs. The ka/ks analyses indicated that the MtMlo and CarMlo genes are subjected to intense purifying selection.

A potent bidirectional promoter from the monocot cereal Eleusine coracana

Ragi bifunctional α-amylase-trypsin inhibitor (RBI) of Eleusine coracana (L.) Gaertn. (finger millet) simultaneously inhibits α-amylase and trypsin. In continuation of previous work on the cloning, expression and characterization of RBI, a bidirectional promoter from finger millet was explored on the basis of experimental observations. Two trypsin inhibitors were identified while purifying RBI from a trypsin-Sepharose column eluent. Using an FPLC gel filtration column, these three inhibitors were purified to homogeneity and subjected to MALDI-TOF-TOF-MS/MS analysis and N-terminal sequencing. Both ragi trypsin inhibitors (RTIs) showed the same N-terminal sequence and considerable sequence similarity to RBI, indicating the presence of a multigene protease inhibitor family in finger millet. To gain insight into the evolution of these genes, the upstream region of RBI was explored by Genome Walking. Interestingly, on sequencing, a genome walking product of ∼1 Kb showed presence of an N-terminal RBI specific primer sequence twice but in opposite directions and leaving an intervening region of ∼0.9 Kb. The intervening region was presumed to represent an E. coracana bidirectional promoter (EcBDP), intuitively having a divergent RBI-RTI gene pair at two sides. For assaying the bidirectionality of promoter activity, a dual reporter GUS-GFP vector construct was made for plant expression containing the reporter genes at two ends of EcBDP, which was used to transform Agrobacterium tumefaciens LBA 4404. Transient plant transformation by recombinant Agrobacterium cells was carried out in onion scale epidermal cells and finger millet seedling leaves. Simultaneous expression of GUS and GFP under EcBDP established it as a potent natural bidirectional promoter from monocot origin, thereby potentially having vast application in cereal gene manipulations. In addition, inducibility of the EcBDP by either abscisic acid or cold treatment, as determined by transient transformation in onion, would substantiate more precise control of gene expression to mitigate the effects of adverse environmental conditions.

Principles for RNA metabolism and alternative transcription initiation within closely spaced promoters

Mammalian transcriptomes are complex and formed by extensive promoter activity. In addition, gene promoters are largely divergent and initiate transcription of reverse-oriented promoter upstream transcripts (PROMPTs). Although PROMPTs are commonly terminated early, influenced by polyadenylation sites, promoters often cluster so that the divergent activity of one might impact another. Here we found that the distance between promoters strongly correlates with the expression, stability and length of their associated PROMPTs. Adjacent promoters driving divergent mRNA transcription support PROMPT formation, but owing to polyadenylation site constraints, these transcripts tend to spread into the neighboring mRNA on the same strand. This mechanism to derive new alternative mRNA transcription start sites (TSSs) is also evident at closely spaced promoters supporting convergent mRNA transcription. We suggest that basic building blocks of divergently transcribed core promoter pairs, in combination with the wealth of TSSs in mammalian genomes, provide a framework with which evolution shapes transcriptomes.

The intergenic region of the maize defensin-like protein genes Def1 and Def2 functions as an embryo-specific asymmetric bidirectional promoter

Bidirectional promoters are identified in diverse organisms with widely varied genome sizes, including bacteria, yeast, mammals, and plants. However, little research has been done on any individual endogenous bidirectional promoter from plants. Here, we describe a promoter positioned in the intergenic region of two defensin-like protein genes, Def1 and Def2 in maize (Zea mays). We examined the expression profiles of Def1 and Def2 in 14 maize tissues by qRT-PCR, and the results showed that this gene pair was expressed abundantly and specifically in seeds. When fused to either green fluorescent protein (GFP) or β-glucuronidase (GUS) reporter genes, P ZmBD1 , P ZmDef1 , and P ZmDef2 were active and reproduced the expression patterns of both Def1 and Def2 genes in transformed immature maize embryos, as well as in developing seeds of transgenic maize. Comparative analysis revealed that PZmBD1 shared most of the expression characteristics of the two polar promoters, but displayed more stringent embryo specificity, delayed expression initiation, and asymmetric promoter activity. Moreover, a truncated promoter study revealed that the core promoters only exhibit basic bidirectional activity, while interacting with necessary cis-elements, which leads to polarity and different strengths. The sophisticated interaction or counteraction between the core promoter and cis-elements may potentially regulate bidirectional promoters.

Regulatory dissection of the CBX5 and hnRNPA1 bi-directional promoter in human breast cancer cells reveals novel transcript variants differentially associated with HP1α down-regulation in metastatic cells

BACKGROUND

The three members of the human heterochromatin protein 1 (HP1) family of proteins, HP1α, HP1β, and HPγ, are involved in chromatin packing and epigenetic gene regulation. HP1α is encoded from the CBX5 gene and is a suppressor of metastasis. CBX5 is down-regulated at the transcriptional and protein level in metastatic compared to non-metastatic breast cancer. CBX5 shares a bi-directional promoter structure with the hnRNPA1 gene. But whereas CBX5 expression is down-regulated in metastatic cells, hnRNAP1 expression is constant. Here, we address the regulation of CBX5 in human breast cancer.

METHODS

Transient transfection and transposon mediated integration of dual-reporter mini-genes containing the bi-directional hnRNPA1 and CBX5 promoter was performed to investigate transcriptional regulation in breast cancer cell lines. Bioinformatics and functional analysis were performed to characterize transcriptional events specifically regulating CBX5 expression. TSA treatment and Chromatin Immunoprecipitation (ChIP) were performed to investigate the chromatin structure along CBX5 in breast cancer cells. Finally, expression of hnRNPA1 and CBX5 mRNA isoforms were measured by quantitative reverse transcriptase PCR (qRT-PCR) in breast cancer tissue samples.

RESULTS

We demonstrate that an hnRNPA1 and CBX5 bi-directional core promoter fragment does not comprise intrinsic capacity for specific CBX5 down-regulation in metastatic cells. Characterization of transcriptional events in the 20 kb CBX5 intron 1 revealed existence of several novel CBX5 transcripts. Two of these encode consensus HP1α protein but used autonomous promoters in intron 1 by which HP1α expression could be de-coupled from the bi-directional promoter. In addition, another CBX5 transcriptional isoform, STET, was discovered. This transcript includes CBX5 exon 1 and part of intron 1 sequences but lacks inclusion of HP1α encoding exons. Inverse correlation between STET and HP1α coding CBX5 mRNA expression was observed in breast cancer cell lines and tissue samples from breast cancer patients.

CONCLUSION

We find that HP1α is down-regulated in a mechanism involving CBX5 promoter downstream sequences and that regulation through alternative polyadenylation and splicing generates a transcript, STET, with potential importance in carcinogenesis.

Beyond the histone tale: HP1α deregulation in breast cancer epigenetics

Heterochromatin protein 1α (HP1α) encoded from the CBX5-gene is an evolutionary conserved protein that binds histone H3 di- or tri-methylated at position lysine 9 (H3K9me2/3), a hallmark for heterochromatin, and has an essential role in forming higher order chromatin structures. HP1α has diverse functions in heterochromatin formation, gene regulation, and mitotic progression, and forms complex networks of gene, RNA, and protein interactions. Emerging evidence has shown that HP1α serves a unique biological role in breast cancer related processes and in particular for epigenetic control mechanisms involved in aberrant cell proliferation and metastasis. However, how HP1α deregulation plays dual mechanistic functions for cancer cell proliferation and metastasis suppression and the underlying cellular mechanisms are not yet comprehensively described. In this paper we provide an overview of the role of HP1α as a new sight of epigenetics in proliferation and metastasis of human breast cancer. This highlights the importance of addressing HP1α in breast cancer diagnostics and therapeutics.

Long non-coding RNAs in corticogenesis: deciphering the non-coding code of the brain

Evidence on the role of long non-coding (lnc) RNAs has been accumulating over decades, but it has been only recently that advances in sequencing technologies have allowed the field to fully appreciate their abundance and diversity. Despite this, only a handful of lncRNAs have been phenotypically or mechanistically studied. Moreover, novel lncRNAs and new classes of RNAs are being discovered at growing pace, suggesting that this class of molecules may have functions as diverse as protein-coding genes. Interestingly, the brain is the organ where lncRNAs have the most peculiar features including the highest number of lncRNAs that are expressed, proportion of tissue-specific lncRNAs and highest signals of evolutionary conservation. In this work, we critically review the current knowledge about the steps that have led to the identification of the non-coding transcriptome including the general features of lncRNAs in different contexts in terms of both their genomic organisation, evolutionary origin, patterns of expression, and function in the developing and adult mammalian brain.

The gene pair PRR11 and SKA2 shares a NF-Y-regulated bidirectional promoter and contributes to lung cancer development

Head-to-head gene pairs represent a unique feature of gene organization in eukaryotes, accounting for >10% of genes in the human genome. Identification and functional analysis of such gene pairs is only in its infancy. Recently, we identified PRR11 as a novel cancer-related gene that is implicated in cell cycle and lung cancer. Here we demonstrate that PRR11 is oriented in a head-to-head configuration with its neighboring gene, SKA2. 5'-RACE assay revealed that the intergenic spacer region between the two genes is <500 bp. Serial luciferase reporter assays demonstrated that a minimal 80-bp intergenic region functions as a core bidirectional promoter to drive basal transcription in both the PRR11 and SKA2 orientations. EMSA and ChIP assays demonstrated that NF-Y binds to and directly transactivates the PRR11-SKA2 bidirectional promoter. SiRNA-mediated NF-Y depletion significantly downregulated PRR11 and SKA2 expression. Expression of both PRR11 and SKA2 was significantly upregulated in lung cancer. Expression of the two genes was highly correlated with each other and with NF-Y expression. Remarkably, high expression of both PRR11 and SKA2 was associated with poorer prognosis in lung cancer patients compared with high expression of one gene or low expression of both genes. Knockdown of PRR11 and/or SKA2 remarkably reduced cell proliferation, migration, and invasion in lung cancer cells. Thus, the PRR11-SKA2 bidirectional transcription unit, which is a novel direct target of NF-Y, is essential for the accelerated proliferation and motility of lung cancer cells and may represent a potential target in the diagnosis and/or treatment of human lung cancer.

Transcriptional and post-transcriptional regulation of transmembrane protein 132A

Transmembrane protein 132A (TMEM132A) was first isolated from rat brain using PCR-selected cDNA subtraction, and it was found to be predominantly expressed in the brain. However, the transcriptional regulation of the TMEM132A gene has not been fully characterized. In this study, we characterized the promoter activity of the 880-bp region upstream of the mouse TMEM132A, identifying several putative sites recognized by transcription factors, which are highly conserved between the mouse and human TMEM132A genes. Using four different mouse cell lines (Neuro2a, NSC-34, NIH3T3, and Raw264.7), we first evaluated the intrinsic levels of TMEM132A mRNA and protein expression. Interestingly, TMEM132A mRNA was expressed in all four cell lines, whereas the protein was negligible in Raw264.7 cells even by transfection of TMEM132A gene. Then, we analyzed the TMEM132A promoter activity using serial deleted constructs, finding it was nearly same pattern in all four cell lines. A mutational analysis of the TMEM132A promoter identified a critical region for its activation just upstream of the transcriptional start site. Finally, we investigated the levels of TMEM132A mRNA and protein after exposure to five different neurotoxic stimuli, including thapsigargin, tunicamycin, serum starvation, homocysteine, and hydrogen peroxide. Treatment with thapsigargin, a calcium modulating agent, markedly attenuated the levels of TMEM132A mRNA and protein in NSC-34 cells. These results give new insight into the mechanisms involved in regulating TMEM132A expression, and suggest that several transcriptional and post-transcriptional pathways regulate TMEM132A expression under developmental and pathophysiological conditions.

Delta-like 4 mRNA is regulated by adjacent natural antisense transcripts

Background

Recent evidence suggests that a majority of RNAs in the genome do not code for proteins. They are located in the sense (S) or antisense (AS) orientation and, to date, the functional significance of these non-coding RNAs (ncRNAs) is poorly understood. Here, we examined the relationship between S and AS transcripts in the regulation of a key angiogenesis gene, Delta-like 4 (Dll4).

Methods

Rapid Amplification of cDNA Ends (RACE) method was used to identify natural antisense transcripts in the Dll4 gene locus in murine and human endothelial cells, referred to as Dll4 Anti-Sense (Dll4-AS). Messenger RNA (mRNA) levels of Dll4 and Dll4-AS were quantified by real-time PCR. The function of Dll4-AS was investigated by overexpression and knocking down of Dll4-AS.

Results

Dll4-AS comprises of three isoforms that map proximal to the Dll4 promoter region. Expression patterns of Dll4-AS isoforms vary among different endothelial cell lines, but are always congruent with those of Dll4. A dual promoter element in the Dll4 locus has been identified that controls the expression of both transcripts. Both Dll4-AS and Dll4 are sensitive to cellular density in that higher cellular density favors their expression. Exogenous Dll4 stimuli such as VEGF, FGF and Notch signaling inhibitor altered both DLL4-AS and DLL4 expression suggesting co-regulation of the transcripts. Also, knocking down of Dll4-AS results in down-regulation of Dll4 expression. As a consequence, endothelial cell proliferation and migration increases in vitro, and sprout formation increases. The regulation of Dll4 by Dll4-AS was also conserved in vivo.

Conclusion

A novel form of non-coding RNA-mediated regulation at the Dll4 locus contributes to vascular developmental processes such as cell proliferation, migration and sprouting.

Electronic supplementary material

The online version of this article (doi:10.1186/s13221-015-0028-9) contains supplementary material, which is available to authorized users.

Neighboring Genes Show Correlated Evolution in Gene Expression

When considering the evolution of a gene’s expression profile, we commonly assume that this is unaffected by its genomic neighborhood. This is, however, in contrast to what we know about the lack of autonomy between neighboring genes in gene expression profiles in extant taxa. Indeed, in all eukaryotic genomes genes of similar expression-profile tend to cluster, reflecting chromatin level dynamics. Does it follow that if a gene increases expression in a particular lineage then the genomic neighbors will also increase in their expression or is gene expression evolution autonomous? To address this here we consider evolution of human gene expression since the human-chimp common ancestor, allowing for both variation in estimation of current expression level and error in Bayesian estimation of the ancestral state. We find that in all tissues and both sexes, the change in gene expression of a focal gene on average predicts the change in gene expression of neighbors. The effect is highly pronounced in the immediate vicinity (<100 kb) but extends much further. Sex-specific expression change is also genomically clustered. As genes increasing their expression in humans tend to avoid nuclear lamina domains and be enriched for the gene activator 5-hydroxymethylcytosine, we conclude that, most probably owing to chromatin level control of gene expression, a change in gene expression of one gene likely affects the expression evolution of neighbors, what we term expression piggybacking, an analog of hitchhiking.

Structural and functional analysis of an asymmetric bidirectional promoter in Arabidopsis thaliana

Bidirectional promoters are relatively abundant in eukaryotic genomes, suggesting that they have an important biological significance. As yet, few of these promoters have been characterized in detail. Here, using a promoter::GUS transgene approach has revealed that the intergenic region of Arabidopsis thaliana divergent genes At1g71850 and At1g71860 is an asymmetric bidirectional promoter, which exhibits an orientation-dependent expression profile. The strength of the forward promoter was greater than that of the reverse promoter, and their tissue specificities were not identical. Deletion analyses revealed that this bidirectional promoter could be divided into three functional regions. The basal level and tissue specificity of the promoter in the reverse orientation were regulated positively by region II and negatively by region III, whereas promoter activity in the forward orientation was regulated negatively by region II and positively by region I. Thus the 52-bp stretch of region II had a dual function, enhancing expression in the reverse orientation and suppressing it in the forward orientation. These results demonstrated that the activity of the At1g71850-At1g71860 bidirectional promoter was modulated by complex interactions between both positive and negative cis-acting elements. These findings will enhance our understanding of the regulatory mechanisms of plant bidirectional promoters.

Long Non-Coding RNA BST2/BISPR is Induced by IFN and Regulates the Expression of the Antiviral Factor Tetherin

Many long non-coding RNAs (lncRNAs) are expressed in cells but only a few have been well characterized. In these cases, lncRNAs have been shown to be key regulators of several cellular processes. Therefore, there is a great need to understand the function of more lncRNAs and their regulation in response to stimuli. Interferon (IFN) is a key molecule in the cellular antiviral response. IFN binding to its receptor activates transcription of several IFN-stimulated genes (ISGs) that function as potent antivirals. In addition, several ISGs are positive or negative regulators of the IFN pathway. This is essential to ensure a strong antiviral response and a later return of the cell to homeostasis. As the ISGs described to date are coding genes, we sought to determine whether IFN also regulates the expression of long non-coding ISGs. To this aim, we used RNA sequencing to analyze the transcriptome of control and HuH7 cells treated with IFNα2. The results show that IFN-treatment regulates the expression of several unknown non-coding transcripts. We have validated two lncRNAs upregulated after treatment with different doses of type I IFNα2 in different cells or with type III IFNλ. These lncRNAs were also induced by influenza and vesicular stomatitis virus mutants unable to block the IFN response, but not by several wild-type lytic viruses tested. These lncRNA genes were named lncISG15 and lncBST2 as they are located close to ISGs ISG15 and BST2, respectively. Interestingly, inhibition experiments showed that lncBST2 is a positive regulator of BST2. Therefore lncBST2 has been renamed BISPR, from BST2 IFN-stimulated positive regulator. Our results may have therapeutic implications as lncBST2/BISPR, but also lncISG15 and their coding neighbors, are increased in cells infected with hepatitis C virus and in the liver of infected patients. These results allow us to hypothesize that several lncRNAs could be activated by IFN to control the potency of the antiviral IFN response.

Uncovering the functional constraints underlying the genomic organization of the odorant-binding protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Identification and functional analysis of long non-coding RNAs in mouse cleavage stage embryonic development based on single cell transcriptome data

Background

Long non-coding RNAs (lncRNAs) regulate embryonic development and cell fate decision in various ways, such as modulation of chromatin modification and post-transcription regulation of gene expression. However, the profiles and roles of lncRNAs in early mammalian development have not yet been demonstrated. Here, we reported a comprehensive analysis of mouse cleavage stage embryonic lncRNA profiles based on public single-cell RNA-seq data.

Results

We reconstructed 50,006 high-confidence transcripts in 22,827 loci, and identified 5563 novel lncRNAs from 3492 loci expressed in mouse cleavage stage embryos. These lncRNAs share similar characteristics with previously reported vertebrate lncRNAs, such as relatively short length, low exon number, low expression level and low sequence conservation. Expression profile analysis revealed that the profiles of lncRNA vary considerably at different stages of cleavage stage embryos, suggesting that many lncRNAs in cleavage stage embryos are stage-specifically expressed. Co-expression network analysis suggested many lncRNAs in cleavage stage embryos are associated with cell cycle regulation, transcription, translation and oxidative phosphorylation to regulate the process of cleavage stage embryonic development.

Conclusions

This study provides the first catalog of lncRNAs expressed in mouse cleavage stage embryos and gives a revealing insight into the molecular mechanism responsible for early embryonic development.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-845) contains supplementary material, which is available to authorized users.

Large-scale phylogenetic classification of fungal chitin synthases and identification of a putative cell-wall metabolism gene cluster in Aspergillus genomes

The cell wall is a protective and versatile structure distributed in all fungi. The component responsible for its rigidity is chitin, a product of chitin synthase (Chsp) enzymes. There are seven classes of chitin synthase genes (CHS) and the amount and type encoded in fungal genomes varies considerably from one species to another. Previous Chsp sequence analyses focused on their study as individual units, regardless of genomic context. The identification of blocks of conserved genes between genomes can provide important clues about the interactions and localization of chitin synthases. On the present study, we carried out an in silico search of all putative Chsp encoded in 54 full fungal genomes, encompassing 21 orders from five phyla. Phylogenetic studies of these Chsp were able to confidently classify 347 out of the 369 Chsp identified (94%). Patterns in the distribution of Chsp related to taxonomy were identified, the most prominent being related to the type of fungal growth. More importantly, a synteny analysis for genomic blocks centered on class IV Chsp (the most abundant and widely distributed Chsp class) identified a putative cell wall metabolism gene cluster in members of the genus Aspergillus, the first such association reported for any fungal genome.

Transcriptional regulation and spatial interactions of head-to-head genes

Background

In eukaryotic genomes, about 10% of genes are arranged in a head-to-head (H2H) orientation, and the distance between the transcription start sites of each gene pair is closer than 1 kb. Two genes in an H2H pair are prone to co-express and co-function. There have been many studies on bidirectional promoters. However, the mechanism by which H2H genes are regulated at the transcriptional level still needs further clarification, especially with regard to the co-regulation of H2H pairs. In this study, we first used the Hi-C data of chromatin linkages to identify spatially interacting H2H pairs, and then integrated ChIP-seq data to compare H2H gene pairs with and without evidence of spatial interactions in terms of their binding transcription factors (TFs). Using ChIP-seq and DNase-seq data, histones and DNase associated with H2H pairs were identified. Furthermore, we looked into the connections between H2H genes in a human co-expression network.

Results

We found that i) Similar to the behaviour of two genes within an H2H pair (intra-H2H pair), a gene pair involving two distinct H2H pairs (inter-H2H pair) which interact with each other spatially, share common transcription factors (TFs); ii) TFs of intra- and inter-H2H pairs are distributed differently. Factors such as HEY1, GABP, Sin3Ak-20, POL2, E2F6, and c-MYC are essential for the bidirectional transcription of intra-H2H pairs; while factors like CTCF, BDP1, GATA2, RAD21, and POL3 play important roles in coherently regulating inter-H2H pairs; iii) H2H gene blocks are enriched with hypersensitive DNase and modified histones, which participate in active transcriptions; and iv) H2H genes tend to be highly connected compared with non-H2H genes in the human co-expression network.

Conclusions

Our findings shed new light on the mechanism of the transcriptional regulation of H2H genes through their linear and spatial interactions. For intra-H2H gene pairs, transcription factors regulate their transcriptions through bidirectional promoters, whereas for inter-H2H gene pairs, transcription factors are likely to regulate their activities depending on the spatial interaction of H2H gene pairs. In this way, two distinctive groups of transcription factors mediate intra- and inter-H2H gene transcriptions respectively, resulting in a highly compact gene regulatory network.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-519) contains supplementary material, which is available to authorized users.

Terminate and make a loop: regulation of transcriptional directionality

•

Transcriptional directionality is controlled by premature transcription termination.

•

Transcriptional directionality is enforced by gene looping.

•

mRNA-specific termination signals and factors are required for gene looping.

Bidirectional promoters are a common feature of many eukaryotic organisms from yeast to humans. RNA Polymerase II that is recruited to this type of promoter can start transcribing in either direction using alternative DNA strands as the template. Such promiscuous transcription can lead to the synthesis of unwanted transcripts that may have negative effects on gene expression. Recent studies have identified transcription termination and gene looping as critical players in the enforcement of promoter directionality. Interestingly, both mechanisms share key components. Here, we focus on recent findings relating to the transcriptional output of bidirectional promoters.

Construction and validation of a novel dual reporter vector for studying mammalian bidirectional promoters

Regulation of gene expression plays important role in cellular functions. With the development of sequencing techniques, more and more genomes are available and genome-wide analyses of genomic structures that may affect gene expression regulation are now possible. Analyses of several genomes have found a class of regulatory regions that contain elements that initiate transcription of two different genes positioned with a head-to-head arrangement in two opposite directions. These regulatory regions are known as bidirectional promoters. Although bidirectional promoters have been known for years, recent genome-scale studies have shown that the regulation of the expression of up to 10% of the genes are controlled by bidirectional promoters. These findings are based mostly on computational work and only a limited number of putative bidirectional promoters have been experimentally validated. Developing methods to study bidirectional promoters will allow researchers to understand how these regions are regulated and the roles that divergent transcription plays in the expression of genes. Here, we have developed a novel dual-fluorescence reporter gene vector to study the transcriptional output of mammalian bidirectional promoters. We demonstrate that this vector is capable of expressing reporter genes under the control of bidirectional promoters, using the known human OSGEP/APEX bidirectional promoter.

Identification and functional characterization of bidirectional gene pairs and their intergenic regions in maize

Background

Bidirectional gene pairs exist as a specific form of gene organization in microorganisms and mammals as well as in model plant species, such as Arabidopsis and rice. Little is known about bidirectional gene pairs in maize, which has a large genome and is one of the most important grain crops.

Results

We conducted a genome-wide search in maize using genome sequencing results from the inbred line B73. In total, 1696 bidirectional transcript pairs were identified using a modified search model. We functionally characterized the promoter activity of the intergenic regions of most of the bidirectional transcript pairs that were expressed in embryos using a maize embryo transient expression system. A comparative study of bidirectional gene pairs performed for three monocot (Zea mays, Sorghum bicolor and Oryza sativa) and two dicot (Arabidopsis thaliana and Glycine max) plant genomes showed that bidirectional gene pairs were abundant in the five plant species. Orthologous bidirectional gene pairs were clearly distinguishable between the monocot and dicot species although the total numbers of orthologous bidirectional genes were similar. Analysis of the gene pairs using the Blast2GO software suite showed that the molecular functions (MF), cellular components (CC), and biological processes (BP) associated with the bidirectional transcripts were similar among the five plant species.

Conclusions

The evolutionary analysis of the function and structure of orthologous bidirectional gene pairs in various plant species revealed a potential pathway of their origin, which may be required for the evolution of a new species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-338) contains supplementary material, which is available to authorized users.