Worldwide Population Structure Research Articles

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History.

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H.vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H.vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

In Silico Analysis of the Biodiversity and Conservation Status of Mitochondrial Cytochrome C Oxidase Subunit 1 (CO1) Gene of Taenia multiceps.

Taenia multiceps resides in the small intestine of carnivores such as dogs, foxes, woles, jackals, while Coenurus cerebralis which is the larval form usually settle in the central nervous system and spinal cord of intermediate hosts like sheep, cattle and goats. The aim of this study was to analyze the haplotype diversity, genetic variation and population structure of the mt-CO1 gene sequences of sheep isolates of T. multiceps had been submitted to GenBank from different countries. A total of 102 sequences from the mt-CO1 gene fragment belong to T. multiceps sheep isolates were used for bioinformatic analyses. Haplotype analysis, phylogenetic analysis and diversity, neutrality, fixation and gene flow analyses were applied to the sequences. As a result, 20 haplotypes together with different multiple nucleotide changes were determined after the sequence analysis. Trimmed fragment length was 337bp hereby 19 polymorphic areas, 12 of which were parsimony informative, were identified, and any insertion-deletion was found. The number of mutations between major haplotypes and the others range from one to nine. The highest (0.72) genetic differentiation (Fst value) was observed between Turkey and Egypt populations while the lowest (- 0.22) was reported from Greece. These findings are important in terms of showing the diversity of nucleotide variation in T. multiceps sheep isolates. This study serves as the basis for future large-scale studies on T. multiceps worldwide epidemiology, bioecology, geographic distribution and population structure.

Soybean aphid biotype 1 genome: Insights into the invasive biology and adaptive evolution of a major agricultural pest

The soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) is a serious pest of the soybean plant, Glycine max, a major world-wide agricultural crop. We assembled a de novo genome sequence of Ap. glycines Biotype 1, from a culture established shortly after this species invaded North America. 20.4% of the Ap. glycines proteome is duplicated. These in-paralogs are enriched with Gene Ontology (GO) categories mostly related to apoptosis, a possible adaptation to plant chemistry and other environmental stressors. Approximately one-third of these genes show parallel duplication in other aphids. But Ap. gossypii, its closest related species, has the lowest number of these duplicated genes. An Illumina GoldenGate assay of 2380 SNPs was used to determine the world-wide population structure of Ap. Glycines. China and South Korean aphids are the closest to those in North America. China is the likely origin of other Asian aphid populations. The most distantly related aphids to those in North America are from Australia. The diversity of Ap. glycines in North America has decreased over time since its arrival. The genetic diversity of Ap. glycines North American population sampled shortly after its first detection in 2001 up to 2012 does not appear to correlate with geography. However, aphids collected on soybean Rag experimental varieties in Minnesota (MN), Iowa (IA), and Wisconsin (WI), closer to high density Rhamnus cathartica stands, appear to have higher capacity to colonize resistant soybean plants than aphids sampled in Ohio (OH), North Dakota (ND), and South Dakota (SD). Samples from the former states have SNP alleles with high FST values and frequencies, that overlap with genes involved in iron metabolism, a crucial metabolic pathway that may be affected by the Rag-associated soybean plant response. The Ap. glycines Biotype 1 genome will provide needed information for future analyses of mechanisms of aphid virulence and pesticide resistance as well as facilitate comparative analyses between aphids with differing natural history and host plant range.

1232. Phylogenomics of Enterococcus faecium From South America: Revisiting Worldwide VRE Population Structure

BackgroundPrevious studies have suggested that the population structure of E. faecium is composed of two main clades; a commensal clade (designated clade B) and a hospital-associated clade (Clade A) that encompass most of the clinical and animal isolates. The phylogenetic analyses leading to these results have been accomplished with the notable absence of isolates from diverse geographical regions (including South America). We aimed to refine the worldwide population structure of E. faecium by including 55 representative genomes from isolates obtained from five Latin American countries recovered between 1998 and 2014.MethodsWe sequenced our 55 representative isolates and selected other 285 genomes, from public databases, obtained across different regions (36 countries), different sources (animal, commensal, and clinical strains) and a wide range of dates of isolation (1946–2017). We characterized the genomes by presence/absence of resistance, virulence and mobile elements, and of CRISPR-cas systems. We analyzed the phylogeny of the entire population, selected the genomes belonging to clade A to examine recombination patterns and performed Bayesian molecular clock analysis excluding recombinant regions.ResultsTwo major clades were identified, as previously reported. However, a higher degree of variation in clade A was found. Indeed, we identified a subclade (subclade I) that diverged ~894 years ago, and clearly distinguished clinical isolates from those of animal origin (distributed among a number of smaller early-branching subclades). A further split within the clinical subclade (subclade II) that diverged around ~371 years ago was also evident. Latin American isolates were distributed within subclades I (48%) and II (42%). Isolates in “animal” branches exhibited an average recombination of 34 Kbp, where it was 5 Kbp and 21 Kbp for subclades I and II, respectively. More resistance determinants were found in subclade II (62%), followed by I (54%) and absence of cas was the norm in the clinical subclades.ConclusionInclusion of E. faecium isolates from diverse geographical region supports a continuous evolution of these organisms causing human infections. Important evolutionary events seem to favor emergence of novel subclades capable to cause important morbidity and mortality.Disclosures J. Munita, Pfizer: Grant Investigator, Research grant. C. Arias, Merck & Co., Inc.: Grant Investigator, Research support. MeMed: Grant Investigator, Research support. Allergan: Grant Investigator, Research support.

Huanglongbing Control: Perhaps the End of the Beginning.

Huanglongbing (HLB) is one of the most destructive citrus plant diseases worldwide. It is associated with the fastidious phloem-limited α-proteobacteria 'Candidatus Liberibacter asiaticus', 'Ca. Liberibacter africanus' and 'Ca. Liberibacter americanus'. In recent years, HLB-associated Liberibacters have extended to North and South America. The causal agents of HLB have been putatively identified, and their transmission pathways and worldwide population structure have been extensively studied. However, very little is known about the epidemiologic relationships of Ca. L. asiaticus, which has limited the scope of HLB research and especially the development of control strategies. HLB-affected plants produce damaged fruits and die within several years. To control the disease, scientists have developed new compounds and screened existing compounds for their antibiotic and antimicrobial activities against the disease. These compounds, however, have very little or even no effect on the disease. The aim of the present review was to compile and compare different methods of HLB disease control with newly developed integrative strategies. In light of recent studies, we also describe how to control the vectors of this disease and the biological control of other citrus plant pathogens. This work could steer the attention of scientists towards integrative control strategies.

Worldwide population structure of the wheat rust fungus Puccinia striiformis in the past

Puccinia striiformis is a basidiomycete causing yellow rust on wheat. The availability of historic samples of this pathogen from the ‘Stubbs collection’ enabled us to investigate past population structure and temporal dynamics on a global scale. A set of 212 single genotype urediniospore isolates, representing samples collected from five continents between 1958 and 1991, were genotyped using 19 polymorphic microsatellite markers. The population genetic analyses revealed the existence of seven genetic groups in the past worldwide P. striiformis population. This genetic grouping generally corresponded with geographical sample origin except for the Middle East, where six of the seven genetic groups were represented. The presence of many genetic groups in the Middle Eastern population reflected a low differentiation from the populations in East Africa (FST=0.052) and in South Asia (FST=0.064). A high diversity and recombinant population structure was observed in China and South Asia, while a clonal population structure was observed in NW Europe, East Africa and the Mediterranean region. The high genetic diversity in the Himalayan region supported recent studies suggesting a putative center of diversity for P. striiformis in this area. Four of the 89 multilocus genotypes detected were resampled in different geographical regions suggesting long-distance migration in the past. Comparison of the past populations with more recent ones, represented by 309 isolates mainly collected between 2001 and 2009, revealed temporal divergence for all populations except for Northwest Europe. Overall, we observed a clear subdivision within the worldwide population structure of P. striiformis and migration in the past.

A comparative study on worldwide genetic diversity and population structure analysis of Plasmodium vivax thrombospondin-related adhesive protein (PvTRAP) and its implications for the vivax vaccine design

Plasmodium vivax thrombospondin-related anonymous protein (PvTRAP) is a promising malaria vaccine candidate; however, it exhibits sequence heterogeneity. Therefore, to design a broadly protective vivax vaccine, it is essential to have adequate information on signatures of selection and geospatial genetic diversity of global PvTRAP. For this purpose, 50 Iranian pvtrap were sequenced and compared with related available global sequences in GenBank. The nucleotide sequence analysis of Iranian pvtrap in comparison with the Sal-1 sequence showed the occurrence of 15 SNPs, and all sites were dimorphic. In total, 12 amino acid substitutions were detected and 2 of which were novel, resulting in 10 haplotypes that 8 of them were not reported in any other geographic regions. In comparison with global population, haplotype and nucleotide diversities were lowest in South Korean populations while higher levels of diversities were observed in Thai and Brazilian P. vivax populations. All 12 amino acid replacements in ectodomain of Iranian PvTRAP were distributed in predicted either B- or T-cells epitope as well as intrinsically unstructured/disordered regions (IURs). The present results revealed that observing the relatively low-level diversity in PvTRAP protein might actually be selected by immune response. In summary, the present analysis in parallel to the limited available published data has shown that genetic diversity in the global pvtrap exhibits low-level diversity and geographic variation. These results are of practical significance for the strategic development and deployment of control measures in particular for development of PvTRAP-based malaria vaccine.

Worldwide Population Structure, Long-Term Demography, and Local Adaptation of Helicobacter pylori.

Helicobacter pylori is an important human pathogen associated with serious gastric diseases. Owing to its medical importance and close relationship with its human host, understanding genomic patterns of global and local adaptation in H. pylori may be of particular significance for both clinical and evolutionary studies. Here we present the first such whole genome analysis of 60 globally distributed strains, from which we inferred worldwide population structure and demographic history and shed light on interesting global and local events of positive selection, with particular emphasis on the evolution of San-associated lineages. Our results indicate a more ancient origin for the association of humans and H. pylori than previously thought. We identify several important perspectives for future clinical research on candidate selected regions that include both previously characterized genes (e.g., transcription elongation factor NusA and tumor necrosis factor alpha-inducing protein Tipα) and hitherto unknown functional genes.

Inferring population structure and demographic history using Y-STR data from worldwide populations.

The Y chromosome is one of the best genetic materials to explore the evolutionary history of human populations. Global analyses of Y chromosomal short tandem repeats (STRs) data can reveal very interesting world population structures and histories. However, previous Y-STR works tended to focus on small geographical ranges or only included limited sample sizes. In this study, we have investigated population structure and demographic history using 17 Y chromosomal STRs data of 979 males from 44 worldwide populations. The largest genetic distances have been observed between pairs of African and non-African populations. American populations with the lowest genetic diversities also showed large genetic distances and coancestry coefficients with other populations, whereas Eurasian populations displayed close genetic affinities. African populations tend to have the oldest time to the most recent common ancestors (TMRCAs), the largest effective population sizes and the earliest expansion times, whereas the American, Siberian, Melanesian, and isolated Atayal populations have the most recent TMRCAs and expansion times, and the smallest effective population sizes. This clear geographic pattern is well consistent with serial founder model for the origin of populations outside Africa. The Y-STR dataset presented here provides the most detailed view of worldwide population structure and human male demographic history, and additionally will be of great benefit to future forensic applications and population genetic studies.

Evidence for genetic diversity of Toxoplasma gondii in selected intermediate hosts in Serbia

To contribute to the insight into the worldwide population structure of Toxoplasma gondii, we genetically characterized a total of eight strains isolated from intermediate hosts including humans, sheep and pigeons in Serbia. Although parasite DNA was detected in 28.2% (60/213) of the human samples from 162 patients serologically suspected of active toxoplasmosis, as well as in 5/7 seropositive pigeons and in 2/12 seropositive sheep examined, multilocus PCR-RFLP genotyping, using SAG1, 5'SAG2, 3'SAG2, GRA6, 5'GRA7 and 3'GRA7 as markers, was successful in only four human isolates (of which one was isolated from both the bronchoalveolar lavage fluid and blood samples of a single patient), one sheep and three pigeons. Of the eight isolates, five were type II (62.5%), one was type III, one was atypical, and one had a type I allele at GRA6 as the single locus genotyped. Although type II, as elsewhere in Europe, predominated, these results may suggest a higher genetic diversity of T. gondii in Serbia, reflecting local environmental contamination and also the geographical position of the country in South-East Europe.

Origin, migration routes and worldwide population genetic structure of the wheat yellow rust pathogen Puccinia striiformis f.sp. tritici.

Analyses of large-scale population structure of pathogens enable the identification of migration patterns, diversity reservoirs or longevity of populations, the understanding of current evolutionary trajectories and the anticipation of future ones. This is particularly important for long-distance migrating fungal pathogens such as Puccinia striiformis f.sp. tritici (PST), capable of rapid spread to new regions and crop varieties. Although a range of recent PST invasions at continental scales are well documented, the worldwide population structure and the center of origin of the pathogen were still unknown. In this study, we used multilocus microsatellite genotyping to infer worldwide population structure of PST and the origin of new invasions based on 409 isolates representative of distribution of the fungus on six continents. Bayesian and multivariate clustering methods partitioned the set of multilocus genotypes into six distinct genetic groups associated with their geographical origin. Analyses of linkage disequilibrium and genotypic diversity indicated a strong regional heterogeneity in levels of recombination, with clear signatures of recombination in the Himalayan (Nepal and Pakistan) and near-Himalayan regions (China) and a predominant clonal population structure in other regions. The higher genotypic diversity, recombinant population structure and high sexual reproduction ability in the Himalayan and neighboring regions suggests this area as the putative center of origin of PST. We used clustering methods and approximate Bayesian computation (ABC) to compare different competing scenarios describing ancestral relationship among ancestral populations and more recently founded populations. Our analyses confirmed the Middle East-East Africa as the most likely source of newly spreading, high-temperature-adapted strains; Europe as the source of South American, North American and Australian populations; and Mediterranean-Central Asian populations as the origin of South African populations. Although most geographic populations are not markedly affected by recent dispersal events, this study emphasizes the influence of human activities on recent long-distance spread of the pathogen.

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over-exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded FST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between-ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17-SNP panel was developed in Atlantic BFT by cross-species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (FST =0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.

Human loci involved in drug biotransformation: worldwide genetic variation, population structure, and pharmacogenetic implications

Understanding the role of inheritance in individual variation in drug response is the focus of pharmacogenetics (PGx). A key part of this understanding is quantifying the role of genetic ancestry in this phenotypic outcome. To provide insight into the relationship between ethnicity and drug response, this study first infers the global distribution of PGx variation and defines its structure. Second, the study evaluates if geographic population structure stems from all PGx loci in general, or if structure is caused by specific genes. Lastly, we identify the genetic variants contributing the greatest proportion of such structure. Our study describes the global genetic structure of PGx loci across the 52 populations of the Human Genome Diversity Cell-Line Panel, the most inclusive set of human populations freely available for studies on human genetic variation. By analysing genetic variation at 1,001 single nucleotide polymorphisms (SNPs) involved in biotransformation of exogenous substances, we describe the between-populations PGx variation, as well geographical groupings of diversity. In addition, with discriminant analysis of principal component (DAPC), we infer how many and which groups of populations are supported by PGx variation, and identify which SNPs actually contribute to the PGx structure between such groups. Our results show that intergenic, synonymous and non-synonymous SNPs show similar levels of genetic variation across the globe. Conversely, loci coding for Cytochrome P450s (mainly metabolizing exogenous substances) show significantly higher levels of genetic diversity between populations than the other gene categories. Overall, genetic variation at PGx loci correlates with geographic distances between populations, and the apportionment of genetic variation is similar to that observed for the rest of the genome. In other words, the pattern of PGx variation has been mainly shaped by the demographic history of our species, as in the case of most of our genes. The population structure defined by PGx loci supports the presence of six genetic clusters reflecting geographic location of samples. In particular, the results of the DAPC analyses show that 27 SNPs substantially contribute to the first three discriminant functions. Among these SNPs, some, such as the intronic rs1403527 of NR1I2 and the non-synonymous rs699 of AGT, are known to be associated with specific drug responses. Their substantial variation between different groups of populations may have important implications for PGx practical applications.

Long-term population size of the North Atlantic humpback whale within the context of worldwide population structure

Once hunted to the brink of extinction, humpback whales (Megaptera novaeangliae) in the North Atlantic have recently been increasing in numbers. However, uncertain information on past abundance makes it difficult to assess the extent of the recovery in this species. While estimates of pre-exploitation abundance based upon catch data suggest the population might be approaching pre-whaling numbers, estimates based on mtDNA genetic diversity suggest they are still only a fraction of their past abundance levels. The difference between the two estimates could be accounted for by inaccuracies in the catch record, by uncertainties surrounding the genetic estimate, or by differences in the timescale to which the two estimates apply. Here we report an estimate of long-term population size based on nuclear gene diversity. We increase the reliability of our genetic estimate by increasing the number of loci, incorporating uncertainty in each parameter and increasing sampling across the geographic range. We report an estimate of long-term population size in the North Atlantic humpback of ~112,000 individuals (95 % CI 45,000–235,000). This value is 2–3 fold higher than estimates based upon catch data. This persistent difference between estimates parallels difficulties encountered by population models in explaining the historical crash of North Atlantic humpback whales. The remaining discrepancy between genetic and catch-record values, and the failure of population models, highlights a need for continued evaluation of whale population growth and shifts over time, and continued caution about changing the conservation status of this population.

Multilocus microsatellite analysis of 'Candidatus Liberibacter asiaticus' associated with citrus Huanglongbing worldwide.

BackgroundHuanglongbing (HLB) is one of the most destructive citrus diseases in the world. The disease is associated with the presence of a fastidious, phloem-limited α- proteobacterium, 'Candidatus Liberibacter asiaticus', 'Ca. Liberibacter africanus' or 'Ca. Liberibacter americanus'. HLB-associated Liberibacters have spread to North America and South America in recent years. While the causal agents of HLB have been putatively identified, information regarding the worldwide population structure and epidemiological relationships for 'Ca. L. asiaticus' is limited. The availability of the 'Ca. L. asiaticus' genome sequence has facilitated development of molecular markers from this bacterium. The objectives of this study were to develop microsatellite markers and conduct genetic analyses of 'Ca. L. asiaticus' from a worldwide collection. Two hundred eighty seven isolates from USA (Florida), Brazil, China, India, Cambodia, Vietnam, Taiwan, Thailand, and Japan were analyzed.ResultsA panel of seven polymorphic microsatellite markers was developed for 'Ca. L. asiaticus'. Microsatellite analyses across the samples showed that the genetic diversity of 'Ca. L. asiaticus' is higher in Asia than Americas. UPGMA and STRUCTURE analyses identified three major genetic groups worldwide. Isolates from India were genetically distinct. East-southeast Asian and Brazilian isolates were generally included in the same group; a few members of this group were found in Florida, but the majority of the isolates from Florida were clustered separately. eBURST analysis predicted three founder haplotypes, which may have given rise to three groups worldwide.ConclusionsOur results identified three major genetic groups of 'Ca. L. asiaticus' worldwide. Isolates from Brazil showed similar genetic makeup with east-southeast Asian dominant group, suggesting the possibility of a common origin. However, most of the isolates recovered from Florida were clustered in a separate group. While the sources of the dominant 'Ca. L. asiaticus' in Florida were not clearly understood, the less-pervasive groups may have been introduced directly from Asia or via Brazil. Notably, the recent outbreak of HLB in Florida probably occurred through multiple introductions. Microsatellite markers developed in this study provide adequate discriminatory power for the identification and differentiation of closely-related isolates, as well as for genetic studies of 'Ca. L. asiaticus'.

Contrasting introduction scenarios among continents in the worldwide invasion of the banana fungal pathogen Mycosphaerella fijiensis

Reconstructing and characterizing introduction routes is a key step towards understanding the ecological and evolutionary factors underlying successful invasions and disease emergence. Here, we aimed to decipher scenarios of introduction and stochastic demographic events associated with the global spread of an emerging disease of bananas caused by the destructive fungal pathogen Mycosphaerella fijiensis. We analysed the worldwide population structure of this fungus using 21 microsatellites and 8 sequence-based markers on 735 individuals from 37 countries. Our analyses designated South-East Asia as the source of the global invasion and supported the location of the centre of origin of M.fijiensis within this area. We confirmed the occurrence of bottlenecks upon introduction into other continents followed by widespread founder events within continents. Furthermore, this study suggested contrasting introduction scenarios of the pathogen between the African and American continents. While potential signatures of admixture resulting from multiple introductions were detected in America, all the African samples examined seem to descend from a single successful founder event. In combination with historical information, our study reveals an original and unprecedented global scenario of invasion for this recently emerging disease caused by a wind-dispersed pathogen.

Worldwide mitochondrial DNA diversity and phylogeography of pilot whales (Globicephala spp.)

Pilot whales (Globicephala spp.) provide an interesting example of recently diverged oceanic species with a complex evolutionary history. The two species have wide but largely non-overlapping ranges. Globicephala melas (long-finned pilot whale; LFPW) has an antitropical distribution and is found in the cold-temperate waters of the North Atlantic and Southern Hemisphere, whereas Globicephala macrorhynchus (short-finned pilot whale; SFPW) has a circumglobal distribution and is found mainly in the tropics and subtropics. To investigate pilot whale evolution and biogeography, we analysed worldwide population structure using mitochondrial DNA (mtDNA) control region sequences (up to 620 bp) from a variety of sources (LFPW = 643; SFPW = 150), including strandings in New Zealand and Tasmania, and whale-meat products purchased on the markets of Japan and Korea. Phylogenetic reconstructions failed to support a reciprocal monophyly of the two species, despite six diagnostic substitutions, possibly because of incomplete lineage sorting or inadequate phylogenetic information. Both species had low haplotype and nucleotide diversity compared to other abundant widespread cetaceans (LFPW, π = 0.35%; SFPW, π = 0.87%) but showed strong mtDNA differentiation between oceanic basins. Strong levels of structuring were also found at the regional level. In LFPW, phylogeographic patterns were suggestive either of a recent demographic expansion or selective sweep acting on the mtDNA. For SFPW, the waters around Japan appear to represent a centre of diversity, with two genetically-distinct forms, as well as a third population of unknown origin. The presence of multiple unique haplotypes among SFPW from South Japan, together with previously documented morphological and ecological differences, suggests that the southern form represents a distinct subspecies and/or evolutionary significant unit. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98, 729–744.

Pathway Analysis of <I>Ceratitis capitata</I> (Diptera: Tephritidae) Using Mitochondrial DNA

The polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method described in a previous publication to document worldwide population structure of Mediterranean fruit flies and analyze probable sources of fly introductions is re-evaluated. Analysis of representatives of the RFLP-based markers using the ND5-tRNA-ND4 mitochondrial locus through DNA sequencing shows that the original markers produced haplotypes that are homoplasic and exclude valuable information. DNA sequences also identified errors in data generated with the protocol because of failure of the XbaI enzyme to efficiently cut PCR product. In all, the PCR-RFLP method is able to document haplotypes present in medfly populations outside of Sub-Saharan Africa and is a cost-effective approach to pathway analysis. New primer pairs for analysis of this locus are reported that reduce the number of PCR steps required to analyze the RFLP marker and the time of each PCR run. For comparison of the PCR-RFLP method, the utility of a DNA sequence-based approach to Mediterranean fruit fly pathway analysis is reported using the ND5-ND4 locus and the proposed "barcode region" of the COI gene. Both gene sequences are more informative than the PCR-RFLP method and document high genetic variation within Sub-Saharan Africa. Application of these methods as tools and their interpretation is discussed.

Worldwide distribution of NAT2 diversity: Implications for NAT2 evolutionary history

BackgroundThe N-acetyltransferase 2 (NAT2) gene plays a crucial role in the metabolism of many drugs and xenobiotics. As it represents a likely target of population-specific selection pressures, we fully sequenced the NAT2 coding region in 97 Mandenka individuals from Senegal, and compared these sequences to extant data on other African populations. The Mandenka data were further included in a worldwide dataset composed of 41 published population samples (6,727 individuals) from four continental regions that were adequately genotyped for all common NAT2 variants so as to provide further insights into the worldwide haplotype diversity and population structure at NAT2.ResultsThe sequencing analysis of the NAT2 gene in the Mandenka sample revealed twelve polymorphic sites in the coding exon (two of which are newly identified mutations, C345T and C638T), defining 16 haplotypes. High diversity and no molecular signal of departure from neutrality were observed in this West African sample. On the basis of the worldwide genotyping survey dataset, we found a strong genetic structure differentiating East Asians from both Europeans and sub-Saharan Africans. This pattern could result from region- or population-specific selective pressures acting at this locus, as further suggested in the HapMap data by extremely high values of FST for a few SNPs positions in the NAT2 coding exon (T341C, C481T and A803G) in comparison to the empirical distribution of FST values accross the whole 400-kb region of the NAT gene family.ConclusionPatterns of sequence variation at NAT2 are consistent with selective neutrality in all sub-Saharan African populations investigated, whereas the high level of population differentiation between Europeans and East Asians inferred from SNPs could suggest population-specific selective pressures acting at this locus, probably caused by differences in diet or exposure to other environmental signals.