Searching for bidirectional promoters in Arabidopsis thaliana

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.

Do DEEPER ROOTING 1 Homologs Regulate the Lateral Root Slope Angle in Cucumber (Cucumis sativus)?

The architecture of the root system is fundamental to plant productivity. The rate of root growth, the density of lateral roots, and the spatial structure of lateral and adventitious roots determine the developmental plasticity of the root system in response to changes in environmental conditions. One of the genes involved in the regulation of the slope angle of lateral roots is DEEPER ROOTING 1 (DRO1). Its orthologs and paralogs have been identified in rice, Arabidopsis, and several other species. However, nothing is known about the formation of the slope angle of lateral roots in species with the initiation of lateral root primordia within the parental root meristem. To address this knowledge gap, we identified orthologs and paralogs of the DRO1 gene in cucumber (Cucumis sativus) using a phylogenetic analysis of IGT protein family members. Differences in the transcriptional response of CsDRO1, CsDRO1-LIKE1 (CsDRO1L1), and CsDRO1-LIKE2 (CsDRO1L2) to exogenous auxin were analyzed. The results showed that only CsDRO1L1 is auxin-responsive. An analysis of promoter-reporter fusions demonstrated that the CsDRO1, CsDRO1L1, and CsDRO1L2 genes were expressed in the meristem in cell files of the central cylinder, endodermis, and cortex; the three genes displayed different expression patterns in cucumber roots with only partial overlap. A knockout of individual CsDRO1, CsDRO1L1, and CsDRO1L2 genes was performed via CRISPR/Cas9 gene editing. Our study suggests that the knockout of individual genes does not affect the slope angle formation during lateral root primordia development in the cucumber parental root.

The Rice ILI2 Locus Is a Bidirectional Target of the African Xanthomonas oryzae pv. oryzae Major Transcription Activator-like Effector TalC but Does Not Contribute to Disease Susceptibility

Xanthomonas oryzae pv. oryzae (Xoo) strains that cause bacterial leaf blight (BLB) limit rice (Oryza sativa) production and require breeding more resistant varieties. Transcription activator-like effectors (TALEs) activate transcription to promote leaf colonization by binding to specific plant host DNA sequences termed effector binding elements (EBEs). Xoo major TALEs universally target susceptibility genes of the SWEET transporter family. TALE-unresponsive alleles of clade III OsSWEET susceptibility gene promoter created with genome editing confer broad resistance on Asian Xoo strains. African Xoo strains rely primarily on the major TALE TalC, which targets OsSWEET14. Although the virulence of a talC mutant strain is severely impaired, abrogating OsSWEET14 induction with genome editing does not confer equivalent resistance on African Xoo. To address this contradiction, we postulated the existence of a TalC target susceptibility gene redundant with OsSWEET14. Bioinformatics analysis identified a rice locus named ATAC composed of the INCREASED LEAF INCLINATION 2 (ILI2) gene and a putative lncRNA that are shown to be bidirectionally upregulated in a TalC-dependent fashion. Gain-of-function approaches with designer TALEs inducing ATAC sequences did not complement the virulence of a Xoo strain defective for SWEET gene activation. While editing the TalC EBE at the ATAC loci compromised TalC-mediated induction, multiplex edited lines with mutations at the OsSWEET14 and ATAC loci remained essentially susceptible to African Xoo strains. Overall, this work indicates that ATAC is a probable TalC off-target locus but nonetheless documents the first example of divergent transcription activation by a native TALE during infection.

Comprehensive construction strategy of bidirectional green tissue-specific synthetic promoters

Bidirectional green tissue-specific promoters have important application prospects in genetic engineering and crop genetic improvement. However, there is no report on the application of them, mainly due to undiscovered natural bidirectional green tissue-specific promoters and the lack of a comprehensive approach for the synthesis of these promoters. In order to compensate for this vacancy, the present study reports a novel strategy for the expression regulatory sequence selection and the bidirectional green tissue-specific synthetic promoter construction. Based on this strategy, seven promoters were synthesized and introduced into rice by agrobacterium-mediated transformation. The functional identification of these synthetic promoters was performed by the expression pattern of GFP and GUS reporter genes in two reverse directions in transgenic rice. The results indicated that all the synthetic promoters possessed bidirectional expression activities in transgenic rice, and four synthetic promoters (BiGSSP2, BiGSSP3, BiGSSP6, BiGSSP7) showed highly bidirectional expression efficiencies specifically in green tissues (leaf, sheath, panicle, stem), which could be widely applied to agricultural biotechnology. Our study provided a feasible strategy for the construction of synthetic promoters, and we successfully created four bidirectional green tissue-specific synthetic promoters. It is the first report on bidirectional green tissue-specific promoters that could be efficiently applied in genetic engineering.

Functional characterization of an inducible bidirectional promoter from Fusarium oxysporum f. sp. cubense

Bidirectional promoters (BDPs) are regulatory DNA sequences (~1000 bp long) intervening two genes arranged on opposite strands with their 5' ends in close proximity. These genes are mostly co-expressed; but, instances of anti-correlation and independent transcription have been observed. In fungal systems, BDPs have shown to provide an improved genetic circuit by assembling and regulating transcription of different genes of a common metabolic pathway. We have identified an intergenic region (1063 bp) from the genome of Fusarium oxysporum f. sp. cubense (Foc), a banana root pathogen. This intergenic region regulates the expression of a gene pair required for the breakdown of hemicellulose. For characterization, it was cloned into pCSN44 vector backbone between two reporter genes, namely β-glucuronidase (GUS) and enhanced green fluorescent protein (EGFP). The newly formed vector was transformed into Foc and tested for its bidirectional expression activity. Using histochemical staining and fluorescence microscopy, the kinetics for both, GUS and EGFP expression were tested under different growth conditions respectively. The activity was differentially regulated by inducers such as xylan, arabinogalactan and pectin. This is the first report on the isolation of the intergenic region with inducible bidirectional promoter activity from Fusarium. Characterization of such BDPs will find applications in genetic engineering, metabolic engineering and synthetic biology using fungal systems.

Dynamic regulation of chromatin topology and transcription by inverted repeat-derived small RNAs in sunflower

Transposable elements (TEs) are extremely abundant in complex plant genomes. siRNAs of 24 nucleotides in length control transposon activity in a process that involves de novo methylation of targeted loci. Usually, these epigenetic modifications trigger nucleosome condensation and a permanent silencing of the affected loci. Here, we show that a TE-derived inverted repeat (IR) element, inserted near the sunflower HaWRKY6 locus, dynamically regulates the expression of the gene by altering chromatin topology. The transcripts of this IR element are processed into 24-nt siRNAs, triggering DNA methylation on its locus. These epigenetic marks stabilize the formation of tissue-specific loops in the chromatin. In leaves, an intragenic loop is formed, blocking HaWRKY6 transcription. While in cotyledons (Cots), formation of an alternative loop, encompassing the whole HaWRKY6 gene, enhances transcription of the gene. The formation of this loop changes the promoter directionality, reducing IR transcription, and ultimately releasing the loop. Our results provide evidence that TEs can act as active and dynamic regulatory elements within coding loci in a mechanism that combines RNA silencing, epigenetic modification, and chromatin remodeling machineries.

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.

From intronization to intron loss: How the interplay between mRNA-associated processes can shape the architecture and the expression of eukaryotic genes

Transcription-coupled processes such as capping, splicing, and cleavage/polyadenylation participate in the journey from genes to proteins. Although they are traditionally thought to serve only as steps in the generation of mature mRNAs, a synthesis of available data indicates that these processes could also act as a driving force for the evolution of eukaryotic genes. A theoretical framework for how mRNA-associated processes may shape gene structure and expression has recently been proposed. Factors that promote splicing and cleavage/polyadenylation in this framework compete for access to overlapping or neighboring signals throughout the transcription cycle. These antagonistic interactions allow mechanisms for intron gain and splice site recognition as well as common trends in eukaryotic gene structure and expression to be coherently integrated. Here, I extend this framework further. Observations that largely (but not exclusively) revolve around the formation of DNA-RNA hybrid structures, called R loops, and promoter directionality are integrated. Additionally, the interplay between splicing factors and cleavage/polyadenylation factors is theorized to also affect the formation of intragenic DNA double-stranded breaks thereby contributing to intron loss. The most notable prediction in this proposition is that RNA molecules can mediate intron loss by serving as a template to repair DNA double-stranded breaks. The framework presented here leverages a vast body of empirical observations, logically extending previous suggestions, and generating verifiable predictions to further substantiate the view that the intracellular environment plays an active role in shaping the structure and the expression of eukaryotic genes.

Processes Underlying a Reproductive Barrier in indica-japonica Rice Hybrids Revealed by Transcriptome Analysis

In rice (Oryza sativa), hybrids between indica and japonica subspecies are usually highly sterile, which provides a model system for studying postzygotic reproductive isolation. A killer-protector system, S5, composed of three adjacent genes (ORF3, ORF4, and ORF5), regulates female gamete fertility of indica-japonica hybrids. To characterize the processes underlying this system, we performed transcriptomic analyses of pistils from rice variety Balilla (BL), Balilla with transformed ORF5+ (BL5+) producing sterile female gametes, and Balilla with transformed ORF3+ and ORF5+ (BL3+5+) producing fertile gametes. RNA sequencing of tissues collected before (MMC), during (MEI), and after (AME) meiosis of the megaspore mother cell detected 19,269 to 20,928 genes as expressed. Comparison between BL5+ and BL showed that ORF5+ induced differential expression of 8,339, 6,278, and 530 genes at MMC, MEI, and AME, respectively. At MMC, large-scale differential expression of cell wall-modifying genes and biotic and abiotic response genes indicated that cell wall integrity damage induced severe biotic and abiotic stresses. The processes continued to MEI and induced endoplasmic reticulum (ER) stress as indicated by differential expression of ER stress-responsive genes, leading to programmed cell death at MEI and AME, resulting in abortive female gametes. In the BL3+5+/BL comparison, 3,986, 749, and 370 genes were differentially expressed at MMC, MEI, and AME, respectively. Large numbers of cell wall modification and biotic and abiotic response genes were also induced at MMC but largely suppressed at MEI without inducing ER stress and programed cell death , producing fertile gametes. These results have general implications for the understanding of biological processes underlying reproductive barriers.

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.

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.

A phosphate starvation-driven bidirectional promoter as a potential tool for crop improvement and in vitro plant biotechnology

Phosphate (Pi)-deficient soils are a major limitant factor for crop production in many regions of the world. Despite that plants have innovated several developmental and biochemical strategies to deal with this stress, there are still massive extensions of land which combine several abiotic stresses, including phosphate starvation, that limit their use for plant growth and food production. In several plant species, a genetic programme underlies the biochemical and developmental responses of the organism to cope with low phosphate (Pi) availability. Both protein- and miRNA-coding genes involved in the adaptative response are transcriptionally activated upon Pi starvation. Several of the responsive genes have been identified as transcriptional targets of PHR1, a transcription factor that binds a conserved cis-element called PHR1-binding site (P1BS). Our group has previously described and characterized a minimal genetic arrangement that includes two P1BS elements, as a phosphate-responsive enhancer (EZ2). Here, we report the engineering and successful use of a phosphate-dependent bidirectional promoter, which has been designed and constructed based on the palindromic sequences of the two P1BS elements present in EZ2. This bidirectional promoter has a potential use in both plant in vitro approaches and in the generation of improved crops adapted to Pi starvation and other abiotic stresses.

A phosphate starvation-driven bidirectional promoter as a potential tool for crop improvement and in vitro plant biotechnology

Phosphate (Pi)-deficient soils are a major limitant factor for crop production in many regions of the world. Despite that plants have innovated several developmental and biochemical strategies to deal with this stress, there are still massive extensions of land which combine several abiotic stresses, including phosphate starvation, that limit their use for plant growth and food production. In several plant species, a genetic programme underlies the biochemical and developmental responses of the organism to cope with low phosphate (Pi) availability. Both protein- and miRNA-coding genes involved in the adaptative response are transcriptionally activated upon Pi starvation. Several of the responsive genes have been identified as transcriptional targets of PHR1, a transcription factor that binds a conserved cis-element called PHR1-binding site (P1BS). Our group has previously described and characterized a minimal genetic arrangement that includes two P1BS elements, as a phosphate-responsive enhancer (EZ2). Here, we report the engineering and successful use of a phosphate-dependent bidirectional promoter, which has been designed and constructed based on the palindromic sequences of the two P1BS elements present in EZ2. This bidirectional promoter has a potential use in both plant in vitro approaches and in the generation of improved crops adapted to Pi starvation and other abiotic stresses.

Dissection of TALE-dependent gene activation reveals that they induce transcription cooperatively and in both orientations

Plant-pathogenic Xanthomonas bacteria inject transcription activator-like effector proteins (TALEs) into host cells to specifically induce transcription of plant genes and enhance susceptibility. Although the DNA-binding mode is well-understood it is still ambiguous how TALEs initiate transcription and whether additional promoter elements are needed to support this. To systematically dissect prerequisites for transcriptional initiation the activity of one TALE was compared on different synthetic Bs4 promoter fragments. In addition, a large collection of artificial TALEs spanning the OsSWEET14 promoter was compared. We show that the presence of a TALE alone is not sufficient to initiate transcription suggesting the requirement of additional supporting promoter elements. At the OsSWEET14 promoter TALEs can initiate transcription from various positions, in a synergistic manner of multiple TALEs binding in parallel to the promoter, and even by binding in reverse orientation. TALEs are known to shift the transcriptional start site, but our data show that this shift depends on the individual position of a TALE within a promoter context. Our results implicate that TALEs function like classical enhancer-binding proteins and initiate transcription in both orientations which has consequences for in planta target gene prediction and design of artificial activators.

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.

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.

Isolation and Functional Characterization of Bidirectional Promoters in Rice

Bidirectional promoters, which show great application potential in genetic improvement of plants, have aroused great research interest recently. However, most bidirectional promoters were cloned individually in the studies of single genes. Here, we initiatively combined RNA-seq data and cDNA microarray data to discover the potential bidirectional promoters in rice genome. Based on the expression level and correlation of each adjacent and oppositely transcribed gene pair, we selected four candidate gene pairs. Then, the intergenic region between each pair was isolated and cloned into a dual reporter vector pDX2181 for functional identification. GUS and GFP assays of the transgenic plants indicated that all the intergenic regions showed bidirectional expression activity in various tissues. Through 5' and 3' deletion analysis on one of the above bidirectional promoters, we identified the enhancing region which sharply increased its bidirectional expression efficiency and the essential regions respectively responsible for its 5' and 3' basic expression activity. The bidirectional arrangement of the four gene pairs in six gramineous plants was also analyzed, showing the conserved characteristics of the four bidirectional promoters identified in our study. In addition, two novel cis-sequences conserved in the four bidirectional promoters were discovered by bioinformatic identification. Our study proposes a feasible method for selecting, cloning, and functionally identifying bidirectional promoters as well as for discovering their bidirectional regulatory regions and conserved sequences in rice.

Genome-wide identification of CCA1 targets uncovers an expanded clock network in Arabidopsis

The circadian clock in Arabidopsis exerts a critical role in timing multiple biological processes and stress responses through the regulation of up to 80% of the transcriptome. As a key component of the clock, the Myb-like transcription factor CIRCADIAN CLOCK ASSOCIATED1 (CCA1) is able to initiate and set the phase of clock-controlled rhythms and has been shown to regulate gene expression by binding directly to the evening element (EE) motif found in target gene promoters. However, the precise molecular mechanisms underlying clock regulation of the rhythmic transcriptome, specifically how clock components connect to clock output pathways, is poorly understood. In this study, using ChIP followed by deep sequencing of CCA1 in constant light (LL) and diel (LD) conditions, more than 1,000 genomic regions occupied by CCA1 were identified. CCA1 targets are enriched for a myriad of biological processes and stress responses, providing direct links to clock-controlled pathways and suggesting that CCA1 plays an important role in regulating a large subset of the rhythmic transcriptome. Although many of these target genes are evening expressed and contain the EE motif, a significant subset is morning phased and enriched for previously unrecognized motifs associated with CCA1 function. Furthermore, this work revealed several CCA1 targets that do not cycle in either LL or LD conditions. Together, our results emphasize an expanded role for the clock in regulating a diverse category of genes and key pathways in Arabidopsis and provide a comprehensive resource for future functional studies.

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.

Analysis of developmentally regulated chorion gene promoter architecture via electroporation of silk moth follicles

In the silk moth Bombyx mori, chorion genes of the same developmental specificity are organized in divergently transcribed α/β gene pairs, sharing a common 5' flanking promoter region. This bidirectional promoter contains a complete set of cis-elements responsible for developmentally accurate gene expression. In the present paper, based on the observation that Bombyx chorion gene promoters contain cis-elements for the same transcription factors without concrete evidence on which of them are essential, we address the question as to how promoter architecture (number, orientation and position of common factor binding sites) facilitates developmentally accurate chorion gene regulation. To this end, we constructed several mutated promoter regions of an early-middle gene pair and cloned them upstream of a reporter gene to introduce these plasmid constructs into silk moth follicle epithelial cells via electroporation as an efficient and quick method for transient expression. This is the first time that an ex vivo method had been applied to test the impact of systematic cis-element mutations on a chorion gene promoter. Our results confirmed the importance of the HMGA factor and the role of the GATA factor as an early repressor, and led to a more detailed understanding of which C/EBP sites participate in the regulation of early-middle chorion gene expression.

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.

Characterization and identification of cis-regulatory elements in Arabidopsis based on single-nucleotide polymorphism information

Identifying regulatory elements and revealing their role in gene expression regulation remains a central goal of plant genome research. We exploited the detailed genomic sequencing information of a large number of Arabidopsis (Arabidopsis thaliana) accessions to characterize known and to identify novel cis-regulatory elements in gene promoter regions of Arabidopsis by relying on conservation as the hallmark signal of functional relevance. Based on the genomic layout and the obtained density profiles of single-nucleotide polymorphisms (SNPs) in sequence regions upstream of transcription start sites, the average length of promoter regions in Arabidopsis could be established at 500 bp. Genes associated with high degrees of variability of their respective upstream regions are preferentially involved in environmental response and signaling processes, while low levels of promoter SNP density are common among housekeeping genes. Known cis-elements were found to exhibit a decreased SNP density than sequence regions not associated with known motifs. For 15 known cis-element motifs, strong positional preferences relative to the transcription start site were detected based on their promoter SNP density profiles. Five novel candidate cis-element motifs were identified as consensus motifs of 17 sequence hexamers exhibiting increased sequence conservation combined with evidence of positional preferences, annotation information, and functional relevance for inducing correlated gene expression. Our study demonstrates that the currently available resolution of SNP data offers novel ways for the identification of functional genomic elements and the characterization of gene promoter sequences.

Bidirectional promoters in seed development and related hormone/stress responses

BACKGROUND

Bidirectional promoters are common in genomes but under-studied experimentally, particularly in plants. We describe a targeted identification and selection of a subset of putative bidirectional promoters to identify genes involved in seed development and to investigate possible coordinated responses of gene pairs to conditions important in seed maturation such as desiccation and ABA-regulation.

RESULTS

We combined a search for 100-600 bp intergenic regions in the Arabidopsis genome with a cis-element based selection for those containing multiple copies of the G-box motif, CACGTG. One of the putative bidirectional promoters identified also contained a CE3 coupling element 5 bp downstream of one G-box and is identical to that characterized previously in the HVA1 promoter of barley. CE3 elements are significantly under-represented and under-studied in Arabidopsis. We further characterized the pair of genes associated with this promoter and uncovered roles for two small, previously uncharacterized, plant-specific proteins in Arabidopsis seed development and stress responses.

CONCLUSIONS

Using bioinformatics we identified putative bidirectional promoters involved in seed development and analysed expression patterns for a pair of plant-specific genes in various tissues and in response to hormones/stress. We also present preliminary functional analysis of these genes that is suggestive of roles in seed development.

An intergenic region shared by At4g35985 and At4g35987 in Arabidopsis thaliana is a tissue specific and stress inducible bidirectional promoter analyzed in transgenic arabidopsis and tobacco plants

On chromosome 4 in the Arabidopsis genome, two neighboring genes (calmodulin methyl transferase At4g35987 and senescence associated gene At4g35985) are located in a head-to-head divergent orientation sharing a putative bidirectional promoter. This 1258 bp intergenic region contains a number of environmental stress responsive and tissue specific cis-regulatory elements. Transcript analysis of At4g35985 and At4g35987 genes by quantitative real time PCR showed tissue specific and stress inducible expression profiles. We tested the bidirectional promoter-function of the intergenic region shared by the divergent genes At4g35985 and At4g35987 using two reporter genes (GFP and GUS) in both orientations in transient tobacco protoplast and Agro-infiltration assays, as well as in stably transformed transgenic Arabidopsis and tobacco plants. In transient assays with GFP and GUS reporter genes the At4g35985 promoter (P85) showed stronger expression (about 3.5 fold) compared to the At4g35987 promoter (P87). The tissue specific as well as stress responsive functional nature of the bidirectional promoter was evaluated in independent transgenic Arabidopsis and tobacco lines. Expression of P85 activity was detected in the midrib of leaves, leaf trichomes, apical meristemic regions, throughout the root, lateral roots and flowers. The expression of P87 was observed in leaf-tip, hydathodes, apical meristem, root tips, emerging lateral root tips, root stele region and in floral tissues. The bidirectional promoter in both orientations shows differential up-regulation (2.5 to 3 fold) under salt stress. Use of such regulatory elements of bidirectional promoters showing spatial and stress inducible promoter-functions in heterologous system might be an important tool for plant biotechnology and gene stacking applications.

Major cis-regulatory elements for rice bidirectional promoter activity reside in the 5'-untranslated regions

Bidirectional promoters are defined as those that regulate adjacent genes organized in a divergent fashion (head to head orientation) and separated by <1 kb. In order to dissect bidirectional promoter activity in a model plant, deletion analysis was performed for seven rice promoters using promoter-reporter gene constructs, which identified three promoters to be bidirectional. Regulatory elements located in or close to the 5'-untranslated regions (UTR) of one of the genes (divergent gene pair) were found to be responsible for their bidirectional activity. DNA footprinting analysis identified unique protein binding sites in these promoters. Deletion/alteration of these motifs resulted in significant loss of expression of the reporter genes on either side of the promoter. Changes in the motifs at both the positions resulted in a remarkable decrease in bidirectional activity of the reporter genes flanking the promoter. Based on our results, we propose a novel mechanism for the bidirectionality of rice bidirectional promoters.

Characterization of a T-DNA promoter trap line of Arabidopsis thaliana uncovers a cryptic bi-directional promoter

Investigation of the transgenic Arabidopsis promoter trap line GFP-868 that showed GFP expression only in anthers revealed the T-DNA insertion at 461bp upstream to the hypothetical gene At4g10596 with the GFP reporter gene in head-to-head orientation to the At4g10596 gene. The expression of the At4g10596 gene in wild type and in GFP-868 plant homozygous for T-DNA insertion was comparable and found in all tissues tested, while the GFP expression was restricted to anthers of the GFP-868 plants suggesting that the 461bp fragment separating the two genes in the GFP-868 line is functioning as bi-directional promoter. This 461bp fragment was cloned upstream to the GUS gene in two orientations to test for bi-directional promoter activity. Transgenic Arabidopsis plants carrying either of these constructs showed GUS activity in anthers indicating that this fragment behaves as bi-directional promoter specific to anthers. These results were also supported by the presence of cis-acting motifs such as TATA box and POLLEN1LELAT52 (AGAAA) within the 461bp sequence in both orientations. However, transcripts corresponding to the upstream sequences beyond -461 nucleotides were not detected in the wild type suggesting that this 461bp fragment is a cryptic promoter. The significance of the promoter trap approach and the usefulness of this type of promoter are discussed.

Development of an efficient bi-directional promoter with tripartite enhancer employing three viral promoters

We have developed a novel bi-directional promoter (FsFfCBD) by placing two heterogeneous core-promoters from the Figwort mosaic virus sub-genomic transcript promoter (FsCP, -69 to +31) and Cauliflower mosaic virus 35S promoter (CCP, -89 to +1) respectively on upstream (5') and downstream (3') ends of a tri-hybrid enhancer (FsEFfECE), in reverse orientation. The FsEFfECE domain encompasses three heterologous enhancer fragments from Figwort mosaic virus sub-genomic transcript promoter (FsE, 101 bp, -70 to -170), Figwort mosaic virus full-length transcript promoter (FfE, 196 bp, -249 to -54) and Cauliflower mosaic virus 35S promoter (CE, 254 bp, -343 to -90). The bi-directional nature of the FsFfCBD promoter (coupled to GFP and GUS) was established both in transient systems (onion epidermal cells and tobacco protoplasts) and transgenic plant (Nicotiana tabacum samsun NN) by monitoring the simultaneous expression of GFP and GUS employing fluorescence (for GFP) and biochemical (for GUS) based assays. In transgenic plants, the FsFfCBD promoter was found to be 6.8 and 2.5 times stronger than two parent promoters; Fs and FfC respectively. The bi-directional compound promoter FsFfCBD, composed of three heterologous enhancers with enhanced activity could become a valuable additional tool for efficient plant metabolic engineering and molecular pharming.

The preservation of bidirectional promoter architecture in eukaryotes: what is the driving force?

Background

The bidirectional gene architecture has been studied in many organisms, and the conservation of bidirectional arrangement has received considerable attention. However, the explanation for the evolutionary conservation about this genomic structure is still insufficient. In this study the large scale identification and pathway enrichment analysis for bidirectional genes were performed in several eukaryotes and the comparative analysis of this arrangement between human and mouse were dissected for the purpose of discovering the driving force of the preservation of this genomic structure.

Results

We identified the bidirectional gene pairs in eight different species and found this structure to be prevalent in eukaryotes. The pathway enrichment analysis indicated the bidirectional genes at the genome level are conserved in certain pathways, such as the DNA repair and some other fundamental cellular pathways. The comparative analysis about the gene expression, function, between human and mouse bidirectional genes were also performed and we observed that the selective force of this architecture doesn't derive from the co-regulation between paired genes, but the functional bias of bidirectional genes at whole genome level is observed strengthened during evolution.

Conclusions

Our result validated the coexpression of bidirectional genes; however failed to support their functional relevance. The conservation of bidirectional promoters seems not the result of functional connection between paired genes, but the functional bias at whole genome level, which imply that the genome-wide functional constraint is important for the conservation of bidirectional structure.

Modeling bidirectional transcription using silkmoth chorion gene promoters

Bidirectional transcription is an interesting feature of eukaryotic genomes; yet not all aspects of its mechanism are understood. Silkmoth choriogenesis is a model system for studying transcriptional regulation at the initiation level. As chorion genes comprise a large group of divergently transcribed gene pairs, we are presented with the possibility of investigating the intricacies of bidirectional transcription. Their well characterized 5' regulatory regions and expression profiles lay the foundation for investigating protein:protein and protein:DNA interactions, and RNA polymerase function during oocyte development. In this article we summarize current knowledge on chorion gene regulation and propose an approach to modeling bidirectional transcription using chorion promoters.

A study on promoter characteristics of head-to-head genes in Saccharomyces cerevisiae

Background

Head-to-head (h2h) genes are prone to have association in expression and in functionality and have been shown conserved in evolution. Currently there are many studies on such h2h gene pairs. We found that the previous studies extremely focused on human genome. Furthermore, they only focused on analyses that require only gene or protein sequences but not conducted a systematic investigation on other promoter features such as the binding evidence of specific transcription factors (TFs). This is mainly because of the incomplete resources of higher organisms, though they are relatively of interest, than model organisms such as Saccharomyces cerevisiae. The authors of this study recently integrated nine promoter features of 6603 genes of S. cerevisiae from six databases and five papers. These resources are suitable to conduct a comprehensive analysis of h2h genes in S. cerevisiae.

Results

This study analyzed various promoter features, including transcription boundaries (TSS, 5'UTR and 3'UTR), TATA box, TF binding evidence, TF regulation evidence, DNA bendability and nucleosome occupancy. The expression profiles and gene ontology (GO) annotations were used to measure if two genes are associated. Based on these promoter features, we found that i) the frequency of h2h genes was close to the expectation, namely they were not relatively frequent in genome; ii) the distance between the TSSs of most h2h genes fell into the range of 0-600 bps and was more centralized in 0-200 bps of the highly associated ones; iii) the number of TFs that regulate both h2h genes influenced the co-expression and co-function of the genes, while the number of TFs that bind both h2h genes influenced only the co-expression of the genes; iv) the association of two h2h genes was influenced by the existence of specific TFs such as STP2; v) the association of h2h genes whose bidirectional promoters have no TATA box was slightly higher than those who have TATA boxes; vi) the association of two h2h genes was not influenced by the DNA bendability and nucleosome occupancy.

Conclusions

This study analyzed h2h genes with various promoter features that have not been used in analyzing h2h genes. The results can be applied to other genomes to confirm if the observations of this study are limited to S. cerevisiae or universal in most organisms.

A bidirectional promoter reporter vector for the analysis of the p53/WDR79 dual regulatory element

Analysis of numerous genomes has identified a class of regulatory regions that contain a head-to-head arrangement (5' to 5') on opposite strands of DNA. Often these regulatory regions have fewer than 1000 base pairs separating their corresponding transcription start sites and have been termed as being "bidirectional". This bidirectional arrangement and the divergent gene pairs under the control of these regulatory regions appear to be a common feature within genomes. Establishing methods to study these bidirectional transcriptional promoters, and understanding how they are regulated will allow researchers to gain more insight into the roles that divergent transcription plays in the expression and maintenance of protein coding genes. Recently, the p53 tumor suppressor gene was shown to have a bidirectional gene partner, WDR79. The transcription start sites (TSSs) of human and murine p53 and WDR79 genes are separated by approximately 800 and 930bp, respectively, in a head-to-head fashion, and fit the criteria of what is designated to be a putative bidirectional regulatory region. However, further testing is needed to demonstrate that the region between these genes contains a functional bidirectional promoter. Here, we have developed a bidirectional reporter vector, termed pLucRLuc, to study the transcriptional output of each promoter. This bidirectional reporter vector will allow researchers to determine the output of transcripts mediated by the bidirectional promoters. By focusing our studies on the transcriptional regulation of p53 and its bidirectional gene partner, WDR79, we hope to elucidate key factors that can control and regulate the expression of the p53 and WDR79 genes. Here, we demonstrate that pLucRLuc is a vector capable of expressing reporter genes under the control of bidirectional promoters in multiple human and murine cell lines and that the regulatory region upstream of the p53 and WDR79 TSSs is a bidirectional promoter controlled by common regulatory factors.

Prediction of plant promoters based on hexamers and random triplet pair analysis

Background

With an increasing number of plant genome sequences, it has become important to develop a robust computational method for detecting plant promoters. Although a wide variety of programs are currently available, prediction accuracy of these still requires further improvement. The limitations of these methods can be addressed by selecting appropriate features for distinguishing promoters and non-promoters.

Methods

In this study, we proposed two feature selection approaches based on hexamer sequences: the Frequency Distribution Analyzed Feature Selection Algorithm (FDAFSA) and the Random Triplet Pair Feature Selecting Genetic Algorithm (RTPFSGA). In FDAFSA, adjacent triplet-pairs (hexamer sequences) were selected based on the difference in the frequency of hexamers between promoters and non-promoters. In RTPFSGA, random triplet-pairs (RTPs) were selected by exploiting a genetic algorithm that distinguishes frequencies of non-adjacent triplet pairs between promoters and non-promoters. Then, a support vector machine (SVM), a nonlinear machine-learning algorithm, was used to classify promoters and non-promoters by combining these two feature selection approaches. We referred to this novel algorithm as PromoBot.

Results

Promoter sequences were collected from the PlantProm database. Non-promoter sequences were collected from plant mRNA, rRNA, and tRNA of PlantGDB and plant miRNA of miRBase. Then, in order to validate the proposed algorithm, we applied a 5-fold cross validation test. Training data sets were used to select features based on FDAFSA and RTPFSGA, and these features were used to train the SVM. We achieved 89% sensitivity and 86% specificity.

Conclusions

We compared our PromoBot algorithm to five other algorithms. It was found that the sensitivity and specificity of PromoBot performed well (or even better) with the algorithms tested. These results show that the two proposed feature selection methods based on hexamer frequencies and random triplet-pair could be successfully incorporated into a supervised machine learning method in promoter classification problem. As such, we expect that PromoBot can be used to help identify new plant promoters. Source codes and analysis results of this work could be provided upon request.

Support for multiple classes of local expression clusters in Drosophila melanogaster, but no evidence for gene order conservation

BACKGROUND

Gene order in eukaryotic genomes is not random, with genes with similar expression profiles tending to cluster. In yeasts, the model taxon for gene order analysis, such syntenic clusters of non-homologous genes tend to be conserved over evolutionary time. Whether similar clusters show gene order conservation in other lineages is, however, undecided. Here, we examine this issue in Drosophila melanogaster using high-resolution chromosome rearrangement data.

RESULTS

We show that D. melanogaster has at least three classes of expression clusters: first, as observed in mammals, large clusters of functionally unrelated housekeeping genes; second, small clusters of functionally related highly co-expressed genes; and finally, as previously defined by Spellman and Rubin, larger domains of co-expressed but functionally unrelated genes. The latter are, however, not independent of the small co-expression clusters and likely reflect a methodological artifact. While the small co-expression and housekeeping/essential gene clusters resemble those observed in yeast, in contrast to yeast, we see no evidence that any of the three cluster types are preserved as synteny blocks. If anything, adjacent co-expressed genes are more likely to become rearranged than expected. Again in contrast to yeast, in D. melanogaster, gene pairs with short intergene distance or in divergent orientations tend to have higher rearrangement rates. These findings are consistent with co-expression being partly due to shared chromatin environment.

CONCLUSIONS

We conclude that, while similar in terms of cluster types, gene order evolution has strikingly different patterns in yeasts and in D. melanogaster, although recombination is associated with gene order rearrangement in both.

Gene clustering pattern, promoter architecture, and gene expression stability in eukaryotic genomes

A balance between gene expression stability and evolvability is essential for the long-term maintenance of a living system. In this paper, we studied whether the genetic and epigenetic properties of the promoter affect gene expression variability. We hypothesized that upstream distance and orientation (head-to-head or head-to-tail) are important for the promoter architecture and gene expression variability. We found that in budding yeast genes with a short upstream distance tend to have low gene expression variability, and their promoter is flanked by strongly positioned nucleosomes and tends to have low nucleosome occupancy. These observations suggest that in vivo positioning of the flanking nucleosomes facilitates stable nucleosome depletion at the core promoter region and enhances gene expression stability. Head-to-head genes have, on average, lower gene expression variability, greater nucleosome depletion at the core promoter region, and more strongly positioned nucleosomes that flank the core promoter than do head-to-tail genes. These observations hold for diverse eukaryotes. In complex organisms such as mammals, only a small fraction of head-to-tail genes have retained a short upstream distance, probably because the promoter may not be flanked by a strongly positioned nucleosome on the upstream side.

Flexible metabolism in Metarhizium anisopliae and Beauveria bassiana: role of the glyoxylate cycle during insect pathogenesis

Insect pathogenic fungi such as Metarhizium anisopliae and Beauveria bassiana have an increasing role in the control of agricultural insect pests and vectors of human diseases. Many of the virulence factors are well studied but less is known of the metabolism of these fungi during the course of insect infection or saprobic growth. Here, we assessed enzyme activity and gene expression in the central carbon metabolic pathway, including isocitrate dehydrogenase, aconitase, citrate synthase, malate synthase (MLS) and isocitrate lyase (ICL), with particular attention to the glyoxylate cycle when M. anisopliae and B. bassiana were grown under various conditions. We observed that ICL and MLS, glyoxylate cycle intermediates, were upregulated during growth on 2-carbon compounds (acetate and ethanol) as well as in insect haemolymph. We fused the promoter of the M. anisopliae ICL gene (Ma-icl) to a marker gene (mCherry) and showed that Ma-icl was upregulated when M. anisopliae was grown in the presence of acetate. Furthermore, Ma-icl was upregulated when fungi were engulfed by insect haemocytes as well as during appressorium formation. Addition of the ICL inhibitor 3-nitroproprionate delayed conidial germination and inhibited appressorium formation. These results show that these insect pathogenic fungi have a flexible metabolism that includes the glyoxylate cycle as an integral part of germination, pathogenesis and saprobic growth.

Co-expression of neighbouring genes in Arabidopsis: separating chromatin effects from direct interactions

Background

In all eukaryotic species examined, genes that are chromosomal neighbours are more similar in their expression than random gene pairs. Currently, it is still unclear how much of this local co-expression is caused by direct transcriptional interactions, and how much is due to shared chromatin environments.

Results

We analysed neighbouring genes in Arabidopsis thaliana. At large intergenic distances (>400 bp), divergently and convergently transcribed gene pairs show very similar levels of co-expression, mediated most likely by shared chromatin environments. At gene distances below 400 bp, co-expression is strongly enhanced only for divergently transcribed gene pairs, indicating bi-directional transcription from a single promoter. Conversely, co-expression is suppressed for short convergently or uni-directionally transcribed pairs. This suppression points to transcriptional interference concentrated at the 3' end, e.g., in the context of transcription termination.

Conclusions

Classifying linked gene pairs by their orientation, we are able to partially tease apart the different levels of regional expression modulation. (i) Regional chromatin characteristics modulate the accessibility for regulation and transcription, regardless of gene orientation; the strength of this chromatin effect can be assessed from divergently or convergently transcribed distant neighbours. (ii) Shared promoter regions up to 400 bp in length enhance the co-expression of close bi-directional neighbours. (iii) Transcriptional interference of close neighbours is concentrated at the 3' ends of genes, and reduces co-expression on average by 40%.

Assessing the influence of adjacent gene orientation on the evolution of gene upstream regions in Arabidopsis thaliana

The orientation of flanking genes may influence the evolution of intergenic regions in which cis-regulatory elements are likely to be located: divergently transcribed genes share their 5' regions, resulting either in smaller "private" spaces or in overlapping regulatory elements. Thus, upstream sequences of divergently transcribed genes (bi-directional upstream regions, or URs) may be more constrained than those of uni-directional gene pairs. We investigated this effect by analyzing nucleotide variation segregating within and between Arabidopsis species. Compared to uni-directional URs, bi-directional URs indeed display lower population mutation rate, as well as more low-frequency polymorphisms. Furthermore, we find that bi-directional regions undergo selection for the maintenance of intergenic distance. Altogether, however, we observe considerable variation in evolutionary rates, with putative signatures of selection on two uni-directional upstream regions.