A microsatellite-based multilocus phylogeny of the Drosophila melanogaster species complex

Genetic structure and diversity of coffee (Coffea) across Africa and the Indian Ocean islands revealed using microsatellites

BACKGROUND AND AIMS

The coffee genus (Coffea) comprises 124 species, and is indigenous to the Old World Tropics. Due to its immense economic importance, Coffea has been the focus of numerous genetic diversity studies, but despite this effort it remains insufficiently studied. In this study the genetic diversity and genetic structure of Coffea across Africa and the Indian Ocean islands is investigated.

METHODS

Genetic data were produced using 13 polymorphic nuclear microsatellite markers (simple sequence repeats, SSRs), including seven expressed sequence tag-SSRs, and the data were analysed using model- and non-model-based methods. The study includes a total of 728 individuals from 60 species.

KEY RESULTS

Across Africa and the Indian Ocean islands Coffea comprises a closely related group of species with an overall pattern of genotypes running from west to east. Genetic structure was identified in accordance with pre-determined geographical regions and phylogenetic groups. There is a good relationship between morpho-taxonomic species delimitations and genetic units. Genetic diversity in African and Indian Ocean Coffea is high in terms of number of alleles detected, and Madagascar appears to represent a place of significant diversification in terms of allelic richness and species diversity.

CONCLUSIONS

Cross-species SSR transferability in African and Indian Ocean islands Coffea was very efficient. On the basis of the number of private alleles, diversification in East Africa and the Indian Ocean islands appears to be more recent than in West and West-Central Africa, although this general trend is complicated in Africa by the position of species belonging to lineages connecting the main geographical regions. The general pattern of phylogeography is not in agreement with an overall east to west (Mascarene, Madagascar, East Africa, West Africa) increase in genome size, the high proportion of shared alleles between the four regions or the high numbers of exclusive shared alleles between pairs or triplets of regions.

Allele frequencies of 6 autosomal STR loci in the Xibo nationality with phylogenetic structure among Chinese populations

In the present study, we investigated the genetic polymorphisms of 6 autosomal STR loci Hum-CSF1PO, D13S317, D5S818, D16S539, TH01, and TPOX in the Xibo population of Liaoning, northeastern China as well as its genetic relationships with other populations in China. No significant deviations from Hardy-Weinberg equilibrium could be found for all loci. Allele frequencies in the Xibo population ranged from 0.001 to 0.507. Among all the 6 loci, D16S539 had the highest polymorphism (PIC=0.8632), whereas TPOX had the lowest (PIC=0.5179). A phylogenic tree was constructed using Poptree 2 software. In the phylogenic tree, Xibo population has a distant relationship with the other populations.

Inter-island divergence within Drosophila mauritiana, a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.

Multiple hybridization events between Drosophila simulans and Drosophila mauritiana are supported by mtDNA introgression

The study of speciation has advanced considerably in the last decades because of the increased application of molecular tools. In particular, the quantification of gene flow between recently diverged species could be addressed. Drosophila simulans and Drosophila mauritiana diverged, probably allopatrically, from a common ancestor approximately 250,000 years ago. However, these species share one mitochondrial DNA (mtDNA) haplotype indicative of a recent episode of introgression. To study the extent of gene flow between these species, we took advantage of a large sample of D. mauritiana and employed a range of different markers, i.e. nuclear and mitochondrial sequences, and microsatellites. This allowed us to detect two new mtDNA haplotypes (MAU3 and MAU4). These haplotypes diverged quite recently from haplotypes of the siII group present in cosmopolitan populations of D. simulans. The mean divergence time of the most diverged haplotype (MAU4) is approximately 127,000 years, which is more than 100,000 years before the assumed speciation time. Interestingly, we also found some evidence for gene flow at the nuclear level because an excess of putatively neutral loci shows significantly reduced differentiation between D. simulans and D. mauritiana. Our results suggest that these species are exchanging genes more frequently than previously thought.

Evolution of endogenous sequences of banana streak virus: what can we learn from banana (Musa sp.) evolution?

Endogenous plant pararetroviruses (EPRVs) are viral sequences of the family Caulimoviridae integrated into the nuclear genome of numerous plant species. The ability of some endogenous sequences of Banana streak viruses (eBSVs) in the genome of banana (Musa sp.) to induce infections just like the virus itself was recently demonstrated (P. Gayral et al., J. Virol. 83:6697-6710, 2008). Although eBSVs probably arose from accidental events, infectious eBSVs constitute an extreme case of parasitism, as well as a newly described strategy for vertical virus transmission in plants. We investigated the early evolutionary stages of infectious eBSV for two distinct BSV species-GF (BSGFV) and Imové (BSImV)-through the study of their distribution, insertion polymorphism, and structure evolution among selected banana genotypes representative of the diversity of 60 wild Musa species and genotypes. To do so, the historical frame of host evolution was analyzed by inferring banana phylogeny from two chloroplast regions-matK and trnL-trnF-as well as from the nuclear genome, using 19 microsatellite loci. We demonstrated that both BSV species integrated recently in banana evolution, circa 640,000 years ago. The two infectious eBSVs were subjected to different selective pressures and showed distinct levels of rearrangement within their final structure. In addition, the molecular phylogenies of integrated and nonintegrated BSVs enabled us to establish the phylogenetic origins of eBSGFV and eBSImV.

Sequence determinants of human microsatellite variability

BACKGROUND

Microsatellite loci are frequently used in genomic studies of DNA sequence repeats and in population studies of genetic variability. To investigate the effect of sequence properties of microsatellites on their level of variability we have analyzed genotypes at 627 microsatellite loci in 1,048 worldwide individuals from the HGDP-CEPH cell line panel together with the DNA sequences of these microsatellites in the human RefSeq database.

RESULTS

Calibrating PCR fragment lengths in individual genotypes by using the RefSeq sequence enabled us to infer repeat number in the HGDP-CEPH dataset and to calculate the mean number of repeats (as opposed to the mean PCR fragment length), under the assumption that differences in PCR fragment length reflect differences in the numbers of repeats in the embedded repeat sequences. We find the mean and maximum numbers of repeats across individuals to be positively correlated with heterozygosity. The size and composition of the repeat unit of a microsatellite are also important factors in predicting heterozygosity, with tetra-nucleotide repeat units high in G/C content leading to higher heterozygosity. Finally, we find that microsatellites containing more separate sets of repeated motifs generally have higher heterozygosity.

CONCLUSIONS

These results suggest that sequence properties of microsatellites have a significant impact in determining the features of human microsatellite variability.

Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup

Background

Sexual dimorphism of body size has been the subject of numerous studies, but few have examined sexual shape dimorphism (SShD) and its evolution. Allometry, the shape change associated with size variation, has been suggested to be a main component of SShD. Yet little is known about the relative importance of the allometric and non-allometric components for the evolution of SShD.

Results

We investigated sexual dimorphism in wing shape in the nine species of the Drosophila melanogaster subgroup. We used geometric morphometrics to characterise wing shape and found significant SShD in all nine species. The amount of shape difference and the diversity of the shape changes evolved across the group. However, mapping the divergence of SShD onto the phylogeny of the Drosophila melanogaster subgroup indicated that there is little phylogenetic signal. Finally, allometry accounted for a substantial part of SShD, but did not explain the bulk of evolutionary divergence in SShD because allometry itself was found to be evolutionarily plastic.

Conclusion

SShD in the Drosophila wing can evolve rapidly and therefore shows only weak phylogenetic structure. The variable contribution of allometric and non-allometric components to the evolutionary divergence of SShD and the evolutionary plasticity of allometry suggest that SShD and allometry are influenced by a complex interaction of processes.

Genetic structure of the polymorphic metrosideros (Myrtaceae) complex in the Hwaiian islands using nuclear microsatellite data

Background

Five species of Metrosideros (Myrtaceae) are recognized in the Hawaiian Islands, including the widespread M. polymorpha, and are characterized by a multitude of distinctive, yet overlapping, habit, ecological, and morphological forms. It remains unclear, despite several previous studies, whether the morphological variation within Hawaiian Metrosideros is due to hybridization, genetic polymorphism, phenotypic plasticity, or some combination of these processes. The Hawaiian Metrosideros complex has become a model system to study ecology and evolution; however this is the first study to use microsatellite data for addressing inter-island patterns of variation from across the Hawaiian Islands.

Methodology/Principal Findings

Ten nuclear microsatellite loci were genotyped from 143 individuals of Metrosideros. We took advantage of the bi-parental inheritance and rapid mutation rate of these data to examine the validity of the current taxonomy and to investigate whether Metrosideros plants from the same island are more genetically similar than plants that are morphologically similar. The Bayesian algorithm of the program structure was used to define genetic groups within Hawaiian Metrosideros and the closely related taxon M. collina from the Marquesas and Austral Islands. Several standard and nested AMOVAs were conducted to test whether the genetic diversity is structured geographically or taxonomically.

Conclusions/Significance

The results suggest that Hawaiian Metrosideros have dynamic gene flow, with genetic and morphological diversity structured not simply by geography or taxonomy, but as a result of parallel evolution on islands following rampant island-island dispersal, in addition to ancient chloroplast capture. Results also suggest that the current taxonomy requires major revisions in order to reflect the genetic structure revealed in the microsatellite data.

Evolution of microsatellite loci in the adaptive radiation of Hawaiian honeycreepers

Previous studies have examined germ-line mutations to infer the processes that generate and maintain variability in microsatellite loci. Few studies, however, have examined patterns to infer processes that act on microsatellite loci over evolutionary time. Here, we examine changes in 8 dinucleotide loci across the adaptive radiation of Hawaiian honeycreepers. The loci were found to be highly variable across the radiation, and we did not detect ascertainment bias with respect to allelic diversity or allele size ranges. In examining patterns at the sequence level, we found that changes in flanking regions, repeat motifs, or repeat interruptions were often shared between closely related species and may be phylogenetically informative. Genetic distance measures based on microsatellites were strongly correlated with those based on mitochondrial DNA (mtDNA) sequences as well as with divergence time up to 3 My. Phylogenetic inferences based on microsatellite genetic distances consistently recovered 2 of the 4 honeycreeper clades observed in a tree based on mtDNA sequences but differed from the mtDNA tree in the relationships among clades. Our results confirm that microsatellite loci may be conserved over evolutionary time, making them useful in population-level studies of species that diverged from the species in which they were characterized as long as 5 Ma. Despite this, we found that their use in phylogenetic inference was limited to closely related honeycreeper species.

The population structure of the Cryptosporidium parvum population in Scotland: a complex picture

We genotyped 297 Scottish C. parvum samples using micro- and minisatellites. Treated as a single population, the population structure was epidemic. When regional populations were analysed, there was evidence of sub-population structure variations. This was dependent upon excluding sub-groups exhibiting significant genetic distance from the main population, implying genetic sub-structuring. We tested the hypothesis that these sub-groups originated outside the UK and demonstrated that one sub-group clustered with Peruvian samples. A geographically comprehensive panel of isolates would fully confirm this result. These data indicate limited sub-structuring within a small geographical area, but substantial sub-structuring over larger geographical distances. Host movement influences parasite diversity and population structure, evidenced by strong correlation (r² = 0.9686) between cattle movements and parasite diversity. Thus, the population structure of C. parvum is complex, with sub-populations differing in structure and being influenced by host movements, including the introduction of novel multilocus genotypes from geographically distinct regions.

Variation in fiber number of a male-specific muscle between Drosophila species: a genetic and developmental analysis

We characterize a newly discovered morphological difference between species of the Drosophila melanogaster subgroup. The muscle of Lawrence (MOL) contains about four to five fibers in D. melanogaster and Drosophila simulans and six to seven fibers in Drosophila mauritiana and Drosophila sechellia. The same number of nuclei per fiber is present in these species but their total number of MOL nuclei differs. This suggests that the number of muscle precursor cells has changed during evolution. Our comparison of MOL development indicates that the species difference appears during metamorphosis. We mapped the quantitative trait loci responsible for the change in muscle fiber number between D. sechellia and D. simulans to two genomic regions on chromosome 2. Our data eliminate the possibility of evolving mutations in the fruitless gene and suggest that a change in the twist might be partly responsible for this evolutionary change.

Morphological evolution through multiple cis-regulatory mutations at a single gene

One central, and yet unsolved, question in evolutionary biology is the relationship between the genetic variants segregating within species and the causes of morphological differences between species. The classic neo-darwinian view postulates that species differences result from the accumulation of small-effect changes at multiple loci. However, many examples support the possible role of larger abrupt changes in the expression of developmental genes in morphological evolution. Although this evidence might be considered a challenge to a neo-darwinian micromutationist view of evolution, there are currently few examples of the actual genes causing morphological differences between species. Here we examine the genetic basis of a trichome pattern difference between Drosophila species, previously shown to result from the evolution of a single gene, shavenbaby (svb), probably through cis-regulatory changes. We first identified three distinct svb enhancers from D. melanogaster driving reporter gene expression in partly overlapping patterns that together recapitulate endogenous svb expression. All three homologous enhancers from D. sechellia drive expression in modified patterns, in a direction consistent with the evolved svb expression pattern. To test the influence of these enhancers on the actual phenotypic difference, we conducted interspecific genetic mapping at a resolution sufficient to recover multiple intragenic recombinants. This functional analysis revealed that independent genetic regions upstream of svb that overlap the three identified enhancers are collectively required to generate the D. sechellia trichome pattern. Our results demonstrate that the accumulation of multiple small-effect changes at a single locus underlies the evolution of a morphological difference between species. These data support the view that alleles of large effect that distinguish species may sometimes reflect the accumulation of multiple mutations of small effect at select genes.

High-resolution quantitative trait locus mapping reveals sign epistasis controlling ovariole number between two Drosophila species

Identifying the genes underlying genetically complex traits is of fundamental importance for medicine, agriculture, and evolutionary biology. However, the level of resolution offered by traditional quantitative trait locus (QTL) mapping is usually coarse. We analyze here a trait closely related to fitness, ovariole number. Our initial interspecific mapping between Drosophila sechellia (8 ovarioles/ovary) and D. simulans (15 ovarioles/ovary) identified a major QTL on chromosome 3 and a minor QTL on chromosome 2. To refine the position of the major QTL, we selected 1038 additional recombinants in the region of interest using flanking morphological markers (selective phenotyping). This effort generated approximately one recombination event per gene and increased the mapping resolution by approximately seven times. Our study thus shows that using visible markers to select for recombinants can efficiently increase the resolution of QTL mapping. We resolved the major QTL into two epistatic QTL, QTL3a and QTL3b. QTL3a shows sign epistasis: it has opposite effects in two different genetic backgrounds, the presence vs. the absence of the QTL3b D. sechellia allele. This property of QTL3a allows us to reconstruct the probable order of fixation of the QTL alleles during evolution.

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.

Refugia, dispersal and divergence in a forest archipelago: a study of Streptocarpus in eastern South Africa

We describe a scenario of plant speciation across a relict forest archipelago in South Africa involving Pleistocene habitat expansion-contraction cycles, dispersal and adaptation to lower temperatures. This is the first population level study using molecular data in South African forests and has significant implications for conservation efforts in this area. Populations of the mesophytic forest floor herbs Streptocarpus primulifolius sensu lato and Streptocarpus rexii were sampled throughout their range in the naturally fragmented forests of eastern South Africa in order to investigate population genetic and phylogenetic patterns within the species complex, using nuclear microsatellites, nuclear ribosomal ITS (internal transcribed spacer) sequences and chloroplast genome sequences. S. primulifolius harbours high levels of genetic diversity at both the nuclear (mean HE = 0.50) and the chloroplast level (each population fixed for a unique haplotype). This is consistent with populations of these coastal species being Pleistocene relicts. In contrast, populations of S. rexii in cooler habitats at higher altitudes and lower latitudes harbour little or no nuclear genetic diversity (mean HE = 0.09) and most share a common chloroplast haplotype. The split of S. rexii from populations intermediate between the two species (S. cf. primulifolius) occurred between 0 and 0.44 million years ago according to the calibrated ITS phylogeny of the taxa. The low genetic diversity and homogeneity of S. rexii is congruent with this species having reached its current range during the Holocene. We found no evidence of monophyly of any of the taxa in this study, which we consider a consequence of recent evolution in a fragmented habitat.

Evolution of the chromosomal location of rDNA genes in two Drosophila species subgroups: ananassae and melanogaster

The evolution of the chromosomal location of ribosomal RNA gene clusters and the organization of heterochromatin in the Drosophila melanogaster group were investigated using fluorescence in situ hybridization and DAPI staining to mitotic chromosomes. The investigation of 18 species (11 of which were being examined for the first time) belonging to the melanogaster and ananassae subgroups suggests that the ancestral configuration consists of one nucleolus organizer (NOR) on each sex chromosome. This pattern, which is conserved throughout the melanogaster subgroup, except in D. simulans and D. sechellia, was observed only in the ercepeae complex within the ananassae subgroup. Both sex-linked NORs must have been lost in the lineage leading to D. varians and in the ananassae and bipectinata complexes, whereas new sites, characterized by intra-species variation in hybridization signal size, appeared on the fourth chromosome related to heterochromatic rearrangements. Nucleolar material is thought to be required for sex chromosome pairing and disjunction in a variety of organisms including Drosophila. Thus, either remnant sequences, possibly intergenic spacer repeats, are still present in the sex chromosomes which have lost their NORs (as observed in D. simulans and D. sechellia), or an alternative mechanism has evolved.

A multilocus microsatellite phylogeny of the Drosophila virilis group

We used a set of 48 polymorphic microsatellites derived from Drosophila virilis to infer phylogenetic relationships in the D. virilis clade. Consistent with previous studies, D. virilis and D. lummei were the most basal species of the group. Within the D. montana phylad, the phylogenetic relationship could not be resolved. Special attention was given to the differentiation between D. americana texana, D. americana americana and D. novamexicana. Significant differences between these three groups were detected by F(ST) analyses. Similarly, a model-based clustering method for multilocus genotype data also provided strong support for the presence of three differentiated groups. This genome-wide differentiation between D. americana texana and D. americana americana contrasts with previous analyses based on DNA sequence data.

Microsatellite variation and evolution in the Peromyscus maniculatus species group

Variation at 12 pure-repeat dinucleotide microsatellites from Peromyscus maniculatus was analyzed for samples of all species in the P. maniculatus species group and P. leucopus. Except for one locus (Pml08) that amplified a product only for P. maniculatus, these microsatellites yielded reliable estimates of variation across these species; per-locus polymorphism and allele-size distribution were not significantly different among or between any of the species sampled from mainland populations. Significantly lower levels of variation and the distribution of alleles in the two populations of the insular-endemic P. sejugis were consistent with the expectation of substantial founder effect and suggest a lack of recent gene flow between these populations. Phenetic analyses of genetic distances based on shared allele frequencies uniformly produced well-supported trees that were entirely concordant with the a priori corroborated relationships within the P. maniculatus species group; this result was not obtained with analyses of a genetic distance computed from differences in allele sizes. The microsatellite data do not support the hypothesis that P. sejugis should be considered conspecific with P. maniculatus but yield a strongly supported sister-group association between P. sejugis and P. keeni.

Genetic diversity, clonality and sexuality in Toxoplasma gondii

The majority of Toxoplasma gondii strains from a variety of human and animal sources have been grouped into three highly clonal but closely related lineages. The low occurrence of nucleotide differences among the three predominant lineages and their unusual dimorphic allelic composition suggest that they have arisen from a recent common ancestry. Less than 1% of the previously studied strains contain unique genotypes and high divergence of DNA sequence, and therefore are considered 'exotic' or 'atypical' strains. The seemingly low genetic diversity in T. gondii may have been underestimated because most parasite strains in previous studies were collected from human patients and domestic animals in North America and Europe. To investigate the genetic diversity of T. gondii, we analysed parasite strains isolated from remote geographical regions by multilocus microsatellite sequencing and phylogenetic analysis. The genetic diversity indices, the molecular analysis of microsatellite genotypes and the constructed phylogram considered together suggest that the global T. gondii population is highly diversified and not characteristic of a clonal organism. The most parsimonious hypothesis is that T. gondii presents a complex population structure with a mix of clonal and sexual propagation as a function of the environmental conditions. The comparison between domestic strains data on one hand and wild strains data on the other hand is in favour of more frequent sexual recombinations in wild environment even though Toxoplasma subpopulation in human and domestic animals is largely clonal.

Mutation and evolution of microsatellite loci in Neurospora

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.

Cytonuclear coadaptation in Drosophila: disruption of cytochrome c oxidase activity in backcross genotypes

The cytochrome c oxidase enzyme (COX) is comprised of 10 nuclear-encoded subunits and three mitochondrial-encoded subunits in close physical association in the inner mitochondrial membrane. COX passes electrons from cytochrome c to molecular oxygen and pumps protons into the inner mitochondrial space for ATP production. Selection on nuclear-mitochondrial interactions within species should lead to coadaptation of the proteins comprising this important enzyme. Under this model, there should be relatively little disruption of COX activity when mitochondrial genomes are crossed among strains within species. A more pronounced disruption of activity is expected when the mitochondrial genome is expressed in the nuclear background of a different species. We test these hypotheses in Drosophila using hybridization and backcrossing among lines of D. simulans and D. mauritiana. Disrupted cytonuclear genotypes were constructed using backcrosses between two lines of D. simulans (siI and siII) that introduced each divergent mitochondrial DNA (mtDNA) into each nuclear background due to maternal inheritance of mtDNA. Similar crosses were used to introduce each D. simulans mtDNA into the D. mauritiana maI nuclear background. Reconstituted cytonuclear control genotypes were constructed by backcrossing the initial F1 females to males of the maternal genotype. COX enzyme activities were compared among these disrupted and reconstituted backcross genotypes within and between species. The disruption effect on COX activity was restricted to males of interspecific genotypes. These data support the coadaptation hypothesis and are consistent with predictions that the evolution of modifiers of male mitochondrial dysfunction is hindered by the maternal inheritance of mtDNA. New sequence data for nuclear encoded subunits of COX identified amino acids that may play a role in the disruption effect.

Identification of plu genes and cis-acting elements of PCNA in the Drosophila genus using conservation of gene order

In Drosophila melanogaster, the cell cycle control gene, plutonium (plu), is located between the PCNA and RpS18 genes at position 56F on chromosome arm 2R. We have used a comparative genomic approach to investigate the evolution of the plu gene and to locate conserved cis-acting elements for plu, RpS18 and PCNA. Using primers within coding regions of PCNA and RpS18, we amplified and sequenced the intervening region from twelve Drosophila species. In each species, this region contains a plu gene resembling the D. melanogaster gene in size and in the number and position of introns. The predicted Plu sequence from the different species demonstrates that the first two ankyrin repeats are conserved. Of the transcriptional control elements of D. melanogaster PCNA, we found that three motifs 5' to the PCNA transcription unit are conserved in Drosophila species.

Molecular approaches to the study of myiasis-causing larvae

Among arthropod diseases affecting animals, larval infections - myiases - of domestic and wild animals have been considered important since ancient times. Besides the significant economic losses to livestock worldwide, myiasis-causing larvae have attracted the attention of scientists because some parasitise humans and are of interest in forensic entomology. In the past two decades, the biology, epidemiology, immunology, immunodiagnosis and control methods of myiasis-causing larvae have been focused on and more recently the number of molecular studies have also begun to increase. The 'new technologies' (i.e. molecular biology) are being used to study taxonomy, phylogenesis, molecular identification, diagnosis (recombinant antigens) and vaccination strategies. In particular, more in depth molecular studies have now been performed on Sarcophagidae, Calliphoridae and flies of the Oestridae sister group. This review discusses the most topical issues and recent studies on myiasis-causing larvae using molecular approaches. In the first part, PCR-based techniques and the genes that have already been analysed, or are potentially useful for the molecular phylogenesis and identification of myiasis-causing larvae, are described. The second section deals with the more recent advances concerning taxonomy, phylogenetics, population studies, molecular identification, diagnosis and vaccination.

Conservation of a dinucleotide simple sequence repeat locus in sharks

Recent studies indicate that the flanking region and repeat motif structure of conserved microsatellite loci are useful for phylogenetic inference. Most comparative studies of microsatellite loci involve relatively closely related species, however, primarily because primers developed for one species often amplify only related species. We describe an analysis of a microsatellite locus in lamniform sharks that we estimate has been conserved for a billion years. Combined analysis of the flanking sequence and repeat motif structure resulted in a gene tree comparable to those reported from similar analyses of other genes. The conservation of the simple sequence repeat (SSR), and of the sequence flanking the SSR, is explained by a low substitution rate in sharks coupled with the possibility that mutations which interrupt perfect repeats are lost by replication slippage.

Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae

The phylogeography of an endangered salmonid, European grayling (Thymallus thymallus), was studied based on analysis of 17 nuclear microsatellite DNA loci. In agreement with earlier mitochondrial DNA (mtDNA) studies, phylogenetic relationships of the populations suggested that northern Europe was colonized from two distinct Pleistocene refugia. Furthermore, microsatellites revealed highly supported grouping of mainland Swedish, Norwegian, Danish, German and Slovenian populations, suggesting that grayling from the northwestern and central Europe have descended from their southern conspecifics. The level of divergence between populations was substantial, even across short geographical distances. Although this was in part due to postglacial colonization patterns and contemporary barriers for gene flow, the high divergence estimates between hydrologically connected sampling locations implied efficient interpopulation reproductive isolation. Microsatellites revealed that the populations exhibited, on average, only 3.5 (+/-2.2) alleles per locus, indicating that T. thymallus has strikingly low levels of intrapopulation genetic diversity as compared with other freshwater fish species. Accordingly, as indicated by analysis of molecular variance (AMOVA), only 49.1-58.0% of the total grayling microsatellite diversity resided within populations. A latitudinal genetic diversity gradient, potentially resulting from glaciation-mediated founder events, was not evident. Alternatively, it is possible that grayling display limited dispersal behaviour/capability, leading to low long-term effective population sizes and, consequently, depauperate intrapopulation polymorphism. These findings have implications for conservation of T. thymallus. Importantly, they exemplify that microsatellites can be highly informative for intraspecific phylogeography studies dealing with substantial divergence scales.

Microsatellite analysis of Drosophila melanogaster populations along a microclimatic contrast at lower Nahel Oren canyon, Mount Carmel, Israel

Drosophila melanogaster populations collected at the south-facing slope (SFS) and north-facing slope (NFS) of lower Nahel Oren canyon, Mount Carmel, Israel display significant differences in survival and longevity at temperature, drought, and starvation stresses. Furthermore, significant assortative mating was previously observed between populations of the two slopes. We used a set of 48 microsatellite markers to analyze patterns of genetic differentiation between D. melanogaster populations from both slopes and D. simulans. Consistent with previous reports, we found D. simulans to be well differentiated from D. melanogaster. Genetic differentiation between SFS and NFS D. melanogaster populations was low (F(ST) = 0.0012). Also a tree of individuals based on the proportion of shared alleles and a model-based clustering method provided no evidence for population substructuring.

Stationary distributions of microsatellite loci between divergent population groups of the European rabbit (Oryctolagus cuniculus)

Previous analysis of mitochondrial DNA polymorphism in the native range of the European rabbit (Oryctolagus cuniculus) demonstrated the occurrence of two highly divergent (2 Myr) maternal lineages with a well-defined geographical distribution. Analysis of both protein and immunoglobulin polymorphisms are highly concordant with this pattern of differentiation. However, the present analysis of nine polymorphic microsatellite loci (with a total of 169 alleles) in 24 wild populations reveals severe allele-size homoplasy which vastly underestimates divergence between the main groups of populations in Iberia. Nonetheless, when applied to more recent historical phenomena, this same data set not only confirms the occurrence of a strong bottleneck associated with the colonization of Mediterranean France but also suggests a two-step dispersal scenario that began with gene flow from northern Spain through the Pyrenean barrier and subsequent range expansion into northern France. The strength and appropriateness of applying microsatellites to more recent evolutionary questions is highlighted by the fact that both mtDNA and protein markers lacked the allelic diversity necessary to properly evaluate the colonization of France. The well-documented natural history of European rabbit populations provides an unusually comprehensive framework within which one can appraise the advantages and limitations of microsatellite markers in revealing patterns of genetic differentiation that have occurred across varying degrees of evolutionary time. The degree of size homoplasy presented in our data should serve as a warning to those drawing conclusions from microsatellite data sets which lack a set of complementary comparative markers, or involve long periods of evolutionary history, even within a single species.

Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.

When genetic distance matters: measuring genetic differentiation at microsatellite loci in whole-genome scans of recent and incipient mosquito species

Genetic distance measurements are an important tool to differentiate field populations of disease vectors such as the mosquito vectors of malaria. Here, we have measured the genetic differentiation between Anopheles arabiensis and Anopheles gambiae, as well as between proposed emerging species of the latter taxon, in whole genome scans by using 23-25 microsatellite loci. In doing so, we have reviewed and evaluated the advantages and disadvantages of standard parameters of genetic distance, F(ST), R(ST), (delta mu)(2), and D. Further, we have introduced new parameters, D' and D(K), which have well defined statistical significance tests and complement the standard parameters to advantage. D' is a modification of D, whereas D(K) is a measure of covariance based on Pearson's correlation coefficient. We find that A. gambiae and A. arabiensis are closely related at most autosomal loci but appear to be distantly related on the basis of X-linked chromosomal loci within the chromosomal Xag inversion. The M and S molecular forms of A. gambiae are practically indistinguishable but differ significantly at two microsatellite loci from the proximal region of the X, outside the Xag inversion. At one of these loci, both M and S molecular forms differ significantly from A. arabiensis, but remarkably, at the other locus, A. arabiensis is indistinguishable from the M molecular form of A. gambiae. These data support the recent proposal of genetically differentiated M and S molecular forms of A. gambiae.

Microsatellite variation in natural Drosophila melanogaster populations from New South Wales (Australia) and Tasmania

Microsatellite variation was studied at 48 microsatellite loci in 10 Drosophila melanogaster populations to investigate the population structure on the Australian east coast. Low, but statistically significant population differentiation was observed among most populations. The populations on the Australian mainland did not show evidence for isolation by distance. We conclude that the population structure of D. melanogaster on the Australian mainland is probably the result of a shared history (recent colonization). The observed differences between local D. melanogaster populations probably reflect variation in effective population sizes rather than patterns of gene flow. Two populations from Tasmania were more differentiated from the Australian mainland than a population from Israel, raising the question whether they are derived from the Australian mainland or colonized from a different source population.

Sex in Drosophila mauritiana: a very high level of amino acid polymorphism in a male reproductive protein gene, Acp26Aa

Many genes pertaining to male reproductive functions have been shown to evolve rapidly between species, and evidence increasingly suggest the influence of positive Darwinian selection. The accessory gland protein gene (Acp26Aa) of Drosophila is one such example. In order to understand the mechanism of selection, it is often helpful to examine the pattern of polymorphism. We report here that the level of amino acid polymorphism in the N-terminal quarter of Acp26Aa is high in Drosophila melanogaster and is unprecedented in its sibling species Drosophila mauritiana. We postulate that (1) this N-terminal segment may play a role in sperm competition, and (2) D. mauritiana may have been under much more intense sexual selection than other species. Both postulates have important ramifications and deserve to be tested rigorously.

Mitochondrial genomics flies high

A new pair of papers reports the complete coding sequences for 28 mitochondrial DNAs in Drosophila. By examining the patterns of polymorphism and divergence among functionally distinct classes of synonymous and nonsynonymous nucleotide sites, Bill Ballard provides a comprehensive whole-genome picture of how mtDNA sequence evolution can depart from the strictly neutral and nearly neutral models of molecular evolution.

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Microsatellite analysis indicates genetic differentiation of the neo-sex chromosomes in Drosophila americana americana

The neo-sex chromosomes of Drosophila americana americana were formed by a centric fusion between the autosomal element B and the X chromosome. Previous work has shown that the neo-Y chromosome is not degenerated genetically and that there is no evidence for genetic differentiation between neo-X and neo-Y chromosomes at the sequence level. To further address the genetic differentiation between the neo-sex chromosomes, microsatellites mapping to the neo-sex chromosome of D. a. americana were isolated. Microsatellite analysis indicated a highly significant differentiation between D. a. americana and D. a. texana as well as between the neo-X and neo-Y chromosomes in D. a. americana. Nevertheless, levels of variability were similar among the neo-sex chromosomes, indicating a very recent origin.

Microsatellite variability differs between dinucleotide repeat motifs-evidence from Drosophila melanogaster

Recently, the use of microsatellites as genetic markers has become very popular. While their evolutionary dynamics are not yet fully understood, the emerging picture is that several factors are influencing microsatellite mutation rates. Recent experiments demonstrated a significant effect of repeat motif length on microsatellite mutation rates. Here, we studied the influence of the base composition of the microsatellite. Forty-two microsatellite loci on the second chromosome with the three most abundant dinucleotide repeat motifs (TC/AG, AT/TA, GT/CA) were characterized for six different Drosophila melanogaster populations. Applying ANOVA to the variance in repeat number, we found a significant influence of repeat motif on microsatellite variability. Calculating relative mutation rates, GT/CA appears to have the highest mutation rate, and AT/TA appears to have the lowest. Similar differences in mutation rates were obtained by an alternative method which estimates microsatellite mutation rates from their genomic length distribution.

A test for concordance between the multilocus genealogies of genes and microsatellites in the pathogenic fungus Coccidioides immitis

Uncovering the correct phylogeny of closely related species requires analysis of multiple gene genealogies or, alternatively, genealogies inferred from the multiple alleles found at highly polymorphic loci, such as microsatellites. However, a concern in using microsatellites is that constraints on allele sizes may occur, resulting in homoplasious distributions of alleles, leading to incorrect phylogenies. Seven microsatellites from the pathogenic fungus Coccidioides immitis were sequenced for 20 clinical isolates chosen to represent the known genetic diversity of the pathogen. An organismal phylogeny for C. immitis was inferred from microsatellite-flanking sequence polymorphisms and other restriction fragment length polymorphism-containing loci. Two microsatellite genetic distances were then used to determine phylogenies for C. immitis, and the trees found by these three methods were compared. Congruence between the organismal and microsatellite phylogenies occurred when microsatellite distances were based on simple allele frequency data. However, complex mutation events at some loci made distances based on stepwise mutation models unreliable. Estimates of times of divergence for the two species of C. immitis based on microsatellites were significantly lower than those calculated from flanking sequence, most likely due to constraints on microsatellite allele sizes. Flanking-sequence insertions/deletions significantly decreased the accuracy of genealogical information inferred from microsatellite loci and caused interspecific length homoplasies at one of the seven loci. Our analysis shows that microsatellites are useful phylogenetic markers, although care should be taken to choose loci with appropriate flanking sequences when they are intended for use in evolutionary studies.

Removal of microsatellite interruptions by DNA replication slippage: phylogenetic evidence from Drosophila

Microsatellites are tandem repetitions of short (1-6 bp) motifs. It is widely assumed that microsatellites degenerate through the accumulation of base substitutions in the repeat array. Using a phylogenetic framework, we studied the evolutionary dynamics of interruptions in three Drosophila microsatellite loci. For all three loci, we show that the interruptions in a microsatellite can be lost, resulting in a longer uninterrupted microsatellite stretch. These results indicate that mutations in the microsatellite array do not necessarily lead to decay but may represent only a transition state during the evolution of a microsatellite. Most likely, this purification of interrupted microsatellites is caused by DNA replication slippage.

The phylogeny of closely related species as revealed by the genealogy of a speciation gene, Odysseus

Molecular differentiation between races or closely related species is often incongruent with the reproductive divergence of the taxa of interest. Shared ancient polymorphism and/or introgression during secondary contact may be responsible for the incongruence. At loci contributing to speciation, these two complications should be minimized (1, 2); hence, their variation may more faithfully reflect the history of the species' reproductive differentiation. In this study, we analyzed DNA polymorphism at the Odysseus (OdsH) locus of hybrid sterility between Drosophila mauritiana and Drosophila simulans and were able to verify such a prediction. Interestingly, DNA variation only a short distance away (1.8 kb) appears not to be influenced by the forces that shape the recent evolution of the OdsH coding region. This locus thus may represent a test case of inferring phylogeny of very closely related species.

Divergence of larval morphology between Drosophila sechellia and its sibling species caused by cis-regulatory evolution of ovo/shaven-baby

We report an extreme morphological difference between Drosophila sechellia and related species of the pattern of hairs on first-instar larvae. On the dorsum of most species, the posterior region of the anterior compartment of most segments is covered by a carpet of fine hairs. In D. sechellia, these hairs have been lost and replaced with naked cuticle. Genetic mapping experiments and interspecific complementation tests indicate that this difference is caused, in its entirety, by evolution at the ovo/shaven-baby locus. The pattern of expression of the ovo/shaven-baby transcript is correlated with this morphological change. The altered dorsal cuticle pattern is probably caused by evolution of the cis-regulatory region of ovo/shaven-baby in the D. sechellia lineage.

Microsatellite analysis reveals substantial genetic differentiation between contemporary New World and Old World Holstein Friesian populations

Genotypic data from 39 microsatellite loci typed in 211 animals were used to assess the genetic differentiation between Old World and New World Holstein Friesian cattle populations. Gene diversities were similar in all five Holstein Friesian populations surveyed, ranging from 0.43 to 0.48. A tree of individuals based on the proportion of shared alleles indicated a clear distinction between Old World and New World Holstein Friesian populations. Similarly, genetic differentiation between populations, as measured by FST, was highly significant. Using the split decomposition method, we were able to visualize the significant introgression of New World Holstein Friesian into European Holstein Friesian populations.

The current state of insect molecular systematics: a thriving Tower of Babel

Insect molecular systematics has undergone remarkable recent growth. Advances in methods of data generation and analysis have led to the accumulation of large amounts of DNA sequence data from most major insect groups. In addition to reviewing theoretical and methodological advances, we have compiled information on the taxa and regions sequenced from all available phylogenetic studies of insects. It is evident that investigators have not usually coordinated their efforts. The genes and regions that have been sequenced differ substantially among studies and the whole of our efforts is thus little greater than the sum of its parts. The cytochrome oxidase I, 16S, 18S, and elongation factor-1 alpha genes have been widely used and are informative across a broad range of divergences in insects. We advocate their use as standards for insect phylogenetics. Insect molecular systematics has complemented and enhanced the value of morphological and ecological data, making substantial contributions to evolutionary biology in the process. A more coordinated approach focused on gathering homologous sequence data will greatly facilitate such efforts.

A quantitative genetic analysis of male sexual traits distinguishing the sibling species Drosophila simulans and D. sechellia

A quantitative trait locus (QTL) genetic analysis of morphological and reproductive traits distinguishing the sibling species Drosophila simulans and D. sechellia was carried out in a backcross design, using 38 markers with an average spacing of 8.4 cM. The direction of QTL effects for the size of the posterior lobe was consistent across the identified QTL, indicating directional selection for this trait. Directional selection also appears to have acted on testis length, indicating that sexual selection may have influenced many reproductive traits, although other forms of directional selection cannot be ruled out. Sex comb tooth number exhibited high levels of variation both within and among isofemale lines and showed no evidence for directional selection and, therefore, may not have been involved in the early speciation process. A database search for genes associated with significant QTL revealed a set of candidate loci for posterior lobe shape and size, sex comb tooth number, testis length, tibia length, and hybrid male fertility. In particular, decapentaplegic (dpp), a gene known to influence the genital arch, was found to be associated with the largest LOD peak for posterior lobe shape and size.