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

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.

Pathogen Resistance Depending on Jacalin-Dirigent Chimeric Proteins Is Common among Poaceae but Absent in the Dicot Arabidopsis as Evidenced by Analysis of Homologous Single-Domain Proteins

MonocotJRLs are Poaceae-specific two-domain proteins that consist of a jacalin-related lectin (JRL) and a dirigent (DIR) domain which participate in multiple developmental processes, including disease resistance. For OsJAC1, a monocotJRL from rice, it has been confirmed that constitutive expression in transgenic rice or barley plants facilitates broad-spectrum disease resistance. In this process, both domains of OsJAC1 act cooperatively, as evidenced from experiments with artificially separated JRL- or DIR-domain-containing proteins. Interestingly, these chimeric proteins did not evolve in dicotyledonous plants. Instead, proteins with a single JRL domain, multiple JRL domains or JRL domains fused to domains other than DIR domains are present. In this study, we wanted to test if the cooperative function of JRL and DIR proteins leading to pathogen resistance was conserved in the dicotyledonous plant Arabidopsis thaliana. In Arabidopsis, we identified 50 JRL and 24 DIR proteins, respectively, from which seven single-domain JRL and two single-domain DIR candidates were selected. A single-cell transient gene expression assay in barley revealed that specific combinations of the Arabidopsis JRL and DIR candidates reduced the penetration success of barley powdery mildew. Strikingly, one of these pairs, AtJAX1 and AtDIR19, is encoded by genes located next to each other on chromosome one. However, when using natural variation and analyzing Arabidopsis ecotypes that express full-length or truncated versions of AtJAX1, the presence/absence of the full-length AtJAX1 protein could not be correlated with resistance to the powdery mildew fungus Golovinomyces orontii. Furthermore, an analysis of the additional JRL and DIR candidates in a bi-fluorescence complementation assay in Nicotiana benthamiana revealed no direct interaction of these JRL/DIR pairs. Since transgenic Arabidopsis plants expressing OsJAC1-GFP also did not show increased resistance to G. orontii, it was concluded that the resistance mediated by the synergistic activities of DIR and JRL proteins is specific for members of the Poaceae, at least regarding the resistance against powdery mildew. Arabidopsis lacks the essential components of the DIR-JRL-dependent resistance pathway.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

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

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

LncRNA-MAP3K4 regulates vascular inflammation through the p38 MAPK signaling pathway and cis-modulation of MAP3K4

Chronic vascular inflammation plays a key role in the pathogenesis of atherosclerosis. Long non-coding RNAs (lncRNAs) have emerged as essential inflammation regulators. We identify a novel lncRNA termed lncRNA-MAP3K4 that is enriched in the vessel wall and regulates vascular inflammation. In the aortic intima, lncRNA-MAP3K4 expression was reduced by 50% during the progression of atherosclerosis (chronic inflammation) and 70% during endotoxemia (acute inflammation). lncRNA-MAP3K4 knockdown reduced the expression of key inflammatory factors (eg, ICAM-1, E-selectin, MCP-1) in endothelial cells or vascular smooth muscle cells and decreased monocytes adhesion to endothelium, as well as reducing TNF-α, IL-1β, COX2 expression in macrophages. Mechanistically, lncRNA-MAP3K4 regulates inflammation through the p38 MAPK signaling pathway. lncRNA-MAP3K4 shares a bidirectional promoter with MAP3K4, an upstream regulator of the MAPK signaling pathway, and regulates its transcription in cis. lncRNA-MAP3K4 and MAP3K4 show coordinated expression in response to inflammation in vivo and in vitro. Similar to lncRNA-MAP3K4, MAP3K4 knockdown reduced the expression of inflammatory factors in several different vascular cells. Furthermore, lncRNA-MAP3K4 and MAP3K4 knockdown showed cooperativity in reducing inflammation in endothelial cells. Collectively, these findings unveil the role of a novel lncRNA in vascular inflammation by cis-regulating MAP3K4 via a p38 MAPK pathway.

Additional Layer of Regulation via Convergent Gene Orientation in Yeasts

Convergent gene pairs can produce transcripts with complementary sequences. We had shown that mRNA duplexes form in vivo in Saccharomyces cerevisiae via interactions of mRNA overlapping 3′-ends and can lead to posttranscriptional regulatory events. Here we show that mRNA duplex formation is restricted to convergent genes separated by short intergenic distance, independently of their 3′-untranslated region (UTR) length. We disclose an enrichment in genes involved in biological processes related to stress among these convergent genes. They are markedly conserved in convergent orientation in budding yeasts, meaning that this mode of posttranscriptional regulation could be shared in these organisms, conferring an additional level for modulating stress response. We thus investigated the mechanistic advantages potentially conferred by 3′-UTR mRNA interactions. Analysis of genome-wide transcriptome data revealed that Pat1 and Lsm1 factors, having 3′-UTR binding preference and participating to the remodeling of messenger ribonucleoprotein particles, bind differently these messenger-interacting mRNAs forming duplexes in comparison to mRNAs that do not interact (solo mRNAs). Functionally, messenger-interacting mRNAs show limited translational repression upon stress. We thus propose that mRNA duplex formation modulates the regulation of mRNA expression by limiting their access to translational repressors. Our results thus show that posttranscriptional regulation is an additional factor that determines the order of coding genes.

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.

Gene Repositioning Is Under Constraints After Evolutionary Conserved Gene Neighborhood Separate

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The Effects of Transcription Directions of Transgenes and the gypsy Insulators on the Transcript Levels of Transgenes in Transgenic Arabidopsis

Manipulation of a single abiotic stress-related gene could improve plant performance under abiotic stress conditions. To simultaneously increase plant tolerance to multiple stresses, it is usually required to overexpress two (or more) genes in transgenic plants. The common strategy is to assemble two or more expression cassettes, where each gene has its own promoter and terminator, within the same T-DNA. Does the arrangement of the two expression cassettes affect expression of the two transgenes? Can we use the Drosophila gypsy insulator sequence to increase the expression of the two transgenes? Answers to these questions would contribute to design better transformation vectors to maximize the effects of multi-gene transformation. Two Arabidopsis genes, PP2A-C5 and AVP1, and the gypsy insulator sequence were used to construct six transformation vectors with or without the gypsy insulator bracketing the two expression cassettes: uni-directional transcription, divergent transcription, and convergent transcription. Total RNAs were isolated for reverse transcription- quantitative real-time polymerase chain reaction (RT-qPCR) assays and a thorough statistical analysis was conducted for the RT-qPCR data. The results showed that the gypsy insulator does promote the expression of two transgenes in transgenic plants. Besides, the plants containing the divergent transcription cassettes tend to have more correlated expression of both genes.

Genome-wide identification of physically clustered genes suggests chromatin-level co-regulation in male reproductive development in Arabidopsis thaliana

Co-expression of physically linked genes occurs surprisingly frequently in eukaryotes. Such chromosomal clustering may confer a selective advantage as it enables coordinated gene regulation at the chromatin level. We studied the chromosomal organization of genes involved in male reproductive development in Arabidopsis thaliana. We developed an in-silico tool to identify physical clusters of co-regulated genes from gene expression data. We identified 17 clusters (96 genes) involved in stamen development and acting downstream of the transcriptional activator MS1 (MALE STERILITY 1), which contains a PHD domain associated with chromatin re-organization. The clusters exhibited little gene homology or promoter element similarity, and largely overlapped with reported repressive histone marks. Experiments on a subset of the clusters suggested a link between expression activation and chromatin conformation: qRT-PCR and mRNA in situ hybridization showed that the clustered genes were up-regulated within 48 h after MS1 induction; out of 14 chromatin-remodeling mutants studied, expression of clustered genes was consistently down-regulated only in hta9/hta11, previously associated with metabolic cluster activation; DNA fluorescence in situ hybridization confirmed that transcriptional activation of the clustered genes was correlated with open chromatin conformation. Stamen development thus appears to involve transcriptional activation of physically clustered genes through chromatin de-condensation.

Physical mapping of immune-related genes in Yesso scallop (Patinopecten yessoensis) using fluorescent in situ hybridization

The innate immune system plays a pivotal role in defending invasion of microorganisms for scallops. Previous studies on immune-related genes in the Yesso scallop, Patinopecten yessoensis (Jay, 1857) have mainly focused on characterization and expression pattern in response to bacterial challenge, no research has been carried out on the cytogenetic level yet. In the present study, eight fosmid clones containing the sequences of key immune-related genes (PyNFkB, PyTRAF2, PyTRAF4, PyTRAF7, PyMyd88-1, PyMyd88-3, PyMKK-7 and PyTNFR) were isolated and seven of them were successfully mapped on chromosomes of Patinopecten yessoensis utilizing fluorescence in situ hybridization. Wherein, PyMyd88-1, PyMyd88-3 and PyMKK-7 located on the same chromosome pair with adjacent positions and the other genes were mapped on four non-homologous chromosome pairs, showing a similar distribution to another five model species. The isolation and mapping of such genes of the Yesso scallop will lay a foundation for studies such as assignment of interested genes to chromosomes, construction cytogenetic maps and so on.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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.

DNA Damage-Induced Transcription of Transposable Elements and Long Non-coding RNAs in Arabidopsis Is Rare and ATM-Dependent

Induction and mobilization of transposable elements (TEs) following DNA damage or other stresses has been reported in prokaryotes and eukaryotes. Recently it was discovered that eukaryotic TEs are frequently associated with long non-coding RNAs (lncRNAs), many of which are also upregulated by stress. Yet, it is unknown whether DNA damage-induced transcriptional activation of TEs and lncRNAs occurs sporadically or is a synchronized, genome-wide response. Here we investigated the transcriptome of Arabidopsis wild-type (WT) and ataxia telangiectasia mutated (atm) mutant plants 3 h after induction of DNA damage. In WT, expression of 5.2% of the protein-coding genes is ≥2-fold changed, whereas in atm plants, only 2.6% of these genes are regulated, and the response of genes associated with DNA repair, replication, and cell cycle is largely lost. In contrast, only less than 0.6% of TEs and lncRNAs respond to DNA damage in WT plants, and the regulation of ≥95% of them is ATM-dependent. The ATM-downstream factors BRCA1, DRM1, JMJ30, AGO2, and the ATM-independent AGO4 participate in the regulation of individual TEs and lncRNAs. Remarkably, protein-coding genes located adjacent to DNA damage-responsive TEs and lncRNAs are frequently coexpressed, which is consistent with the hypothesis that TEs and lncRNAs located close to genes commonly function as controlling elements.

A low-temperature-responsive element involved in the regulation of the Arabidopsis thaliana At1g71850/At1g71860 divergent gene pair

The bidirectional promoter of the Arabidopsis thaliana gene pair At1g71850/At1g71860 harbors low-temperature-responsive elements, which participate in anti-correlated transcription regulation of the driving genes in response to environmental low temperature. A divergent gene pair is defined as two adjacent genes organized head to head in opposite orientation, sharing a common promoter region. Divergent gene pairs are mainly coexpressed, but some display opposite regulation. The mechanistic basis of such anti-correlated regulation is not well understood. Here, the regulation of the Arabidopsis thaliana gene pair At1g71850/At1g71860 was investigated. Semi-quantitative RT-PCR and Genevestigator analyses showed that while one of the pair was upregulated by exposure to low temperature, the same treatment downregulated the other. Promoter::GUS fusion transgenes were used to show that this behavior was driven by a bidirectional promoter, which harbored an as-1 motif, associated with the low-temperature response; mutation of this sequence produced a significant decrease in cold-responsive expression. With regard to the as-1 motif in the native orientation repressing the promoter's low-temperature responsiveness, the same as-1 motif introduced in the reverse direction showed a slight enhancement in the promoter's responsiveness to low-temperature exposure, indicating that the orientation of the motif was important for the promoter's activity. These findings provide new insights into the complex transcriptional regulation of bidirectional gene pairs as well as plant stress response.

A Novel Candidate Gene for Temperature-Dependent Sex Determination in the Common Snapping Turtle

Temperature-dependent sex determination (TSD) was described nearly 50 years ago. Researchers have since identified many genes that display differential expression at male- vs. female-producing temperatures. Yet, it is unclear whether these genes (1) are involved in sex determination per se, (2) are downstream effectors involved in differentiation of ovaries and testes, or (3) are thermo-sensitive but unrelated to gonad development. Here we present multiple lines of evidence linking CIRBP to sex determination in the snapping turtle, Chelydra serpentina We demonstrate significant associations between a single nucleotide polymorphism (SNP) (c63A > C) in CIRBP, transcript levels in embryonic gonads during specification of gonad fate, and sex in hatchlings from a thermal regime that produces mixed sex ratios. The A allele was induced in embryos exposed to a female-producing temperature, while expression of the C allele did not differ between female- and male-producing temperatures. In accord with this pattern of temperature-dependent, allele-specific expression, AA homozygotes were more likely to develop ovaries than AC heterozygotes, which, in turn, were more likely to develop ovaries than CC homozygotes. Multiple regression using SNPs in CIRBP and adjacent loci suggests that c63A > C may be the causal variant or closely linked to it. Differences in CIRBP allele frequencies among turtles from northern Minnesota, southern Minnesota, and Texas reflect small and large-scale latitudinal differences in TSD pattern. Finally, analysis of CIRBP protein localization reveals that CIRBP is in a position to mediate temperature effects on the developing gonads. Together, these studies strongly suggest that CIRBP is involved in determining the fate of the bipotential gonad.

Integrated mapping and characterization of the gene underlying the okra leaf trait in Gossypium hirsutum L

Diverse leaf morphology has been observed among accessions of Gossypium hirsutum, including okra leaf, which has advantages and disadvantages in cotton production. The okra leaf locus has been mapped to chromosome 15 of the Dt subgenome, but the underlying gene has yet to be identified. In this study, we used a combination of targeted association analysis, F2 population-based fine mapping, and comparative sequencing of orthologues to identify a candidate gene underlying the okra leaf trait in G. hirsutum. The okra leaf gene identified, GhOKRA, encoded a homeodomain leucine-zipper class I protein, whose closely related genes in several other plant species have been shown to be involved in regulating leaf morphology. The transcript levels of GhOKRA in shoot apices were positively correlated with the phenotypic expression of the okra leaf trait. Of the multiple sequence variations observed in the coding region among GrOKRA of Gossypium raimondii and GhOKRA-Dt of normal and okra/superokra leaf G. hirsutum accessions, a non-synonymous substitution near the N terminus and the variable protein sequences at the C terminus may be related to the leaf shape difference. Our results suggest that both transcription and protein activity of GhOKRA may be involved in regulating leaf shape. Furthermore, we found that non-reciprocal homoeologous recombination, or gene conversion, may have played a role in the origin of the okra leaf allele. Our results provided tools for further investigating and understanding the fundamental biological processes that are responsible for the cotton leaf shape variation and will help in the design of cotton plants with an ideal leaf shape for enhanced cotton production.

A comparative gene analysis with rice identified orthologous group II HKT genes and their association with Na(+) concentration in bread wheat

BACKGROUND

Although the HKT transporter genes ascertain some of the key determinants of crop salt tolerance mechanisms, the diversity and functional role of group II HKT genes are not clearly understood in bread wheat. The advanced knowledge on rice HKT and whole genome sequence was, therefore, used in comparative gene analysis to identify orthologous wheat group II HKT genes and their role in trait variation under different saline environments.

RESULTS

The four group II HKTs in rice identified two orthologous gene families from bread wheat, including the known TaHKT2;1 gene family and a new distinctly different gene family designated as TaHKT2;2. A single copy of TaHKT2;2 was found on each homeologous chromosome arm 7AL, 7BL and 7DL and each gene was expressed in leaf blade, sheath and root tissues under non-stressed and at 200 mM salt stressed conditions. The proteins encoded by genes of the TaHKT2;2 family revealed more than 93% amino acid sequence identity but ≤52% amino acid identity compared to the proteins encoded by TaHKT2;1 family. Specifically, variations in known critical domains predicted functional differences between the two protein families. Similar to orthologous rice genes on chromosome 6L, TaHKT2;1 and TaHKT2;2 genes were located approximately 3 kb apart on wheat chromosomes 7AL, 7BL and 7DL, forming a static syntenic block in the two species. The chromosomal region on 7AL containing TaHKT2;1 7AL-1 co-located with QTL for shoot Na(+) concentration and yield in some saline environments.

CONCLUSION

The differences in copy number, genes sequences and encoded proteins between TaHKT2;2 homeologous genes and other group II HKT gene families within and across species likely reflect functional diversity for ion selectivity and transport in plants. Evidence indicated that neither TaHKT2;2 nor TaHKT2;1 were associated with primary root Na(+) uptake but TaHKT2;1 may be associated with trait variation for Na(+) exclusion and yield in some but not all saline environments.

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.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

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

Transcriptome analysis of a petal anthocyanin polymorphism in the arctic mustard, Parrya nudicaulis

Angiosperms are renown for their diversity of flower colors. Often considered adaptations to pollinators, the most common underlying pigments, anthocyanins, are also involved in plants’ stress response. Although the anthocyanin biosynthetic pathway is well characterized across many angiosperms and is composed of a few candidate genes, the consequences of blocking this pathway and producing white flowers has not been investigated at the transcriptome scale. We take a transcriptome-wide approach to compare expression differences between purple and white petal buds in the arctic mustard, Parrya nudicaulis, to determine which genes’ expression are consistently correlated with flower color. Using mRNA-Seq and de novo transcriptome assembly, we assembled an average of 722 bp per gene (49.81% coding sequence based on the A. thaliana homolog) for 12,795 genes from the petal buds of a pair of purple and white samples. Our results correlate strongly with qRT-PCR analysis of nine candidate genes in the anthocyanin biosynthetic pathway where chalcone synthase has the greatest difference in expression between color morphs (P/W = ∼7×). Among the most consistently differentially expressed genes between purple and white samples, we found 3× more genes with higher expression in white petals than in purple petals. These include four unknown genes, two drought-response genes (CDSP32, ERD5), a cold-response gene (GR-RBP2), and a pathogen defense gene (DND1). Gene ontology analysis of the top 2% of genes with greater expression in white relative to purple petals revealed enrichment in genes associated with stress responses including cold, drought and pathogen defense. Unlike the uniform downregulation of chalcone synthase that may be directly involved in the loss of petal anthocyanins, the variable expression of several genes with greater expression in white petals suggest that the physiological and ecological consequences of having white petals may be microenvironment-dependent.

Neighboring genes show interchromosomal colocalization after their separation

The order of genes on eukaryotic chromosomes is nonrandom. Some neighboring genes show order conservation among species, while some neighboring genes separate during evolution. Here, we investigated whether neighboring genes show interactions after their separation. We found that neighboring gene pairs tend to show interchromosomal colocalization (i.e., nuclear colocalization) in the species in which they separate. These nuclear colocalized separated neighboring gene pairs 1) show neighborhood conservation in more species, 2) tend to be regulated by the same transcription factor, and 3) tend to be regulated by the same histone modification. These results suggest a mechanism by which neighboring genes could retain nuclear proximity after their separation.

Suppression of expression between adjacent genes within heterologous modules in yeast

Recent studies have shown that proximal arrangement of multiple genes can have complex effects on gene expression. For example, in the case of heterologous gene expression modules, certain arrangements of the selection marker and the gene expression cassette may have unintended consequences that limit the predictability and interpretability of module behaviors. The relationship between arrangement and expression has not been systematically characterized within heterologous modules to date. In this study, we quantitatively measured gene expression patterns of the selection marker (KlURA3 driven by the promoter, pKlURA) and the gene expression cassette (GFP driven by the galactose-inducible GAL1 promoter, pGAL1) in all their possible relative arrangements in Saccharomyces cerevisiae. First, we observed that pKlURA activity depends strongly on the relative arrangement and the activity of pGAL1. Most notably, we observed transcriptional suppression in the case of divergent arrangements: pKlURA activity was reduced when pGAL1 was inactive. Based on our nucleosome occupancy data, we attribute the observed transcriptional reduction to nucleosome repositioning. Second, we observed that pGAL1 activity also depends on the relative arrangement of pKlURA. In particular, strains with divergent promoters showed significantly different pGAL1 activation patterns from other strains, but only when their growth was compromised by lack of uracil. We reasoned that this difference in pGAL1 activation patterns arises from arrangement-dependent pKlURA activity that can affect the overall cell physiology (i.e., cell growth and survival in the uracil-depleted condition). Our results underscore the necessity to consider ramifications of promoter arrangement when using synthetic gene expression modules.

Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana

Short noncoding RNAs have been demonstrated to play important roles in regulation of gene expression and stress responses, but the repertoire and functions of long noncoding RNAs (lncRNAs) remain largely unexplored, particularly in plants. To explore the role of lncRNAs in disease resistance, we used a strand-specific RNA-sequencing approach to identify lncRNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana. Antisense transcription was found in c. 20% of the annotated A. thaliana genes. Several noncoding natural antisense transcripts responsive to F. oxysporum infection were found in genes implicated in disease defense. While the majority of the novel transcriptionally active regions (TARs) were adjacent to annotated genes and could be an extension of the annotated transcripts, 159 novel intergenic TARs, including 20 F. oxysporum-responsive lncTARs, were identified. Ten F. oxysporum-induced lncTARs were functionally characterized using T-DNA insertion or RNA-interference knockdown lines, and five were demonstrated to be related to disease development. Promoter analysis suggests that some of the F. oxysporum-induced lncTARs are direct targets of transcription factor(s) responsive to pathogen attack. Our results demonstrated that strand-specific RNA sequencing is a powerful tool for uncovering hidden levels of transcriptome and that IncRNAs are important components of the antifungal networks in A. thaliana.

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.

Comparative aspects of kisspeptin gene regulation

Kisspeptin plays an important role in the onset of puberty through stimulation of gonadotropin-releasing hormone (GnRH), a master molecule of reproduction. Furthermore, the existence of multiple kisspeptins is evident in most vertebrate species. Therefore, elucidating the regulatory mechanisms of the kisspeptin genes is important to understand the functions of multiple kisspeptin forms in the brain. This review focuses on the comparative aspects of kisspeptin gene regulation with an emphasis on the role of environmental signals including gonadal steroids, photoperiods and metabolic signals. These environmental signals differently regulate the kisspeptin genes distinctively in each species. In addition, photoperiodic regulation of the kisspeptin genes alters during sexual maturational, suggesting interactions between the gonadal hormone pathway and the photoperiod pathway. Further studies of the regulatory mechanisms of kisspeptin genes especially in teleosts which possess multiple kisspeptin/kisspeptin receptor systems will help to understand the precise role of multiple kisspeptin forms in different species.

Analysis of cellulose synthase genes from domesticated apple identifies collinear genes WDR53 and CesA8A: partial co-expression, bicistronic mRNA, and alternative splicing of CESA8A

Cellulose synthase (CesA) genes constitute a complex multigene family with six major phylogenetic clades in angiosperms. The recently sequenced genome of domestic apple, Malus×domestica, was mined for CesA genes, by blasting full-length cellulose synthase protein (CESA) sequences annotated in the apple genome against protein databases from the plant models Arabidopsis thaliana and Populus trichocarpa. Thirteen genes belonging to the six angiosperm CesA clades and coding for proteins with conserved residues typical of processive glycosyltransferases from family 2 were detected. Based on their phylogenetic relationship to Arabidopsis CESAs, as well as expression patterns, a nomenclature is proposed to facilitate further studies. Examination of their genomic organization revealed that MdCesA8-A is closely linked and co-oriented with WDR53, a gene coding for a WD40 repeat protein. The WDR53 and CesA8 genes display conserved collinearity in dicots and are partially co-expressed in the apple xylem. Interestingly, the presence of a bicistronic WDR53-CesA8A transcript was detected in phytoplasma-infected phloem tissues of apple. The bicistronic transcript contains a spliced intergenic sequence that is predicted to fold into hairpin structures typical of internal ribosome entry sites, suggesting its potential cap-independent translation. Surprisingly, the CesA8A cistron is alternatively spliced and lacks the zinc-binding domain. The possible roles of WDR53 and the alternatively spliced CESA8 variant during cellulose biosynthesis in M.×domestica are discussed.

Long noncoding RNAs in C. elegans

Thousands of long noncoding RNAs (lncRNAs) have been found in vertebrate animals, a few of which have known biological roles. To better understand the genomics and features of lncRNAs in invertebrates, we used available RNA-seq, poly(A)-site, and ribosome-mapping data to identify lncRNAs of Caenorhabditis elegans. We found 170 long intervening ncRNAs (lincRNAs), which had single- or multiexonic structures that did not overlap protein-coding transcripts, and about sixty antisense lncRNAs (ancRNAs), which were complementary to protein-coding transcripts. Compared to protein-coding genes, the lncRNA genes tended to be expressed in a stage-dependent manner. Approximately 25% of the newly identified lincRNAs showed little signal for sequence conservation and mapped antisense to clusters of endogenous siRNAs, as would be expected if they serve as templates and targets for these siRNAs. The other 75% tended to be more conserved and included lincRNAs with intriguing expression and sequence features associating them with processes such as dauer formation, male identity, sperm formation, and interaction with sperm-specific mRNAs. Our study provides a glimpse into the lncRNA content of a nonvertebrate animal and a resource for future studies of lncRNA function.

Prediction of operon-like gene clusters in the Arabidopsis thaliana genome based on co-expression analysis of neighboring genes

Operon-like arrangements of genes occur in eukaryotes ranging from yeasts and filamentous fungi to nematodes, plants, and mammals. In plants, several examples of operon-like gene clusters involved in metabolic pathways have recently been characterized, e.g. the cyclic hydroxamic acid pathways in maize, the avenacin biosynthesis gene clusters in oat, the thalianol pathway in Arabidopsis thaliana, and the diterpenoid momilactone cluster in rice. Such operon-like gene clusters are defined by their co-regulation or neighboring positions within immediate vicinity of chromosomal regions. A comprehensive analysis of the expression of neighboring genes therefore accounts a crucial step to reveal the complete set of operon-like gene clusters within a genome. Genome-wide prediction of operon-like gene clusters should contribute to functional annotation efforts and provide novel insight into evolutionary aspects acquiring certain biological functions as well. We predicted co-expressed gene clusters by comparing the Pearson correlation coefficient of neighboring genes and randomly selected gene pairs, based on a statistical method that takes false discovery rate (FDR) into consideration for 1469 microarray gene expression datasets of A. thaliana. We estimated that A. thaliana contains 100 operon-like gene clusters in total. We predicted 34 statistically significant gene clusters consisting of 3 to 22 genes each, based on a stringent FDR threshold of 0.1. Functional relationships among genes in individual clusters were estimated by sequence similarity and functional annotation of genes. Duplicated gene pairs (determined based on BLAST with a cutoff of E<10(-5)) are included in 27 clusters. Five clusters are associated with metabolism, containing P450 genes restricted to the Brassica family and predicted to be involved in secondary metabolism. Operon-like clusters tend to include genes encoding bio-machinery associated with ribosomes, the ubiquitin/proteasome system, secondary metabolic pathways, lipid and fatty-acid metabolism, and the lipid transfer system.

A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat

To improve our understanding of the organization and regulation of the wheat (Triticum aestivum) gene space, we established a transcription map of a wheat chromosome (3B) by hybridizing a newly developed wheat expression microarray with bacterial artificial chromosome pools from a new version of the 3B physical map as well as with cDNA probes derived from 15 RNA samples. Mapping data for almost 3,000 genes showed that the gene space spans the whole chromosome 3B with a 2-fold increase of gene density toward the telomeres due to an increase in the number of genes in islands. Comparative analyses with rice (Oryza sativa) and Brachypodium distachyon revealed that these gene islands are composed mainly of genes likely originating from interchromosomal gene duplications. Gene Ontology and expression profile analyses for the 3,000 genes located along the chromosome revealed that the gene islands are enriched significantly in genes sharing the same function or expression profile, thereby suggesting that genes in islands acquired shared regulation during evolution. Only a small fraction of these clusters of cofunctional and coexpressed genes was conserved with rice and B. distachyon, indicating a recent origin. Finally, genes with the same expression profiles in remote islands (coregulation islands) were identified suggesting long-distance regulation of gene expression along the chromosomes in wheat.

Gene expression: sizing it all up

Genomic architecture appears to be a largely unexplored component of gene expression. That architecture can be related to chromatin domains, transposable element neighborhoods, epigenetic modifications of the genome, and more. Although surely not the end of the story, we are learning that when it comes to gene expression, size is also important. We have been surprised to find that certain patterns of expression, tissue specific versus constitutive, or high expression versus low expression, are often associated with physical attributes of the gene and genome. Multiple studies have shown an inverse relationship between gene expression patterns and various physical parameters of the genome such as intron size, exon size, intron number, and size of intergenic regions. An increase in expression level and breadth often correlates with a decrease in the size of physical attributes of the gene. Three models have been proposed to explain these relationships. Contradictory results were found in several organisms when expression level and expression breadth were analyzed independently. However, when both factors were combined in a single study a novel relationship was revealed. At low levels of expression, an increase in expression breadth correlated with an increase in genic, intergenic, and intragenic sizes. Contrastingly, at high levels of expression, an increase in expression breadth inversely correlated with the size of the gene. In this article we explore the several hypotheses regarding genome physical parameters and gene expression.

Comparative expression profiles for KiSS-1 and REN genes in preeclamptic and healthy placental tissues

OBJECTIVE

The aim of the present work was to look at differences in the placental tissue expression of KiSS-1 and REN genes from preeclamptic and healthy pregnant women, that could account for a possible synergistic function for both genes in the pathogenesis of preeclampsia.

STUDY DESIGN

This case-control study involved 27 preeclamptic women and 27 normoevolutive pregnant women. cDNA was obtained from placental tissue to carry out qPCR for both KiSS-1 and REN genes in order to compare mRNA expression levels in the studied groups. Statistical analysis showed expression differences that correlate with clinical and/or biochemical variables.

RESULTS

Higher expression for KiSS-1 in PEE vs. control woman (p=0.001) was observed, whereas no difference was observed for REN expression (p=0.300) when all the subjects were included. However, REN expression was significant higher when the samples were stratified according to preeclampsia severity. For 18 mild preeclamptic patients the p-value was p=0.001 compared to their controls, while for the remaining nine with severe preeclampsia the expression became significant (p=0.001).

CONCLUSION

Our results suggest that the high KiSS-1 expression seen in preeclamptic patients is in accordance with its role as an inhibitor of trophoblast invasiveness and maintained until the end of gestation. On the other hand, aggressive therapeutic management and/or severity status of patients have a direct effect on placental REN expression levels, masking the natural high expression of this gene on preeclamptic placental tissue. Therefore it was not possible to establish a real concordant expression profile for KiSS-1 and REN genes.

Coexpression of linked gene pairs persists long after their separation

In many eukaryotes, physically linked gene pairs tend to be coexpressed. However, it is still controversial to what extent this neighbor coexpression is maintained by selection and to what extent it is nonselective, purely mechanistic "leaky expression." Here, we analyze expression patterns of gene pairs that have lost their linkage in the evolution of Saccharomyces cerevisiae since its last common ancestor with Kluyveromyces waltii or that were never linked in the S. cerevisiae lineage but became neighbors in a related yeast. We demonstrate that coexpression of many linked genes is retained long after their separation and is thus likely to be functionally important. In addition, unlinked gene pairs that recently became neighbors in other yeast species tend to be coexpressed in S. cerevisiae. This suggests that natural selection often favors chromosomal rearrangements in which coexpressed genes become neighbors. Contrary to previous suggestions, selectively favorable coexpression appears not to be restricted to bidirectional promoters.

Minimal regulatory spaces in yeast genomes

BACKGROUND

The regulatory information encoded in the DNA of promoter regions usually enforces a minimal, non-zero distance between the coding regions of neighboring genes. However, the size of this minimal regulatory space is not generally known. In particular, it is unclear if minimal promoter size differs between species and between uni- and bi-directionally acting regulatory regions.

RESULTS

Analyzing the genomes of 11 yeasts, we show that the lower size limit on promoter-containing regions is species-specific within a relatively narrow range (80-255 bp). This size limit applies equally to regions that initiate transcription on one or both strands, indicating that bi-directional promoters and uni-directional promoters are constrained similarly. We further find that young, species-specific regions are on average much longer than older regions, suggesting either a bias towards deletions or selection for genome compactness in yeasts. While the length evolution of promoter-less intergenic regions is well described by a simplistic, purely neutral model, regions containing promoters typically show an excess of unusually long regions. Regions flanked by divergently transcribed genes have a bi-modal length distribution, with short lengths found preferentially among older regions. These old, short regions likely harbor evolutionarily conserved bi-directionally active promoters. Surprisingly, some of the evolutionarily youngest regions in two of the eleven species (S. cerevisiae and K. waltii) are shorter than the lower limit observed in older regions.

CONCLUSIONS

The minimal chromosomal space required for transcriptional regulation appears to be relatively similar across yeast species, and is the same for uni-directional and bi-directional promoters. New intergenic regions created by genome rearrangements tend to evolve towards the more narrow size distribution found among older regions.

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.

More than the sum of its parts--how to achieve a specific transcriptional response to abiotic stress

A rapid and appropriate response to stress is key to survival. A major part of plant adaptation to abiotic stresses is regulated at the level of gene expression. The regulatory steps involved in accurate expression of stress related genes need to be tailored to the specific stress for optimal plant performance. Accumulating evidence suggests that there are several processes contributing to signalling specificity: post-translational activation and selective nuclear import of transcription factors, regulation of DNA accessibility by chromatin modifying and remodelling enzymes, and cooperation between two or more response elements in a stress-responsive promoter. These mechanisms should not be viewed as independent events, instead the nuclear DNA is in a complex landscape where many proteins interact, compete, and regulate each other. Hence future studies should consider an integrated view of gene regulation composed of numerous chromatin associated proteins in addition to transcription factors. Although most studies have focused on a single regulatory mechanism, it is more likely the combined actions of several mechanisms that provide a stress specific output. In this review recent progress in abiotic stress signalling is discussed with emphasis on possible mechanisms for generating specific responses.

The transcriptome landscape of Arabidopsis male meiocytes from high-throughput sequencing: the complexity and evolution of the meiotic process

Meiosis is essential for eukaryotic sexual reproduction, with two consecutive rounds of nuclear divisions, allowing production of haploid gametes. Information regarding the meiotic transcriptome should provide valuable clues about global expression patterns and detailed gene activities. Here we used RNA sequencing to explore the transcriptome of a single plant cell type, the Arabidopsis male meiocyte, detecting the expression of approximately 20 000 genes. Transcription of introns of >400 genes was observed, suggesting previously unannotated exons. More than 800 genes may be preferentially expressed in meiocytes, including known meiotic genes. Of the 3378 Pfam gene families in the Arabidopsis genome, 3265 matched meiocyte-expressed genes, and 18 gene families were over-represented in male meiocytes, including transcription factor and other regulatory gene families. Expression was detected for many genes thought to encode meiosis-related proteins, including MutS homologs (MSHs), kinesins and ATPases. We identified more than 1000 orthologous gene clusters that are also expressed in meiotic cells of mouse and fission yeast, including 503 single-copy genes across the three organisms, with a greater number of gene clusters shared between Arabidopsis and mouse than either share with yeast. Interestingly, approximately 5% transposable element genes were apparently transcribed in male meiocytes, with a positive correlation to the transcription of neighboring genes. In summary, our RNA-Seq transcriptome data provide an overview of gene expression in male meiocytes and invaluable information for future functional studies.