The genomics of gene expression

Extreme context specificity in differential allelic expression

Variability in cis-regulation of gene expression has been implicated in the phenotypic manifestation of complex traits including common, multifactorial diseases. The differential expression of alleles due to polymorphism in cis-regulatory elements is common in the human genome, but there is a paucity of information about the context specificity of these control elements. In this study, we examined the differential allelic expression (DAE) of BMP5 in human mesenchymal tissues obtained from 16 donors undergoing joint replacement for treatment of osteoarthritis. We observed significant differences in BMP5 allelic output, with allelic ratios greater than 4:1 (P < 10(-20)) in the tissues of some donors. We also discovered a significant variability in allelic expression within the different tissues of donors. For 12 of our donors, we examined the allelic expression of BMP5 in two different regions of cartilage: cartilage adjacent to the site of the osteoarthritic lesion and cartilage distal from the lesion. Five of these 12 donors demonstrated highly significant differences (P < or = 10(-8)) in allelic expression between the different regions of their cartilage. Using DAE as a phenotype, we attempted to map tissue-specific cis-regulatory polymorphisms, and we identified a single nucleotide polymorphism located downstream of BMP5, which was significantly associated with DAE in some but not all of the examined tissues. These findings suggest that allelic expression can be highly context specific and that when interrogating the cis-regulatory control of a particular gene, one cannot necessarily assume that allelic expression is conserved across different tissues or even across different regions of the same tissue.

Identification of putative cis-regulatory elements in Cryptosporidium parvum by de novo pattern finding

Background

Cryptosporidium parvum is a unicellular eukaryote in the phylum Apicomplexa. It is an obligate intracellular parasite that causes diarrhea and is a significant AIDS-related pathogen. Cryptosporidium parvum is not amenable to long-term laboratory cultivation or classical molecular genetic analysis. The parasite exhibits a complex life cycle, a broad host range, and fundamental mechanisms of gene regulation remain unknown. We have used data from the recently sequenced genome of this organism to uncover clues about gene regulation in C. parvum. We have applied two pattern finding algorithms MEME and AlignACE to identify conserved, over-represented motifs in the 5' upstream regions of genes in C. parvum. To support our findings, we have established comparative real-time -PCR expression profiles for the groups of genes examined computationally.

Results

We find that groups of genes that share a function or belong to a common pathway share upstream motifs. Different motifs are conserved upstream of different groups of genes. Comparative real-time PCR studies show co-expression of genes within each group (in sub-sets) during the life cycle of the parasite, suggesting co-regulation of these genes may be driven by the use of conserved upstream motifs.

Conclusion

This is one of the first attempts to characterize cis-regulatory elements in the absence of any previously characterized elements and with very limited expression data (seven genes only). Using de novo pattern finding algorithms, we have identified specific DNA motifs that are conserved upstream of genes belonging to the same metabolic pathway or gene family. We have demonstrated the co-expression of these genes (often in subsets) using comparative real-time-PCR experiments thus establishing evidence for these conserved motifs as putative cis-regulatory elements. Given the lack of prior information concerning expression patterns and organization of promoters in C. parvum we present one of the first investigations of gene regulation in this important human pathogen.

Protein metabolism in marine animals: the underlying mechanism of growth

Growth is a fundamental process within all marine organisms. In soft tissues, growth is primarily achieved by the synthesis and retention of proteins as protein growth. The protein pool (all the protein within the organism) is highly dynamic, with proteins constantly entering the pool via protein synthesis or being removed from the pool via protein degradation. Any net change in the size of the protein pool, positive or negative, is termed protein growth. The three inter-related processes of protein synthesis, degradation and growth are together termed protein metabolism. Measurement of protein metabolism is vital in helping us understand how biotic and abiotic factors affect growth and growth efficiency in marine animals. Recently, the developing fields of transcriptomics and proteomics have started to offer us a means of greatly increasing our knowledge of the underlying molecular control of protein metabolism. Transcriptomics may also allow us to detect subtle changes in gene expression associated with protein synthesis and degradation, which cannot be detected using classical methods. A large literature exists on protein metabolism in animals; however, this chapter concentrates on what we know of marine ectotherms; data from non-marine ectotherms and endotherms are only discussed when the data are of particular relevance. We first consider the techniques available to measure protein metabolism, their problems and what validation is required. Protein metabolism in marine organisms is highly sensitive to a wide variety of factors, including temperature, pollution, seasonality, nutrition, developmental stage, genetics, sexual maturation and moulting. We examine how these abiotic and biotic factors affect protein metabolism at the level of whole-animal (adult and larval), tissue and cellular protein metabolism. Available gene expression data, which help us understand the underlying control of protein metabolism, are also discussed. As protein metabolism appears to comprise a significant proportion of overall metabolic costs in marine organisms, accurate estimates of the energetic cost per unit of synthesised protein are important. Measured costs of protein metabolism are reviewed, and the very high variability in reported costs highlighted. Two major determinants of protein synthesis rates are the tissue concentration of RNA, often expressed as the RNA to protein ratio, and the RNA activity (k(RNA)). The effects of temperature, nutrition and developmental stage on RNA concentration and activity are considered. This chapter highlights our complete lack of knowledge of protein metabolism in many groups of marine organisms, and the fact we currently have only limited data for animals held under a narrow range of experimental conditions. The potential assistance that genomic methods may provide in increasing our understanding of protein metabolism is described.

Free energy of DNA duplex formation on short oligonucleotide microarrays

DNA/DNA duplex formation is the basic mechanism that is used in genome tiling arrays and SNP arrays manufactured by Affymetrix. However, detailed knowledge of the physical process is still lacking. In this study, we show a free energy analysis of DNA/DNA duplex formation these arrays based on the positional-dependent nearest-neighbor (PDNN) model, which was developed previously for describing DNA/RNA duplex formation on expression microarrays. Our results showed that the two ends of a probe contribute less to the stability of the duplexes and that there is a microarray surface effect on binding affinities. We also showed that free energy cost of a single mismatch depends on the bases adjacent to the mismatch site and obtained a comprehensive table of the cost of a single mismatch under all possible combination of adjacent bases. The mismatch costs were found to be correlated with those determined in aqueous solution. We further demonstrate that the DNA copy number estimated from the SNP array correlates negatively with the target length; this is presumably caused by inefficient PCR amplification for long fragments. These results provide important insights into the molecular mechanisms of microarray technology and have implications for microarray design and the interpretation of observed data.

A change of expression in the conserved signaling gene MKK7 is associated with a selective sweep in the western house mouse Mus musculus domesticus

Changes in gene expression are known to occur between closely related species, but it is not yet clear how many of these are due to random fixation of allelic variants or due to adaptive events. In a microarray survey between subspecies of the Mus musculus complex, we identified the mitogen-activated protein-kinase-kinase MKK7 as a candidate for change in gene expression. Quantitative PCR experiments with multiple individuals from each subspecies confirmed a specific and significant up-regulation in the testis of M. m. domesticus. Northern blot analysis shows that this is due to a new transcript that is not found in other tissues, nor in M. m. musculus. A cis-trans test via allele specific expression analysis of the MKK7 gene in F1 hybrids between domesticus and musculus shows that the expression change is mainly caused by a mutation located in cis. Nucleotide diversity was found to be significantly reduced in a window of at least 20 kb around the MKK7 locus in domesticus, indicative of a selective sweep. Because the MKK7 gene is involved in modulating a kinase signalling cascade in a stress response pathway, it seems a plausible target for adaptive differences between subspecies, although the functional role of the new testis-specific transcripts will need to be further studied.

Heritability and tissue specificity of expression quantitative trait loci

Variation in gene expression is heritable and has been mapped to the genome in humans and model organisms as expression quantitative trait loci (eQTLs). We applied integrated genome-wide expression profiling and linkage analysis to the regulation of gene expression in fat, kidney, adrenal, and heart tissues using the BXH/HXB panel of rat recombinant inbred strains. Here, we report the influence of heritability and allelic effect of the quantitative trait locus on detection of cis- and trans-acting eQTLs and discuss how these factors operate in a tissue-specific context. We identified several hundred major eQTLs in each tissue and found that cis-acting eQTLs are highly heritable and easier to detect than trans-eQTLs. The proportion of heritable expression traits was similar in all tissues; however, heritability alone was not a reliable predictor of whether an eQTL will be detected. We empirically show how the use of heritability as a filter reduces the ability to discover trans-eQTLs, particularly for eQTLs with small effects. Only 3% of cis- and trans-eQTLs exhibited large allelic effects, explaining more than 40% of the phenotypic variance, suggestive of a highly polygenic control of gene expression. Power calculations indicated that, across tissues, minor differences in genetic effects are expected to have a significant impact on detection of trans-eQTLs. Trans-eQTLs generally show smaller effects than cis-eQTLs and have a higher false discovery rate, particularly in more heterogeneous tissues, suggesting that small biological variability, likely relating to tissue composition, may influence detection of trans-eQTLs in this system. We delineate the effects of genetic architecture on variation in gene expression and show the sensitivity of this experimental design to tissue sampling variability in large-scale eQTL studies.

Sodium channel gene expression associated with pyrethroid resistant house flies and German cockroaches

The voltage-gated sodium channels of the insect nervous system are the primary target of DDT and pyrethroid insecticides. The loss of target site sensitivity to insecticides resulting from a substitution of leucine to phenylalanine, termed the L-to-F kdr mutation, in the sodium channel of the insect nervous system is known to be important in insecticide resistance. Yet, little is known about the molecular basis underlying the genotype and kdr-mediated resistance phenotype relationship. Here we report a systematic study of resistance-associated kdr allelic expression within and among resistant and susceptible house fly and German cockroach populations. We compared genomic DNA and RNA sequences within the same individuals from different insect strains, finding no correlation for the kdr allele at the genomic DNA level with levels of susceptibility or resistance to insecticide. However, there was a strong correlation between kdr allele expression and the levels of insecticide resistance. This correlation is probably regulated through RNA variation and RNA editing. These results suggest a role for posttranscriptional regulation in the connection of the sodium channel genotype and its mutation-mediated resistance phenotype.

The challenge of Down syndrome

Down syndrome (DS) has been recognized as a clinical entity for about 150 years, but it is only recently that there has been hope for the possibility to understand its pathogenesis and to use this information to devise approaches for the prevention and treatment of its numerous features. The earlier pessimism was due to several reasons, including: (i) the nature of the genetic defect that leads to the syndrome; (ii) the multiplicity of systems involved; and (iii) the high degree of variability of the phenotype. However, science has now caught up with the problem, and recent developments, especially in genetics, genomics, developmental biology and neuroscience, suggest that these potential impediments might not be as arduous as once appeared. As a result, basic research on DS is now rapidly accelerating, and there is hope that the findings will be translatable into benefit for people with DS.

Amplification of low quantity bacterial RNA for microarray studies: time-course analysis of Leptospirillum ferrooxidans under nitrogen-fixing conditions

We have developed a method for the amplification of low quantity total bacterial RNA for DNA microarrays analysis. Current methods are based on the linear amplification by the in vitro transcription from the T7 promoter, similar to that used for eukaryotic mRNA amplification. For the incorporation of T7 promoter, the prokaryotic RNA must be enzymatically modified for the incorporation of a polyA tail at the 3' end to emulate the eukaryotic mRNA. The method we describe and validate herein avoids this step by the direct and random incorporation of the T7 promoter. From 500 ng of total bacterial RNA, we obtained 130-150 microg of antisense RNA, such products being good substrate for fluorescent labelling and DNA microarray analysis. The method was validated with bacterial samples from which it is very difficult to obtain sufficient amounts and quality of total RNA for global gene expression analysis. This is critical for low cell density growing microorganisms, environmental samples, or many extremophiles where the composition of the cultural media severely affects the RNA yield, like in the case of the acidophile and iron oxidizer Gram-negative bacterium Leptospirillum ferrooxidans. We further validated our amplification method in parallel experiments with non-amplified RNA by following the expression of the L. ferrooxidans nif regulon along the time-course of growth.

Variation within and among species in gene expression: raw material for evolution

Heritable variation in regulatory or coding regions is the raw material for evolutionary processes. The advent of microarrays has recently promoted examination of the extent of variation in gene expression within and among taxa and examination of the evolutionary processes affecting variation. This review examines these issues. We find: (i) microarray-based measures of gene expression are precise given appropriate experimental design; (ii) there is large inter-individual variation, which is composed of a minor nongenetic component and a large heritable component; (iii) variation among populations and species appears to be affected primarily by neutral drift and stabilizing selection, and to a lesser degree by directional selection; and (iv) neutral evolutionary divergence in gene expression becomes nonlinear with greater divergence times due to functional constraint. Evolutionary analyses of gene expression reviewed here provide unique insights into partitioning of regulatory variation in nature. However, common limitations of these studies include the tendency to assume a linear relationship between expression divergence and species divergence, and failure to test explicit hypotheses that involve the ecological context of evolutionary divergence.

Genetic variation in human gene expression

Gene expression variation has been the focus of many studies in the past few years. The relevance of gene regulation and gene expression to disease and the development of the technologies used to screen large numbers of genes simultaneously have allowed this rapid development. In this review we discuss issues relating to the biological information one obtains from such studies and the biological significance and use of signals from mapping of gene expression variation.

Genetical genomics: combining genetics with gene expression analysis

The biological mechanisms that link genetic variation and its phenotypic outcome stand as a central puzzle in biology. Geneticists have usually approached this problem by trying to identify genetic variants that underlie the trait in question. Ten years ago, microarray technology opened a second front by making it possible to compare expression levels for most active genes under a variety of genetic and environmental conditions. A typical study reveals up- or down-regulation of genes or pathways associated with a phenotype (case/control) or condition (treated/untreated). In the past few years, a number of groups have started to combine gene expression studies with genetic linkage analysis, leading to a new synergy between these approaches. In this strategy, expression levels are treated as quantitative phenotypes and genetic variants that influence gene expression are sought. Several studies have shown that mRNA levels for many genes are heritable, thus amenable to genetic analysis. Quantitative trait loci mapping efforts have led to the initial characterization of genetic regulation in 'cis' probably because of variants in the gene's own regulatory regions, as well as in 'trans', i.e. by loci elsewhere in the genome. The existence of some 'master regulators' that each affects expression levels of hundreds of genes is an important finding that will surely enrich our understanding of regulatory networks. Although this novel field is still developing, understanding the genetic basis of molecular phenotypes such as gene expression is expected to shed light on the intermediate processes that connect genotype to cellular and organismal traits and represents a critical step towards true systems biology.

Kdr allelic variation in pyrethroid resistant mosquitoes, Culex quinquefasciatus (S.)

The loss of target site sensitivity to insecticides resulting from a substitution of leucine to phenylalanine, termed the kdr mutation, in the voltage-gated sodium channel of the insect nervous system is known to be important in insecticide resistance. However, little is known about the molecular basis of the genotype and kdr-mediated resistance phenotype relationship. This study investigated whether the functional polymorphism of the L-to-F kdr mutation that determines resistance phenotype undergoes DNA variation or goes through transcriptional regulatory variation. We detected no correlation for the kdr allele at the genomic DNA level with levels of susceptibility and resistance to insecticide. However, we find a strong correlation between the kdr allelic expression and levels of insecticide resistance and susceptibility through RNA allelic variation and RNA editing. These findings shed new light on the role of transcriptional regulation in the kdr-mediated resistance in mosquitoes and its connection with the genotype-resistance phenotype relationship.

Epigenetics and phenotypic variation in mammals

What causes phenotypic variation? By now it is clear that phenotype is a result of the interaction between genotype and environment, in addition to variation not readily attributable to either. Epigenetic phenomena associated with phenotypic variation at the biochemical, cellular, tissue, and organism level are now well recognized and are likely to contribute to the "intangible variation" alluded to. While it is clear that epigenetic modifications are mitotically heritable, the fidelity of this process is not well understood. Inheritance through more than one generation of meioses is even less well studied. So it remains unclear to what extent epigenetic changes contribute to phenotypic variation in natural populations. How might such evidence be obtained? What are the features of phenotypes that might suggest an epigenetic component? How much of the epigenetic component is truly independent of genetic changes? The answers to such questions must come from studies designed specifically to detect subtle, stochastically determined phenotypic variation in suitable animal models.

Neutral and adaptive variation in gene expression

Variation among populations in gene expression should be related to the accumulation of random-neutral changes and evolution by natural selection. The following evolutionary analysis has general applicability to biological and medical science because it accounts for genetic relatedness and identifies patterns of expression variation that are affected by natural selection. To identify genes evolving by natural selection, we allocate the maximum among-population variation to genetic distance and then examine the remaining variation relative to a hypothesized important ecological parameter (temperature). These analyses measure the expression of metabolic genes in common-gardened populations of the fish Fundulus heteroclitus whose habitat is distributed along a steep thermal gradient. Although much of the variation in gene expression fits a null model of neutral drift, the variation in expression for 22% of the genes that regress with habitat temperature was far greater than could be accounted for by genetic distance alone. The most parsimonious explanation for among-population variation for these genes is evolution by natural selection. In addition, many metabolic genes have patterns of variation incongruent with neutral evolution: They have too much or too little variation. These patterns of biological variation in expression may reflect important physiological or ecological functions.

Ubiquitous cancer genes: Multipurpose molecules for protein micro-arrays

Multipurpose genes in the human genome which are over-expressed in a large variety of different cancers have been identified. Forty-two of the 19,016 human genes annotated to date (0.2%) are ubiquitously over-expressed in half or more of the 36 investigated human cancers. Of these genes, 15 are involved in protein biosynthesis and folding, six of them in glycolysis. A group of 13 solid tumours over-express almost all (39-42 of 42) ubiquitous cancer genes, suggesting a common mechanism underlying these cancers. Others, such as endocrine cancers, have only a few over-expressed ubiquitous cancer genes. The proteins for which these genes code or the corresponding antibodies are candidates for small protein microarrays aiming at maximum information with only a limited number of proteins. Since the over-expression pattern varies from cancer to cancer, distinction between different cancer classes is possible using one single set of protein or antibody molecules.

The transcriptome of Toxoplasma gondii

BACKGROUND

Toxoplasma gondii gives rise to toxoplasmosis, among the most prevalent parasitic diseases of animals and man. Transformation of the tachzyoite stage into the latent bradyzoite-cyst form underlies chronic disease and leads to a lifetime risk of recrudescence in individuals whose immune system becomes compromised. Given the importance of tissue cyst formation, there has been intensive focus on the development of methods to study bradyzoite differentiation, although the molecular basis for the developmental switch is still largely unknown.

RESULTS

We have used serial analysis of gene expression (SAGE) to define the Toxoplasma gondii transcriptome of the intermediate-host life cycle that leads to the formation of the bradyzoite/tissue cyst. A broad view of gene expression is provided by >4-fold coverage from nine distinct libraries (approximately 300,000 SAGE tags) representing key developmental transitions in primary parasite populations and in laboratory strains representing the three canonical genotypes. SAGE tags, and their corresponding mRNAs, were analyzed with respect to abundance, uniqueness, and antisense/sense polarity and chromosome distribution and developmental specificity.

CONCLUSION

This study demonstrates that phenotypic transitions during parasite development were marked by unique stage-specific mRNAs that accounted for 18% of the total SAGE tags and varied from 1-5% of the tags in each developmental stage. We have also found that Toxoplasma mRNA pools have a unique parasite-specific composition with 1 in 5 transcripts encoding Apicomplexa-specific genes functioning in parasite invasion and transmission. Developmentally co-regulated genes were dispersed across all Toxoplasma chromosomes, as were tags representing each abundance class, and a variety of biochemical pathways indicating that trans-acting mechanisms likely control gene expression in this parasite. We observed distinct similarities in the specificity and expression levels of mRNAs in primary populations (Day-6 post-sporozoite infection) that occur prior to the onset of bradyzoite development that were uniquely shared with the virulent Type I-RH laboratory strain suggesting that development of RH may be arrested. By contrast, strains from Type II-Me49B7 and Type III-VEGmsj contain SAGE tags corresponding to bradyzoite genes, which suggests that priming of developmental expression likely plays a role in the greater capacity of these strains to complete bradyzoite development.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

The quantitative genetics of transcription

Quantitative geneticists have become interested in the heritability of transcription and detection of expression quantitative trait loci (eQTLs). Linkage mapping methods have identified major-effect eQTLs for some transcripts and have shown that regulatory polymorphisms in cis and in trans affect expression. It is also clear that these mapping strategies have little power to detect polygenic factors, and some new statistical approaches are emerging that paint a more complex picture of transcriptional heritability. Several studies imply pervasive non-additivity of transcription, transgressive segregation and epistasis, and future studies will soon document the extent of genotype-environment interaction and population structure at the transcriptional level. The implications of these findings for genotype-phenotype mapping and modeling the evolution of transcription are discussed.

Allelic Expression Imbalance of Human mu Opioid Receptor (OPRM1) Caused by Variant A118G

As a primary target for opioid drugs and peptides, the mu opioid receptor (OPRM1) plays a key role in pain perception and addiction. Genetic variants of OPRM1 have been implicated in predisposition to drug addiction, in particular the single nucleotide polymorphism A118G, leading to an N40D substitution, with an allele frequency of 10-32%, and uncertain functions. We have measured allele-specific mRNA expression of OPRM1 in human autopsy brain tissues, using A118G as a marker. In 8 heterozygous samples measured, the A118 mRNA allele was 1.5-2.5-fold more abundant than the G118 allele. Transfection into Chinese hamster ovary cells of a cDNA representing only the coding region of OPRM1, carrying adenosine, guanosine, cytidine, and thymidine in position 118, resulted in 1.5-fold lower mRNA levels only for OPRM1-G118, and more than 10-fold lower OPRM1 protein levels, measured by Western blotting and receptor binding assay. After transfection and inhibition of transcription with actinomycin D, analysis of mRNA turnover failed to reveal differences in mRNA stability between A118 and G118 alleles, indicating a defect in transcription or mRNA maturation. These results indicate that OPRM1-G118 is a functional variant with deleterious effects on both mRNA and protein yield. Clarifying the functional relevance of polymorphisms associated with susceptibility to a complex disorder such as drug addiction provides a foundation for clinical association studies.

Challenges for identifying functionally important genetic variation: the promise of combining complementary research strategies

Strategies for the identification of functional genetic variation underlying phenotypic traits of ecological and evolutionary importance have received considerable attention in the literature recently. This paper aims to bring together and compare the relative strengths and limitations of various potentially useful research strategies for dissecting functionally important genetic variation in a wide range of organisms. We briefly explore the relative strengths and limitations of traditional and emerging approaches and evaluate their potential use in free-living populations. While it is likely that much of the progress in functional genetic analyses will rely on progress in traditional model species, it is clear that with prudent choices of methods and appropriate sampling designs, much headway can be also made in a diverse range of species. We suggest that combining research approaches targeting different functional and biological levels can potentially increase understanding the genetic basis of ecological and evolutionary processes both in model and non-model organisms.

Reverse transcriptional profiling: non-correspondence of transcript level variation and proximal promoter polymorphism

BACKGROUND

Variation in gene expression between two Drosophila melanogaster strains, as revealed by transcriptional profiling, seldom corresponded to variation in proximal promoter sequence for 34 genes analyzed. Two sets of protein-coding genes were selected from pre-existing microarray data: (1) those whose expression varied significantly and reproducibly between strains, and (2) those whose transcript levels did not vary. Only genes whose regulation of expression was uncharacterized were chosen. At least one kB of the proximal promoters of 15-19 genes in each set was sequenced and compared between strains (Oregon R and Russian 2b).

RESULTS

Of the many promoter polymorphisms, 89.6% were SNPs and 10.4% were indels, including homopolymer tracts, microsatellite repeats, and putative transposable element footprints. More than half of the SNPs were changes within a nucleotide class. Hypothetically, genes differing in expression between the two strains should have more proximal promoter polymorphisms than those whose expression is similar. The number, frequency, and type of polymorphism, however, were the same in both sets of genes. In fact, the promoters of six genes with significantly different mRNA expression were identical in sequence.

CONCLUSION

For these genes, sequences external to the proximal promoter, such as enhancers or in trans, must play a greater role than the proximal promoter in transcriptomic variation between D. melanogaster strains.