Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia

Despite its morphological similarity to the other species in the Drosophila melanogaster species complex, D. sechellia has evolved distinct physiological and behavioral adaptations to its host plant Morinda citrifolia, commonly known as Tahitian Noni. The odor of the ripe fruit of M. citrifolia originates from hexanoic and octanoic acid. D. sechellia is attracted to these two fatty acids, whereas the other species in the complex are repelled. Here, using interspecies hybrids between D. melanogaster deficiency mutants and D. sechellia, we showed that the Odorant-binding protein 57e (Obp57e) gene is involved in the behavioral difference between the species. D. melanogaster knock-out flies for Obp57e and Obp57d showed altered behavioral responses to hexanoic acid and octanoic acid. Furthermore, the introduction of Obp57d and Obp57e from D. simulans and D. sechellia shifted the oviposition site preference of D. melanogaster Obp57d/e^KO flies to that of the original species, confirming the contribution of these genes to D. sechellia's specialization to M. citrifolia. Our finding of the genes involved in host-plant determination may lead to further understanding of mechanisms underlying taste perception, evolution of plant–herbivore interactions, and speciation.     

A screen for recessive speciation genes expressed in the gametes of F1 hybrid yeast

Diploid hybrids of Saccharomyces cerevisiae and its closest relative, Saccharomyces paradoxus, are viable, but the sexual gametes they produce are not. One of several possible causes of this gamete inviability is incompatibility between genes from different species—such incompatible genes are usually called “speciation genes.” In diploid F1 hybrids, which contain a complete haploid genome from each species, the presence of compatible alleles can mask the effects of (recessive) incompatible speciation genes. But in the haploid gametes produced by F1 hybrids, recessive speciation genes may be exposed, killing the gametes and thus preventing F1 hybrids from reproducing sexually. Here I present the results of an experiment to detect incompatibilities that kill hybrid gametes. I transferred nine of the 16 S. paradoxus chromosomes individually into S. cerevisiae gametes and tested the ability of each to replace its S. cerevisiae homeolog. All nine chromosomes were compatible, producing nine viable haploid strains, each with 15 S. cerevisiae chromosomes and one S. paradoxus chromosome. Thus, none of these chromosomes contain speciation genes that were capable of killing the hybrid gametes that received them. This is a surprising result that suggests that such speciation genes do not play a major role in yeast speciation.     

Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus

Interspecific hybrid lethality and sterility are a consequence of divergent evolution between species and serve to maintain the discrete identities of species. The evolution of hybrid incompatibilities has been described in widely accepted models by Dobzhansky and Muller where lineage-specific functional divergence is the essential characteristic of hybrid incompatibility genes. Experimentally tractable models are required to identify and test candidate hybrid incompatibility genes. Several Drosophila melanogaster genes involved in hybrid incompatibility have been identified but none has yet been shown to have functionally diverged in accordance with the Dobzhansky-Muller model. By introducing transgenic copies of the X-linked Hybrid male rescue (Hmr) gene into D. melanogaster from its sibling species D. simulans and D. mauritiana, we demonstrate that Hmr has functionally diverged to cause F1 hybrid incompatibility between these species. Consistent with the Dobzhansky-Muller model, we find that Hmr has diverged extensively in the D. melanogaster lineage, but we also find extensive divergence in the sibling-species lineage. Together, these findings implicate over 13% of the amino acids encoded by Hmr as candidates for causing hybrid incompatibility. The exceptional level of divergence at Hmr cannot be explained by neutral processes because we use phylogenetic methods and population genetic analyses to show that the elevated amino-acid divergence in both lineages is due to positive selection in the distant past—at least one million generations ago. Our findings suggest that multiple substitutions driven by natural selection may be a general phenomenon required to generate hybrid incompatibility alleles.     

The parasexual cycle in Candida albicans provides an alternative pathway to meiosis for the formation of recombinant strains

Candida albicans has an elaborate, yet efficient, mating system that promotes conjugation between diploid a and α strains. The product of mating is a tetraploid a/α cell that must undergo a reductional division to return to the diploid state. Despite the presence of several “meiosis-specific” genes in the C. albicans genome, a meiotic program has not been observed. Instead, tetraploid products of mating can be induced to undergo efficient, random chromosome loss, often producing strains that are diploid, or close to diploid, in ploidy. Using SNP and comparative genome hybridization arrays we have now analyzed the genotypes of products from the C. albicans parasexual cycle. We show that the parasexual cycle generates progeny strains with shuffled combinations of the eight C. albicans chromosomes. In addition, several isolates had undergone extensive genetic recombination between homologous chromosomes, including multiple gene conversion events. Progeny strains exhibited altered colony morphologies on laboratory media, demonstrating that the parasexual cycle generates phenotypic variants of C. albicans. In several fungi, including Saccharomyces cerevisiae and Schizosaccharomyces pombe, the conserved Spo11 protein is integral to meiotic recombination, where it is required for the formation of DNA double-strand breaks. We show that deletion of SPO11 prevented genetic recombination between homologous chromosomes during the C. albicans parasexual cycle. These findings suggest that at least one meiosis-specific gene has been re-programmed to mediate genetic recombination during the alternative parasexual life cycle of C. albicans. We discuss, in light of the long association of C. albicans with warm-blooded animals, the potential advantages of a parasexual cycle over a conventional sexual cycle.     

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

Large-scale trends in the evolution of gene structures within 11 animal genomes

We have used the annotations of six animal genomes (Homo sapiens, Mus musculus, Ciona intestinalis, Drosophila melanogaster, Anopheles gambiae, and Caenorhabditis elegans) together with the sequences of five unannotated Drosophila genomes to survey changes in protein sequence and gene structure over a variety of timescales—from the less than 5 million years since the divergence of D. simulans and D. melanogaster to the more than 500 million years that have elapsed since the Cambrian explosion. To do so, we have developed a new open-source software library called CGL (for “Comparative Genomics Library”). Our results demonstrate that change in intron–exon structure is gradual, clock-like, and largely independent of coding-sequence evolution. This means that genome annotations can be used in new ways to inform, corroborate, and test conclusions drawn from comparative genomics analyses that are based upon protein and nucleotide sequence similarities.     

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Conserved regulation of MAP kinase expression by PUF RNA-binding proteins

Mitogen-activated protein kinase (MAPK) and PUF (for Pumilio and FBF [fem-3 binding factor]) RNA-binding proteins control many cellular processes critical for animal development and tissue homeostasis. In the present work, we report that PUF proteins act directly on MAPK/ERK-encoding mRNAs to downregulate their expression in both the Caenorhabditis elegans germline and human embryonic stem cells. In C. elegans, FBF/PUF binds regulatory elements in the mpk-1 3′ untranslated region (3′ UTR) and coprecipitates with mpk-1 mRNA; moreover, mpk-1 expression increases dramatically in FBF mutants. In human embryonic stem cells, PUM2/PUF binds 3′UTR elements in both Erk2 and p38α mRNAs, and PUM2 represses reporter constructs carrying either Erk2 or p38α 3′ UTRs. Therefore, the PUF control of MAPK expression is conserved. Its biological function was explored in nematodes, where FBF promotes the self-renewal of germline stem cells, and MPK-1 promotes oocyte maturation and germ cell apoptosis. We found that FBF acts redundantly with LIP-1, the C. elegans homolog of MAPK phosphatase (MKP), to restrict MAPK activity and prevent apoptosis. In mammals, activated MAPK can promote apoptosis of cancer cells and restrict stem cell self-renewal, and MKP is upregulated in cancer cells. We propose that the dual negative regulation of MAPK by both PUF repression and MKP inhibition may be a conserved mechanism that influences both stem cell maintenance and tumor progression.     

Substitution as a mechanism for genetic robustness: the duplicated deacetylases Hst1p and Sir2p in Saccharomyces cerevisiae

How duplicate genes provide genetic robustness remains an unresolved question. We have examined the duplicated histone deacetylases Sir2p and Hst1p in Saccharomyces cerevisiae and find that these paralogs with non-overlapping functions can provide genetic robustness against null mutations through a substitution mechanism. Hst1p is an NAD⁺-dependent histone deacetylase that acts with Sum1p to repress a subset of midsporulation genes. However, hst1Δ mutants show much weaker derepression of target loci than sum1Δ mutants. We show that this modest derepression of target loci in hst1Δ strains occurs in part because Sir2p substitutes for Hst1p. Sir2p contributes to repression of the midsporulation genes only in the absence of Hst1p and is recruited to target promoters by a physical interaction with the Sum1 complex. Furthermore, when Sir2p associates with the Sum1 complex, the complex continues to repress in a promoter-specific manner and does not spread. Our results imply that after the duplication, SIR2 and HST1 subfunctionalized. The single SIR2/HST1 gene from Kluyveromyces lactis, a closely related species that diverged prior to the duplication, can suppress an hst1Δ mutation in S. cerevisiae as well as interact with Sir4p in S. cerevisiae. In addition, the existence of two distinct protein interaction domains for the Sir and Sum1 complexes was revealed through the analysis of a chimeric Sir2–Hst1 molecule. Therefore, the ability of Sir2p to substitute for Hst1p probably results from a retained but reduced affinity for the Sum1 complex that is a consequence of subfunctionalization via the duplication, degeneration, and complementation mechanism. These results suggest that the evolutionary path of duplicate gene preservation may be an important indicator for the ability of duplicated genes to contribute to genetic robustness.     

Complex genetic interactions in a quantitative trait locus

Whether in natural populations or between two unrelated members of a species, most phenotypic variation is quantitative. To analyze such quantitative traits, one must first map the underlying quantitative trait loci. Next, and far more difficult, one must identify the quantitative trait genes (QTGs), characterize QTG interactions, and identify the phenotypically relevant polymorphisms to determine how QTGs contribute to phenotype. In this work, we analyzed three Saccharomyces cerevisiae high-temperature growth (Htg) QTGs (MKT1, END3, and RHO2). We observed a high level of genetic interactions among QTGs and strain background. Interestingly, while the MKT1 and END3 coding polymorphisms contribute to phenotype, it is the RHO2 3′UTR polymorphisms that are phenotypically relevant. Reciprocal hemizygosity analysis of the Htg QTGs in hybrids between S288c and ten unrelated S. cerevisiae strains reveals that the contributions of the Htg QTGs are not conserved in nine other hybrids, which has implications for QTG identification by marker-trait association. Our findings demonstrate the variety and complexity of QTG contributions to phenotype, the impact of genetic background, and the value of quantitative genetic studies in S. cerevisiae.     

Asymmetrical reinforcement and Wolbachia infection in Drosophila

Reinforcement refers to the evolution of increased mating discrimination against heterospecific individuals in zones of geographic overlap and can be considered a final stage in the speciation process. One the factors that may affect reinforcement is the degree to which hybrid matings result in the permanent loss of genes from a species' gene pool. Matings between females of Drosophila subquinaria and males of D. recens result in high levels of offspring mortality, due to interspecific cytoplasmic incompatibility caused by Wolbachia infection of D. recens. Such hybrid inviability is not manifested in matings between D. recens females and D. subquinaria males. Here we ask whether the asymmetrical hybrid inviability is associated with a corresponding asymmetry in the level of reinforcement. The geographic ranges of D. recens and D. subquinaria were found to overlap across a broad belt of boreal forest in central Canada. Females of D. subquinaria from the zone of sympatry exhibit much stronger levels of discrimination against males of D. recens than do females from allopatric populations. In contrast, such reproductive character displacement is not evident in D. recens, consistent with the expected effects of unidirectional cytoplasmic incompatibility. Furthermore, there is substantial behavioral isolation within D. subquinaria, because females from populations sympatric with D. recens discriminate against allopatric conspecific males, whereas females from populations allopatric with D. recens show no discrimination against any conspecific males. Patterns of general genetic differentiation among populations are not consistent with patterns of behavioral discrimination, which suggests that the behavioral isolation within D. subquinaria results from selection against mating with Wolbachia-infected D. recens. Interspecific cytoplasmic incompatibility may contribute not only to post-mating isolation, an effect already widely recognized, but also to reinforcement, particularly in the uninfected species. The resulting reproductive character displacement not only increases behavioral isolation from the Wolbachia-infected species, but may also lead to behavioral isolation between populations of the uninfected species. Given the widespread occurrence of Wolbachia among insects, it thus appears that there are multiple ways by which these endosymbionts may directly and indirectly contribute to reproductive isolation and speciation.     

The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism

A central question in evolutionary biology concerns the developmental processes by which new phenotypes arise. An exceptional example of evolutionary innovation is the single-celled seed trichome in Gossypium (“cotton fiber”). We have used fiber development in Gossypium as a system to understand how morphology can rapidly evolve. Fiber has undergone considerable morphological changes between the short, tightly adherent fibers of G. longicalyx and the derived long, spinnable fibers of its closest relative, G. herbaceum, which facilitated cotton domestication. We conducted comparative gene expression profiling across a developmental time-course of fibers from G. longicalyx and G. herbaceum using microarrays with ~22,000 genes. Expression changes between stages were temporally protracted in G. herbaceum relative to G. longicalyx, reflecting a prolongation of the ancestral developmental program. Gene expression and GO analyses showed that many genes involved with stress responses were upregulated early in G. longicalyx fiber development. Several candidate genes upregulated in G. herbaceum have been implicated in regulating redox levels and cell elongation processes. Three genes previously shown to modulate hydrogen peroxide levels were consistently expressed in domesticated and wild cotton species with long fibers, but expression was not detected by quantitative real time-PCR in wild species with short fibers. Hydrogen peroxide is important for cell elongation, but at high concentrations it becomes toxic, activating stress processes that may lead to early onset of secondary cell wall synthesis and the end of cell elongation. These observations suggest that the evolution of long spinnable fibers in cotton was accompanied by novel expression of genes assisting in the regulation of reactive oxygen species levels. Our data suggest a model for the evolutionary origin of a novel morphology through differential gene regulation causing prolongation of an ancestral developmental program.     

Natural biocombinatorics in the polyketide synthase genes of the actinobacterium Streptomyces avermitilis

Modular polyketide synthases (PKSs) of bacteria provide an enormous reservoir of natural chemical diversity. Studying natural biocombinatorics may aid in the development of concepts for experimental design of genes for the biosynthesis of new bioactive compounds. Here we address the question of how the modularity of biosynthetic enzymes and the prevalence of multiple gene clusters in Streptomyces drive the evolution of metabolic diversity. The phylogeny of ketosynthase (KS) domains of Streptomyces PKSs revealed that the majority of modules involved in the biosynthesis of a single compound evolved by duplication of a single ancestor module. Using Streptomyces avermitilis as a model organism, we have reconstructed the evolutionary relationships of different domain types. This analysis suggests that 65% of the modules were altered by recombinational replacements that occurred within and between biosynthetic gene clusters. The natural reprogramming of the biosynthetic pathways was unambiguously confined to domains that account for the structural diversity of the polyketide products and never observed for the KS domains. We provide examples for natural acyltransferase (AT), ketoreductase (KR), and dehydratase (DH)–KR domain replacements. Potential sites of homologous recombination could be identified in interdomain regions and within domains. Our results indicate that homologous recombination facilitated by the modularity of PKS architecture is the most important mechanism underlying polyketide diversity in bacteria.     

Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting

The phylogenetic relationship of the now fully sequenced species Drosophila erecta and D. yakuba with respect to the D. melanogaster species complex has been a subject of controversy. All three possible groupings of the species have been reported in the past, though recent multi-gene studies suggest that D. erecta and D. yakuba are sister species. Using the whole genomes of each of these species as well as the four other fully sequenced species in the subgenus Sophophora, we set out to investigate the placement of D. erecta and D. yakuba in the D. melanogaster species group and to understand the cause of the past incongruence. Though we find that the phylogeny grouping D. erecta and D. yakuba together is the best supported, we also find widespread incongruence in nucleotide and amino acid substitutions, insertions and deletions, and gene trees. The time inferred to span the two key speciation events is short enough that under the coalescent model, the incongruence could be the result of incomplete lineage sorting. Consistent with the lineage-sorting hypothesis, substitutions supporting the same tree were spatially clustered. Support for the different trees was found to be linked to recombination such that adjacent genes support the same tree most often in regions of low recombination and substitutions supporting the same tree are most enriched roughly on the same scale as linkage disequilibrium, also consistent with lineage sorting. The incongruence was found to be statistically significant and robust to model and species choice. No systematic biases were found. We conclude that phylogenetic incongruence in the D. melanogaster species complex is the result, at least in part, of incomplete lineage sorting. Incomplete lineage sorting will likely cause phylogenetic incongruence in many comparative genomics datasets. Methods to infer the correct species tree, the history of every base in the genome, and comparative methods that control for and/or utilize this information will be valuable advancements for the field of comparative genomics.     

Using likelihood-free inference to compare evolutionary dynamics of the protein networks of H. pylori and P. falciparum

Gene duplication with subsequent interaction divergence is one of the primary driving forces in the evolution of genetic systems. Yet little is known about the precise mechanisms and the role of duplication divergence in the evolution of protein networks from the prokaryote and eukaryote domains. We developed a novel, model-based approach for Bayesian inference on biological network data that centres on approximate Bayesian computation, or likelihood-free inference. Instead of computing the intractable likelihood of the protein network topology, our method summarizes key features of the network and, based on these, uses a MCMC algorithm to approximate the posterior distribution of the model parameters. This allowed us to reliably fit a flexible mixture model that captures hallmarks of evolution by gene duplication and subfunctionalization to protein interaction network data of Helicobacter pylori and Plasmodium falciparum. The 80% credible intervals for the duplication–divergence component are [0.64, 0.98] for H. pylori and [0.87, 0.99] for P. falciparum. The remaining parameter estimates are not inconsistent with sequence data. An extensive sensitivity analysis showed that incompleteness of PIN data does not largely affect the analysis of models of protein network evolution, and that the degree sequence alone barely captures the evolutionary footprints of protein networks relative to other statistics. Our likelihood-free inference approach enables a fully Bayesian analysis of a complex and highly stochastic system that is otherwise intractable at present. Modelling the evolutionary history of PIN data, it transpires that only the simultaneous analysis of several global aspects of protein networks enables credible and consistent inference to be made from available datasets. Our results indicate that gene duplication has played a larger part in the network evolution of the eukaryote than in the prokaryote, and suggests that single gene duplications with immediate divergence alone may explain more than 60% of biological network data in both domains.     

The genetic system controlling recombination in the silkworm

The nature of recombination modifiers was investigated in Bombyx mori lines selected for high (H) and low (L) recombination rates between the p^S and Y loci in chromosome 2. Since the mean recombination rates for the H x L and L x H F₁ crosses were approximately intermediate between those of high and low lines, the cytoplasmic maternal effect and difference in the activity of recombination modifiers between marked and unmarked second chromosomes were not detected. The H x (L x H), H x (H x L), L x (L x H) and L x (H x L) backcrosses indicated the presence of additive and dominance effects of marked and unmarked second chromosomes and the remaining chromosomes.——Recombination rates between the p^S and Y loci in chromosome 2 and half-nonrecombination rates between the pe and re loci in chromosome 5 of high and low lines indicated that these recombination modifiers caused changes in the recombination frequency between p^S and Y in chromosome 2, but not between pe and re in chromosome 5.——There were no differences in viability between individuals having the second chromosomes of the recombinant types [p^S +, p^Y (H); p^S +, + Y (L)] and those of the nonrecombinant types [p^S Y, p + (H); p^S Y, + + (L)] in both high and low lines. Mean recombination rates measured in cis [p^S Y/p + (H); p^S Y/+ + (L)] and trans [p^S +/p Y (H); p^S +/+ Y (L)] males were the same in the high but not in the low line. No segregation of a single recombination modifier was indicated by the distribution of recombination rates measured in trans males [p^S +/p Y (H); p^S +/+ Y (L)] of high and low lines. Accordingly, the recombination modifiers distributed on chromosome 2 in the heterozygous condition were not gross chromosomal aberrations, but polygenic factors in the low line.     

Evolutionary Changes in Gene Expression,Coding Sequence and Copy-Number at the Cyp6g1 Locus Contribute to Resistance to Multiple Insecticides in Drosophila

Widespread use of insecticides has led to insecticide resistance in many populations of insects. In some populations, resistance has evolved to multiple pesticides. In Drosophila melanogaster, resistance to multiple classes of insecticide is due to the overexpression of a single cytochrome P450 gene, Cyp6g1. Overexpression of Cyp6g1 appears to have evolved in parallel in Drosophila simulans, a sibling species of D. melanogaster, where it is also associated with insecticide resistance. However, it is not known whether the ability of the CYP6G1 enzyme to provide resistance to multiple insecticides evolved recently in D. melanogaster or if this function is present in all Drosophila species. Here we show that duplication of the Cyp6g1 gene occurred at least four times during the evolution of different Drosophila species, and the ability of CYP6G1 to confer resistance to multiple insecticides exists in D. melanogaster and D. simulans but not in Drosophila willistoni or Drosophila virilis. In D. virilis, which has multiple copies of Cyp6g1, one copy confers resistance to DDT and another to nitenpyram, suggesting that the divergence of protein sequence between copies subsequent to the duplication affected the activity of the enzyme. All orthologs tested conferred resistance to one or more insecticides, suggesting that CYP6G1 had the capacity to provide resistance to anthropogenic chemicals before they existed. Finally, we show that expression of Cyp6g1 in the Malpighian tubules, which contributes to DDT resistance in D. melanogaster, is specific to the D. melanogaster–D. simulans lineage. Our results suggest that a combination of gene duplication, regulatory changes and protein coding changes has taken place at the Cyp6g1 locus during evolution and this locus may play a role in providing resistance to different environmental toxins in different Drosophila species.     

Triplex forming oligonucleotide targeted to 3'UTR downregulates the expression of the bcl-2 proto-oncogene in HeLa cells

The bcl-2 proto-oncogene is overexpressed in a variety of human cancers and plays an important role in programmed cell death. Recent reports implied that the 3′-untranslated region (3′UTR) functions effectively in the regulation of gene expression. Here, we attempt to assay the ability of triplex forming oligonucleotides (TFOs) to inhibit expression of a target gene in vivo and to examine the potential of the 3′UTR of the bcl-2 proto-oncogene in the regulation of bcl-2 gene expression. To do this, we have developed a novel cellular system that involves transfection of a Doxycyclin inducible expression plasmid containing the bcl-2 ORF and the 3′UTR together with a TFO targeted to the 3′UTR of the bcl-2 proto-oncogene. Phosphorothioate-modified TFO targeted to the 3′UTR of the bcl-2 gene significantly downregulated the expression of the bcl-2 gene in HeLa cells as demonstrated by western blotting. Our results indicate that blocking the functions of the 3′UTR using the TFO can downregulate the expression of the targeted gene, and suggest that triplex strategy is a promising approach for oligonucleotide-based gene therapy. In addition, triplex-based sequence targeting may provide a useful tool for studying the regulation of gene expression.     

Distribution and Expression of Elicitin Genes in the Interspecific Hybrid Oomycete Phytophthora alni

Phytophthora alni subsp. alni, P. alni subsp. multiformis, and P. alni subsp. uniformis are responsible for alder disease in Europe. Class I and II elicitin gene patterns of P. alni subsp. alni, P. alni subsp. multiformis, P. alni subsp. uniformis, and the phylogenetically close species P. cambivora and P. fragariae were studied through mRNA sequencing and 3′ untranslated region (3′UTR)-specific PCRs and sequencing. The occurrence of multiple 3′UTR sequences in association with identical elicitin-encoding sequences in P. alni subsp. alni indicated duplication/recombination events. The mRNA pattern displayed by P. alni subsp. alni demonstrated that elicitin genes from all the parental genomes are actually expressed in this allopolyploid taxon. The complementary elicitin patterns resolved confirmed the possible involvement of P. alni subsp. multiformis and P. alni subsp. uniformis in the genesis of the hybrid species P. alni subsp. alni. The occurrence of multiple and common elicitin gene sequences throughout P. cambivora, P. fragariae, and P. alni sensu lato, not observed in other Phytophthora species, suggests that duplication of these genes occurred before the radiation of these species.