Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “synonymous mutation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31Linked to original sources
PubMed2001

Evolution of Mhc class II B genes in Darwin's finches and their closest relatives: birth of a new gene.

The 15 extant species of Darwin's finches on the Galápagos and Cocos Islands are the products of an unfinished adaptive radiation from a founder flock of birds related to the South American species Tiaris obscura. Molecular characterization of their major histocompatibility complex ( Mhc) class II B genes has revealed the existence of several related groups of sequences (presumably encoded in distinct loci) from which one (group 5) stands out because of its low divergence over extended time periods. Analysis of group 5 exon 2 and intron 2 sequences has revealed that the encoding locus apparently arose 2-3 million years ago in the Tiaris group of South and Central American Thraupini. The locus shows no evidence of inactivation, but displays a very low degree of polymorphism, both in terms of number of alleles and genetic distances between alleles. Some of the polymorphism, however, appears to be trans-specific. All the observed intergenic differences can be explained by point mutations and most of the exon 2 changes represent non-synonymous substitutions, although the rate of non-synonymous and synonymous substitutions appears to be the same. The origin of the new locus is explained by the birth-and-death model of Mhc evolution with two important extensions. First, the ancestor of the group 5 genes may have arisen without new gene duplication and second, the birth of the new group may have been brought about by a switch from balancing to directional selection. The ancestor of the group 5 genes may have been a classical class II B allele (one of many) which directional selection fixed in the ancestral population and drove into the category of nonclassical genes.

Animals↗
PubMed2002

Molecular cloning and phylogenetic analysis of new human endogenous retrovirus HERV-W family in cancer cells.

A human endogenous retroviral family (HERV-W) has recently been described that is related to multiple sclerosis-associated retrovirus (MSRV) sequences. By using the PCR approach with human genomic DNA derived from cancer cell lines (HepG2, Jurkat, MCF7, UO-31), five env fragments of HERV-W family were newly identified and analyzed. They showed a high degree of nucleotide sequence similarity (94-99%) with that of the HERV-W. Translation of the env fragments showed no frameshift and termination codon by deletion/insertion or point mutation in clones HepG2-1 and JUR-3. The ratio of synonymous to non-synonymous substitutions indicated that negative selective pressure is acting on HepG2-1 and JUR-3 sequences. These env gene sequences could be associated with an active provirus in human cancer cells (HepG2 and Jurkat). The HepG2-1 and JUR-3 showed sister relationship with the HERV-W and W-7-1 derived from human Chromosome (Chr) 7. Phylogenetic analysis from the HERV-W family indicated close relationships of the env gene sequences across human chromosomes.

Amino Acid Substitution↗
PubMed2006

Detection of feline calicivirus (FCV) from vaccinated cats and phylogenetic analysis of its capsid genes.

We analysed genogroups of four feline calcivirus (FCV) isolates (FCV-S, H10, Ao198-1 and ML89) obtained from cats that experienced FCV infection after having been vaccinated against FCV. New PCR primer sets (8F/8R, Ao-S/Ao-A, cp-S/cp-A) were also designed, since the conventional Seal primer failed to amplify the target sequences in two samples. The genogroups of the four isolates as well as eight global and 17 domestic strains were determined by phylogenetic analysis of their amino acid sequences. One out of the four strains (25%) isolated in this study, H10, was grouped into genogroup I, along with the vaccine strains F9 and FCV-255. The other three isolates (75%) belonged to genogroup II. Thus, there were more isolates in genogroup II than in genogroup I. However, the antibody values of the four isolates against cat anti-F9 antisera were significantly decreased. There may be no relationship between the neutralizing antibody titre and genogroup. Amino acid sequence alignment of the four isolates showed that only a single amino acid in region C, which is involved in neutralization epitopes, was different in ML89 strain from that of F9. The other three strains, H10, Ao198-1 and FCV-B, shared the same amino acid sequence with F9. Alignment of amino acids for linear epitopes in the F9 strain, which are located at regions D and E, showed variations in 5' hypervariable region (HVR) of E, whereas D and conE had only synonymous substitutions i.e. no change in the amino acid sequence. This mutation in 5' HVR of region E suggested a vaccine breakdown, as the region is known to be essential for antigenicity. The genogroup II FCV is likely to be the cause of the FCV infection in this study, while the vaccine strains belong to genogroup I. Thus, the existing vaccine may need reevaluation for its effectiveness.

Amino Acid Sequence↗
PubMedDate not supplied

Genotypic and phenotypic heterogeneity of African Americans with primary iron overload.

Primary iron overload may be relatively common in African Americans, but its cause is incompletely understood. Thus, we evaluated genotype and phenotype characteristics of unselected African American index patients with primary iron overload who reside in central Alabama. All had hepatic iron concentration > or =30 micromol/g dry wt or > or =2.0 g of iron mobilized by phlebotomy to achieve iron depletion. Genotype analyses were performed in African American control subjects from the same region. There were 23 patients (19 men, 4 women); mean age at diagnosis was 52 +/- 12 years (1 SD) (range 32-69 years). Nine (39.1%) reported that they consumed > or =45 g of ethanol daily; five had chronic hepatitis C. Eight had some form of hemoglobinopathy or thalassemia. Mean serum transferrin saturation was 56 +/- 28% (range 15-100%). The geometric mean serum ferritin at diagnosis was 1076 ng/mL [95% confidence interval 297-3473 ng/mL]. Increased stainable liver iron was observed in hepatocytes only in 4 patients, in macrophages only in 8 patients, and in hepatocytes and macrophages in 8 patients. The mean quantity of iron mobilized by phlebotomy (corrected for iron absorbed during treatment) was 5.3 +/- 2.0 g (range 4.0-8.4 g). Iron removed by phlebotomy was greater in patients with hemoglobinopathy or thalassemia than in those without these forms of anemia (6.6 +/- 1.3 g vs 3.9 +/- 1.6 g, respectively; P = 0.0144). Daily consumption of > or =45 g of ethanol or chronic hepatitis C was not associated with an increased or decreased amount of phlebotomy-mobilized iron, on the average. The percentage of index patients positive for HFE C282Y was greater than that of controls (P = 0.0058). The respective percentages of phenotype positivity for HFE H63D, D6S105(8), and HLA-A*03 were similar in patients and controls. HFE S65C, I105T, and G93R were not detected in index or control subjects. Two of 13 patients were heterozygous for the ferroportin allele nt 744 G-->T (Q248H), although the phenotype frequency of this allele was similar in patients and 39 controls. Synonymous ferroportin alleles were also detected in some patients. The ceruloplasmin mutation nt 1099C-->T (exon 6; Arg367Cys) was detected in 1 of 2 patients tested. Abnormal alleles of beta-2 microglobulin, Nramp2, TFR2, hepcidin, or IRP2 alleles were not detected in either of the 2 patients so tested. We conclude that primary iron overload in African Americans is not the result of the mutation of a single gene. HFE C282Y, ferroportin 744 G-->T, and common forms of heritable anemia appear to account for increased iron absorption or retention in some patients.

Adult↗
PubMed2006

Genetic diversity and function in the human cytosolic sulfotransferases.

Amino-acid substitutions, which result from common nonsynonymous (NS) polymorphisms, may dramatically alter the function of the encoded protein. Gaining insight into how these substitutions alter function is a step toward acquiring predictability. In this study, we incorporated gene resequencing, functional genomics, amino-acid characterization and crystal structure analysis for the cytosolic sulfotransferases (SULTs) to attempt to gain predictability regarding the function of variant allozymes. Previously, four SULT genes were resequenced in 118 DNA samples. With additional resequencing of the remaining eight SULT family members in the same DNA samples, a total of 217 polymorphisms were revealed. Of 64 polymorphisms identified within 8785 bp of coding regions from SULT genes examined, 25 were synonymous and 39 were NS. Overall, the proportion of synonymous changes was greater than expected from a random distribution of mutations, suggesting the presence of a selective pressure against amino-acid substitutions. Functional data for common variants of five SULT genes have been previously published. These data, together with the SULT1A1 variant allozyme data presented in this paper, showed that the major mechanism by which amino acid changes altered function in a transient expression system was through decreases in immunoreactive protein rather than changes in enzyme kinetics. Additional insight with regard to mechanisms by which NS single nucleotide polymorphisms alter function was sought by analysis of evolutionary conservation, physicochemical properties of the amino-acid substitutions and crystal structure analysis. Neither individual amino-acid characteristics nor structural models were able to accurately and reliably predict the function of variant allozymes. These results suggest that common amino-acid substitutions may not dramatically alter the protein structure, but affect interactions with the cellular environment that are currently not well understood.

Amino Acid Substitution↗
PubMed2006

A different function for a member of an ancient and highly conserved cytochrome P450 family: from essential sterols to plant defense.

CYP51 sterol demethylases are the only cytochrome P450 enzymes with a conserved function across the animal, fungal, and plant kingdoms (in the synthesis of essential sterols). These highly conserved enzymes, which are important targets for cholesterol-lowering drugs, antifungal agents, and herbicides, are regarded as the most ancient member cytochrome P450 family. Here we present a report of a CYP51 enzyme that has acquired a different function. We show that the plant enzyme AsCYP51H10 is dispensable for synthesis of essential sterols and has been recruited for the production of antimicrobial compounds (avenacins) that confer disease resistance in oats. The AsCyp51H10 gene is synonymous with Sad2, a gene that we previously had defined by mutation as being required for avenacin synthesis. In earlier work, we showed that Sad1, the gene encoding the first committed enzyme in the avenacin pathway (beta-amyrin synthase), had arisen by duplication and divergence of a cycloartenol synthase-like gene. Together these data indicate an intimate evolutionary connection between the sterol and avenacin pathways. Sad1 and Sad2 lie within 70 kb of each other and are expressed specifically in the epidermal cells of the root tip, the site of accumulation of avenacins. These findings raise intriguing questions about the recruitment, coevolution, and regulation of the components of this specialized defense-related metabolic pathway.

Amino Acid Sequence↗
PubMed2003

Natural allelic variants of breast cancer resistance protein (BCRP) and their relationship to BCRP expression in human intestine.

The aim of this study was to identify the extent of genetic variability in breast cancer resistance protein (BCRP) in humans. We first analysed the sequence of BCRP cDNA from human livers and from human intestines phenotyped for expression of intestinal BCRP. We then determined the frequency of all known coding single nucleotide polymorphisms (cSNPs) using DNA from individuals representing 11 different ethnic populations. Nine SNPs including four non-synonymous and three synonymous cSNPs and two intronic SNPs were identified. Of the missense mutations, exon 2 SNP (G34A) resulted in a V12M change; exon 5 SNP (C421A) resulted in a Q141K substitution; exon 6 SNP (A616C) resulted in an I206L amino acid substitution; and exon 15 SNP (A1768T) resulted in a N590Y change in the BCRP protein. The two most frequent polymorphisms identified in the human population studied were the G34A and C421A transitions. There was marked variation in BCRP genotypes and allele frequencies in the different populations. BCRP mRNA was phenotyped in human small bowel intestinal samples by real-time polymerase chain reaction and BCRP protein was analysed on immunoblots of tissue from the same individuals. There was a 78-fold variation in expression of BCRP mRNA and significant variation in BCRP protein expression in human intestine. Expression of intestinal BCRP mRNA and protein was not different between persons expressing the common Gln141 allele compared to the Lys141 allele. Thus, common natural allelic variants of BCRP have been identified, and did not influence interindividual variation in expression of BCRP mRNA in human intestine, but remain to be tested for their effect on BCRP function.

ATP Binding Cassette Transporter, Subfamily G, Mem↗
PubMed2001

Comparison of genomic and predicted amino acid sequences of respiratory and enteric bovine coronaviruses isolated from the same animal with fatal shipping pneumonia.

The complete genome sequences are reported here of two field isolates of bovine coronavirus (BCoV), which were isolated from respiratory and intestinal samples of the same animal experiencing fatal pneumonia during a bovine shipping fever epizootic. Both genomes contained 31028 nucleotides and included 13 open reading frames (ORFs) flanked by 5'- and 3'-untranslated regions (UTRs). ORF1a and ORF1b encode replicative polyproteins pp1a and pp1ab, respectively, that contain all of the putative functional domains documented previously for the closest relative, mouse hepatitis virus. The genomes of the BCoV isolates differed in 107 positions, scattered throughout the genome except the 5'-UTR. Differences in 25 positions were non-synonymous and were located in all proteins except pp1b. Six replicase mutations were identified within or immediately downstream of the predicted largest pp1a-derived protein, p195/p210. Single amino acid changes within p195/p210 as well as within the S glycoprotein might contribute to the different phenotypes of the BCoV isolates.

Amino Acid Sequence↗
PubMed2004

A domain-specific usherin/collagen IV interaction may be required for stable integration into the basement membrane superstructure.

Usherin is a basement membrane protein encoded by the gene associated with Usher syndrome type IIa, the most common deaf/blind disorder. This report demonstrates a specific interaction between type IV collagen and usherin in the basement membrane, with a 1:1 stoichiometry for binding. Genetic and biochemical approaches were used to explore the role of type IV collagen binding in usherin function. We demonstrate binding occurs between the LE domain of usherin and the 7S domain of type IV collagen. A purified fusion peptide comprising the first four LE modules was shown to compete with full-length recombinant usherin for type IV collagen binding. However, synonymous fusion peptides with single amino acid substitutions resulting from missense mutations that were known to cause Usher syndrome type IIa in humans, failed to compete. Only mutations in loop b of the LE domain abolished binding activity. Co-immunoprecipitation and western blot analysis of testicular basement membranes from the Alport mouse model show a 70% reduction in type IV collagen is associated with a similar reduction in usherin, suggesting the usherin/collagen (IV) interaction stabilizes usherin in the basement membrane. Thus, the domain-specific interaction between usherin and type IV collagen appears essential to usherin stability in vivo, and loss of this interaction may result in Usher pathology in humans.

Animals↗
PubMed2002

DNA variation in the acidic chitinase locus (ChiA) region in Arabis gemmifera and its related species.

We analyzed DNA variation at the acidic chitinase (ChiA) locus of Arabis gemmifera and among its eight related species. Nucleotide diversity (pi) of the entire locus in A. gemmifera was 0.0032, which was one third that of A. thaliana. In A. gemmifera, an excess of unique polymorphisms yielded significantly negative results with the tests of Tajima and Fu and Li. The McDonald and Kreitman test revealed that the ratio of nucleotide replacement to synonymous changes in A. gemmifera was significantly greater than those between A. gemmifera and A. glabra, A. gemmifera and A. griffithiana, A. gemmifera and A. korshinskyi, A. gemmifera and A. wallichii, and A. gemmifera and A. himalaica. These results indicated that the neutrality assumption, the equilibrium population assumption, or both, could not be applied to the ChiA locus of A. gemmifera. The small size and relative isolation of local subpopulations of A. gemmifera could explain the excess of unique polymorphisms and the high proportion of replacement changes. The specific sampling scheme of this study, where one strain each was sampled from each local subpopulation might also have led to an excess of singletons. Interspecific comparison among Arabidopsis, Arabis and Cardaminopsis species showed that Ka was always lower than Ks, providing evidence against the adaptive evolution of ChiA. However, Ka/Ks was greater between closely related species than between more distant related species. ChiA had a higher level of replacement divergence and a lower level of synonymous divergence compared than did Adh. We suggest that both the mutation rate at the nucleotide level and the selective constraints at the protein level are lower in ChiA than in Adh.

Analysis of Variance↗
PubMed1989

Anomalous and selective DNA mutations of the Old World monkey alpha-globin genes.

It has been a widely accepted hypothesis that the molecular clock slows down during evolution of higher primates. By molecular cloning and nucleotide sequence comparison of a rhesus macaque alpha-globin gene to its homologs in man, orangutan, olive baboon, and other mammals, we demonstrate a burst of evolution of the baboon alpha-globin gene since its separation from the rhesus macaque. This mutation burst has occurred only at the nonsynonymous sites but not the synonymous sites. Its magnitude is at least 10-fold higher than the synonymous substitution rates in higher primates and as high as the synonymous substitution rates of the rodent lineage. On the contrary, the rate of synonymous site substitutions in the alpha-globin genes of either the rhesus macaque or the olive baboon is several times slower than that of human. Our data demonstrate an anomalous exception to the slow rates of molecular evolution in higher primates and provide strong evidence for a recently accelerated evolution of a primate globin gene under an as yet unknown selective force(s).

Amino Acid Sequence↗
PubMed2006

Absence of coding mutations in the X-linked genes neuroligin 3 and neuroligin 4 in individuals with autism from the IMGSAC collection.

Neuroligin abnormalities have been recently implicated in the aetiology of autism spectrum disorders (ASD), given the finding of point mutations in the two X-linked genes NLGN3 and NLGN4X and the important role of neuroligins in synaptogenesis. To enquire on the relevance and frequency of neuroligin mutations in ASD, we performed a mutation screening of NLGN3 and NLGN4X in a sample of 124 autism probands from the International Molecular Genetic Study of Autism Consortium (IMGSAC). We identified a new non-synonymous variant in NLGN3 (Thr632Ala), which is likely to be a rare polymorphism. Our data indicate that coding mutations in these genes are very rarely associated to ASD.

Autistic Disorder↗
PubMed1988

Molecular evolution of the human and simian immunodeficiency viruses.

Molecular evolution and phylogeny of different human immunodeficiency virus type 1 (HIV1) strains, of a type 2 (HIV2) strain, and of two simian immunodeficiency viruses (SIVAGM and SIVMAC) have been studied by comparing the nucleotide sequences of the two regions of their pol genes which encode the reverse transcriptase (RT) and endonuclease/integrase (EN). The analyses show that the different HIV 1s form one cluster (HIV1 group) and that the SIVs and HIV2 form another (HIV2 group). When the entire genomes of a HIV1, a HIV2, and the two SIVs were compared, the SIVAGM showed a unique pattern of mutation accumulations; that is, the SIVAGM has accumulated more nonsynonymous changes than synonymous changes in the RT and EN regions after its recent divergence from SIVMAC-142, and, furthermore, it has a deletion of approximately 350 bp in the region between the pol and env genes. The SIVAGM was apparently derived from cell cultures infected with a macaque isolate, SIVMAC-251. The contamination provides an opportunity to measure the maximum rate of evolution in the SIVAGM by comparing its DNA sequence to those of SIVMAC-251 and SIVMAC-142. The analysis shows that the rates are given approximately by (1.95 +/- 1.37) x 10(-3)/site/year for one SIVAGM sequence and (5.18 +/- 2.25) x 10(-3)/site/year for another.

Base Sequence↗
PubMed2003

Impact of selection, mutation rate and genetic drift on human genetic variation.

The accumulation of genome-wide information on single nucleotide polymorphisms in humans provides an unprecedented opportunity to detect the evolutionary forces responsible for heterogeneity of the level of genetic variability across loci. Previous studies have shown that history of recombination events has produced long haplotype blocks in the human genome, which contribute to this heterogeneity. Other factors, however, such as natural selection or the heterogeneity of mutation rates across loci, may also lead to heterogeneity of genetic variability. We compared synonymous and non-synonymous variability within human genes with their divergence from murine orthologs. We separately analyzed the non-synonymous variants predicted to damage protein structure or function and the variants predicted to be functionally benign. The predictions were based on comparative sequence analysis and, in some cases, on the analysis of protein structure. A strong correlation between non-synonymous, benign variability and non-synonymous human-mouse divergence suggests that selection played an important role in shaping the pattern of variability in coding regions of human genes. However, the lack of correlation between deleterious variability and evolutionary divergence shows that a substantial proportion of the observed non-synonymous single-nucleotide polymorphisms reduces fitness and never reaches fixation. Evolutionary and medical implications of the impact of selection on human polymorphisms are discussed.

Animals↗
PubMed2002

Low nucleotide diversity at the pal1 locus in the widely distributed Pinus sylvestris.

Nucleotide polymorphism in Scots pine (Pinus sylvestris) was studied in the gene encoding phenylalanine ammonia-lyase (Pal, EC 4.3.1.5). Scots pine, like many other pine species, has a large current population size. The observed levels of inbreeding depression suggest that Scots pine may have a high mutation rate to deleterious alleles. Many Scots pine markers such as isozymes, RFLPs, and microsatellites are highly variable. These observations suggest that the levels of nucleotide variation should be higher than those in other plant species. A 2,045-bp fragment of the pal1 locus was sequenced from five megagametophytes each from a different individual from each of four populations, from northern and southern Finland, central Russia, and northern Spain. There were 12 segregating sites in the locus. The synonymous site overall nucleotide diversity was only 0.0049. In order to compare pal1 with other pine genes, sequence was obtained from two alleles of 11 other loci (total length 4,606 bp). For these, the synonymous nucleotide diversity was 0.0056. These estimates are lower than those from other plants. This is most likely because of a low mutation rate, as estimated from between-pine species synonymous site divergence. In other respects, Scots pine has the characteristics of a species with a large effective population. There was no linkage disequilibrium even between closely linked sites. This resulted in high haplotype diversity (14 different haplotypes among 20 sequences). This could also give rise to high per locus diversity at the protein level. Divergence between populations in the main range was low, whereas an isolated Spanish population had slightly lower diversity and higher divergence than the remaining populations.

ATP-Binding Cassette Transporters↗
PubMed1994

Contrasting rates of nucleotide substitution in the X-linked and Y-linked zinc finger genes.

We have sequenced the entire exon (approximately 1.180 bp) encoding the zinc finger domain of the X-linked and Y-linked zinc finger genes (ZFX and ZFY, respectively) in the orangutan, the baboon, the squirrel monkey, and the rat; a total of 9,442 bp were sequenced. The ratio of the rates of synonymous substitution in the ZFY and ZFX genes is estimated to be 2.1 in primates. This is close to the ratio of 2.3 estimated from primate ZFY and ZFX intron sequences and supports the view that the male-to-female ratio of mutation rate in humans in considerably higher than 1 but not extremely large. The ratio of synonymous substitution rates in ZFY and ZFX is estimated to be 1.3 in the rat lineage but 4.2 in the mouse lineage. The former is close to the estimate (1.4) from introns. The much higher ratio in the mouse lineage (not statistically significant) might have arisen from relaxation of selective constraints. The synonymous divergence between mouse and rat ZFX is considerably lower than that between mouse and rat autosomal genes, agreeing with previous observations and providing some evidence for stronger selective constraints on synonymous changes in X-linked genes than in autosomal genes. At the protein level ZFX has been highly conserved in all placental mammals studied while ZFY has been well conserved in primates and foxes but has evolved rapidly in mice and rats, possibly due to relaxation of functional constraints as a result of the development of X-inactivation of ZFX in rodents. The long persistence of the ZFY-ZFX gene pair in mammals provides some insight into the process of degeneration of Y-linked genes.

Animals↗
PubMed1983

Noise immunity of the genetic code.

Error detection and correction properties are fundamental for informative codes. Hamming's distance allows us to study this noise resistance. We present codes characterized by the resistance optimization to nonsense mutational effects. The calculation of the cumulated Hamming's distance allowing to determine the number of optimal codes and their structure can be detailed. The principle of these laws of optimization of resistance consists of choosing constituent codons connected by mutational neighbouring in such a way that random application of mutations on such a code minimize the occurrence of nonsense n-uplets or terminators. New coding symmetries are then described and screened using Galois's polynomials properties and Baudot's code. Such a study can be applied to any length of the codons. Here we present the principles of this optimization for the most simple doublet codes. Another constraint is discussed: the distribution of optimal subcodes for synonymity and the frequencies of utilization of the different codons. We compare these results to those of the present genetic code, and we observe that all coded amino acids (except the particular case of SER) are using optimal sub-codes of synonymity. This work suggests that the appearance of the genetic code was provoked by mutations while optimizing on several levels its resistance to their effects. Thus genetic coding would have been the best automata that could be produced in prebiotic conditions.

Base Sequence↗
PubMed2000

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.

Animals↗