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

The non-canonically organized members of MIR395 gene family in Brassica juncea are associated with developmentally regulated, sulfate-stress responsive bidirectional promoters that exhibit orientation-dependent differential transcriptional activity

Several MICRORNA genes belonging to same family or different families are often found in homologous or non-homologous clusters. Among the various classes, head-to-head arranged genes form one of the largest categories of non-canonically organized genes. Such head-to-head arranged, non-canonically organized genes possibly share cis-regulatory region with the intergenic sequence having the potential to function as bi-directional promoter (BDP). The transcriptional regulation of head-to-head arranged genes, especially with bidirectional promoters, remains an enigma. In the past, bidirectional promoters have been characterized for a small set of protein-coding gene pairs in plants; however, to the best of our knowledge, no such study has been carried so far for MICRORNA genes. The present study thus functionally characterizes bidirectional promoters associated with members of MIR395 family, which is evolutionary conserved and is most frequently occurring cluster across plant kingdom. In Arabidopsis thaliana, the MIR395 gene family contains six members with two head-to-head arranged gene pairs- MIR395A-B and MIR395E-F. This organization was found to be conserved at seven loci for MIR395A-B, and eleven loci for MIR395E-F in five Brassica sps. Sequence analysis of the putative bidirectional promoters revealed variation in length, GC content and distribution of strict TATA-box. Comparatively higher level of conservation at both the ends of the bidirectional promoters, corresponding to ca. 250 bp upstream of 5'end of the respective MIRNA precursor, was observed. These conserved regions harbour several abiotic stress (nutrient, salt, drought) and hormone (ABA, ethylene) responsive cis-motifs. Functional characterization of putative bidirectional promoters associated with MIR395A-B and MIR395E-F from Arabidopsis and their respective orthologs from Brassica juncea (Bj_A08 MIR395A-B, Bj_B03 MIR395A-B, Bj_A07.1 MIR395E-F and Bj_A07.2 MIR395E-F) was carried out using a dual-reporter vector with β-glucuronidase (GUS) and Green Fluorescent Protein (GFP). Analysis of transcriptional regulation of the two reporter genes - GUS and GFP during developmental stages confirmed their bidirectional nature. Orientation-dependent differential reporter activity indicated asymmetric nature of the promoters. Comparison of the reporter activity amongst orthologs, paralogs and homeologs revealed regulatory diversification, an outcome expected in polyploid genomes. Interestingly, reporter gene activities driven by selected bidirectional promoters were also observed in anther and siliques apart vegetative tissues indicating role of miR395 in anther and fruit development. Finally, we evaluated the activity of reporter genes driven under transcriptional regulation of bidirectional promoters under normal and sulfate-deprived conditions which revealed asymmetric inducibility under sulfate-starvation, in agreement with the known role of miR395 in sulfate homeostasis.

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

BACKGROUND

In research to improve the quality of transgenic crops, it is often necessary to introduce multiple functionally related genes into recipient plants simultaneously to improve crop genetic traits effectively. Compared with unidirectional promoters, bidirectional promoters simultaneously regulate the expression of multiple genes and improve the efficiency of biotechnology. Therefore, in this study, bidirectional gene pairs were systematically analyzed in Gossypium hirsutum TM-1, and the structure, function and evolutionary relationships of the bidirectional genes were analyzed. The endogenous bidirectional promoters of cotton were mined, and their specific regulatory elements and biological functions were explored to provide useful promoter resources and a theoretical basis for cultivating new cotton germplasms with excellent fiber quality.

RESULTS

Using an improved search model, a total of 1,383 bidirectional transcript pairs were identified in the Gossypium hirsutum TM-1 genome, and their gene structure and functional annotations were systematically analyzed. Thirty bidirectional intergenic sequences were randomly screened for promoter activity analysis via a transient expression system, and 25 intergenic sequences were found to have bidirectional promoter activity. Comparative analysis of the bidirectional gene profiles of the four cotton subspecies revealed that these subspecies presented abundant bidirectional gene pairs with high homology and that the bidirectional genes in the cotton subspecies were more similar in terms of their molecular functions, cellular components and biological processes. In addition, parallel analysis of bidirectional genes in dicotyledons and monocotyledons revealed that abundant bidirectional gene pairs exist in different species. Although the total number of orthologous bidirectional genes was similar, there was a significant difference in the number of orthologous bidirectional gene pairs between dicotyledons and monocotyledons. This evolutionary analysis of the function and structure of homologous bidirectional gene pairs in different varieties and different subspecies of the same species revealed potential pathways by which these gene pairs originated, which may be necessary for the evolution of a new species.

CONCLUSION

In this study, many bidirectional gene pairs in Gossypium hirsutum TM-1 were identified using computer programming, and systematic analysis was conducted to explore their functions and evolutionary relationships. In addition, the promoter activity of the bidirectional intergenic sequences was verified. The combination of computer programming screening, experimental validation and other methods is expected to provide preferred bidirectional promoters for transgenic breeding work via multigene cotransformation methods, and this information is valuable for genetic engineering research and applications.

Structural and Functional Analysis of a Bidirectional Promoter from Gossypium hirsutum in Arabidopsis

Stacked traits have become an important trend in the current development of genomically modified crops. The bidirectional promoter can not only prevent the co-suppression of multigene expression, but also increase the efficiency of the cultivation of transgenic plants with multigenes. In Gossypium hirsutum, Ghrack1 and Ghuhrf1 are head-to-head gene pairs located on chromosome D09. We cloned the 1429-bp intergenic region between the Ghrack1 and Ghuhrf1 genes from Gossypium hirsutum. The cloned DNA fragment GhZU had the characteristics of a bidirectional promoter, with 38.7% G+C content, three CpG islands and no TATA-box. Using gfp and gus as reporter genes, a series of expression vectors were constructed into young leaves of tobacco. The histochemical GUS (Beta-glucuronidase) assay and GFP (green fluorescence protein) detection results indicated that GhZU could drive the expression of the reporter genes gus and gfp simultaneously in both orientations. Furthermore, we transformed the expression vectors into Arabidopsis and found that GUS was concentrated at vigorous growth sites, such as the leaf tip, the base of the leaves and pod, and the stigma. GFP was also mainly expressed in the epidermis of young leaves. In summary, we determined that the intergenic region GhZU was an orientation-dependent bidirectional promoter, and this is the first report on the bidirectional promoter from Gossypium hirsutum. Our findings in this study are likely to enhance understanding on the regulatory mechanisms of plant bidirectional promoters.

Synthesis of Seed-Specific Bidirectional Promoters for Metabolic Engineering of Anthocyanin-Rich Maize

Tissue-specific promoters play an important role in plant molecular farming. Here, we describe a strategy to modify the tissue specificity of a maize embryo-specific bidirectional promoter PZmBD1. Six types of cis-elements, i.e. RY repeats (R), GCN4 (G), the prolamin box (P), Skn-1 (S), and the ACGT and AACA (A) motifs, were collected and fused to PZmBD1 to generate eight chimeric putative bidirectional promoters. Qualitative and quantitative analysis of reporter genes driven by the promoters showed that two promoters exhibited high seed-specific bidirectional activity in maize transient and stable transformed systems. The stronger one was chosen and fused to the intergenic region of two gene clusters consisting of four anthocyanin biosynthesis-related genes (ZmBz1, ZmBz2, ZmC1 and ZmR2) and seven reporter genes, resulting in the first embryo and endosperm anthocyanin-rich purple maize. Anthocyanin analysis showed that the total anthocyanin content reaches 2,910 mg kg-1 DW in transgenic maize and cyanidin is the major anthocyanin in transgenic maize, as in natural varieties. The expression profile analysis of endogenous genes showed that the anthocyanin biosynthesis pathway was activated by two transgenic transcription factor genes ZmC1 and ZmR2. Our results indicate that both the modification strategy and these functionally characterized tissue-specific bidirectional promoters generated could be used for genetic research and development of plant biotechnology products. The anthocyanin-rich purple maize could provide economic natural colorants for the food and beverage industry, and valuable germplasm for developing anthocyanin-rich fresh corn.

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.

Analysis of Promoters of Arabidopsis thaliana Divergent Gene Pair SERAT3;2 and IDH-III Shows SERAT3;2 Promoter is Nested Within the IDH-III Promoter

Intergenic regions of divergent gene pairs show bidirectional promoter activity but whether regulatory sequences for gene expression in opposite directions are shared is not established. In this study, promoters of divergently arranged gene pair At4g35640-At4g35650 (SERAT3;2-IDH-III) of Arabidopsis thaliana were analyzed to identify overlapping regulatory regions. Both genes showed the highest expression in flower buds and flowers. 5' RACE experiments extended the intergenic region from 161 bp shown in TAIR annotation to 512 bp. GUS analysis of transgenic A. thaliana plants carrying the 691 bp fragment (512 bp intergenic region plus 5' UTR of both the genes) linked to uidA gene revealed that SERAT3;2 promoter drives gene expression in the tapetum, whereas IDH-III promoter functions specifically in microspores/pollen. Serial 5' deletion of the 691 bp fragment showed SERAT3;2 promoter extends up to -355 position, whereas IDH-III promoter encompasses the 512 bp intergenic region. In transgenics, uidA transcript levels were lower than native SERAT3;2 and IDH-III transcripts indicating presence of additional cis regulatory elements beyond the 691 bp fragment. The present study demonstrated for the first time occurrence of a nested promoter in plants and identified a novel bidirectional promoter capable of driving gene expression in tapetum and microspores/pollen.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

Functional characterization of open chromatin in bidirectional promoters of rice

Bidirectional gene pairs tend to be highly coregulated and function in similar biological processes in eukaryotic genomes. Structural features and functional consequences of bidirectional promoters (BDPs) have received considerable attention among diverse species. However, the underlying mechanisms responsible for the bidirectional transcription and coexpression of BDPs remain poorly understood in plants. In this study, we integrated DNase-seq, RNA-seq, ChIP-seq and MNase-seq data and investigated the effect of physical DNase I hypersensitive site (DHS) positions on the transcription of rice BDPs. We found that the physical position of a DHS relative to the TSS of bidirectional gene pairs can affect the expression of the corresponding genes: the closer a DHS is to the TSS, the higher is the expression level of the genes. Most importantly, we observed that the distribution of DHSs plays a significant role in the regulation of transcription and the coexpression of gene pairs, which are possibly mediated by orchestrating the positioning of histone marks and canonical nucleosomes around BDPs. Our results demonstrate that the combined actions of chromatin structures with DHSs, which contain functional cis-elements for interaction with transcriptional machinery, may play an important role in the regulation of the bidirectional transcription or coexpression in rice BDPs. Our findings may help to enhance the understanding of DHSs in the regulation of bidirectional gene pairs.

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.

DLGP: A database for lineage-conserved and lineage-specific gene pairs in animal and plant genomes

The conservation of gene organization in the genome with lineage-specificity is an invaluable resource to decipher their potential functionality with diverse selective constraints, especially in higher animals and plants. Gene pairs appear to be the minimal structure for such kind of gene clusters that tend to reside in their preferred locations, representing the distinctive genomic characteristics in single species or a given lineage. Despite gene families having been investigated in a widespread manner, the definition of gene pair families in various taxa still lacks adequate attention. To address this issue, we report DLGP (http://lcgbase.big.ac.cn/DLGP/) that stores the pre-calculated lineage-based gene pairs in currently available 134 animal and plant genomes and inspect them under the same analytical framework, bringing out a set of innovational features. First, the taxonomy or lineage has been classified into four levels such as Kingdom, Phylum, Class and Order. It adopts all-to-all comparison strategy to identify the possible conserved gene pairs in all species for each gene pair in certain species and reckon those that are conserved in over a significant proportion of species in a given lineage (e.g. Primates, Diptera or Poales) as the lineage-conserved gene pairs. Furthermore, it predicts the lineage-specific gene pairs by retaining the above-mentioned lineage-conserved gene pairs that are not conserved in any other lineages. Second, it carries out pairwise comparison for the gene pairs between two compared species and creates the table including all the conserved gene pairs and the image elucidating the conservation degree of gene pairs in chromosomal level. Third, it supplies gene order browser to extend gene pairs to gene clusters, allowing users to view the evolution dynamics in the gene context in an intuitive manner. This database will be able to facilitate the particular comparison between animals and plants, between vertebrates and arthropods, and between monocots and eudicots, accounting for the significant contribution of gene pairs to speciation and diversification in specific lineages.