Genetic diversity of populations of Old World screwworm fly, Chrysomya bezziana, causing traumatic myiasis of livestock in the Gulf region and implications for control by sterile insect technique

Abstract.  Fly larvae were collected from 181 cases of traumatic myiasis in livestock in 10 regions of four countries in the Middle East Gulf region: Iran, Iraq, Saudi Arabia and Oman. The predominant fly species responsible for cases was the Old World screwworm (OWS) fly, Chrysomya bezziana Villeneuve (Diptera: Calliphoridae). In cases from Iran and Oman, which included non-OWS fly species, OWS fly was found solely responsible for 67.6% of cases and jointly with other fly species for a further 12.7% of cases. The major hosts were sheep and goats, together comprising 84.6% of the total, which reflects their predominance among the livestock of these Gulf countries. The major site of wounding on sheep and goats was the tail (40.3%), followed by female genitalia (14.0%). The 3' terminal 715 nucleotides of the mitochondrial cytochrome b gene were sequenced for 178 larvae of OWS. Five haplotypes were identified: three had been recorded previously in the region (two were common throughout and one was unique to Oman), and two were newly identified, one from southern Iraq and the other from Saudi Arabia, both in regions sampled for the first time. The haplotypes varied from one another only at one or two nucleotide sites, equivalent to an intraspecific difference of 0.14–0.28% across the entire 715-bp fragment. There was a single statistically significant association between host species and haplotype in Saudi Arabia, a first such record for OWS fly. The small degree of genetic diversity between geographical populations of OWS fly within the Gulf region suggests that a single Gulf colony could be used to implement the sterile insect technique within an integrated control programme.     

Old World screwworm fly, Chrysomya bezziana, occurs as two geographical races

A morphological and molecular analysis was undertaken with the objective of identifying markers for geographical populations of Old World screwworm flies, Chrysomya bezziana Villeneuve (Diptera: Calliphoridae). The morphological analysis involved 192 adult flies from 14 countries, and the molecular analysis involved 45 larvae or adults from 14 populations in 11 countries. Principal components and cluster analysis of 10 morphological characters indicated that flies from Papua New Guinea (PNG) were a distinct group and most similar to flies from nearby Asian islands (Java, Sabah). There was poor resolution of other geographical regions, but some support for clustering of flies from Africa or India. Cladistic analysis of mitochondrial DNA sequences gave strong support for recognizing two races of Old World screwworm, one from sub-Saharan Africa and the other from the Gulf region and Asia. This latter race could be further divided into two lineages, i.e. one from mainland Asia (from Iraq to the Malay Peninsula) and the other from two islands of PNG.     

Mean life expectancy of Old World screwworm fly,Chrysomya bezziana,inferred from the reproductive age-structure of native females caught on swormlure-baited sticky traps

Abstract. The reproductive status of native (wild) screwworm fly, Chrysomya bezziana, caught on swormlure-baited sticky traps in Papua New Guinea is described. A total of 1122 females and 25 males were trapped. Of these females 595 were scored for insemination status and stage of ovarian development (on a scale of 2–10) of which 20% were in suitable condition for assignment to first, second and third ovarian cycles. Of the nulliparous females, only 17% were inseminated at stage 3 of ovarian development, 70% at stage 4, 93-97% at stages 5 and 6, and all of stages 7–10 (gravids). All parous females were inseminated. More than half of the captured females were parous (58%) and only 7% of the total were gravid. Proportions of females in ovarian cycles 1, 2 and 3 were 41%, 50% and 9% respectively. Survival of female Ch. bezziana in the laboratory was adequately described by lognormal and Gompertz survival functions, for both of which the mortality rate is an increasing function of reproductive age. Analysis of the reproductive age distribution of native females estimated their mean life-expectancy at 9 days under the prevailing mean field temperature of 26.5^oC. This equates to completion of 1.7 ovarian cycles and an estimated mean lifetime fecundity of 146 female progeny. The survival models, which also allowed responsiveness of females to swormlure-baited traps (female trappability) to vary according to their stage of ovarian development, indicated significant age-dependent trapping bias. These findings are compared with similar data for the New World screwworm fly, Cochliomyia hominivorax.     

Specific detection of the Old World screwworm fly,Chrysomya bezziana,in bulk fly trap catches using real‐time PCR

The Old World screwworm fly (OWS), Chrysomya bezziana Villeneuve (Diptera: Calliphoridae), is a myiasis‐causing blowfly of major concern for both animals and humans. Surveillance traps are used in several countries for early detection of incursions and to monitor control strategies. Examination of surveillance trap catches is time‐consuming and is complicated by the presence of morphologically similar flies that are difficult to differentiate from Ch. bezziana, especially when the condition of specimens is poor. A molecular‐based method to confirm or refute the presence of Ch. bezziana in trap catches would greatly simplify monitoring programmes. A species‐specific real‐time polymerase chain reaction (PCR) assay was designed to target the ribosomal DNA internal transcribed spacer 1 (rDNA ITS1) of Ch. bezziana. The assay uses both species‐specific primers and an OWS‐specific Taqman^® MGB probe. Specificity was confirmed against morphologically similar and related Chrysomya and Cochliomyia species. An optimal extraction protocol was developed to process trap catches of up to 1000 flies and the assay is sensitive enough to detect one Ch. bezziana in a sample of 1000 non‐target species. Blind testing of 29 trap catches from Australia and Malaysia detected Ch. bezziana with 100% accuracy. The probability of detecting OWS in a trap catch of 50 000 flies when the OWS population prevalence is low (one in 1000 flies) is 63.6% for one extraction. For three extractions (3000 flies), the probability of detection increases to 95.5%. The real‐time PCR assay, used in conjunction with morphology, will greatly increase screening capabilities in surveillance areas where OWS prevalence is low.     

Seasonality of Old World screwworm myiasis in the Mesopotamia valley in Iraq

Following the first recorded introduction of the Old World screwworm fly (OWS), Chrysomya bezziana Villeneuve (Diptera: Calliphoridae), into the Mesopotamia valley in Iraq in September 1996, cases of livestock myiasis caused by OWS developed a distinctly seasonal pattern. The annual cycle of clinical OWS cases is explained here on the basis of environmental variables that affect the different life-cycle stages of C. bezziana. This analysis suggests that low temperatures restricted pupal development during the winter, whereas the dispersal of adult flies was constrained by hot/dry summer conditions. A restricted number of OWS foci persisted throughout the year. In these foci, pupal development was fastest during the autumn months. In autumn, rapid multiplication, lasting several OWS generations, allowed subsequent adult fly dispersal across the valley floor during the winter. Hence, the monthly incidence of clinical OWS cases in livestock peaked during December-January and was lowest during July-August. In addition to temperature and humidity, vegetation cover played a role in OWS distribution. Hence the majority of OWS cases were clustered in the medium density type of vegetation [normalized difference vegetation index (NDVI) values of 0.2-0.4] along the main watercourses in the marshy Mesopotamia valley. Although sheep were the host most commonly infested by C. bezziana, local sheep density was not found to be a major factor in disease spread. Satellite imagery and the application of Geographical Information System (GIS) tools were found to be valuable in understanding the distribution of OWS in relation to vegetation and watercourses. The presence of screwworm in Iraq, at the perimeter of the intercontinental OWS distribution, may give rise to major seasonal flare-ups.     

The potential geographical distribution of the Old World screw-worm fly, Chrysomya bezziana

The potential geographical distribution and relative abundance of the Old World screw-worm fly, Chrysomya bezziana Villeneuve (Diptera: Calliphoridae) as determined by climate, was assessed using CLIMEX, a computer program for matching climates. CLIMEX describes the relative growth and persistence of animal populations in relation to climate. The observed global distribution of C.bezziana was compared with the potential distribution predicted by CLIMEX. The differences in the two distributions indicate the areas at risk of colonization, with particular reference to Australia and the Americas. According to the model, the potential area of permanent colonization in Australia extends south to the mid-coast of New South Wales. Comparison of areas suitable for permanent establishment with the potential summer distribution indicates that large additional areas, carrying most of the continent's livestock, could be colonized in the summer months. Seasonal population growth indices are presented for three ports in Australia at which screw-worm fly specimens have been collected by quarantine authorities. They indicate the relative risk associated with introductions at different places in different seasons and so provide valuable planning information for quarantine authorities. The CLIMEX predictions for C.bezziana in North America are shown to be similar to the recorded distribution limits of the New World screw-worm fly, Cochliomyia hominivorax (Coquerel). The fly could also colonize South America, as far south as southern Brazil and midway through Argentina.     

Phylogenetic analysis of New World screwworm fly, Cochliomyia hominivorax, suggests genetic isolation of some Caribbean island populations following colonization from South America

Abstract.  Larval infestations of the New World screwworm (NWS) fly, Cochliomyia hominivorax , cause considerable economic losses through the direct mortality and reduced production of livestock. Since the 1950s, NWS populations in North and Central America have been the target of virtually continuous eradication attempts by sterile insect technique (SIT). Nevertheless, in some areas, such as Jamaica, SIT-based control programmes have failed. Reasons for the failure of SIT-based programmes in some locations are unknown, but it is hypothesized that failure may be related to the mating incompatibility between sterile and wild flies or to the existence of sexually incompatible cryptic species. Accordingly, the current research investigates intraspecific phylogenetic relationships and associated biogeographic patterns between NWS populations from the Caribbean and South America, which represent those populations involved in, or earmarked for, forthcoming SIT programmes. Uniquely, this study also includes analyses of two North American samples, collected in Texas in 1933 and 1953 prior to initiation of the SIT-based eradication programme. The study utilizes three nucleotide datasets: elongation factor-1α (nuclear); cytochrome oxidase subunit 1 (mitochondrial), and 12S rRNA (mitochondrial). Phylogenetic analysis of these data, representing populations from across the Caribbean, South America and Texas, indicates sub-structuring of fly populations on several of the larger Caribbean islands, suggesting a period of isolation and/or founder effects following colonization from South America; significantly, our findings do not support a North American origin for Cuban flies. The importance of these findings in the light of proposed SIT programmes in the region is discussed.     

基于细胞色素b序列的鲹科分子系统发育

测定了中国鲹科8属9种鱼的细胞色素b基因的全序列(1141bp),结合来自GenBank中分布于美国、安哥拉、希腊以及巴拿马的鲹科4属14种鱼的相应同源序列生成供系统发育分析的序列矩阵,用最大简约法和邻接法构建分子系统树。结果显示:(1)支持科下设四个亚科(鲹亚科,亚科,鲳鲹亚科,鰆鲹科)阶元的分类系统;(2)亚科属下不宜设亚属分类阶元;(3)及达副叶鲹与丽叶鲹亲缘关系较近,应同属于副叶鲹属;(4)我国传统的鱼类检索系统将狮鼻鲳鲹误鉴定为卵形鲳鲹,建议予以修正。     

Evolutionary distinctness of grayling {Thymallus thymallus) inhabiting the Adriatic river system, as based on mtDNA variation

The northern part of the peri-Adriatic river system is inhabited by a native grayling population. In comparison with the Danubian population most of the Adriatic population is regarded as morphologically distinct and therefore recognized as racially differentiated. DNA sequence analysis of the 394-bp fragment of the mtDNA control region and 363-bp fragment of the cytochrome b gene region from grayling specimens of the Adriatic, Danubian and Atlantic river systems revealed nine composite mtDNA haplotypes. The genetic divergence of the Adriatic and the Danubian-Atlantic haplotypes ranged from 2.52 to 3.10% and from 3.38 to 3.92% in the cytochrome b gene and control region, respectively. On the basis of sequence data an attempt was made to elucidate the phylogeography of the Adriatic population. The molecular data indicate an ancient monophyletic origin and also confirm its distinct position within the Thymallus thymallus species, in spite of the fact that stocking with the Danubian type occurred in the Adriatic river system. Therefore the Adriatic population deserves special attention. As a first step in its preservation, its formal distinction by using the name Adriatic grayling is proposed.     

Phylogeographic history of Patagonian lizards of the Liolaemus elongatus complex (Iguania: Liolaemini) based on mitochondrial and nuclear DNA sequences

In this study, we present a phylogeographic analysis of a group of lizards distributed in north‐western Patagonia, the Liolaemus elongatus complex. We sequenced 581 individuals for one mitochondrial gene (cytochrome‐b), and for a subset, we sequenced another mitochondrial gene (12S rRNA) and two nuclear genes: kinesin family member 24 (KIF24) and the anonymous nuclear locus LDAB1D. We estimated gene trees, mitochondrial and nuclear haploytpe networks, standard molecular diversity indices, genetic distances between lineages and Bayesian skyline plots. Our results provide evidence for recognition of seven species previously described within the L. elongatus complex: Liolaemus antumalguen, Liolaemus chillanensis, Liolaemus carlosgarini, Liolaemus burmeisteri, Liolaemus smaug, Liolaemus elongatus and Liolaemus crandalli, but we did not find sufficient evidence to support Liolaemus choique, Liolaemus shitan or Liolaemus sp. 6 as distinct species. We identified four candidate species (Liolaemus sp. 1, Liolaemus sp. 2, Liolaemus sp. 3 and Liolaemus sp. 7), and we discuss evolutionary processes that may have contributed to the origin of these lineages and their taxonomic and conservation implications.     

Phylogeny and distribution of crested gibbons (genus Nomascus) based on mitochondrial cytochrome b gene sequence data

Crested gibbons, genus Nomascus, are endemic to the Indochinese bioregion and occur only in Vietnam, Laos, Cambodia, and southern China. However, knowledge about the number of species to be recognized and their exact geographical distributions is still limited. To further elucidate the evolutionary history of crested gibbon species and to settle their distribution ranges, we analyzed the complete mitochondrial cytochrome b gene from 79 crested gibbon individuals from known locations. Based on our findings, crested gibbons should be classified into seven species. Within N. concolor, we recognize two subspecies, N. concolor concolor and N. concolor lu. Phylogenetic reconstructions indicate that the northernmost species, N. hainanus, N. nasutus, and N. concolor branched off first, suggesting that the genus originated in the north and successively migrated to the south. The most recent splits within Nomascus occurred between N. leucogenys and N. siki, and between Nomascus sp. and N. gabriellae. Based on our data, the currently postulated distributions of the latter four species have to be revised. Our study shows that molecular methods are a useful tool to elucidate phylogenetic relationships among crested gibbons and to determine species boundaries. Am. J. Primatol. 72:1047–1054, 2010. © 2010 Wiley‐Liss, Inc.     

Phylogeography of Coreoperca whiteheadi (Perciformes: Coreoperca) in China based on mitochondrial and nuclear gene sequences

The phylogeography of Coreoperca whiteheadi was studied using sequences from the mitochondrial cytochrome b (cytb) gene and first intron of the S7 ribosomal protein gene. The reciprocal monophyly mitochondrial cytb and nuclear sequences strongly suggest that C. whiteheadi has significantly reduced gene flow between populations and may be in the process of speciation. Phylogenetic analysis revealed phylogeographical structuring into two major clades (highly supported by bootstrap analysis) corresponding to the two regions divided by Nanling-Wuyi Mountain Range. Subsequent nested clade analysis confirmed this structure. The results suggest that Nanling-Wuyi Mountain Range represents a major phylogeographic barrier for C. whiteheadi. Six evolutionarily significant units of C. whiteheadi were designated from the genetic analysis for conservation and management.     

Genetic diversity and demographic history of the Old World Bollworm,Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae), in Ethiopia inferred from mitochondrial gene sequences

The Old World bollworm, Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae), is a globally distributed agricultural and horticultural insect pest. Despite the economic importance of this insect in Ethiopia, its genetic diversity and demographic history are poorly understood. We examined the nucleotide variation of the mitochondrial cytochrome c oxidase subunit I (COI) gene fragment of 74 H. armigera individuals from six collection sites in Ethiopia. We recorded 15 COI haplotypes in H. armigera, ten globally shared and five exclusive to Ethiopia (HaET15, HaET14, HaET10, HaET7, and HaET4). Haplotype HaET1 was the most widely geographically distributed and frequent (71.62%). Analysis of molecular variance (AMOVA) revealed a high and significant level of variation within H. armigera populations (θ_ST  = −0.0135). Negative values of the neutrality test and nonsignificant index of mismatch distribution supported the demographic expansion of H. armigera populations in Ethiopia; furthermore, this was also supported by the nonsignificant values of the sum of squared deviations (SSD) and raggedness index (r). The high genetic variation and population expansion of H. armigera have immense implications for devising locally adapted management strategies in area‐wide integrated pest management IPM programs. However, a comprehensive study of H. armigera genetic diversity and population structure using various molecular markers is needed for future confirmation.     

Analysis of genetic diversity of domestic water buffaloes and anoas based on variations in the mitochondrial gene for cytochrome b

There are two major groups of domestic water buffaloes in East and Southeast Asia: swamp buffaloes and river buffaloes. Genetic diversity among swamp and river buffaloes was studied by DNA sequence analysis of the mitochondrial gene for cytochrome b. The results showed that each of the two groups has mitochondrial DNA (mtDNA) with a specific cytochrome b haplotype. The pairwise nucleotide sequence divergence was calculated to be 2·67% between swamp and river buffaloes, suggesting that they might have diverged from the ancestral populations of Asiatic domestic water buffaloes, approximately 1 million years ago. In addition, the sequences of the same gene from three subspecies of anoa (lowland, mountain and quarles anoa) were determined and compared with that of a domestic water buffalo. The sequence divergence was 1·2% for mountain anoa vs quarles anoa, 3·6% for mountain anoa vs lowland anoa and 3·3% for quarles anoa vs lowland anoa. Moreover, the sequence divergence between water buffaloes and anoas was found to be approximately 3·33%. Our results provide molecular evidence to support the taxonomic classification, namely, that Asiatic buffaloes may be classified into four lineages, swamp buffalo, river buffalo, lowland anoa and mountain plus quarles anoa. However, the sequence divergence values among these four groups were lower than the sequence divergence values found in the genus and subgenus levels within the subfamily Bovinae. In particular, in contrast to some proposed taxonomic classifications, our results indicated that mtDNA in the water buffaloes and anoas did not diverge at the genus level.     

Evidence for three major clades within the snapping shrimp genus Alpheus inferred from nuclear and mitochondrial gene sequence data

The snapping shrimp genus Alpheus is among the most diverse of caridean shrimps, and analyses of taxa separated by the Isthmus of Panama have been used to estimate rates of molecular evolution. Although seven morphological groups have been informally suggested, no formal phylogenetic analysis of the genus has been previously attempted. Here we infer the phylogenetic relationships within Alpheus using sequence data from two nuclear genes, glucose-6-phosphate isomerase and elongation factor-1alpha, and from the mitochondrial gene cytochrome oxidase I. Three major clades corresponding to previously noted morphological features were identified. Discrepancies between earlier informal morphological groupings and molecular analyses largely consisted of species whose morphologies were not entirely typical of the group to which they had been assigned. The traditional placements of shrimp with highly sessile lifestyles and consequently simplified morphologies were also not supported by molecular analyses. Phylogenies for Alpheus suggest that specialized ecological requirements (e.g., symbiotic associations and estuarine habitats) and modified claw morphologies have evolved independently several times. These new analyses also support the sister species status of transisthmian pairs analyzed previously, although very similar pairs were not always resolved with the more slowly evolving nuclear loci. In addition, six new cryptic species were identified in the course of these studies plus a seventh whose status remains to be determined.     

Taxonomy, evolutionary history and biogeography of the broad-toothed field mouse (Apodemus mystacinus) in the eastern Mediterranean area based on mitochondrial and nuclear genes

The broad-toothed field mouse ( Apodemus mystacinus ) is distributed throughout the Balkan Peninsula, Asia Minor and the Middle East. It is generally split into two different specific entities: Apodemus epimelas occurs on the Balkan Peninsula and A. mystacinus inhabits Asia Minor and the Middle East. This analysis, based on two mitochondrial regions (cytochrome b and the D-loop) and the interstitial retinol binding protein (IRBP) nuclear gene, confirms an important level of genetic divergence between the animals from these regions and their separation from each other at least 4.2–5.1 Mya, which is in favour of a distinct specific status. Finally, the broad-toothed field mice from south-western Turkey appear to be closely related to the animals from Crete but highly distinct from the populations of the other Oriental regions. This supports a distinct subspecific level ( A. m. rhodius ) for the insular animals and also for those from south-western Turkey. From a biogeographical point of view, it can be assumed that either late Pliocene or early Pleistocene cooling led to the isolation of two main groups of A. mystacinus , one in the Balkan region and the other one in Turkey and the Near East (Syria and Israel). In this region, it is suggested that a more recent event appeared during the Quaternary period, isolating broad-toothed field mice in Crete and leading to the appearance of two well-differentiated genetic groups: one in Crete and south-western Turkey, and the other widespread in northern and eastern Turkey as well as in the Near East.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 53–63.     

基于Cytb和COⅠ基因部分序列的15种麻蝇之间的系统发育关系(双翅目：麻蝇科)

本研究通过测序Cytb基因和COⅠ基因的部分序列来推定15种麻蝇之间的系统发育关系。在世界麻蝇名录中,本研究的15种麻蝇能够代表麻蝇属Sarcophaga的6个亚属。连接序列（972 bp）被用于系统发育分析;分析方法包括了了最大简约法、最大似然法以及贝叶斯法。我们的结果提示了亚麻蝇亚属Parasarcophaga、别麻蝇亚属Boettcherisca以及红麻蝇亚属Liopygia的单系性,同时也表明蛇麻蝇亚属Liosarcophaga和德麻蝇亚属Pandelleisca并不是单源的。不过,目前的研究并不能分辨野德麻蝇S. (Pandelleisca) similis和峨眉叉麻蝇S. (Robineauella) coei的系统发育位置。此外,最大简约分析和似然功能分析在scopariiformis-iwuensis进化枝和polystylata-hui进化枝的关系上产生了不一致的系统发育推断。因此,后续研究不仅需要其他的分子标记,也需要更多的分类取样。     

Phylogenetic relationships between Tuber pseudoexcavatum, a Chinese truffle, and other Tuber species based on parsimony and distance analysis of four different gene sequences

Phylogenetic relationships between Tuber pseudoexcavatum and other Tuber species were investigated by studying the sequences of four genes: 5.8S-ITS2, beta-tubulin, protein kinase C and elongation factor 1alpha. The four phylogenetic trees allowed to differentiate the black truffle clade, composed of two subclades, one comprising the Asian black truffles (T. indicum, T. sinense, T. himalayense) and the Perigord black truffle (T. melanosporum), the second comprising T. pseudoexcavatum and T. brumale. These two subclades diverged relatively early. We propose a common ancestor, located between Europe and China, to all the black truffles. The T. brumale/pseudoexcavatum subclade would have started to diverge and migrate first, T. brumale towards Europe through a northern route and T. pseudoexcavatum towards China. Later the T. melanosporum subclade would have started to migrate through the same route, T. melanosporum towards Europe and T. indicum towards China, leading to vicariant species.     

基于线粒体ND1和16S rRNA基因序列探讨中国蝴蝶12科的系统发育关系(鳞翅目: 双孔次亚目: 蝶类)

对中国12科共32种代表蝶类的ND1基因和16S rRNA 基因进行了序列测定（包括新测30种ND1基因和9种16S rRNA基因）和比较分析, 同时采用邻接法、最大似然法和贝叶斯法构建了12科蝶类的系统发育树, 探讨了其高级分类群的系统发育关系。序列分析的结果显示: 经比对处理后的两个基因总长度为869 bp, 其中保守位点373个, 可变位点496个, 简约信息位点375个; A+T的平均含量为80.2％, 明显高于C+G的平均含量19.8％。分子系统树表明: 蛱蝶科不是单系群; 珍蝶类、斑蝶类和喙蝶类位于蛱蝶科内; 粉蝶科和凤蝶科具有共同祖先。据此建议: 绢蝶科应归入凤蝶科; 蚬蝶科归入灰蝶科; 珍蝶类、斑蝶类和喙蝶类作为蛱蝶科中的亚科, 眼蝶类从蛱蝶科中分离出来独立成科。另外, 环蝶类的系统分类地位还有待于进一步研究。     

Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences

We sequenced mitochondrial (cytochrome b, 12S rRNA) and nuclear (IRBP, RAG1) genes for 17 species of the Old World murine genus Mus, drawn primarily from the Eurasian subgenus Mus. Phylogenetic analysis of the newly and previously available sequences support recognition of four subgenera within Mus (Mus, Coelomys, Nannomys, and Pyromys), with an unresolved basal polytomy. Our data further indicate that the subgenus Mus contains three distinct 'species groups': (1) a Mus booduga Species Group, also including Mus terricolor and Mus fragilicauda (probably also Mus famulus); (2) a Mus cervicolor Species Group, also including Mus caroli and Mus cookii; and (3) a Mus musculus Species Group, also including Mus macedonicus, Mus spicilegus, and Mus spretus. Species diversity in Eurasian Mus is probably explicable in terms of several phases of range expansion and vicariance, and by a propensity within the group to undergo biotope transitions. IRBP and RAG1 molecular clocks for Mus date the origin of subgenera to around 5-6 mya and the origin of Species Groups within subgenus Mus to around 2-3 mya. The temporal pattern of evolution among Eurasian Mus is more complex than that within the Eurasian temperate genus Apodemus.