Large clusters of co-expressed genes in the Drosophila genome

A novel bioinformatics approach identifies candidate genes for the synthesis and feruloylation of arabinoxylan

Arabinoxylans (AXs) are major components of graminaceous plant cell walls, including those in the grain and straw of economically important cereals. Despite some recent advances in identifying the genes encoding biosynthetic enzymes for a number of other plant cell wall polysaccharides, the genes encoding enzymes of the final stages of AX synthesis have not been identified. We have therefore adopted a novel bioinformatics approach based on estimation of differential expression of orthologous genes between taxonomic divisions of species. Over 3 million public domain cereal and dicot expressed sequence tags were mapped onto the complete sets of rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) genes, respectively. It was assumed that genes in cereals involved in AX biosynthesis would be expressed at high levels and that their orthologs in dicotyledonous plants would be expressed at much lower levels. Considering all rice genes encoding putative glycosyl transferases (GTs) predicted to be integral membrane proteins, genes in the GT43, GT47, and GT61 families emerged as much the strongest candidates. When the search was widened to all other rice or Arabidopsis genes predicted to encode integral membrane proteins, cereal genes in Pfam family PF02458 emerged as candidates for the feruloylation of AX. Our analysis, known activities, and recent findings elsewhere are most consistent with genes in the GT43 families encoding beta-1,4-xylan synthases, genes in the GT47 family encoding xylan alpha-1,2- or alpha-1,3-arabinosyl transferases, and genes in the GT61 family encoding feruloyl-AX beta-1,2-xylosyl transferases.

Evidence for de novo evolution of testis-expressed genes in the Drosophila yakuba/Drosophila erecta clade

The mutational origin and subsequent evolution of de novo genes, which are hypothesized to be genes of recent origin that are not obviously related to ancestral coding sequence, are poorly understood. However, accumulating evidence suggests that such genes may often function in male reproduction. Here we use testis-derived expressed sequence tags (ESTs) from Drosophila yakuba to identify genes that have likely arisen either in D. yakuba or in the D. yakuba/D. erecta ancestor. We found several such genes, which show testis-biased expression and are often X-linked. Comparative data indicate that three of these genes have very short open reading frames, which suggests the possibility that a significant number of testis-biased de novo genes in the D. yakuba/D. erecta clade may be noncoding RNA genes. These data, along with previously published data from D. melanogaster, support the idea that many de novo Drosophila genes function in male reproduction and that a small region of the X chromosome in the melanogaster subgroup may be a hotspot for the evolution of novel testis-biased genes.

DI. Up-regulation of the Ku heterodimer in Drosophila testicular cyst cells

In Drosophila, developing germline cysts in testis are enveloped by two somatic cyst cells essential for germline development and male reproduction. The cyst cells continue development along with the germline. However, the mechanisms of somatic gene expression in testes are poorly understood. We report transcriptional up-regulation of the Ku heterodimer in cyst cells. The initial up-regulation is independent of germline, and transcription is further augmented during spermatogenesis. Abundance of Ku in the cyst cell cytoplasm suggests the role for Ku subunits in the regulation of sperm individualization.

The role of translocation and selection in the emergence of genetic clusters and modules

Biomolecular studies point increasingly to the importance of modularity in the organization of the genome. Processes such as the maintenance of metabolism are controlled by suites of genes that act as distinct, self-contained units, or modules. One effect is to promote stability of inherited characters. Despite the obvious importance of genetic modules, the mechanisms by which they form and persist are not understood. One clue is that functionally related genes tend to cluster together. Here we show that genetic translocation, recombination, and natural selection play a central role in this process. We distill the question of emerging genetic modularity into three simulation experiments that show: (1) a tendency, under natural selection, for essential genes to co-locate on the same chromosome and to settle in fixed loci; (2) that genes associated with a particular function tend to form functional clusters; and (3) that genes within a functional cluster tend to become arranged in transcription order. The results also imply that high proportions of junk DNA are essential to the process.

Structure and expression of a novel gene family showing male germline specific expression in Drosophila melanogaster

We report the characterization of two novel genes of Drosophila melanogaster, named mst36Fa and mst36Fb. They define a novel gene family, showing identical time and tissue-specificity limited to male germ cells where their transcription starts during meiotic prophase. These two genes encode for two slightly basic proteins highly homologous to each other and fairly rich in leucine and glutamic acid. Although strictly clustered, these genes utilize different promoter regions as revealed by examination of transgenic flies bearing mst36F-promoter-lacZ reporter constructs and by reverse transcription-polymerase chain reaction assays. Our data suggest that at least one gene (mst36Fa) of the cluster is under translational repression until spermiogenesis suggesting a putative role in the spermatides differentiation. The present study is aimed at the structural analysis of these genes.

Genomic and functional evolution of the Drosophila melanogaster sperm proteome

In addition to delivering a haploid genome to the egg, sperm have additional critical functions, including egg activation, origination of the zygote centrosome and delivery of paternal factors. Despite this, existing knowledge of the molecular basis of sperm form and function is limited. We used whole-sperm mass spectrometry to identify 381 proteins of the Drosophila melanogaster sperm proteome (DmSP). This approach identified mitochondrial, metabolic and cytoskeletal proteins, in addition to several new functional categories. We also observed nonrandom genomic clustering of sperm genes and underrepresentation on the X chromosome. Identification of widespread functional constraint on the proteome indicates that sexual selection has had a limited role in the overall evolution of D. melanogaster sperm. The relevance of the DmSP to the study of mammalian sperm function and fertilization mechanisms is demonstrated by the identification of substantial homology between the DmSP and proteins of the mouse axoneme accessory structure.

Genome-wide patterns of expression in Drosophila pure species and hybrid males. II. Examination of multiple-species hybridizations, platforms, and life cycle stages

Species often produce sterile hybrids early in their evolutionary divergence, and some evidence suggests that hybrid sterility may be associated with deviations or disruptions in gene expression. In support of this idea, many studies have shown that a high proportion of male-biased genes are underexpressed, compared with non-sex-biased genes, in sterile F1 male hybrids of Drosophila species. In this study, we examined and compared patterns of misexpression in sterile F1 male hybrids of Drosophila simulans and 2 of its sibling species, Drosophila mauritiana and Drosophila sechellia, at both the larval and adult life stages. We analyzed hybrids using both commercial Drosophila melanogaster microarrays and arrays we developed from reverse transcriptase-polymerase chain reactions of spermatogenesis and reproduction-related transcripts from these species (sperm array). Although the majority of misexpressed transcripts were underexpressed, a disproportionate number of the overexpressed transcripts were located on the X chromosome. We detected a high overlap in the genes misexpressed between the 2 species pairs, and our sperm array was better at detecting such misexpression than the D. melanogaster array, suggesting possible weaknesses in the use of an array designed from another species. We found only minimal misexpression in the larval samples with the sperm array, suggesting that disruptions in spermatogenesis occur after this life stage. Further study of these misexpressed loci may allow us to identify precisely where disruptions in the spermatogenesis pathway occur.

A global genomic transcriptional code associated with CNS-expressed genes

Highly conserved non-coding DNA regions (HCNR) occur frequently in vertebrate genomes, but their functional roles remain unclear. Here, we provide evidence that a large portion of HCNRs are enriched for binding sites for Sox, POU and Homeodomain transcription factors, and such HCNRs can act as cis-regulatory regions active in neural stem cells. Strikingly, these HCNRs are linked to several hundreds of genes expressed in the developing CNS and they may exert locus-wide regulatory effects on multiple genes flanking their genomic location. Moreover, these data imply a unifying transcriptional logic for a large set of CNS-expressed genes in which Sox and POU proteins act as generic promoters of transcription while Homeodomain proteins control the spatial expression of genes through active repression.

Shared long-range regulatory elements coordinate expression of a gene cluster encoding nicotinic receptor heteromeric subtypes

The nicotinic acetylcholine receptor (nAChR) beta4/alpha3/alpha5 gene cluster encodes several heteromeric transmitter receptor subtypes that are essential for cholinergic synaptic transmission in adrenal gland, autonomic ganglia, pineal gland, and several nuclei in the central nervous system. However, the transcriptional mechanisms coordinating expression of these subunit genes in different cell populations are unknown. Here, we used transgenic methods to investigate long-range transcriptional control of the cluster. A 132-kb P1-derived artificial chromosome (PAC) encoding the rat cluster recapitulated the neurally- and endocrine-restricted expression patterns of the endogenous beta4/alpha3/alpha5 genes. Mutation of ETS factor binding sites in an enhancer, beta43', embedded in the beta4 3'-untranslated exon resulted in greatly diminished beta4, alpha3, and alpha5 expression in adrenal gland and to a lesser extent in the superior cervical ganglion (SCG) but not in other tissues. Phylogenetic sequence analyses revealed several conserved noncoding regions (CNRs) upstream of beta4 and alpha5. Deletion of one of them (CNR4) located 20 kb upstream of beta4 resulted in a dramatic decrease in beta4 and alpha3 expression in the pineal gland and SCG. CNR4 was sufficient to direct LacZ transgene expression to SCG neurons, which express the endogenous beta4alpha3alpha5 subunits, and pineal cells, which express the endogenous beta4alpha3 combination. Finally, CNR4 was able to direct transgene expression to major sites of expression of the endogenous cluster in the brain. Together, our findings support a model in which cell type-specific shared long-range regulatory elements are required for coordinate expression of clustered nAChR genes.

Identification of genomic sites that bind the Drosophila suppressor of Hairy-wing insulator protein

Eukaryotic genomes are divided into independent transcriptional domains by DNA elements known as insulators. The gypsy insulator, a 350-bp element isolated from the Drosophila gypsy retrovirus, contains twelve degenerate binding sites for the Suppressor of Hairy-wing [Su(Hw)] protein. Su(Hw) associates with over 500 non-gypsy genomic sites, the functions of which are largely unknown. Using a bioinformatics approach, we identified 37 putative Su(Hw) insulators (pSIs) that represent regions containing clustered matches to the gypsy insulator Su(Hw) consensus binding sequence. The majority of these pSIs contain fewer than four Su(Hw) binding sites, with only seven showing in vivo Su(Hw) association, as demonstrated by chromatin immunoprecipitation. To understand the properties of the pSIs, these elements were tested for enhancer-blocking capabilities using a transgene assay system. In a complementary set of experiments, effects of the pSIs on transcriptional regulation of genes at the natural genomic location were determined. Our data suggest that pSIs have complex genomic functions and, in some cases, establish insulators. These studies provide the first direct evidence that the Su(Hw) protein contributes to the regulation of gene expression in the Drosophila genome through the establishment of endogenous insulators.

Transcriptome profiling of human hepatocytes treated with Aroclor 1254 reveals transcription factor regulatory networks and clusters of regulated genes

Background

Aroclor 1254 is a well-known hepatotoxin and consists of a complex mixture of polychlorinated biphenyls (PCBs), some of which have the ability to activate the aryl hydrocarbon receptor (AhR) and other transcription factors (TFs). Altered transcription factor expression enables activation of promoters of many genes, thereby inducing a regulatory gene network. In the past, computational approaches were not applied to understand the combinatorial interplay of TFs acting in concert after treatment of human hepatocyte cultures with Aroclor 1254. We were particularly interested in interrogating promoters for transcription factor binding sites of regulated genes.

Results

Here, we present a framework for studying a gene regulatory network and the large-scale regulation of transcription on the level of chromatin structure. For that purpose, we employed cDNA and oligomicroarrays to investigate transcript signatures in human hepatocyte cultures treated with Aroclor 1254 and found 910 genes to be regulated, 52 of which code for TFs and 47 of which are involved in cell cycle and apoptosis. We identified regulatory elements proximal to AhR binding sites, and this included recognition sites for the transcription factors ETS, SP1, CREB, EGR, NF-kB, NKXH, and ZBP. Notably, ECAT and TBP binding sites were identified for Aroclor 1254-induced and E2F, MAZ, HOX, and WHZ for Aroclor 1254-repressed genes. We further examined the chromosomal distribution of regulated genes and observed a statistically significant high number of gene pairs within a distance of 200 kb. Genes regulated by Aroclor 1254, are much closer located to each other than genes distributed randomly all over the genome. 37 regulated gene pairs are even found to be directly neighbored. Within these directly neighbored gene pairs, not all genes were bona fide targets for AhR (primary effect). Upon further analyses many were targets for other transcription factors whose expression was regulated by Aroclor 1254 (secondary effect).

Conclusion

We observed coordinate events in transcript regulation upon treatment of human hepatocytes with Aroclor 1254 and identified a regulatory gene network of different TFs acting in concert. We determined molecular rules for transcriptional regulation to explain, in part, the pleiotropic effect seen in animals and humans upon exposure to Aroclor 1254.

Does lack of recombination enhance asymmetric evolution among duplicate genes? Insights from the Drosophila melanogaster genome

Gene duplication has different outcomes: pseudogenization (death of one of the two copies), gene amplification (both copies remain the same), sub-functionalization (both copies are required to perform the ancestral function) and neo-functionalization (one copy acquires a new function). Asymmetric evolution (one copy evolves faster than the other) is usually seen as a signature of neo-functionalization. However, it has been proposed that sub-functionalization could also generate asymmetric evolution among duplicate genes when they experience different local recombination rates. Indeed, the low recombination copy is expected to evolve faster because of Hill-Robertson effects. Here we tested this idea with about 100 pairs of young duplicates from the Drosophila melanogaster genome. Looking only at young duplicates allowed us to compare recombination rates and evolutionary rates on a similar time-scale contrary to previous work. We found that dispersed pairs tend to evolve more asymmetrically than tandem ones. Among dispersed copies, the low recombination copy tends to be the fast-evolving one. We also tested the possibility that all this was explained by a confounding factor (expression level) but found no evidence for it. In conclusion, our results do support the idea that asymmetric evolution among duplicates is enhanced by restricted recombination. However, further work is needed to clearly distinguish between sub-functionalization and neo-functionalization for the asymmetrically-evolving duplicate pairs that we found.

Characterization of the Drosophila melanogaster genome at the nuclear lamina

The nuclear lamina binds chromatin in vitro and is thought to function in its organization, but genes that interact with it are unknown. Using an in vivo approach, we identified approximately 500 Drosophila melanogaster genes that interact with B-type lamin (Lam). These genes are transcriptionally silent and late replicating, lack active histone marks and are widely spaced. These factors collectively predict lamin binding behavior, indicating that the nuclear lamina integrates variant and invariant chromatin features. Consistently, proximity of genomic regions to the nuclear lamina is partly conserved between cell types, and induction of gene expression or active histone marks reduces Lam binding. Lam target genes cluster in the genome, and these clusters are coordinately expressed during development. This genome-wide analysis gives clear insight into the nature and dynamic behavior of the genome at the nuclear lamina, and implies that intergenic DNA functions in the global organization of chromatin in the nucleus.

Molecular network and chromosomal clustering of genes involved in synaptic plasticity in the hippocampus

Gene transcription is required for establishing and maintaining the enduring form of long term potentiation (LTP). However, the transcriptome and its associated molecular programs that support LTP are not well understood. The purpose of this study was to identify activity-regulated genes (ARGs) and their molecular pathways that are modulated by LTP induction and to investigate the genomic mechanism for coordinating the transcription of ARGs. We performed time course DNA microarray analyses on the mouse dentate gyrus to determine the temporal genomic expression profiles of ARGs in response to LTP-inducing tetanic stimulation. Our studies uncovered ARGs that regulate various cellular processes, including the structure and function of the synapse, and offered an overview of the dynamic molecular programs that are probably important for LTP. Surprisingly, we found that ARGs are clustered on chromosomes, and ARG clusters are conserved during evolution. Although ARGs in the same cluster have apparently different molecular properties, they are functionally correlated by regulating LTP. In addition, ARGs in specific clusters are co-regulated by the cAMP-response element-binding protein. We propose that chromosomal clustering provides a genomic mechanism for coordinating the transcription of ARGs involved in LTP.

Identification of FGF-dependent genes in the Drosophila tracheal system

The embryonic development of the tracheal system of the fruit fly Drosophila provides a paradigm for genetic studies of branching morphogenesis. Efforts of many laboratories have identified Branchless (Bnl, a fibroblast growth factor homologue) and Breathless (Btl, the receptor homologue) as crucial factors at many stages of tracheal system development. The downstream targets of the Bnl/Btl signalling cascade, however, remain mostly unknown. Misexpression of the bnl gene results in specific tracheal phenotypes that lead to larval death. We characterised the transcriptional profiles of targeted over-expression of bnl in the embryonic trachea and of loss-of-function bnl(P1) mutant embryos. Gene expression data was mapped to high-throughput in situ hybridisation based ImaGO-annotation. Thus, we identified and confirmed by quantitative PCR 13 Bnl-dependent genes that are expressed in cells within and outside of the tracheal system.

Regional activation of chromosomal arm 7q with and without gene amplification in taxane-selected human ovarian cancer cell lines

Taxanes are important drugs in the treatment of ovarian and other cancers, but their efficacy is limited by intrinsic and acquired drug resistance. Expression of the multidrug transporter P-glycoprotein, encoded by the MDR1 (ABCB1) gene, is one of the causes of clinical drug resistance to taxanes. To study the mechanisms of MDR1 activation related to taxanes, we established 11 multidrug-resistant variants from six ovarian cancer cell lines by continuous exposure to either paclitaxel or docetaxel. We profiled gene expression and gene copy number alterations in these cell lines using cDNA microarrays and identified a cluster of genes coactivated with MDR1 in 7q21.11-13. Regional activation was evident in nine resistant variants displaying a coexpression pattern of up to 22 genes over an 8-Mb area, including SRI, MGC4175, CLDN12, CROT, and CDK6. In six of these variants, regional activation was driven by gene copy number alterations, with low-level gains or high-level amplifications spanning the involved region. However, three variants displayed regional increases in gene expression even without concomitant gene copy number changes. These results suggest that regional gene activation may be a fundamental mechanism for acquired drug resistance, with or without changes in gene dosage. In addition to numerical and structural chromosomal changes driven by genome instability in cancer cells, other mechanisms might be involved in MDR1 regional activation, such as chromatin remodeling and DNA or histone modifications of the 7q21 region.

Novel genes derived from noncoding DNA in Drosophila melanogaster are frequently X-linked and exhibit testis-biased expression

Descriptions of recently evolved genes suggest several mechanisms of origin including exon shuffling, gene fission/fusion, retrotransposition, duplication-divergence, and lateral gene transfer, all of which involve recruitment of preexisting genes or genetic elements into new function. The importance of noncoding DNA in the origin of novel genes remains an open question. We used the well annotated genome of the genetic model system Drosophila melanogaster and genome sequences of related species to carry out a whole-genome search for new D. melanogaster genes that are derived from noncoding DNA. Here, we describe five such genes, four of which are X-linked. Our RT-PCR experiments show that all five putative novel genes are expressed predominantly in testes. These data support the idea that these novel genes are derived from ancestral noncoding sequence and that new, favored genes are likely to invade populations under selective pressures relating to male reproduction.

A dicistronic gene pair within a cluster of "EF-hand" protein genes in the genomes of Drosophila species

Androcam is a Drosophila melanogaster calmodulin-related protein that functions specifically in the testis. We show that the Acam gene is part of a cluster of three intronless genes arranged in a head-to-tail manner. The additional genes also encode calmodulin-related proteins with testis-specific transcription. Acam and the 5'-most gene (gene1) generate monocistronic transcripts. Surprisingly, the central gene (gene2) is transcribed only as a dicistronic transcript with Acam. A similar cluster is found in D. yakuba. In D. pseudoobscura, the cluster contains four genes: two Acam-type genes downstream of a single gene related to both gene1 and gene2 and a fourth weakly related gene. Nevertheless, the D. pseudoobscura cluster also generates a dicistronic transcript from a gene pair analogous to the gene2-Acam pair. A cotranscribed gene1/2-Acam gene pair may be the founding feature of this locus. Although Acam protein is present in D. melanogaster and D. pseudoobscura testes, cognate proteins for the gene1/2-type ORFs are not detectable by immunoblotting and mass spectrometry techniques.

Transcriptome coexpression map of human embryonic stem cells

Background

Human embryonic stem (ES) cells hold great promise for medicine and science. The transcriptome of human ES cells has been studied in detail in recent years. However, no systematic analysis has yet addressed whether gene expression in human ES cells may be regulated in chromosomal domains, and no chromosomal domains of coexpression have been identified.

Results

We report the first transcriptome coexpression map of the human ES cell and the earliest stage of ES differentiation, the embryoid body (EB), for the analysis of how transcriptional regulation interacts with genomic structure during ES self-renewal and differentiation. We determined the gene expression profiles from multiple ES and EB samples and identified chromosomal domains showing coexpression of adjacent genes on the genome. The coexpression domains were not random, with significant enrichment in chromosomes 8, 11, 16, 17, 19, and Y in the ES state, and 6, 11, 17, 19 and 20 in the EB state. The domains were significantly associated with Giemsa-negative bands in EB, yet showed little correlation with known cytogenetic structures in ES cells. Different patterns of coexpression were revealed by comparative transcriptome mapping between ES and EB.

Conclusion

The findings and methods reported in this investigation advance our understanding of how genome organization affects gene expression in human ES cells and help to identify new mechanisms and pathways controlling ES self-renewal or differentiation.

A method for similarity search of genomic positional expression using CAGE

With the advancement of genome research, it is becoming clear that genes are not distributed on the genome in random order. Clusters of genes distributed at localized genome positions have been reported in several eukaryotes. Various correlations have been observed between the expressions of genes in adjacent or nearby positions along the chromosomes depending on tissue type and developmental stage. Moreover, in several cases, their transcripts, which control epigenetic transcription via processes such as transcriptional interference and genomic imprinting, occur in clusters. It is reasonable that genomic regions that have similar mechanisms show similar expression patterns and that the characteristics of expression in the same genomic regions differ depending on tissue type and developmental stage. In this study, we analyzed gene expression patterns using the cap analysis gene expression (CAGE) method for exploring systematic views of the mouse transcriptome. Counting the number of mapped CAGE tags for fixed-length regions allowed us to determine genomic expression levels. These expression levels were normalized, quantified, and converted into four types of descriptors, allowing the expression patterns along the genome to be represented by character strings. We analyzed them using dynamic programming in the same manner as for sequence analysis. We have developed a novel algorithm that provides a novel view of the genome from the perspective of genomic positional expression. In a similarity search of expression patterns across chromosomes and tissues, we found regions that had clusters of genes that showed expression patterns similar to each other depending on tissue type. Our results suggest the possibility that the regions that have sense–antisense transcription show similar expression patterns between forward and reverse strands.

Through the advancement of genome research, it is becoming clear that genes are not distributed on the genome in random order. Clusters of genes distributed at localized genome positions have been reported in several eukaryotes. Various correlations have been observed between the expressions of genes in adjacent or nearby positions along the chromosomes depending on tissue type and developmental stage. It is reasonable that genomic regions that have similar mechanisms show similar expression patterns. In this study, the authors analyzed gene expression patterns using the computational algorithm of similarity search for exploring systematic views of the mouse transcriptome. They found regions that had clusters of highly expressed genes in certain tissue types whose expression patterns showed strong similarity to each other. This work aims to provide additional insight into genome-wide mechanisms of transcription.

Operon information improves gene expression estimation for cDNA microarrays

BACKGROUND

In prokaryotic genomes, genes are organized in operons, and the genes within an operon tend to have similar levels of expression. Because of co-transcription of genes within an operon, borrowing information from other genes within the same operon can improve the estimation of relative transcript levels; the estimation of relative levels of transcript abundances is one of the most challenging tasks in experimental genomics due to the high noise level in microarray data. Therefore, techniques that can improve such estimations, and moreover are based on sound biological premises, are expected to benefit the field of microarray data analysis

RESULTS

In this paper, we propose a hierarchical Bayesian model, which relies on borrowing information from other genes within the same operon, to improve the estimation of gene expression levels and, hence, the detection of differentially expressed genes. The simulation studies and the analysis of experiential data demonstrated that the proposed method outperformed other techniques that are routinely used to estimate transcript levels and detect differentially expressed genes, including the sample mean and SAM t statistics. The improvement became more significant as the noise level in microarray data increases.

CONCLUSION

By borrowing information about transcriptional activity of genes within classified operons, we improved the estimation of gene expression levels and the detection of differentially expressed genes.

The SuUR gene influences the distribution of heterochromatic proteins HP1 and SU(VAR)3–9 on nurse cell polytene chromosomes of Drosophila melanogaster

We have investigated the distribution of three heterochromatic proteins [SUppressor of UnderReplication (SUUR), heterochromatin protein 1 (HP1), and SU(VAR)3-9] in chromosomes of nurse cells (NCs) and have compared the data obtained with the distribution of the same proteins in salivary gland (SG) chromosomes. In NC chromosomes, the SU(VAR)3-9 protein was found in pericentric heterochromatin and at 223 sites on euchromatic arms, while in SG chromosomes, it was mainly restricted to the chromocenter. In NC chromosomes, the HP1 and SUUR proteins bind to 331 and 256 sites, respectively, which are almost twice the number of sites in SG chromosomes. The distribution of the HP1 and SU(VAR)3-9 proteins depends on the SuUR gene. A mutation in this gene results in a dramatic decrease in the amount of SU(VAR)3-9 binding sites in autosomes. In the X chromosome, these sites are relocated in comparison to the SuUR (+), and their total number only varies slightly. HP1 binding sites are redistributed in chromosomes of SuUR mutants, and their overall number did not change as considerably as SU(VAR)3-9. These data together point to an interaction of these three proteins in Drosophila NC chromosomes.

Remarkable site specificity of local transposition into the Hsp70 promoter of Drosophila melanogaster

Heat-shock genes have numerous features that ought to predispose them to insertional mutagenesis via transposition. To elucidate the evolvability of heat-shock genes via transposition, we have exploited a local transposition technique and Drosophila melanogaster strains with EPgy2 insertions near the Hsp70 gene cluster at 87A7 to produce numerous novel EPgy2 insertions into these Hsp70 genes. More than 50% of 45 independent insertions were made into two adjacent nucleotides in the proximal promoter at positions -96 and -97, and no insertions were into a coding or 3'-flanking sequence. All inserted transposons were in inverse orientation to the starting transposon. The frequent insertion into nucleotides -96 and -97 is consistent with the DNase hypersensitivity, absence of nucleosomes, flanking GAGA-factor-binding sites, and nucleotide sequence of this region. These experimental insertions recapitulated many of the phenotypes of natural transposition into Hsp70: reduced mRNA expression, less Hsp70 protein, and decreased inducible thermotolerance. The results suggest that the distinctive features of heat-shock promoters, which underlie the massive and rapid expression of heat-shock genes upon heat shock, also are a source of evolutionary variation on which natural selection can act.

Clusters of adjacent and similarly expressed genes across normal human tissues complicate comparative transcriptomic discovery

Transcriptomic techniques are valuable tools with which to validate genetic and biological hypotheses and are now widely available for research. However, with the exception of tumor biology, comparative genomics analyses have been difficult to use as discovery engines to describe biologically relevant expression changes. We propose that physical proximity of human genes correlates with similar mRNA expression, so that increased expression might include a disease-relevant gene and many other genes in the adjacent region. To increase the efficiency of combining susceptibility gene mapping and interpretation of transcriptomics, we developed a method to identify clusters of adjacent and similarly expressed genes. Gene expression profiles for 28,945 genes across 101 normal human tissues were obtained from the Gene Logic BioExpress system. The expression similarity for genes in sliding-windows was measured using average pair-wise Pearson correlation coefficients. We identified 187 clusters (p < 10e-4) of co-regulated genes, including 2648 genes, or 9.1% of all genes considered and termed these "clusters of adjacent and similarly expressed genes" (CASEGs). Genes in 15 (8.2%) of these clusters demonstrate a significant co-expression enrichment (p < 10e-10). This study demonstrates the coordinate expression of neighboring genes and provides a comprehensive view of expression-based compartmentalization of the human genome, which can be overlaid on genetic susceptibility gene maps.

Drosophila melanogaster: a case study of a model genomic sequence and its consequences

The sequencing and annotation of the Drosophila melanogaster genome, first published in 2000 through collaboration between Celera Genomics and the Drosophila Genome Projects, has provided a number of important contributions to genome research. By demonstrating the utility of methods such as whole-genome shotgun sequencing and genome annotation by a community "jamboree," the Drosophila genome established the precedents for the current paradigm used by most genome projects. Subsequent releases of the initial genome sequence have been improved by the Berkeley Drosophila Genome Project and annotated by FlyBase, the Drosophila community database, providing one of the highest-quality genome sequences and annotations for any organism. We discuss the impact of the growing number of genome sequences now available in the genus on current Drosophila research, and some of the biological questions that these resources will enable to be solved in the future.

Developmental timing in Dictyostelium is regulated by the Set1 histone methyltransferase

Histone-modifying enzymes have enormous potential as regulators of the large-scale changes in gene expression occurring during differentiation. It is unclear how different combinations of histone modification coordinate regimes of transcription during development. We show that different methylation states of lysine 4 of histone H3 (H3K4) mark distinct developmental phases of the simple eukaryote, Dictyostelium. We demonstrate that the enzyme responsible for all mono, di and tri-methylation of H3K4 is the Dictyostelium homolog of the Set1 histone methyltransferase. In the absence of Set1, cells display unusually rapid development, characterized by precocious aggregation of amoebae into multicellular aggregates. Early differentiation markers are abundantly expressed in growing set1 cells, indicating the differentiation program is ectopically activated during growth. This phenotype is caused specifically by the loss of Set1 catalytic activity. Set1 mutants induce premature differentiation in wild-type cells, indicating Set1 regulates production of an extra-cellular factor required for the correct perception of growth conditions. Microarray analysis of the set1 mutants reveals genomic clustering of mis-expressed genes, suggesting a requirement for Set1 in the regulation of chromatin-mediated events at gene clusters.

Molecular characterisation of the t(1;15)(p22;q22) translocation in the prostate cancer cell line LNCaP

Although chromosome translocations are well-documented recurrent events in hematological malignancies and soft tissue sarcomas, their significance in carcinomas is less clear. We report here the molecular characterization of the reciprocal translocation t(1;15)(p22;q22) in the prostate carcinoma cell line, LNCaP. The chromosome 1 breakpoint was localized to a single BAC clone, RP11-290M5, by sequential FISH analysis of clones selected from the NCBI chromosome 1 map. This was further refined to a 580-bp region by Southern blot analysis. A 2.85-kb fragment spanning the der(1) breakpoint was amplified by long-range inverse PCR. The breakpoint on chromosome 1 was shown to lie between the CYR61 and the DDAH1 genes with the der(1) junctional sequence linking the CYR61 gene to the TSPAN3 (TM4SF8) gene on chromosome 15. Confirmatory PCR and FISH mapping of the der(15) showed loss of chromosome material proximal to the breakpoint on chromosome 15, containing the PSTPIP1 and RCN2 genes. On the available evidence we conclude that this translocation does not result in an in-frame gene fusion. Comparative expressed sequence hybridization (CESH) and comparative genomic hybridization (CGH) analysis, showed relative down-regulation of gene expression surrounding the breakpoint, but no gross change in genomic copy number. Real-time quantitative RT-PCR for genes around the breakpoint supported the CESH data. Therefore, here we may have revealed a gene down-regulation mechanism associated with a chromosome translocation, either through small deletion at the breakpoint or through another means of chromosome domain related gene regulation.

Comprehensive analysis of pathway or functionally related gene expression in the National Cancer Institute's anticancer screen

We have analyzed the level of gene coregulation, using gene expression patterns measured across the National Cancer Institute's 60 tumor cell panels (NCI(60)), in the context of predefined pathways or functional categories annotated by KEGG (Kyoto Encyclopedia of Genes and Genomes), BioCarta, and GO (Gene Ontology). Statistical methods were used to evaluate the level of gene expression coherence (coordinated expression) by comparing intra- and interpathway gene-gene correlations. Our results show that gene expression in pathways, or groups of functionally related genes, has a significantly higher level of coherence than that of a randomly selected set of genes. Transcriptional-level gene regulation appears to be on a "need to be" basis, such that pathways comprising genes encoding closely interacting proteins and pathways responsible for vital cellular processes or processes that are related to growth or proliferation, specifically in cancer cells, such as those engaged in genetic information processing, cell cycle, energy metabolism, and nucleotide metabolism, tend to be more modular (lower degree of gene sharing) and to have genes significantly more coherently expressed than most signaling and regular metabolic pathways. Hierarchical clustering of pathways based on their differential gene expression in the NCI(60) further revealed interesting interpathway communications or interactions indicative of a higher level of pathway regulation. The knowledge of the nature of gene expression regulation and biological pathways can be applied to understanding the mechanism by which small drug molecules interfere with biological systems.

Chromosomal clustering and GATA transcriptional regulation of intestine-expressed genes in C. elegans

We used mRNA tagging to identify genes expressed in the intestine of C. elegans. Animals expressing an epitope-tagged protein that binds the poly-A tail of mRNAs (FLAG::PAB-1) from an intestine-specific promoter (ges-1) were used to immunoprecipitate FLAG::PAB-1/mRNA complexes from the intestine. A total of 1938 intestine-expressed genes (P<0.001) were identified using DNA microarrays. First, we compared the intestine-expressed genes with those expressed in the muscle and germline, and identified 510 genes enriched in all three tissues and 624 intestine-, 230 muscle- and 1135 germ line-enriched genes. Second, we showed that the 1938 intestine-expressed genes were physically clustered on the chromosomes, suggesting that the order of genes in the genome is influenced by the effect of chromatin domains on gene expression. Furthermore, the commonly expressed genes showed more chromosomal clustering than the tissue-enriched genes, suggesting that chromatin domains may influence housekeeping genes more than tissue-specific genes. Third, in order to gain further insight into the regulation of intestinal gene expression, we searched for regulatory motifs. This analysis found that the promoters of the intestine genes were enriched for the GATA transcription factor consensus binding sequence. We experimentally verified these results by showing that the GATA motif is required in cis and that GATA transcription factors are required in trans for expression of these intestinal genes.

Unusual linkage patterns of ligands and their cognate receptors indicate a novel reason for non-random gene order in the human genome

BACKGROUND

Prior to the sequencing of the human genome it was typically assumed that, tandem duplication aside, gene order is for the most part random. Numerous observers, however, highlighted instances in which a ligand was linked to one of its cognate receptors, with some authors suggesting that this may be a general and/or functionally important pattern, possibly associated with recombination modification between epistatically interacting loci. Here we ask whether ligands are more closely linked to their receptors than expected by chance.

RESULTS

We find no evidence that ligands are linked to their receptors more closely than expected by chance. However, in the human genome there are approximately twice as many co-occurrences of ligand and receptor on the same human chromosome as expected by chance. Although a weak effect, the latter might be consistent with a past history of block duplication. Successful duplication of some ligands, we hypothesise, is more likely if the cognate receptor is duplicated at the same time, so ensuring appropriate titres of the two products.

CONCLUSION

While there is an excess of ligands and their receptors on the same human chromosome, this cannot be accounted for by classical models of non-random gene order, as the linkage of ligands/receptors is no closer than expected by chance. Alternative hypotheses for non-random gene order are hence worth considering.

The pattern of chromosome folding in interphase is outlined by the linear gene density profile

Spatial organization of chromatin in the interphase nucleus plays a role in gene expression and inheritance. Although it appears not to be random, the principles of this organization are largely unknown. In this work, we show an explicit relationship between the intranuclear localization of various chromosome segments and the pattern of gene distribution along the genome sequence. Using a 7-megabase-long region of the Drosophila melanogaster chromosome 2 as a model, we observed that the six gene-poor chromosome segments identified in the region interact with components of the nuclear matrix to form a compact stable cluster. The six gene-rich segments form a spatially segregated unstable cluster dependent on nonmatrix nuclear proteins. The resulting composite structure formed by clusters of gene-rich and gene-poor regions is reproducible between the nuclei. We suggest that certain aspects of chromosome folding in interphase are predetermined and can be inferred through in silico analysis of chromosome sequence, using gene density profile as a manifestation of "folding code."

Uncovering evolutionary patterns of gene expression using microarrays

The advent of microarray technology is providing new insights into fundamental questions in evolutionary biology. Here, we review the recent literature on the use of microarrays to study the evolution of genome-wide patterns of gene expression within and between species. Large levels of variation in gene expression patterns have been observed at the intra and interspecific level, and a substantial fraction of transcriptional variation has a genetic component that is contributed by changes in both cis-acting and trans-acting regulatory elements. We argue that there is solid evidence to show that the temporal dynamics of transcriptional variation is largely determined by natural selection, with the fraction of the transcriptome more closely related to sex and reproduction evolving more rapidly.

Chromosomal clustering of a human transcriptome reveals regulatory background

Background

There has been much evidence recently for a link between transcriptional regulation and chromosomal gene order, but the relationship between genomic organization, regulation and gene function in higher eukaryotes remains to be precisely defined.

Results

Here, we present evidence for organization of a large proportion of a human transcriptome into gene clusters throughout the genome, which are partly regulated by the same transcription factors, share biological functions and are characterized by non-housekeeping genes. This analysis was based on the cardiac transcriptome identified by our genome-wide array analysis of 55 human heart samples. We found 37% of these genes to be arranged mainly in adjacent pairs or triplets. A significant number of pairs of adjacent genes are putatively regulated by common transcription factors (p = 0.02). Furthermore, these gene pairs share a significant number of GO functional classification terms. We show that the human cardiac transcriptome is organized into many small clusters across the whole genome, rather than being concentrated in a few larger clusters.

Conclusion

Our findings suggest that genes expressed in concert are organized in a linear arrangement for coordinated regulation. Determining the relationship between gene arrangement, regulation and nuclear organization as well as gene function will have broad biological implications.

Cross-species comparison of Drosophila male accessory gland protein genes

Drosophila melanogaster males transfer seminal fluid proteins along with sperm during mating. Among these proteins, ACPs (Accessory gland proteins) from the male's accessory gland induce behavioral, physiological, and life span reduction in mated females and mediate sperm storage and utilization. A previous evolutionary EST screen in D. simulans identified partial cDNAs for 57 new candidate ACPs. Here we report the annotation and confirmation of the corresponding Acp genes in D. melanogaster. Of 57 new candidate Acp genes previously reported in D. melanogaster, 34 conform to our more stringent criteria for encoding putative male accessory gland extracellular proteins, thus bringing the total number of ACPs identified to 52 (34 plus 18 previously identified). This comprehensive set of Acp genes allows us to dissect the patterns of evolutionary change in a suite of proteins from a single male-specific reproductive tissue. We used sequence-based analysis to examine codon bias, gene duplications, and levels of divergence (via dN/dS values and ortholog detection) of the 52 D. melanogaster ACPs in D. simulans, D. yakuba, and D. pseudoobscura. We show that 58% of the 52 D. melanogaster Acp genes are detectable in D. pseudoobscura. Sequence comparisons of ACPs shared and not shared between D. melanogaster and D. pseudoobscura show that there are separate classes undergoing distinctly dissimilar evolutionary dynamics.

WebGestalt: an integrated system for exploring gene sets in various biological contexts

High-throughput technologies have led to the rapid generation of large-scale datasets about genes and gene products. These technologies have also shifted our research focus from 'single genes' to 'gene sets'. We have developed a web-based integrated data mining system, WebGestalt (http://genereg.ornl.gov/webgestalt/), to help biologists in exploring large sets of genes. WebGestalt is composed of four modules: gene set management, information retrieval, organization/visualization, and statistics. The management module uploads, saves, retrieves and deletes gene sets, as well as performs Boolean operations to generate the unions, intersections or differences between different gene sets. The information retrieval module currently retrieves information for up to 20 attributes for all genes in a gene set. The organization/visualization module organizes and visualizes gene sets in various biological contexts, including Gene Ontology, tissue expression pattern, chromosome distribution, metabolic and signaling pathways, protein domain information and publications. The statistics module recommends and performs statistical tests to suggest biological areas that are important to a gene set and warrant further investigation. In order to demonstrate the use of WebGestalt, we have generated 48 gene sets with genes over-represented in various human tissue types. Exploration of all the 48 gene sets using WebGestalt is available for the public at http://genereg.ornl.gov/webgestalt/wg_enrich.php.

Why repetitive DNA is essential to genome function

There are clear theoretical reasons and many well-documented examples which show that repetitive, DNA is essential for genome function. Generic repeated signals in the DNA are necessary to format expression of unique coding sequence files and to organise additional functions essential for genome replication and accurate transmission to progeny cells. Repetitive DNA sequence elements are also fundamental to the cooperative molecular interactions forming nucleoprotein complexes. Here, we review the surprising abundance of repetitive DNA in many genomes, describe its structural diversity, and discuss dozens of cases where the functional importance of repetitive elements has been studied in molecular detail. In particular, the fact that repeat elements serve either as initiators or boundaries for heterochromatin domains and provide a significant fraction of scaffolding/matrix attachment regions (S/MARs) suggests that the repetitive component of the genome plays a major architectonic role in higher order physical structuring. Employing an information science model, the 'functionalist' perspective on repetitive DNA leads to new ways of thinking about the systemic organisation of cellular genomes and provides several novel possibilities involving repeat elements in evolutionarily significant genome reorganisation. These ideas may facilitate the interpretation of comparisons between sequenced genomes, where the repetitive DNA component is often greater than the coding sequence component.

Global analysis of IL-2 target genes: identification of chromosomal clusters of expressed genes

T lymphocytes play a central role in controlling adaptive immune responses. IL-2 critically regulates both T cell growth and death and is involved in maintaining peripheral tolerance, but the molecules involved in these and other IL-2 actions are only partially known. We now provide a comprehensive compendium of the genes expressed in T cells and of those regulated by IL-2 based on a combination of DNA microarrays and serial analysis of gene expression (SAGE). The newly identified IL-2 target genes include many genes previously linked to apoptosis in other cellular systems that may contribute to IL-2-dependent survival functions. We also studied the mRNA expression of known regulators of signaling pathways for their induction in response to IL-2 in order to identify potential novel positive and/or negative feedback regulators of IL-2 signaling. We show that IL-2 regulates only a limited number of these genes. These include suppressors of cytokine signaling (SOCS) 1, SOCS2, dual-specificity phosphatase (DUSP) 5, DUSP6 and non-receptor type phosphatase-7 (PTPN7). Additionally, we provide evidence that many genes expressed in T cells locate in chromosomal clusters, and that select IL-2-regulated genes are located in at least two clusters, one at 5q31, a known cytokine gene cluster, and the other at 6p21.3, a region that contains genes encoding the tumor necrosis factor (TNF) superfamily members TNF, LT-alpha and LT-beta.

Genomic analysis of Drosophila chromosome underreplication reveals a link between replication control and transcriptional territories

In Drosophila polytene chromosomes, most late-replicating regions remain underreplicated. A loss-of-function mutant of the suppressor of underreplication [Su(UR)] gene suppresses underreplication (UR), whereas extra copies of this gene enhance the level and number of regions showing UR. By combining DNA microarray analysis with manipulation of the number of Su(UR) gene copies, we achieved genomic-scale molecular identification of 1,036 genes that are arranged in clusters located in 52 UR chromosomal regions. These regions overlap extensively (96%) but are not completely identical with late-replicating regions of mitotically dividing Kc cells in culture. Reanalysis of published gene expression profiles revealed that genomic regions defined by replication properties include clusters of coordinately expressed genes. Genomic regions that are UR in polytene chromosomes and late replicated in Kc cell chromosomes show a particularly common association with transcriptional territories that are expressed in testis/males but not ovary/females or embryos. An attractive hypothesis for future testing is that factors involved in replication control, such as SU(UR), may interact physically with those involved in epigenetic silencing of transcription territories.

Chromosomal clustering of genes controlled by the aire transcription factor

Autoimmune regulator (aire) is a transcription factor that controls the self-reactivity of the T cell repertoire. Although previous results indicate that it exerts this function in part by promoting ectopic expression of a battery of peripheral-tissue antigens in epithelial cells of the thymic medulla, recent data argue for additional roles in negative selection of thymocytes by medullary cells. As one approach to exploring such roles, we performed computational analyses of microarray data on medullary RNA transcripts from aire-deficient versus wild-type mice, focusing on the genomic localization of aire-controlled genes. Our results highlight this molecule's transcriptional activating and silencing roles and reveal a significant degree of clustering of its target genes. On a local scale, aire-regulated clusters appeared punctate, with aire-controlled and aire-independent genes often being interspersed. This pattern suggests that aire's action may not be a simple reflection of the wide action of a chromatin remodeling enzyme. Analysis of the identity of certain of the clustered genes was evocative of aire's potential roles in antigen presentation and the coordination of intrathymic cell migration: for example, major histocompatibility complex class I and class II gene products and certain chemokine genes are targets of aire-regulated transcription.

Protein evolution in the context of Drosophila development

The tempo at which a protein evolves depends not only on the rate at which mutations arise but also on the selective effects that those mutations have at the organismal level. It is intuitive that proteins functioning during different stages of development may be predisposed to having mutations of different selective effects. For example, it has been hypothesized that changes to proteins expressed during early development should have larger phenotypic consequences because later stages depend on them. Conversely, changes to proteins expressed much later in development should have smaller consequences at the organismal level. Here we assess whether proteins expressed at different times during Drosophila development vary systematically in their rates of evolution. We find that proteins expressed early in development and particularly during mid-late embryonic development evolve unusually slowly. In addition, proteins expressed in adult males show an elevated evolutionary rate. These two trends are independent of each other and cannot be explained by peculiar rates of mutation or levels of codon bias. Moreover, the observed patterns appear to hold across several functional classes of genes, although the exact developmental time of the slowest protein evolution differs among each class. We discuss our results in connection with data on the evolution of development.

Interspecific chromosome-wide transcription profiles reveal the existence of mammalian-specific and species-specific chromosome domains

A long-range exploration of expression levels through wide chromosome territories was carried out in three species (pig, cattle, and chicken) by aligning EST counts against the human genome. This strategy made it possible to produce expression profiles that were very similar between pig and cattle and that were significantly correlated with chicken levels of expression. In parallel with these alignments, we developed a statistical approach enabling us to screen genomic regions for both underexpression and overexpression at the chromosome level within a given species, as well as interspecifically. The observed correlations are indicative of the existence of interspecifically conserved domains of gene expression, not only for housekeeping genes (which are highly expressed), but also for regions where genes are significantly underexpressed. Furthermore, our strategy made it possible to point out regions that are differentially regulated between species. These expression data were crossed with available comparative mapping information for pigs and cattle, suggesting that coregulated regions are syntenic in various mammals.

Regulated chromatin domain comprising cluster of co-expressed genes in Drosophila melanogaster

Recently, the phenomenon of clustering of co-expressed genes on chromosomes was discovered in eukaryotes. To explore the hypothesis that genes within clusters occupy shared chromatin domains, we performed a detailed analysis of transcription pattern and chromatin structure of a cluster of co-expressed genes. We found that five non-homologous genes (Crtp, Yu, CK2betates, Pros28.1B and CG13581) are expressed exclusively in Drosophila melanogaster male germ-line and form a non-interrupted cluster in the 15 kb region of chromosome 2. The cluster is surrounded by genes with broader transcription patterns. Analysis of DNase I sensitivity revealed 'open' chromatin conformation in the cluster and adjacent regions in the male germ-line cells, where all studied genes are transcribed. In contrast, in somatic tissues where the cluster genes are silent, the domain of repressed chromatin encompassed four out of five cluster genes and an adjacent non-cluster gene CG13589 that is also silent in analyzed somatic tissues. The fifth cluster gene (CG13581) appears to be excluded from the chromatin domain occupied by the other four genes. Our results suggest that extensive clustering of co-expressed genes in eukaryotic genomes does in general reflect the domain organization of chromatin, although domain borders may not exactly correspond to the margins of gene clusters.

Global transcriptome analysis of the C57BL/6J mouse testis by SAGE: evidence for nonrandom gene order

BACKGROUND

We generated the gene expression profile of the total testis from the adult C57BL/6J male mice using serial analysis of gene expression (SAGE). Two high-quality SAGE libraries containing a total of 76 854 tags were constructed. An extensive bioinformatic analysis and comparison of SAGE transcriptomes of the total testis, testicular somatic cells and other mouse tissues was performed and the theory of male-biased gene accumulation on the X chromosome was tested.

RESULTS

We sorted out 829 genes predominantly expressed from the germinal part and 944 genes from the somatic part of the testis. The genes preferentially and specifically expressed in total testis and testicular somatic cells were identified by comparing the testis SAGE transcriptomes to the available transcriptomes of seven non-testis tissues. We uncovered chromosomal clusters of adjacent genes with preferential expression in total testis and testicular somatic cells by a genome-wide search and found that the clusters encompassed a significantly higher number of genes than expected by chance. We observed a significant 3.2-fold enrichment of the proportion of X-linked genes specific for testicular somatic cells, while the proportions of X-linked genes specific for total testis and for other tissues were comparable. In contrast to the tissue-specific genes, an under-representation of X-linked genes in the total testis transcriptome but not in the transcriptomes of testicular somatic cells and other tissues was detected.

CONCLUSION

Our results provide new evidence in favor of the theory of male-biased genes accumulation on the X chromosome in testicular somatic cells and indicate the opposite action of the meiotic X-inactivation in testicular germ cells.

Comparison of the oxidative phosphorylation (OXPHOS) nuclear genes in the genomes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

An analysis of nuclear-encoded oxidative phosphorylation genes in Drosophila and Anopheles reveals that pairs of duplicated genes have strikingly different expression patterns.

Background

In eukaryotic cells, oxidative phosphorylation (OXPHOS) uses the products of both nuclear and mitochondrial genes to generate cellular ATP. Interspecies comparative analysis of these genes, which appear to be under strong functional constraints, may shed light on the evolutionary mechanisms that act on a set of genes correlated by function and subcellular localization of their products.

Results

We have identified and annotated the Drosophila melanogaster, D. pseudoobscura and Anopheles gambiae orthologs of 78 nuclear genes encoding mitochondrial proteins involved in oxidative phosphorylation by a comparative analysis of their genomic sequences and organization. We have also identified 47 genes in these three dipteran species each of which shares significant sequence homology with one of the above-mentioned OXPHOS orthologs, and which are likely to have originated by duplication during evolution. Gene structure and intron length are essentially conserved in the three species, although gain or loss of introns is common in A. gambiae. In most tissues of D. melanogaster and A. gambiae the expression level of the duplicate gene is much lower than that of the original gene, and in D. melanogaster at least, its expression is almost always strongly testis-biased, in contrast to the soma-biased expression of the parent gene.

Conclusions

Quickly achieving an expression pattern different from the parent genes may be required for new OXPHOS gene duplicates to be maintained in the genome. This may be a general evolutionary mechanism for originating phenotypic changes that could lead to species differentiation.

Quantitative evolutionary genomics: differential gene expression and male reproductive success in Drosophila melanogaster

We combined traditional quantitative genetics and oligonucleotide microarrays to examine within-population genetic variation in a trait closely related to fitness. The trait, male reproductive success under competitive conditions (MCRS), is of central importance to both life-history and sexual-selection theory. We identified 27 candidate genes whose expression levels were associated with within-population variation in MCRS. "High" MCRS was associated with low expression of a cytochrome P450 that causes pesticide resistance, suggesting a fitness cost to resistance. Two groups of metabolic proteins (glutathione transferases and phosphatases) were significantly over-represented, and a large portion of the candidates are genes involved in oxidative stress resistance, energy acquisition or energy storage. Genes expressed in accessory glands and testes were not over-represented among differentially expressed genes, but testis-expressed genes were significantly more likely to be upregulated in high MCRS genotypes. Finally, nine candidate genes that we identified had no previous functional annotation, and this experiment suggests that they play a role in male reproductive success.