Analysis and recognition of the GAGA transcription factor binding sites in Drosophila genes

Developmental and Housekeeping Genes: Two Types of Genetic Organization in the Drosophila Genome

We developed a procedure for locating genes on Drosophila melanogaster polytene chromosomes and described three types of chromosome structures (gray bands, black bands, and interbands), which differed markedly in morphological and genetic properties. This was reached through the use of our original methods of molecular and genetic analysis, electron microscopy, and bioinformatics data processing. Analysis of the genome-wide distribution of these properties led us to a bioinformatics model of the Drosophila genome organization, in which the genome was divided into two groups of genes. One was constituted by 65, in which the genome was divided into two groups, 62 genes that are expressed in most cell types during life cycle and perform basic cellular functions (the so-called "housekeeping genes"). The other one was made up of 3162 genes that are expressed only at particular stages of development ("developmental genes"). These two groups of genes are so different that we may state that the genome has two types of genetic organization. Different are the timings of their expression, chromatin packaging levels, the composition of activating and deactivating proteins, the sizes of these genes, the lengths of their introns, the organization of the promoter regions of the genes, the locations of origin recognition complexes (ORCs), and DNA replication timings.

The misregulation of mitochondria-associated genes caused by GAGA-factor lack promotes autophagic germ cell death in Drosophila testes

The Drosophila GAGA-factor encoded by the Trithorax-like (Trl) gene is DNA-binding protein with unusually wide range of applications in diverse cell contexts. In Drosophila spermatogenesis, reduced GAGA expression caused by Trl mutations induces mass autophagy leading to germ cell death. In this work, we investigated the contribution of mitochondrial abnormalities to autophagic germ cell death in Trl gene mutants. Using a cytological approach, in combination with an analysis of high-throughput RNA sequencing (RNA-seq) data, we demonstrated that the GAGA deficiency led to considerable defects in mitochondrial ultrastructure, by causing misregulation of GAGA target genes encoding essential components of mitochondrial molecular machinery. Mitochondrial anomalies induced excessive production of reactive oxygen species and their release into the cytoplasm, thereby provoking oxidative stress. Changes in transcription levels of some GAGA-independent genes in the Trl mutants indicated that testis cells experience ATP deficiency and metabolic aberrations, that may trigger extensive autophagy progressing to cell death.

The effects of maternal care on the developmental transcriptome and metatranscriptome of a wild bee

Maternal care acts as a strong environmental stimulus that can induce phenotypic plasticity in animals and may also alter their microbial communities through development. Here, we characterize the developmental metatranscriptome of the small carpenter bee, Ceratina calcarata, across developmental stages and in the presence or absence of mothers. Maternal care had the most influence during early development, with the greatest number and magnitude of differentially expressed genes between maternal care treatments, and enrichment for transcription factors regulating immune response in motherless early larvae. Metatranscriptomic data revealed fungi to be the most abundant group in the microbiome, with Aspergillus the most abundant in early larvae raised without mothers. Finally, integrative analysis between host transcriptome and metatranscriptome highlights several fungi correlating with developmental and immunity genes. Our results provide characterizations of the influence of maternal care on gene expression and the microbiome through development in a wild bee.

The complex role of transcription factor GAGA in germline death during Drosophila spermatogenesis: transcriptomic and bioinformatic analyses

The GAGA protein (also known as GAF) is a transcription factor encoded by the Trl gene in D. melanogaster. GAGA is involved in the regulation of transcription of many genes at all stages of fly development and life. Recently, we investigated the participation of GAGA in spermatogenesis and discovered that Trl mutants experience massive degradation of germline cells in the testes. Trl underexpression induces autophagic death of spermatocytes, thereby leading to reduced testis size. Here, we aimed to determine the role of the transcription factor GAGA in the regulation of ectopic germline cell death. We investigated how Trl underexpression affects gene expression in the testes. We identified 15,993 genes in three biological replicates of our RNA-seq analysis and compared transcript levels between hypomorphic Trl ^R85/Trl ³⁶² and Oregon testes. A total of 2,437 differentially expressed genes were found, including 1,686 upregulated and 751 downregulated genes. At the transcriptional level, we detected the development of cellular stress in the Trl-mutant testes: downregulation of the genes normally expressed in the testes (indicating slowed or abrogated spermatocyte differentiation) and increased expression of metabolic and proteolysis-related genes, including stress response long noncoding RNAs. Nonetheless, in the Flybase Gene Ontology lists of genes related to cell death, autophagy, or stress, there was no enrichment with GAGA-binding sites. Furthermore, we did not identify any specific GAGA-dependent cell death pathway that could regulate spermatocyte death. Thus, our data suggest that GAGA deficiency in male germline cells leads to an imbalance of metabolic processes, impaired mitochondrial function, and cell death due to cellular stress.

Motif models proposing independent and interdependent impacts of nucleotides are related to high and low affinity transcription factor binding sites in Arabidopsis

Position weight matrix (PWM) is the traditional motif model representing the transcription factor (TF) binding sites. It proposes that the positions contribute independently to TFs binding affinity, although this hypothesis does not fit the data perfectly. This explains why PWM hits are missing in a substantial fraction of ChIP-seq peaks. To study various modes of the direct binding of plant TFs, we compiled the benchmark collection of 111 ChIP-seq datasets for Arabidopsis thaliana, and applied the traditional PWM, and two alternative motif models BaMM and SiteGA, proposing the dependencies of the positions. The variation in the stringency of the recognition thresholds for the models proposed that the hits of PWM, BaMM, and SiteGA models are associated with the sites of high/medium, any, and low affinity, respectively. At the medium recognition threshold, about 60% of ChIP-seq peaks contain PWM hits consisting of conserved core consensuses, while BaMM and SiteGA provide hits for an additional 15% of peaks in which a weaker core consensus is compensated through intra-motif dependencies. The presence/absence of these dependencies in the motifs of alternative/traditional models was confirmed by the dependency logo DepLogo visualizing the position-wise partitioning of the alignments of predicted sites. We exemplify the detailed analysis of ChIP-seq profiles for plant TFs CCA1, MYC2, and SEP3. Gene ontology (GO) enrichment analysis revealed that among the three motif models, the SiteGA had the highest portions of genes with the significantly enriched GO terms among all predicted genes. We showed that both alternative motif models provide for traditional PWM greater extensions in predicted sites for TFs MYC2/SEP3 with condition/tissue specific functions, compared to those for TF CCA1 with housekeeping functions. Overall, the combined application of standard and alternative motif models is beneficial to detect various modes of the direct TF-DNA interactions in the maximal portion of ChIP-seq loci.

[Lack of GAGA protein in Trl mutants causes massive cell death in Drosophila spermatogenesis and oogenesis]

Белок дрозофилы GAGA (GAF) является фактором эпигенетической регуляции транскрипции большой группы генов с широким разнообразием клеточных функций. GAF кодируется геном Trithorax-like (Trl), который экспрессируется в различных органах и тканях на всех стадиях онтогенеза дрозофилы. Мутации этого гена вызывают множественные нарушения развития. В предыдущих работах мы показали, что этот белок необходим для развития половой системы как самцов, так и самок дрозофилы. Снижение экспрессии гена Trl приводило ко множественным нарушениям спермато- и оогенеза. Одно из значительных нарушений было связано с массовой деградацией и потерей клеток зародышевого пути, что позволило предположить, что этот белок вовлечен в регуляцию клеточной гибели. В представленной работе мы провели более детальное цитологическое исследование, чтобы определить, какой тип гибели клеток зародышевого пути характерен для Trl-мутантов, и происходят ли нарушения или изменения этого процесса по сравнению с нормой. Полученные результаты показали, что недостаток белка GAF вызывает массовую гибель клеток зародышевого пути как у самок, так и самцов дрозофилы, но проявляется эта гибель в зависимости от пола по-разному. У самок, мутантных по гену Trl, фенотипически этот процесс не отличается от нормы и в гибнущих яйцевых камерах выявлены признаки апоптоза и аутофагии клеток зародышевого пути. У самцов, мутантных по гену Trl, в отличие от самок, не обнаружены признаки апоптоза. У самцов мутации Trl индуцируют массовую гибель клеток через аутофагию, что не характерно для сперматогенеза дрозофилы и не описано ранее ни в норме, ни у мутаций по другим генам. Таким образом, недостаток GAF у мутантов Trl приводит к усилению апоптотической и аутофагической гибели клеток зародышевого пути. Эктопическая клеточная гибель и атрофия зародышевой линии, вероятно, связаны с нарушением экспрессии генов-мишеней GAGAфактора, среди которых есть гены, регулирующие как апоптоз, так и аутофагию.

Migration of primordial germline cells is negatively regulated by surrounding somatic cells during early embryogenesis in Drosophila melanogaster

Cell migration is an important morphogenetic process necessary at different stages of individual development and body functioning. The initiation and maintenance of the cell movement state requires the activation of many factors involved in the regulation of transcription, signal transduction, adhesive interactions, modulation of membranes and the cytoskeleton. However, cell movement depends on the status of both migrating and surrounding cells, interacting with each other during movement. The surrounding cells or cell matrix not only form a substrate for movement, but can also participate in the spatio-temporal regulation of the migration. At present, there is no exact understanding of the genetic mechanisms of this regulation. To determine the role of the cell environment in the regulation of individual cell migration, we studied the migration of primordial germline cells (PGC) during early embryogenesis in Drosophila melanogaster. Normally, PGC are formed at the 3rd stage of embryogenesis at the posterior pole of the embryo. During gastrulation (stages 6-7), PGC as a consolidated cell group passively transfers into the midgut primordium. Further, PGC are individualized, acquire an amoeboid form, and actively move through the midgut epithelium and migrate to the 5-6 abdominal segment of the embryo, where they form paired embryonic gonads. We screened for genes expressed in the epithelium surrounding PGC during early embryogenesis and affecting their migration. We identified the myc, Hph, stat92E, Tre-1, and hop genes, whose RNA interference leads to premature active PGC migration at stages 4-7 of embryogenesis. These genes can be divided into two groups: 1) modulators of JAK/STAT pathway activity inducing PGC migration (stat92E, Tre-1, hop), and 2) myc and Hph involved in epithelial morphogenesis and polarization, i. e. modifying the permeability of the epithelial barrier. Since a depletion of each of these gene products resulted in premature PGC migration, we can conclude that, normally, the somatic environment negatively regulates PGC migration during early Drosophila embryogenesis.

GAGA Regulates Border Cell Migration in Drosophila

Collective cell migration is a complex process that happens during normal development of many multicellular organisms, as well as during oncological transformations. In Drosophila oogenesis, a small set of follicle cells originally located at the anterior tip of each egg chamber become motile and migrate as a cluster through nurse cells toward the oocyte. These specialized cells are referred to as border cells (BCs) and provide a simple and convenient model system to study collective cell migration. The process is known to be complexly regulated at different levels and the product of the slow border cells (slbo) gene, the C/EBP transcription factor, is one of the key elements in this process. However, little is known about the regulation of slbo expression. On the other hand, the ubiquitously expressed transcription factor GAGA, which is encoded by the Trithorax-like (Trl) gene was previously demonstrated to be important for Drosophila oogenesis. Here, we found that Trl mutations cause substantial defects in BC migration. Partially, these defects are explained by the reduced level of slbo expression in BCs. Additionally, a strong genetic interaction between Trl and slbo mutants, along with the presence of putative GAGA binding sites within the slbo promoter and enhancer, suggests the direct regulation of this gene by GAGA. This idea is supported by the reduction in the slbo-Gal4-driven GFP expression within BC clusters in Trl mutant background. However, the inability of slbo overexpression to compensate defects in BC migration caused by Trl mutations suggests that there are other GAGA target genes contributing to this process. Taken together, the results define GAGA as another important regulator of BC migration in Drosophila oogenesis.

Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes

Microsatellite repeat DNA is best known for its length mutability, which is implicated in several neurological diseases and cancers, and often exploited as a genetic marker. Less well-known is the body of work exploring the widespread and surprisingly diverse functional roles of microsatellites. Recently, emerging evidence includes the finding that normal microsatellite polymorphism contributes substantially to the heritability of human gene expression on a genome-wide scale, calling attention to the task of elucidating the mechanisms involved. At present, these are underexplored, but several themes have emerged. I review evidence demonstrating roles for microsatellites in modulation of transcription factor binding, spacing between promoter elements, enhancers, cytosine methylation, alternative splicing, mRNA stability, selection of transcription start and termination sites, unusual structural conformations, nucleosome positioning and modification, higher order chromatin structure, noncoding RNA, and meiotic recombination hot spots.

Hidden heterogeneity of transcription factor binding sites: A case study of SF-1

Steroidogenic factor 1 (SF-1) belongs to a small group of the transcription factors that bind DNA only as a monomer. Three different approaches-Sitecon, SiteGA, and oPWM-constructed using the same training sample of experimentally confirmed SF-1 binding sites have been used to recognize these sites. The appropriate prediction thresholds for recognition models have been selected. Namely, the thresholds concordant by false positive or negative rates for various methods were used to optimize the discrimination of steroidogenic gene promoters from the datasets of non-specific promoters. After experimental verification, the models were used to analyze the ChIP-seq data for SF-1. It has been shown that the sets of sites recognized by different models overlap only partially and that an integration of these models allows for identification of SF-1 sites in up to 80% of the ChIP-seq loci. The structures of the sites detected using the three recognition models in the ChIP-seq peaks falling within the [-5000, +5000] region relative to the transcription start sites (TSS) extracted from the FANTOM5 project have been analyzed. The MATLIGN classified the frequency matrices for the sites predicted by oPWM, Sitecon, and SiteGA into two groups. The first group is described by oPWM/Sitecon and the second, by SiteGA. Gene ontology (GO) analysis has been used to clarify the differences between the sets of genes carrying different variants of SF-1 binding sites. Although this analysis in general revealed a considerable overlap in GO terms for the genes carrying the binding sites predicted by oPWM, Sitecon, or SiteGA, only the last method elicited notable trend to terms related to negative regulation and apoptosis. The results suggest that the SF-1 binding sites are different in both their structure and the functional annotation of the set of target genes correspond to the predictions by oPWM+Sitecon and SiteGA. Further application of Homer software for de novo identification of enriched motifs in ChIP-Seq data for SF-1ChIP-seq dataset gave the data similar to oPWM+Sitecon.

GAGA factor maintains nucleosome-free regions and has a role in RNA polymerase II recruitment to promoters

Previous studies have shown that GAGA Factor (GAF) is enriched on promoters with paused RNA Polymerase II (Pol II), but its genome-wide function and mechanism of action remain largely uncharacterized. We assayed the levels of transcriptionally-engaged polymerase using global run-on sequencing (GRO-seq) in control and GAF-RNAi Drosophila S2 cells and found promoter-proximal polymerase was significantly reduced on a large subset of paused promoters where GAF occupancy was reduced by knock down. These promoters show a dramatic increase in nucleosome occupancy upon GAF depletion. These results, in conjunction with previous studies showing that GAF directly interacts with nucleosome remodelers, strongly support a model where GAF directs nucleosome displacement at the promoter and thereby allows the entry Pol II to the promoter and pause sites. This action of GAF on nucleosomes is at least partially independent of paused Pol II because intergenic GAF binding sites with little or no Pol II also show GAF-dependent nucleosome displacement. In addition, the insulator factor BEAF, the BEAF-interacting protein Chriz, and the transcription factor M1BP are strikingly enriched on those GAF-associated genes where pausing is unaffected by knock down, suggesting insulators or the alternative promoter-associated factor M1BP protect a subset of GAF-bound paused genes from GAF knock-down effects. Thus, GAF binding at promoters can lead to the local displacement of nucleosomes, but this activity can be restricted or compensated for when insulator protein or M1BP complexes also reside at GAF bound promoters.

Transcriptional regulation is critical for proper gene expression in response to environmental changes and developmental programs. Eukaryotes have evolved multiple mechanisms by which transcription factors regulate transcription. One mechanism is the reorganization of chromatin to allow Pol II recruitment. Another is the release of promoter-proximal paused Pol II, where Pol II transcription that is halted 20–60 bases downstream of the transcription start site (TSS) is allowed to enter into productive elongation through the gene body. The Drosophila transcription factor GAF binds to genes that undergo pausing and interacts with nucleosome remodelers and the pausing factor NELF. Thus, GAF can regulate multiple points necessary for transcription, but its mechanistic role is not fully understood genome-wide. We depleted GAF from cells and examined the genome-wide changes in Pol II and nucleosome distributions across genes. We found that GAF depletion reduces polymerase density at genes where GAF binds just upstream of the TSS, and results in nucleosomes moving into the promoter region. Our results show that GAF is important for maintaining the promoter accessibility, allowing Pol II to be recruited to promoters and enter the pause sites downstream of the TSS. Thus, GAF is critical for providing the chromatin environment necessary for the proper control of gene expression.

The polymorphisms in the promoter of HSP90 gene and their association with heat tolerance of bay scallop

The heat shock protein 90 (HSP90) is a highly abundant and ubiquitous molecular chaperone which plays essential roles in many cellular processes. In the present study, the messenger RNA (mRNA) expressions of HSP90 after acute heat stress were investigated in two bay scallop populations (Argopecten irradians irradians and Argopecten irradians concentricus). The heat-resistant scallop A. i. concentricus, which is distributed in Zhanjiang, China, exhibited significantly higher induction of HSP90 compared with that of the heat-sensitive scallop A. i. irradians, which is distributed in Qinhuangdao, China. The promoter sequence of HSP90 gene from bay scallop (AiHSP90) was cloned, and the polymorphisms within this region were investigated by sequencing to analyze their association with heat tolerance. A total of six single nucleotide polymorphisms (SNPs), including -1167 T-C, -1023 A-C, -799 C-T, -774 A-G, -686 C-T, and -682 A-C, were identified in the amplified promoter region, and most of them affected the putative transcription factor binding sites except for locus -1167. All the six SNP sites were found to be associated with heat tolerance after Hardy-Weinberg equilibrium (HWE) and association analysis. Moreover, haplotypes CACACC and TCTATC were also found to be associated with heat tolerance based on the result of linkage disequilibrium and association analysis. The results provided insights into the molecular mechanisms underlying the thermal adaptation of different congener endemic bay scallops, which suggested that the increased heat tolerance of A. i. concentricus (compared with A. i. irradians) was associated with the higher expression of AiHSP90. Meanwhile, the six genotypes (-1167 TT, -1023 CC, -799 TT, -774 GG, -686 CC, and -682 AA) and two haplotypes (CACACC and TCTATC) could be used as potential markers for scallop selection breeding with higher heat tolerance.

Activity of heat shock genes' promoters in thermally contrasting animal species

Heat shock gene promoters represent a highly conserved and universal system for the rapid induction of transcription after various stressful stimuli. We chose pairs of mammalian and insect species that significantly differ in their thermoresistance and constitutive levels of Hsp70 to compare hsp promoter strength under normal conditions and after heat shock (HS). The first pair includes the HSPA1 gene promoter of camel (Camelus dromedarius) and humans. It was demonstrated that the camel HSPA1A and HSPA1L promoters function normally in vitro in human cell cultures and exceed the strength of orthologous human promoters under basal conditions. We used the same in vitro assay for Drosophila melanogaster Schneider-2 (S2) cells to compare the activity of the hsp70 and hsp83 promoters of the second species pair represented by Diptera, i.e., Stratiomys singularior and D. melanogaster, which dramatically differ in thermoresistance and the pattern of Hsp70 accumulation. Promoter strength was also monitored in vivo in D. melanogaster strains transformed with constructs containing the S. singularior hsp70 ORF driven either by its own promoter or an orthologous promoter from the D. melanogaster hsp70Aa gene. Analysis revealed low S. singularior hsp70 promoter activity in vitro and in vivo under basal conditions and after HS in comparison with the endogenous promoter in D. melanogaster cells, which correlates with the absence of canonical GAGA elements in the promoters of the former species. Indeed, the insertion of GAGA elements into the S. singularior hsp70 regulatory region resulted in a dramatic increase in promoter activity in vitro but only modestly enhanced the promoter strength in the larvae of the transformed strains. In contrast with hsp70 promoters, hsp83 promoters from both of the studied Diptera species demonstrated high conservation and universality.

Computational analysis of auxin responsive elements in the Arabidopsis thaliana L. genome

Auxin responsive elements (AuxRE) were found in upstream regions of target genes for ARFs (Auxin response factors). While Chip-seq data for most of ARFs are still unavailable, prediction of potential AuxRE is restricted by consensus models that detect too many false positive sites. Using sequence analysis of experimentally proven AuxREs, we revealed both an extended nucleotide context pattern for AuxRE itself and three distinct types of its coupling motifs (Y-patch, AuxRE-like, and ABRE-like), which together with AuxRE may form the composite elements. Computational analysis of the genome-wide distribution of the predicted AuxREs and their impact on auxin responsive gene expression allowed us to conclude that: (1) AuxREs are enriched around the transcription start site with the maximum density in 5'UTR; (2) AuxREs mediate auxin responsive up-regulation, not down-regulation. (3) Directly oriented single AuxREs and reverse multiple AuxREs are mostly associated with auxin responsiveness. In the composite AuxRE elements associated with auxin response, ABRE-like and Y-patch are 5'-flanking or overlapping AuxRE, whereas AuxRE-like motif is 3'-flanking. The specificity in location and orientation of the coupling elements suggests them as potential binding sites for ARFs partners.

Genetic basis of olfactory cognition: extremely high level of DNA sequence polymorphism in promoter regions of the human olfactory receptor genes revealed using the 1000 Genomes Project dataset

The molecular mechanism of olfactory cognition is very complicated. Olfactory cognition is initiated by olfactory receptor proteins (odorant receptors), which are activated by olfactory stimuli (ligands). Olfactory receptors are the initial player in the signal transduction cascade producing a nerve impulse, which is transmitted to the brain. The sensitivity to a particular ligand depends on the expression level of multiple proteins involved in the process of olfactory cognition: olfactory receptor proteins, proteins that participate in signal transduction cascade, etc. The expression level of each gene is controlled by its regulatory regions, and especially, by the promoter [a region of DNA about 100–1000 base pairs long located upstream of the transcription start site (TSS)]. We analyzed single nucleotide polymorphisms using human whole-genome data from the 1000 Genomes Project and revealed an extremely high level of single nucleotide polymorphisms in promoter regions of olfactory receptor genes and HLA genes. We hypothesized that the high level of polymorphisms in olfactory receptor promoters was responsible for the diversity in regulatory mechanisms controlling the expression levels of olfactory receptor proteins. Such diversity of regulatory mechanisms may cause the great variability of olfactory cognition of numerous environmental olfactory stimuli perceived by human beings (air pollutants, human body odors, odors in culinary etc.). In turn, this variability may provide a wide range of emotional and behavioral reactions related to the vast variety of olfactory stimuli.

Temporal regulation of Dpp signaling output in the Drosophila wing

BACKGROUND

The Decapentaplegic (Dpp) signaling pathway is used in many developmental and homeostatic contexts, each time resulting in cellular responses particular to that biological niche. The flexibility of Dpp signaling is clearly evident in epithelial cells of the Drosophila wing imaginal disc. During larval stages of development, Dpp functions as a morphogen, patterning the wing developmental field and stimulating tissue growth. A short time later, however, as wing-epithelial cells exit the cell cycle and begin to differentiate, Dpp is a critical determinant of vein-cell fate. It is likely that the Dpp signaling pathway regulates different sets of target genes at these two developmental time points.

RESULTS

To identify mechanisms that temporally control the transcriptional output of Dpp signaling in this system, we have taken a gene expression profiling approach. We identified genes affected by Dpp signaling at late larval or early pupal developmental time points, thereby identifying patterning- and differentiation-specific downstream targets, respectively.

CONCLUSIONS

Analysis of target genes and transcription factor binding sites associated with these groups of genes revealed potential mechanisms by which target-gene specificity of the Dpp signaling pathway is temporally regulated. In addition, this approach revealed novel mechanisms by which Dpp affects the cellular differentiation of wing-veins.

Novel arrangement and comparative analysis of hsp90 family genes in three thermotolerant species of Stratiomyidae (Diptera)

The heat shock proteins belonging to the Hsp90 family (Hsp83 in Diptera) play a crucial role in the protection of cells due to their chaperoning functions. We sequenced hsp90 genes from three species of the family Stratiomyidae (Diptera) living in thermally different habitats and characterized by extraordinarily high thermotolerance. The sequence variation and structure of the hsp90 family genes were compared with previously described features of hsp70 copies isolated from the same species. Two functional hsp83 genes were found in the species studied, that are arranged in tandem orientation at least in one of them. This organization was not previously described. Stratiomyidae hsp83 genes share a high level of identity with hsp83 of Drosophila, and the deduced protein possesses five conserved amino acid sequence motifs characteristic of the Hsp90 family as well as the C-terminus MEEVD sequence characteristic of the cytosolic isoform. A comparison of the hsp83 promoters of two Stratiomyidae species from thermally contrasting habitats demonstrated that while both species contain canonical heat shock elements in the same position, only one of the species contains functional GAF-binding elements. Our data indicate that in the same species, hsp83 family genes show a higher evolution rate than the hsp70 family.

Identification of functionally significant elements in the second intron of the Drosophila melanogaster Trithorax-like gene

It is known that a lot of genes having a distinct expression pattern require the complex system of transcription regulation. The regulatory regions of such genes can include not only the 5'-flanking regions, but also other regions, particularly their intron sequences. The Drosophila melanogaster Trithorax-like (Trl) gene, encoding the GAGA protein, is one of the genes with complex expression pattern. GAGA is one of a few transcription factors that can regulate gene expression at multiple levels. The GAGA-mediated modulation of expression seems to be linked with modifications of the chromatin structure. Nowadays, the regulatory potential of the Trl 5'-flanking region that contains multiple GAGA binding sites has been analyzed, but the presence of the functionally significant elements in other Trl regions has not been examined. We found DNase I hypersensitive sites, evolutionary-conserved sequences and numerous GAGA binding sites in the second intron of the Trl gene. Interestingly, these sequences localize in two main regions of the intron in immediate proximity to preferred regions of transposon insertions. Additionally, we revealed that deletion of the intron fragment in the Trl(1-72) mutants caused an alteration of the Trl expression pattern. These results allow us to conclude that the second intron of the Trl gene contains functionally significant elements.

Regulatory elements of Caenorhabditis elegans ribosomal protein genes

Background

Ribosomal protein genes (RPGs) are essential, tightly regulated, and highly expressed during embryonic development and cell growth. Even though their protein sequences are strongly conserved, their mechanism of regulation is not conserved across yeast, Drosophila, and vertebrates. A recent investigation of genomic sequences conserved across both nematode species and associated with different gene groups indicated the existence of several elements in the upstream regions of C. elegans RPGs, providing a new insight regarding the regulation of these genes in C. elegans.

Results

In this study, we performed an in-depth examination of C. elegans RPG regulation and found nine highly conserved motifs in the upstream regions of C. elegans RPGs using the motif discovery algorithm DME. Four motifs were partially similar to transcription factor binding sites from C. elegans, Drosophila, yeast, and human. One pair of these motifs was found to co-occur in the upstream regions of 250 transcripts including 22 RPGs. The distance between the two motifs displayed a complex frequency pattern that was related to their relative orientation.

We tested the impact of three of these motifs on the expression of rpl-2 using a series of reporter gene constructs and showed that all three motifs are necessary to maintain the high natural expression level of this gene. One of the motifs was similar to the binding site of an orthologue of POP-1, and we showed that RNAi knockdown of pop-1 impacts the expression of rpl-2. We further determined the transcription start site of rpl-2 by 5’ RACE and found that the motifs lie 40–90 bases upstream of the start site. We also found evidence that a noncoding RNA, contained within the outron of rpl-2, is co-transcribed with rpl-2 and cleaved during trans-splicing.

Conclusions

Our results indicate that C. elegans RPGs are regulated by a complex novel series of regulatory elements that is evolutionarily distinct from those of all other species examined up until now.