PCA and Clustering Reveal Alternate mtDNA Phylogeny of N and M Clades

Phylogenetic trees based on mtDNA polymorphisms are often used to infer the history of recent human migrations. However, there is no consensus on which method to use. Most methods make strong assumptions which may bias the choice of polymorphisms and result in computational complexity which limits the analysis to a few samples/polymorphisms. For example, parsimony minimizes the number of mutations, which biases the results to minimizing homoplasy events. Such biases may miss the global structure of the polymorphisms altogether, with the risk of identifying a "common" polymorphism as ancient without an internal check on whether it either is homoplasic or is identified as ancient because of sampling bias (from oversampling the population with the polymorphism). A signature of this problem is that different methods applied to the same data or the same method applied to different datasets results in different tree topologies. When the results of such analyses are combined, the consensus trees have a low internal branch consensus. We determine human mtDNA phylogeny from 1737 complete sequences using a new, direct method based on principal component analysis (PCA) and unsupervised consensus ensemble clustering. PCA identifies polymorphisms representing robust variations in the data and consensus ensemble clustering creates stable haplogroup clusters. The tree is obtained from the bifurcating network obtained when the data are split into k = 2,3,4,...,kmax clusters, with equal sampling from each haplogroup. Our method assumes only that the data can be clustered into groups based on mutations, is fast, is stable to sample perturbation, uses all significant polymorphisms in the data, works for arbitrary sample sizes, and avoids sample choice and haplogroup size bias. The internal branches of our tree have a 90% consensus accuracy. In conclusion, our tree recreates the standard phylogeny of the N, M, L0/L1, L2, and L3 clades, confirming the African origin of modern humans and showing that the M and N clades arose in almost coincident migrations. However, the N clade haplogroups split along an East-West geographic divide, with a "European R clade" containing the haplogroups H, V, H/V, J, T, and U and a "Eurasian N subclade" including haplogroups B, R5, F, A, N9, I, W, and X. The haplogroup pairs (N9a, N9b) and (M7a, M7b) within N and M are placed in nonnearest locations in agreement with their expected large TMRCA from studies of their migrations into Japan. For comparison, we also construct consensus maximum likelihood, parsimony, neighbor joining, and UPGMA-based trees using the same polymorphisms and show that these methods give consistent results only for the clade tree. For recent branches, the consensus accuracy for these methods is in the range of 1-20%. From a comparison of our haplogroups to two chimp and one bonobo sequences, and assuming a chimp-human coalescent time of 5 million years before present, we find a human mtDNA TMRCA of 206,000 +/- 14,000 years before present.     

Molecular Dissection of the Basal Clades in the Human Y Chromosome Phylogenetic Tree

One hundred and forty-six previously detected mutations were more precisely positioned in the human Y chromosome phylogeny by the analysis of 51 representative Y chromosome haplogroups and the use of 59 mutations from literature. Twenty-two new mutations were also described and incorporated in the revised phylogeny. This analysis made it possible to identify new haplogroups and to resolve a deep trifurcation within haplogroup B2. Our data provide a highly resolved branching in the African-specific portion of the Y tree and support the hypothesis of an origin in the north-western quadrant of the African continent for the human MSY diversity.     

MtDNA and Y-Chromosome Lineages in the Yakut Population

The structure of female (mtDNA) and male (Y-chromosome haplotypes) lineages in the Yakut population was examined. To determine mtDNA haplotypes, sequencing of hypervariable segment I and typing of haplotype-specific point substitutions in the other parts of the mtDNA molecule were performed. Y haplogroups were identified through typing of biallelic polymorphisms in the nonrecombining part of the chromosome. Haplotypes within haplogroups were analyzed with seven microsatellite loci. Mitochondrial gene pool of Yakuts is mainly represented by the lineages of eastern Eurasian origin (haplogroups A, B, C, D, G, and F). In Yakuts haplogroups C and D showing the total frequency of almost 80% and consisting of 12 and 10 different haplopypes, respectively, were the most frequent and diverse. The total part of the lineages of western Eurasian origin (Caucasoid) was about 6% (4 haplotypes, haplogroups H, J, and U). Most of Y chromosomes in the Yakut population (87%) belonged to haplogroup N3 (HG16), delineated by the T–C substitution at the Tat locus. Chromosomes of haplogroup N3 displayed the presence of 19 microsatellite haplotypes, the most frequent of which encompassed 54% chromosomes of this haplogroup. Median network of haplogroup N3 in Yakuts demonstrated distinct starlike phylogeny. Male lineages of Yakuts were shown to be closest to those of Eastern Evenks.     

Genetic Structure and Antimicrobial Resistance of Escherichia coli and Cryptic Clades in Birds with Diverse Human Associations

The manner and extent to which birds associate with humans may influence the genetic attributes and antimicrobial resistance of their commensal Escherichia communities through strain transmission and altered selection pressures. In this study, we determined whether the distribution of the different Escherichia coli phylogenetic groups and cryptic clades, the occurrence of 49 virulence associated genes, and/or the prevalence of resistance to 12 antimicrobials differed between four groups of birds from Australia with contrasting types of human association. We found that birds sampled in suburban and wilderness areas had similar Escherichia communities. The Escherichia communities of backyard domestic poultry were phylogenetically distinct from the Escherichia communities sourced from all other birds, with a large proportion (46%) of poultry strains belonging to phylogenetic group A and a significant minority (17%) belonging to the cryptic clades. Wild birds sampled from veterinary and wildlife rehabilitation centers (in-care birds) carried Escherichia isolates that possessed particular virulence-associated genes more often than Escherichia isolates from birds sampled in suburban and wilderness areas. The Escherichia isolates from both the backyard poultry and in-care birds were more likely to be multidrug resistant than the Escherichia isolates from wild birds. We also detected a multidrug-resistant E. coli strain circulating in a wildlife rehabilitation center, reinforcing the importance of adequate hygiene practices when handling and caring for wildlife. We suggest that the relatively high frequency of antimicrobial resistance in the in-care birds and backyard poultry is due primarily to the use of antimicrobials in these animals, and we recommend that the treatment protocols used for these birds be reviewed.     

Tracing the Route of Modern Humans out of Africa by Using 225 Human Genome Sequences from Ethiopians and Egyptians

The predominantly African origin of all modern human populations is well established, but the route taken out of Africa is still unclear. Two alternative routes, via Egypt and Sinai or across the Bab el Mandeb strait into Arabia, have traditionally been proposed as feasible gateways in light of geographic, paleoclimatic, archaeological, and genetic evidence. Distinguishing among these alternatives has been difficult. We generated 225 whole-genome sequences (225 at 8× depth, of which 8 were increased to 30×; Illumina HiSeq 2000) from six modern Northeast African populations (100 Egyptians and five Ethiopian populations each represented by 25 individuals). West Eurasian components were masked out, and the remaining African haplotypes were compared with a panel of sub-Saharan African and non-African genomes. We showed that masked Northeast African haplotypes overall were more similar to non-African haplotypes and more frequently present outside Africa than were any sets of haplotypes derived from a West African population. Furthermore, the masked Egyptian haplotypes showed these properties more markedly than the masked Ethiopian haplotypes, pointing to Egypt as the more likely gateway in the exodus to the rest of the world. Using five Ethiopian and three Egyptian high-coverage masked genomes and the multiple sequentially Markovian coalescent (MSMC) approach, we estimated the genetic split times of Egyptians and Ethiopians from non-African populations at 55,000 and 65,000 years ago, respectively, whereas that of West Africans was estimated to be 75,000 years ago. Both the haplotype and MSMC analyses thus suggest a predominant northern route out of Africa via Egypt.     

Phylogenetic Framework and Molecular Signatures for the Main Clades of the Phylum Actinobacteria

Summary: The phylum Actinobacteria harbors many important human pathogens and also provides one of the richest sources of natural products, including numerous antibiotics and other compounds of biotechnological interest. Thus, a reliable phylogeny of this large phylum and the means to accurately identify its different constituent groups are of much interest. Detailed phylogenetic and comparative analyses of >150 actinobacterial genomes reported here form the basis for achieving these objectives. In phylogenetic trees based upon 35 conserved proteins, most of the main groups of Actinobacteria as well as a number of their superageneric clades are resolved. We also describe large numbers of molecular markers consisting of conserved signature indels in protein sequences and whole proteins that are specific for either all Actinobacteria or their different clades (viz., orders, families, genera, and subgenera) at various taxonomic levels. These signatures independently support the existence of different phylogenetic clades, and based upon them, it is now possible to delimit the phylum Actinobacteria (excluding Coriobacteriia) and most of its major groups in clear molecular terms. The species distribution patterns of these markers also provide important information regarding the interrelationships among different main orders of Actinobacteria. The identified molecular markers, in addition to enabling the development of a stable and reliable phylogenetic framework for this phylum, also provide novel and powerful means for the identification of different groups of Actinobacteria in diverse environments. Genetic and biochemical studies on these Actinobacteria-specific markers should lead to the discovery of novel biochemical and/or other properties that are unique to different groups of Actinobacteria.     

Molecular Evidence for the Major Clades of Placental Mammals

Higher-level relationships among placental mammals, as well as the historical biogeography of this group against the backdrop of continental fragmentation and reassembly, remain poorly understood. Here, we analyze two independent molecular data sets that represent all placental orders. The first data set includes six genes (A2AB, IRBP, vWF, 12S rRNA, tRNA valine, 16S rRNA; total = 5.71 kb) for 26 placental taxa and two marsupials; the second data set includes 2.95 kb of exon 11 of the BRCA1 gene for 51 placental taxa and four marsupials. We also analyzed a concatenation of these data sets (8.66 kb) for 26 placentals and one marsupial. Unrooted and rooted analyses were performed with parsimony, distance methods, maximum likelihood, and a Bayesian approach. Unrooted analyses provide convincing support for a fundamental separation of placental orders into groups with southern and northern hemispheric origins according to the current fossil record. On rooted trees, one or both of these groups are monophyletic depending on the position of the root. Maximum likelihood and Bayesian analyses with the BRCA1 and combined 8.66 kb data sets provide strong support for the monophyly of the northern hemisphere group (Boreoeutheria). Boreoeutheria is divided into Laurasiatheria (Carnivora + Cetartiodactyla + Chiroptera + Eulipotyphla + Perissodactyla + Pholidota) and Euarchonta (Dermoptera + Primates + Scandentia) + Glires (Lagomorpha + Rodentia). The southern hemisphere group is either monophyletic or paraphyletic, depending on the method of analysis used. Within this group, Afrotheria (Proboscidea + Sirenia + Hyracoidea + Tubulidentata + Macroscelidea + Afrosoricida) is monophyletic. A unique nine base-pair deletion in exon 11 of the BRCA1 gene also supports Afrotheria monophyly. Given molecular dates that suggest that the southern hemisphere group and Boreoeutheria diverged in the Early Cretaceous, a single trans-hemispheric dispersal event may have been of fundamental importance in the early history of crown-group Eutheria. Parallel adaptive radiations have subsequently occurred in the four major groups: Laurasiatheria, Euarchonta + Glires, Afrotheria, and Xenarthra.     

Mycobacterium tuberculosis Phylogeny Reconstruction Based on Combined Numerical Analysis with IS1081, IS6110, VNTR, and DR-Based Spoligotyping Suggests the Existence of Two New Phylogeographical Clades

This paper deals with phylogenetic relationships among a set of 90 clinical strains representative of the worldwide diversity of the Mycobacterium tuberculosis complex (Kremer et al. 1999) using eight independent genetic markers: IS6110, IS1081, the direct repeat (DR) locus, and five variable number of tandem DNA repeat loci (VNTR). In a preliminary experiment, phylogenetic trees based on single markers were constructed that led to the detection of some similarities between the VNTR-based and the spoligotyping-based phylogenetic trees. In the second step, a more global phenetic approach based on pairwise comparison of strains within each typing system was used, followed by calculations of mean genetic distances based on all the eight loci and the use of the neighbor-joining algorithm for tree reconstruction. This analysis confirmed our preliminary observations and suggested the existence of at least two new phylogeographical clades of M. tuberculosis, one defined as the ``East African–Indian family' (EA-I), which may find its origin on the African or Asian continents, and the other as the ``Latin American and Mediterranean' (LA-M) family. The existence of these two families was also validated by an independent phylogenetic analysis of spoligotyping on a larger set of shared types (n= 252) and further corroborated by VNTR and katG–gyrA results. The potential origin of these families of bacilli is discussed based on cattle domestication and human migration history. In conclusion, the information contained in insertion sequence and repetitive DNAs may serve as a model for the phylogenetic reconstruction of the M. tuberculosis complex.     

Introducing the Algerian Mitochondrial DNA and Y-Chromosome Profiles into the North African Landscape

North Africa is considered a distinct geographic and ethnic entity within Africa. Although modern humans originated in this Continent, studies of mitochondrial DNA (mtDNA) and Y-chromosome genealogical markers provide evidence that the North African gene pool has been shaped by the back-migration of several Eurasian lineages in Paleolithic and Neolithic times. More recent influences from sub-Saharan Africa and Mediterranean Europe are also evident. The presence of East-West and North-South haplogroup frequency gradients strongly reinforces the genetic complexity of this region. However, this genetic scenario is beset with a notable gap, which is the lack of consistent information for Algeria, the largest country in the Maghreb. To fill this gap, we analyzed a sample of 240 unrelated subjects from a northwest Algeria cosmopolitan population using mtDNA sequences and Y-chromosome biallelic polymorphisms, focusing on the fine dissection of haplogroups E and R, which are the most prevalent in North Africa and Europe respectively. The Eurasian component in Algeria reached 80% for mtDNA and 90% for Y-chromosome. However, within them, the North African genetic component for mtDNA (U6 and M1; 20%) is significantly smaller than the paternal (E-M81 and E-V65; 70%). The unexpected presence of the European-derived Y-chromosome lineages R-M412, R-S116, R-U152 and R-M529 in Algeria and the rest of the Maghreb could be the counterparts of the mtDNA H1, H3 and V subgroups, pointing to direct maritime contacts between the European and North African sides of the western Mediterranean. Female influx of sub-Saharan Africans into Algeria (20%) is also significantly greater than the male (10%). In spite of these sexual asymmetries, the Algerian uniparental profiles faithfully correlate between each other and with the geography.     

Analysis of the Trojan Y-Chromosome eradication strategy for an invasive species

The Trojan Y-Chromosome (TYC) strategy, an autocidal genetic biocontrol method, has been proposed to eliminate invasive alien species. In this work, we analyze the dynamical system model of the TYC strategy, with the aim of studying the viability of the TYC eradication and control strategy of an invasive species. In particular, because the constant introduction of sex-reversed trojan females for all time is not possible in practice, there arises the question: What happens if this injection is stopped after some time? Can the invasive species recover? To answer that question, we perform a rigorous bifurcation analysis and study the basin of attraction of the recovery state and the extinction state in both the full model and a certain reduced model. In particular, we find a theoretical condition for the eradication strategy to work. Additionally, the consideration of an Allee effect and the possibility of a Turing instability are also studied in this work. Our results show that: (1) with the inclusion of an Allee effect, the number of the invasive females is not required to be very low when the introduction of the sex-reversed trojan females is stopped, and the remaining Trojan Y-Chromosome population is sufficient to induce extinction of the invasive females; (2) incorporating diffusive spatial spread does not produce a Turing instability, which would have suggested that the TYC eradication strategy might be only partially effective, leaving a patchy distribution of the invasive species.     

流式细胞仪分离精子法的研究进展

以X精子和Y精子的DNA含量差别为基础的流式细胞仪分离精子技术是有效的哺乳动物性别控制方法。目前,流式细胞仪分离精子的速度与上世纪90年代初期相比已提高了30多倍,以3000个/秒的常规速度分离X精子或Y精子的准确率可达到85%～90%。迄今,用分离精子人工授精已产下了数万头的动物后代,而分离人的精子也开始应用到了避免X染色体连锁遗传疾病的生殖医学工程中。     

Identification and Phylogeny of the First T Cell Epitope Identified from a Human Gut Bacteroides Species

Host T cell reactivity toward gut bacterial epitopes has been recognized as part of disease pathogenesis. However, the specificity of T cells that recognize this vast number of epitopes has not yet been well described. After colonizing a C57BL/6J germ-free mouse with the human gut symbiotic bacteria Bacteroides thetaiotaomicron, we isolated a T cell that recognized these bacteria in vitro. Using this T cell, we mapped the first known non-carbohydrate T cell epitope within the phylum Bacteroidetes. The T cell also reacted to two other additional Bacteroides species. We identified the peptide that stimulated the T cell by using a genetic approach. Genomic data from the epitope-positive and epitope-negative bacteria explain the cross-reactivity of the T cell to multiple species. This epitope degeneracy should shape our understanding of the T cell repertoire stimulated by the complex microbiome residing in the gastrointestinal tract in both healthy and disease states.     

Phylogeny of the Dictyoptera

Abstract. Relationships among six key dictyopteran taxa (Mantodea; Blattodea (excluding Cryptocercidae); Cryptocercidae; Mastotermes darwiniensis , Termopsidae and Kalotermitidae [Isoptera]) are analysed based on seventy morphological, developmental and behavioural characters. The fossil record and the 'living fossil' genera Cryptocercus, Mastotermes and Archotermopsis are discussed in detail. Exact analysis of the character state matrix by implicit enumeration (Hennig86) resulted in one cladogram, strongly supporting Blattodea + Cryptocercidae as a sister group to Mantodea, with the Isoptera as a sister group to that complex. Arrangements within the termites are equivocal, with Termopsidae and Mastotermes darwiniensis possible as the relatively most primitive element of Isoptera.     

Phylogeny and Biogeography of the Stapeliads

The stapeliads (Apocynaceae — Asclepiadoideae — Ceropegieae) form a group of highly succulent, practically leafless plants with a very diverse floral morphology. They are found in the drier parts of the Old World from southern Africa eastwards to Burma. 55 characters are selected on which a cladistic analysis is based and these are discussed in detail where necessary. From this a phylogeny is produced for the stapeliads. It is shown that they are monophyletic and that they are nested within a larger group among the tribe Ceropegieae of the subfamily Asclepiadoideae. Several of the important innovations that have led to the diversification of the stapeliads in semidesertic regions are discussed. The stapeliads appear to be mainly fly-pollinated. The methods of attraction of flies are mentioned and it is suggested that there are four basic floral types among the stapeliads.     

Gene Pool Structure of Eastern Ukrainians as Inferred from the Y-Chromosome Haplogroups

Y chromosomes from representative sample of Eastern Ukrainians (94 individuals) were analyzed for composition and frequencies of haplogroups, defined by 11 biallelic loci located in non-recombining part of the chromosome (SRY1532, YAP, 92R7, DYF155S2, 12f2, Tat, M9, M17, M25,M89, andM56). In the Ukrainian gene, pool six haplogroups were revealed: E, F (including G and I), J, N3, P, and R1a1. These haplogroups were earlier detected in a study of Y-chromosome diversity on the territory of Europe as a whole. The major haplogroup in the Ukrainian gene pool, haplogroup R1a1 (earlier designated HG3), accounted for about 44% of all Y chromosomes in the sample examined. This haplogroup is thought to mark the migration patterns of the early Indo-Europeans and is associated with the distribution of the Kurgan archaeological culture. The second major haplogroup is haplogroup F (21.3%), which is a combination of the lineages differing by the time of appearance. Haplogroup P found with the frequency of 9.6%, represents the genetic contribution of the population originating from the ancient autochthonous population of Europe. Haplogroups J and E (11.7 and 4.2%, respectively) mark the migration patterns of the Middle-Eastern agriculturists during the Neolithic. The presence of the N3 lineage (9.6%) is likely explained by a contribution of the assimilated Finno–Ugric tribes. The data on the composition and frequencies of Y-chromosome haplogroups in the sample studied substantially supplement the existing picture of the male lineage distribution in the Eastern Slav population.     

Deep Phylogenetic Analysis of Haplogroup G1 Provides Estimates of SNP and STR Mutation Rates on the Human Y-Chromosome and Reveals Migrations of Iranic Speakers

Y-chromosomal haplogroup G1 is a minor component of the overall gene pool of South-West and Central Asia but reaches up to 80% frequency in some populations scattered within this area. We have genotyped the G1-defining marker M285 in 27 Eurasian populations (n= 5,346), analyzed 367 M285-positive samples using 17 Y-STRs, and sequenced ~11 Mb of the Y-chromosome in 20 of these samples to an average coverage of 67X. This allowed detailed phylogenetic reconstruction. We identified five branches, all with high geographical specificity: G1-L1323 in Kazakhs, the closely related G1-GG1 in Mongols, G1-GG265 in Armenians and its distant brother clade G1-GG162 in Bashkirs, and G1-GG362 in West Indians. The haplotype diversity, which decreased from West Iran to Central Asia, allows us to hypothesize that this rare haplogroup could have been carried by the expansion of Iranic speakers northwards to the Eurasian steppe and via founder effects became a predominant genetic component of some populations, including the Argyn tribe of the Kazakhs. The remarkable agreement between genetic and genealogical trees of Argyns allowed us to calibrate the molecular clock using a historical date (1405 AD) of the most recent common genealogical ancestor. The mutation rate for Y-chromosomal sequence data obtained was 0.78×10^-9 per bp per year, falling within the range of published rates. The mutation rate for Y-chromosomal STRs was 0.0022 per locus per generation, very close to the so-called genealogical rate. The “clan-based” approach to estimating the mutation rate provides a third, middle way between direct farther-to-son comparisons and using archeologically known migrations, whose dates are subject to revision and of uncertain relationship to genetic events.