A screen for fast evolving genes from Drosophila

In an attempt to quantify the rates of protein sequence divergence in Drosophila, we have devised a screen to differentiate between slow and fast evolving genes. We find that over one-third of randomly drawn cDNAs from a Drosophila melanogaster library do not cross-hybridize with Drosophila virilis DNA, indicating that they evolve with a very high rate. To determine the evolutionary characteristics of such protein sequences, we sequenced their homologs from a more closely related species (Drosophila yakuba). The amino acid substitution rates among these cDNAs are among the fastest known and several are only about 2-fold lower than the corresponding values for silent substitutions. An analysis of within-species polymorphisms for one of these sequences reveals an exceptionally high number of polymorphic amino acid positions, indicating that the protein is not under strong negative selection. We conclude that the Drosophila genome harbors a substantial proportion of genes with a very high divergence rate.

Establishment of a high-efficiency SNP-based framework marker set for Arabidopsis

The major goal of this project was the establishment of a tool for rapid mapping of new mutations and genotyping in Arabidopsis consisting of at least 100 evenly spaced framework markers. We assembled a single nucleotide polymorphism (SNP)-based marker set consisting of 112 polymorphic sites with average spacing of 1.15 Mbp derived from an SNP database that we recently developed. This information was used to set up efficient SNP detection reactions based on multiplexed primer extension assays. The 112 Columbia (Col-0)/C24 framework markers were used to assemble 18 multiplexed SNaPshot assays with which up to eight separate loci can be genotyped in a single-tube/single-capillary format. In addition, for 110 framework markers matrix-assisted laser desorption/ionization time of flight (MALDI-ToF) assays have been established for high throughput analyses. We demonstrated the usefulness and the robustness of both procedures of this tool by genotyping 48 BC3F1 individuals created between the accessions Col-0 and C24. Subsets of 10-62 of the established markers discriminate between various combinations of the accessions Col-0, C24, Landsberg erecta (Ler), Cape Verdi Islands (Cvi) and Niederzenz (Nd). Using a subset of 17 evenly distributed and established SNP markers that are also polymorphic between Ler and Col-0, we were able to rapidly map a mutant gene (tbr1) to an interval of 2.3 Mbp in an Ler (tbr1) x Col-0 cross.

The evolutionary analysis of "orphans" from the Drosophila genome identifies rapidly diverging and incorrectly annotated genes

In genome projects of eukaryotic model organisms, a large number of novel genes of unknown function and evolutionary history ("orphans") are being identified. Since many orphans have no known homologs in distant species, it is unclear whether they are restricted to certain taxa or evolve rapidly, either because of a lack of constraints or positive Darwinian selection. Here we use three criteria for the selection of putatively rapidly evolving genes from a single sequence of Drosophila melanogaster. Thirteen candidate genes were chosen from the Adh region on the second chromosome and 1 from the tip of the X chromosome. We succeeded in obtaining sequence from 6 of these in the closely related species D. simulans and D. yakuba. Only 1 of the 6 genes showed a large number of amino acid replacements and in-frame insertions/deletions. A population survey of this gene suggests that its rapid evolution is due to the fixation of many neutral or nearly neutral mutations. Two other genes showed "normal" levels of divergence between species. Four genes had insertions/deletions that destroy the putative reading frame within exons, suggesting that these exons have been incorrectly annotated. The evolutionary analysis of orphan genes in closely related species is useful for the identification of both rapidly evolving and incorrectly annotated genes.

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.

From genes to individuals: developmental genes and the generation of the phenotype

The success of the genetic approach to developmental biology has provided us with a suite of genes that are involved in the regulation of ontogenetic pathways. It is therefore time to ask whether and how such genes might be involved in the generation of adaptive phenotypes. Unfortunately, the current results do not provide a clear answer. Most of the genes that have been studied by developmental biologists affect early embryonic traits with significant effects on the whole organism. These genes are often highly conserved which allows us to do comparative studies even across phyla. However, whether the same genes are also involved in short-term ecological adaptations remains unclear. The suggestion that early acting ontogenetic genes may also affect late phenotypes comes from the genetic analysis of quantitative traits like bristle numbers in Drosophila. A rough mapping of the major loci affecting these traits shows that these loci might correspond to well known early acting genes. On the other hand, there are also many minor effect loci that are as yet uncharacterized. We suggest that these minor loci might correspond to a different class of genes. In comparative studies of randomly drawn cDNAs from Drosophila we find that there is a large group of genes that evolve fast and that are significantly under-represented in normal genetic screens. We speculate that these genes might provide a large, as yet poorly understood, reservoir of genes that might be involved in the evolution of quantitative traits and short-term adaptations.

Phenotypic landscapes: phenological patterns in wild and cultivated barley

Phenotypic variation in natural populations is the outcome of the joint effects of environmentally induced adaptations and neutral processes on the genetic architecture of quantitative traits. In this study, we examined the role of adaptation in shaping wild barley phenotypic variation along different environmental gradients. Detailed phenotyping of 164 wild barley (Hordeum spontaneum) accessions from Israel (of the Barley1K collection) and 18 cultivated barley (H. vulgare) varieties was conducted in common garden field trials. Cluster analysis based on phenotypic data indicated that wild barley in this region can be differentiated into three ecotypes in accordance with their ecogeographical distribution: north, coast and desert. Population differentiation (Q(ST) ) for each trait was estimated using a hierarchical Bayesian model and compared to neutral differentiation (F(ST) ) based on 42 microsatellite markers. This analysis indicated that the three clusters diverged in morphological but not in reproductive characteristics. To address the issue of phenotypic variation along environmental gradients, climatic and soil gradients were compared with each of the measured traits given the geographical distance between sampling sites using a partial Mantel test. Flowering time and plant growth were found to be differentially correlated with climatic and soil characteristic gradients, respectively. The H. vulgare varieties were superior to the H. spontaneum accessions in yield components, yet resembled the Mediterranean types in vegetative characteristics and flowering time, which may indicate the geographical origin of domesticated barley.

DNA sequence variation of wild barley Hordeum spontaneum (L.) across environmental gradients in Israel

Wild barley Hordeum spontaneum (L.) shows a wide geographic distribution and ecological diversity. A key question concerns the spatial scale at which genetic differentiation occurs and to what extent it is driven by natural selection. The Levant region exhibits a strong ecological gradient along the North-South axis, with numerous small canyons in an East-West direction and with small-scale environmental gradients on the opposing North- and South-facing slopes. We sequenced 34 short genomic regions in 54 accessions of wild barley collected throughout Israel and from the opposing slopes of two canyons. The nucleotide diversity of the total sample is 0.0042, which is about two-thirds of a sample from the whole species range (0.0060). Thirty accessions collected at 'Evolution Canyon' (EC) at Nahal Oren, close to Haifa, have a nucleotide diversity of 0.0036, and therefore harbor a large proportion of the genetic diversity. There is a high level of genetic clustering throughout Israel and within EC, which roughly differentiates the slopes. Accessions from the hot and dry South-facing slope have significantly reduced genetic diversity and are genetically more distinct from accessions from the North-facing slope, which are more similar to accessions from other regions in Northern Israel. Statistical population models indicate that wild barley within the EC consist of three separate genetic clusters with substantial gene flow. The data indicate a high level of population structure at large and small geographic scales that shows isolation-by-distance, and is also consistent with ongoing natural selection contributing to genetic differentiation at a small geographic scale.

A comparison of homologous developmental genes from Drosophila and Tribolium reveals major differences in length and trinucleotide repeat content

The flour beetle Tribolium castaneum has become an important model organism for comparative studies of insect development. Many developmentally important genes have now been cloned from both Tribolium and Drosophila and their expression characteristics were studied. We analyze here the complete coding sequences of 17 homologous gene pairs from D. melanogaster and T. castaneum, most of which encode transcription factors. We find that the Tribolium genes are on average 30% shorter than their Drosophila homologues. This appears to be due largely to the almost-complete absence of trinucleotide repeats in the coding sequences of Tribolium as well as the generally lower degree of internal repetitiveness. Clusters of polar and other amino acids such as glutamine, proline, and serine, which are often considered to be important for transcriptional activation domains in Drosophila, are almost completely absent in Tribolium. Codon usage is generally less biased in Tribolium, although we find a similar tendency for the preference of G- or C-ending codons and a higher bias in conserved subregions of the proteins as in Drosophila. Most of the aminoacid substitutions in the DNA-binding domains of the transcription factors occur at residues that do not make a specific contact to DNA, suggesting that the recognition sequences are likely to be conserved between the two species.

The Drosophila gene 2A5 complements the defect in mitochondrial F1-ATPase assembly in yeast lacking the molecular chaperone Atp11p

Assembly of mitochondrial F1-ATPase in Saccharomyces cerevisiae requires the molecular chaperone, Atp11p. Database searches have identified protein sequences from Schizosaccharomyces pombe and two species of Drosophila that are homologous to S. cerevisiae Atp11p. A cDNA encoding the putative Atp11p from Drosophila yakuba was shown to complement the respiratory deficient phenotype of yeast harboring an atp11::HIS3 disruption allele. Furthermore, the product of this Drosophila gene was shown to interact with the S. cerevisiae F1 beta subunit in the yeast two-hybrid assay. These results indicate that Atp11p function is conserved in higher eukaryotes.

Sequence and expression of DmMKLP1, a homolog of the human MKLP1 kinesin-like protein from Drosophila melanogaster

We have isolated the Drosophila gene DmMKLP1, which has a high similarity to members of the mitotic kinesin-like subfamily of kinesin proteins. DmMKLP1 has no known close relatives in the Drosophila genome and can therefore be assumed to be the ortholog of human MKLP1 and hamster CHOI kinesin-like proteins. In situ hybridization reveals a homogeneous maternal expression in the early embryo and a terminally restricted expression pattern at blastoderm stage. Later, the expression becomes increasingly restricted to the developing central nervous system, where it remains expressed at least until the end of embryogenesis.