Co-expression of neighbouring genes in Arabidopsis: separating chromatin effects from direct interactions
© Chen et al; licensee BioMed Central Ltd. 2010
Received: 25 November 2009
Accepted: 16 March 2010
Published: 16 March 2010
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.
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.
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%.
Eukaryotic gene order is not random with respect to gene functions or expression patterns: clusters of co-expressed genes are routinely observed in a wide range of species . For example, 'housekeeping' genes expressed in many tissues tend to cluster in humans ; directly neighbouring gene pairs are more likely to act in the metabolic pathway ; and neighbouring genes that are transcribed divergently (←→) are unusually frequent in the human genome and show particularly strong co-expression . Co-expression is not restricted to close neighbours; in yeast, significant co-expression is observed over distances covering dozens of genes .
Several mechanisms have been proposed to explain this local co-expression of neighbouring gene pairs, including shared promoters or transcription factor binding sites, transcriptional read-through, and chromatin remodelling [1, 6]. In the yeast genome, for example, much local co-expression may stem from chromatin domains that as a whole switch between euchromatin and heterochromatin states . It is not clear to what extent this observation can be extrapolated to other model species: in mammalian genomes, there appears to be no corresponding compartmentalization of tissue-specific genes into co-regulated regions . Thus, the relative contributions of different local co-expression mechanisms seem to vary across species . This may not be surprising: fine details of gene regulation are under strong, species-specific selective pressures, and are affected by multiple internal and external cellular factors [8, 9].
Due to the compactness of its genome , the thale cress Arabidopsis thaliana is a multi-cellular eukaryote particularly well suited to the study of local co-expression. Short intergenic regions mean that neighbouring genes are physically very close to each other. As all proposed mechanisms for local co-expression are likely to act in a distance-dependent way , this should enhance their signal compared to genomes with larger and more variable gene distances. Several previous studies have examined local co-expression in A. thaliana.
Similar to previous results in fungi and animals [reviewed in ], neighbouring genes in A. thaliana are on average more co-expressed than expected by chance . 5-10% of A. thaliana genes are located in highly co-expressed clusters of 2-4 neighbouring genes [12, 13]; however, these are mostly not in the same Gene Ontology (GO) categories, and there is no detectable micro-syntenty of co-expressed clusters to corresponding clusters in rice . Evolutionary conservation of gene neighbourhood was stronger for highly co-expressed gene pairs, as well as for gene pairs sharing at least one GO category. The genes in 60% of metabolic pathways show significant clustering on A. thaliana chromosomes, a feature also found in other eukaryotes [3, 11]. While co-expresion of tandemly duplicated genes provides some of the signal, neighbour co-expression remains significant even when excluding duplicated genes [3, 11, 13].
Several previous studies have examined the role of relative gene orientation in determining co-expression levels of neighbouring genes, and have found that divergent transcription can account for some but not all of the co-expression effect [11–13, 15]. Divergent transcription of A. thaliana is preferentially associated with certain GO categories , and co-expression is increased between genes of related functions . Thus, while the possible sharing of promoter regions explains part of the co-expression pattern, it is highly likely that chromatin effects also play an important role in A. thaliana. This is supported by the finding that duplicated genes which are parts of larger duplications retain strongly correlated expression patterns, as long as their neighbourhood is conserved .
In this study, we dissect co-expression patterns of neighbouring gene pairs in the model plant species A. thaliana, paying particular attention to the relative orientation of the genes and to intergenic distances. We find the strongest co-expression among uni-directional gene pairs; this contrasts with results in humans and yeast, where bi-directional gene pairs are the most strongly co-expressed. This contradiction is resolved when examining only closely spaced gene pairs (<400 bp), where bi-directional pairs are the most strongly co-expressed also in A. thaliana. In line with expectations, this indicates that sharing of regulatory elements enhances co-expression. The opposite effect is seen for shared 3' regions in convergent gene pairs: these show decreased co-expression levels. The latter observation appears consistent with a model of transcription termination in which the polymerase overshoots the polyadenylation addition site [18–20], and may hence interfere with the transcription of closely positioned downstream genes.
Results and discussion
Co-expression extends over large distances
Why does co-expression drop so fast? Consider three consecutive genes A-B-C. If both A-B and B-C are correlated in their expression, shouldn't it follow that A-C are correlated similarly? Assume that the corresponding expression patterns (a, b, c) are correlated such that variation in a explains a fraction R ab 2 of variation in b, and variation in b explains a fraction R bc 2 of variation in c. If there is no direct interaction of a and c, we would expect a correlation R ac 2 = R ab 2 × R bc 2 between a and c. If we insert the numbers read from Figure 1, we thus expect a co-expression value of approx. 0.01 for pairs with 1 intervening gene, which is actually substantially lower than the observed value. This indicates that there are indeed local factors that enhance co-expression for pairs that are not direct neighbours.
Different orientations show different co-expression levels
By far the strongest co-expression is seen for direct neighbours, on which we focused next. To ensure a priori independence of the genes, we removed all genes that overlap with other (protein coding or non-coding) genes in the genome according to the TAIR7 annotation [21, 22]. The remaining 24,638 pairs of non-overlapping protein coding gene pairs were further classified according to their relative orientations, resulting in 6,077 bi-directional (or divergent, ←→), 13,574 uni-directional (→→ or ←←), and 4,987 convergent pairs (→←).
Co-expression of directly neighbouring genes
P** (compared to random pairs)
Among the three relative orientations, uni-directional pairs showed the highest average expression correlation, as was also observed previously . This contrasts with findings in the human genome, where bi-directional pairs are more highly co-expressed . Due to the vast evolutionary distance and different ecology, direct comparisons are of course difficult; however, one possible reason for this difference between humans and plants is the existence of bicistronic and fused monocistronic transcripts in the A. thaliana genome . Bicistronic transcripts contain two separate open reading frames (ORFs), leading to perfect co-expression of the ORFs; fused monocistronic transcripts arise by the fusion of two adjacent transcripts in particular tissues or environments, leading to elevated co-expression. To test if the higher average co-expression of uni-directional pairs might be attributed to bicistronic and fused monocistronic transcripts, we used a previously determined set of such genes (58 bicistronic and 30 fused monocistronic transcripts ). We found an average co-expression of 0.108, only slightly higher than the value obtained for all uni-directional transcripts. After excluding these bicistronic and fused monocistronic transcripts, uni-directional gene pairs are still more co-expressed than the two alternative relative orientations (P = 3.2 × 10-10 and 0.016 for convergent and bi-directional pairs respectively, Wilcoxon rank sum tests).
Long intergene distances do not reduce co-expression
Co-expression of gene pairs with short and long intergenic distances
Mean co-expression *
< 400 bp
> 400 bp
Direct interactions of the polymerases transcribing two neighbouring genes will be very different for convergent and divergent pairs. However, at distances >400 bp we find no difference between the co-expression of convergently and divergently transcribed gene pairs. This observation is strong evidence that not direct interactions, but shared chromatin regions are in fact responsible for co-expression at these larger distances. Thus, we conclude that co-expression of distant genes (>400 bp) is mostly due to chromatin effects. The cutoff 400 bp was chosen here based on Figure 2; a lower threshold of 300 bp leads to qualitatively very similar results.
Short bi-directional gene pairs are strongly co-expressed
As seen from Figure 2, the three orientations exhibit very different co-expression patterns at short distances (<400 bp). Bi-directional gene pairs closer than 400 bp are more co-expressed than either more distant bi-directional pairs (P = 1.3 × 10-7) or close neighbours in other orientations (P = 6.8 × 10-12 and P = 3.1 × 10-5 for convergent and uni-directional gene pairs, respectively). This is consistent with the existence of shared cis-regulatory motifs or promoters in some overlapping promoter regions. At least in yeast [25, 26] and mammals [27, 28], approximately half of all promoters initiate transcription in both directions, and it is conceivable that the same applies in Arabidopsis.
Transcriptional interference between short uni-directional and convergent gene pairs
As seen in Figure 2, short convergent as well as short uni-directional gene pairs show reduced co-expression levels compared to more distant pairs of the same orientation (P = 0.044 and P = 0.0012, respectively). This reduced co-expression suggests some kind of transcriptional interference between the neighbouring genes.
What could be responsible? That this reduced co-expression is found only when at least one of the two genes flankes the intergenic region with its 3' end suggests an involvement of transcription termination. Indeed, transcription by RNA-polymerase II  (and probably also polymerase I ) usually extends beyond the poly(A) site. In the 'torpedo' model of transcription termination , the mRNA is released by cleavage while RNA is still being synthesised. The remaining RNA trails out of the active polymerase, and is chewed up by an exonuclease. Once the 'exonuclease-torpedo' catches up, it dislodges the polymerase, thereby terminating transcription. Thus, the polymerase of an upstream gene may transcribe into the downstream neighbour, interfering with the neighbour's transcription (e.g., by blocking the processing of transcription termination, or by direct collision with another polymerase transcribing the downstream gene).
Overlapping gene pairs on opposing strands
P(compared to random pairs)c
Overlapping gene pairs on different strands (groups 1-3) actually form sense-antisense (SA) pairs. Thus, the lower co-expression of these genes could potentially be explained by RNA interference, which would reduce both the expression abundance and the co-expression of SA pairs. As shown in Table 3, partially overlapping genes are in fact more abundant (group 1) or equally abundant (group 2) as randomly picked genes. This indicates that at least for partially overlapping genes, expression is not strongly influenced by RNA interference (This result remained when we compared overlapping genes to subsets of the full genome with similar gene length, exon number or GO functional category distributions; results not shown).
Direct neighbours are more often co-expressed than are more distant genes. Beyond approximately 400 bp, the level of co-expression is independent of the intergenic distance, and is the same for divergent and convergent pairs. These results suggest that beyond 400 bp, co-expression is due to shared chromatin domains rather than to direct transcriptional interactions. This chromatin-mediated effect leads to an increase in the expression correlation coefficient of about 0.04 across our set of 1,436 microarray experiments (averaging over the co-expression of convergent and divergent pairs >400 bp apart, and subtracting the co-expression expected for random pairs; uni-directional pairs were excluded here as they are potentially affected by transcriptional read-through even at larger distances).
Using the strength of this effect as a reference, we can get an approximate estimate of the contributions of chromatin-mediated effects and direct interactions. We estimate that for close bi-directional pairs, shared promoter regions are on average responsible for roughly 70% of the co-expression signal ((0.15 - 0.04)/0.15), while the shared chromatin environment accounts for the remaining 30%. Further, for close convergent gene pairs, the co-expression mediated by a shared chromatin environment is reduced by roughly 40% through transcriptional interference ((0.040 - 0.024)/0.040).
As long as there are no intervening genes, chromatin-mediated co-expression does not seem to diminish with increasing intergenic distance. This observation suggests that chromatin domain establishment is regulated such that domains extend as far as needed, regardless of intergenic distance. This conclusion appears consistent with recent results on the role of nucleosome organization and 3-D chromatin structure in gene regulation .
At short distances, direct transcriptional interactions become important. For divergent gene pairs, the most important direct interactions are probably the sharing of transcriptional regulatory sites, and bi-directional initiation of RNA-polymerase II transcription from a shared promoter; such bi-directional promoters appear to be the rule rather than an exception in yeast and mammals . That we found elevated co-expression of bi-directional gene pairs only below ~400 bp strongly supports this notion, as this range is similar to the distance found between the two divergent peaks of transcription in other eukaryotes . Beyond promoting co-expression, bi-directional transcription from a single promoter may also serve as a positive feedback loop to support continued expression.
A rather different direct transcriptional interaction is likely responsible for the reduced co-expression found for those pairs where the intergenic region is flanked by at least one 3' gene end (uni-directional and convergent pairs). In transcription termination, the polymerase overshoots the poly(A)-site before being dislodged by a following exonuclease [19, 20]; this may result in continued transcription through the intergenic region and into the neighbouring gene, interfering with the neighbour's transcription.
In contrast to previous results in humans and yeast [4, 7], we found that uni-directional gene pairs were on average more co-expressed than bi-directional pairs. This is presumably due to the more compact spacing of A. thaliana genes, with a higher fraction of closely-spaced uni-directional gene pairs. When comparing the co-expression of the different orientations at fixed intergenic distances, our results agree with those in other eukaryotes.
We used release 7 of the A. thaliana genome annotation from TAIR [21, 22]. Neighbouring gene pairs were defined as protein-coding genes that are direct chromosomal neighbours, with no intervening genes (coding or non-coding). We restricted our dataset to genes that do not overlap with other genes. We also compiled a separate set of overlapping protein-coding gene pairs on opposing strands. We identified putative tandem duplicate genes based on the protocol developed in Ref. : performing pairwise blast searches between all neighbours, we removed those pairs with E-value < 0.2 from further analysis. The conservative cutoff 0.2 has been shown previously to lead to good sensitivity and specificity in the identification of even ancient gene duplicates .
We downloaded pre-processed expression data from the TAIR database ftp://ftp.arabidopsis.org/home/tair/Microarrays/analyzed_data/affy_data_1436_10132005.zip. The dataset contains 1,436 hybridization experiments using the Affymetrix A. thaliana ATH1 (25K) array, which contains 22,810 probe sets. The data were previously normalized using robust multi-array average (RMA) method, according to ftp://ftp.arabidopsis.org/home/tair/Microarrays/analyzed_data/README. Pre-calculated Log 2 values of signal densities were used to calculate co-expression values (Pearson correlation coefficients across all experiments) and expression abundances.
We considered all gene pairs on a given chromosome, i.e., each gene was paired with every other gene on the chromosome. If either of the two genes overlapped with another transcript, the pair was removed from the analysis. Distances were defined (i) as the number of intervening annotated transcripts, or (ii) as the number of base pairs between the transcripts (i.e., between the respective start and end points of transcription).
Co-expression values for individual pairs of genes were defined as Pearson's correlation coefficient between the two expression vectors across normalised hybridisation experiments [1, 2, 5–7, 12–14, 16]. To assess general trends of co-expression, we averaged co-expression values over all gene pairs at a given distance (and sometimes with a given orientation). Averaging over signed correlation coefficients does not allow to distinguish lack of co-expression from a superposition of equal amounts of positive and negative co-expression; however, this is justified as the overall distribution of co-expression values is approximately Gaussian (data not shown).
To compare our results to expectations, we constructed a data set of 100,000 pairs of randomly chosen protein-coding genes, for which we also calculated co-expression values. Mean co-expression of this dataset was 0.014. This positive 'random' co-expression is possibly due to the fact that several large gene sets are functionally correlated (such as 'growth-related' or 'stress-related' genes), and hence 'random' pairs usually contain pairs that are positively correlated (while the sets are not necessarily negatively correlated between them).
The statistical significance of differences in mean co-expression values was assessed using Wilcoxon rank sum tests. Throughout the manuscript, we use a threshold α = 0.05 to determine statistical significance. When multiple comparsions were performed (e.g., among the three orientations), we sometimes summarize results by reporting only P < the largest P-value.
We thank Guang-Zhong Wang for helpful discussions.
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nature reviews. 2004, 5 (4): 299-310. 10.1038/nrg1319.PubMedView ArticleGoogle Scholar
- Lercher MJ, Urrutia AO, Hurst LD: Clustering of housekeeping genes provides a unified model of gene order in the human genome. Nature genetics. 2002, 31 (2): 180-183. 10.1038/ng887.PubMedView ArticleGoogle Scholar
- Lee JM, Sonnhammer EL: Genomic gene clustering analysis of pathways in eukaryotes. Genome Res. 2003, 13 (5): 875-882. 10.1101/gr.737703.PubMed CentralPubMedView ArticleGoogle Scholar
- Trinklein ND, Aldred SF, Hartman SJ, Schroeder DI, Otillar RP, Myers RM: An abundance of bidirectional promoters in the human genome. Genome Res. 2004, 14 (1): 62-66. 10.1101/gr.1982804.PubMed CentralPubMedView ArticleGoogle Scholar
- Lercher MJ, Hurst LD: Co-expressed Yeast Genes Cluster Over a Long Range but are not Regularly Spaced. Journal of Molecular Biology. 2006, 359 (3): 825-831. 10.1016/j.jmb.2006.03.051.PubMedView ArticleGoogle Scholar
- Lercher MJ, Blumenthal T, Hurst LD: Coexpression of neighboring genes in Caenorhabditis elegans is mostly due to operons and duplicate genes. Genome Res. 2003, 13 (2): 238-243. 10.1101/gr.553803.PubMed CentralPubMedView ArticleGoogle Scholar
- Batada NN, Urrutia AO, Hurst LD: Chromatin remodelling is a major source of coexpression of linked genes in yeast. Trends in Genetics. 2007, 23 (10): 480-484. 10.1016/j.tig.2007.08.003.PubMedView ArticleGoogle Scholar
- Gilbert SF: Mechanisms for the environmental regulation of gene expression: ecological aspects of animal development. Journal of biosciences. 2005, 30 (1): 65-74. 10.1007/BF02705151.PubMedView ArticleGoogle Scholar
- Lanctot C, Cheutin T, Cremer M, Cavalli G, Cremer T: Dynamic genome architecture in the nuclear space: regulation of gene expression in three dimensions. Nature reviews. 2007, 8 (2): 104-115. 10.1038/nrg2041.PubMedView ArticleGoogle Scholar
- The Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.View ArticleGoogle Scholar
- Williams EJ, Bowles DJ: Coexpression of neighboring genes in the genome of Arabidopsis thaliana. Genome Res. 2004, 14 (6): 1060-1067. 10.1101/gr.2131104.PubMed CentralPubMedView ArticleGoogle Scholar
- Ren XY, Fiers M, Stiekema WJ, Nap JP: Local coexpression domains of two to four genes in the genome of Arabidopsis. Plant Physiology. 2005, 138 (2): 923-934. 10.1104/pp.104.055673.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhan S, Horrocks J, Lukens LN: Islands of co-expressed neighbouring genes in Arabidopsis thaliana suggest higher-order chromosome domains. Plant Journal. 2006, 45 (3): 347-357. 10.1111/j.1365-313X.2005.02619.x.PubMedView ArticleGoogle Scholar
- Ren XY, Stiekema WJ, Nap JP: Local coexpression domains in the genome of rice show no microsynteny with Arabidopsis domains. Plant Molecular Biology. 2007, 65 (1-2): 205-217. 10.1007/s11103-007-9209-0.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Q, Wan L, Li DY, Zhu LH, Qian MP, Deng MH: Searching for bidirectional promoters in Arabidopsis thaliana. Bmc Bioinformatics. 2009, 10 (Suppl 1): S29-10.1186/1471-2105-10-S1-S29.PubMed CentralPubMedView ArticleGoogle Scholar
- Krom N, Ramakrishna W: Comparative analysis of divergent and convergent gene pairs and their expression patterns in rice, Arabidopsis, and Populus. Plant Physiology. 2008, 147 (4): 1763-1773. 10.1104/pp.108.122416.PubMed CentralPubMedView ArticleGoogle Scholar
- Casneuf T, De Bodt S, Raes J, Maere S, Peer Van de Y: Nonrandom divergence of gene expression following gene and genome duplications in the flowering plant Arabidopsis thaliana. Genome Biology. 2006, 7 (2): 10.1186/gb-2006-7-2-r13.
- Kawauchi J, Mischo H, Braglia P, Rondon A, Proudfoot NJ: Budding yeast RNA polymerases I and II employ parallel mechanisms of transcriptional termination. Genes & Development. 2008, 22 (8): 1082-1092.View ArticleGoogle Scholar
- Tollervey D: Molecular biology: termination by torpedo. Nature. 2004, 432 (7016): 456-457. 10.1038/432456a.PubMedView ArticleGoogle Scholar
- Luo W, Johnson AW, Bentley DL: The role of Rat1 in coupling mRNA 3'-end processing to transcription termination: implications for a unified allosteric-torpedo model. Genes Dev. 2006, 20 (8): 954-965. 10.1101/gad.1409106.PubMed CentralPubMedView ArticleGoogle Scholar
- Poole RL: The TAIR Database. Methods Mol Biol. 2007, 406: 179-212. full_text.PubMedGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008, D1009-1014. 36 Database
- Thimmapuram J, Duan H, Liu L, Schuler MA: Bicistronic and fused monocistronic transcripts are derived from adjacent loci in the Arabidopsis genome. RNA. 2005, 11 (2): 128-138. 10.1261/rna.7114505.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu X, Han B: Evolutionary conservation of neighbouring gene pairs in plants. Gene. 2009, 437 (1-2): 71-79. 10.1016/j.gene.2009.02.012.PubMedView ArticleGoogle Scholar
- Neil H, Malabat C, d'Aubenton-Carafa Y, Xu Z, Steinmetz LM, Jacquier A: Widespread bidirectional promoters are the major source of cryptic transcripts in yeast. Nature. 2009, 457 (7232): 1038-1042. 10.1038/nature07747.PubMedView ArticleGoogle Scholar
- Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Munster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM: Bidirectional promoters generate pervasive transcription in yeast. Nature. 2009, 457 (7232): 1033-1037. 10.1038/nature07728.PubMed CentralPubMedView ArticleGoogle Scholar
- Core LJ, Waterfall JJ, Lis JT: Nascent RNA sequencing reveals widespread pausing and divergent initiation at human promoters. Science. 2008, 322 (5909): 1845-1848. 10.1126/science.1162228.PubMed CentralPubMedView ArticleGoogle Scholar
- Seila AC, Calabrese JM, Levine SS, Yeo GW, Rahl PB, Flynn RA, Young RA, Sharp PA: Divergent transcription from active promoters. Science. 2008, 322 (5909): 1849-1851. 10.1126/science.1162253.PubMed CentralPubMedView ArticleGoogle Scholar
- Babu MM, Janga SC, de Santiago I, Pombo A: Eukaryotic gene regulation in three dimensions and its impact on genome evolution. Current Opinion in Genetics & Development. 2008, 18 (6): 571-582.View ArticleGoogle Scholar
- Seila AC, Core LJ, Lis JT, Sharp PA: Divergent transcription: a new feature of active promoters. Cell Cycle. 2009, 8 (16): 2557-2564.PubMedView ArticleGoogle Scholar