Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration

Background

The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history.

Results

We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3–4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi.

Conclusions

This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

     

Sonication of mouse sperm membranes reveals distinct protein domains.

Molecular interactions between sperm and zona pellucida (ZP) during mammalian fertilization are not well characterized. To begin to characterize sperm components that are involved in sperm-ZP interactions, we isolated and density fractionated sperm membranes. The membrane fractions recovered from a density fractionation protocol were characterized, and sonication was compared with vortexing for preparation of sperm membranes by examining the distribution of proteins in the membrane fractions obtained from these 2 protocols. Biochemical and microscopic analyses were used to determine the composition of the sonicated membrane fractions, and immunoblotting was used to identify fractions containing some of the previously suggested ZP3 receptors. Transmission electron microscopy revealed that bands 1-3 contained membrane vesicles and band 4 contained axonemal and midpiece fragments. SDS-PAGE revealed that bands 1 and 2 shared many proteins, but band 3 contained a number of unique proteins. Surface labeling with 125I demonstrated that bands 1 and 2 contained the majority of the sperm surface protein markers, whereas band 3 contained minor amounts of surface markers. Lectin-binding characteristics of sperm membrane glycoproteins were used to compare the relative distribution of glycosylated proteins in vortexed or sonicated membrane preparations. These characterizations indicate that sonication enhanced the differential distribution of sperm membrane proteins among the density fractions and suggests that this method is preferable for preparation of membrane fractions to be used for identification of proteins that mediate sperm-egg interactions.     

Y-chromosome evidence of southern origin of the East Asian-specific haplogroup O3-M122

The prehistoric peopling of East Asia by modern humans remains controversial with respect to early population migrations. Here, we present a systematic sampling and genetic screening of an East Asian-specific Y-chromosome haplogroup (O3-M122) in 2,332 individuals from diverse East Asian populations. Our results indicate that the O3-M122 lineage is dominant in East Asian populations, with an average frequency of 44.3%. The microsatellite data show that the O3-M122 haplotypes in southern East Asia are more diverse than those in northern East Asia, suggesting a southern origin of the O3-M122 mutation. It was estimated that the early northward migration of the O3-M122 lineages in East Asia occurred approximately 25,000-30,000 years ago, consistent with the fossil records of modern humans in East Asia.     

Phylogeography of lions (Panthera leo ssp.) reveals three distinct taxa and a late Pleistocene reduction in genetic diversity

Lions were the most widespread carnivores in the late Pleistocene, ranging from southern Africa to the southern USA, but little is known about the evolutionary relationships among these Pleistocene populations or the dynamics that led to their extinction. Using ancient DNA techniques, we obtained mitochondrial sequences from 52 individuals sampled across the present and former range of lions. Phylogenetic analysis revealed three distinct clusters: (i) modern lions, Panthera leo ; (ii) extinct Pleistocene cave lions, which formed a homogeneous population extending from Europe across Beringia (Siberia, Alaska and western Canada); and (iii) extinct American lions, which formed a separate population south of the Pleistocene ice sheets. The American lion appears to have become genetically isolated around 340 000 years ago, despite the apparent lack of significant barriers to gene flow with Beringian populations through much of the late Pleistocene. We found potential evidence of a severe population bottleneck in the cave lion during the previous interstadial, sometime after 48 000 years, adding to evidence from bison, mammoths, horses and brown bears that megafaunal populations underwent major genetic alterations throughout the last interstadial, potentially presaging the processes involved in the subsequent end-Pleistocene mass extinctions.     

Site-directed mutagenesis of the yeast PRP20/1SRM1 gene reveals distinct activity domains in the protein product

Prp20/Srm1, a homolog of the mammalian protein RCC1 in Saccharomyces cerevisiae, binds to double-stranded DNA (dsDNA) through a multicomponent complex in vitro. This dsDNA-binding capability of the Prp20 complex has been shown to be cell-cycle dependent; affinity for dsDNA is lost during DNA replication. By analyzing a number of temperature sensitive (ts) prp20 alleles produced in vivo and in vitro, as well as site-directed mutations in highly conserved positions in the imperfect repeats that make up the protein, we have determined a relationship between the residues at these positions, cell viability, and the dsDNA-binding abilities of the Prp20 complex. These data reveal that the essential residues for Prp20 function are located mainly in the second and the third repeats at the amino-terminus and the last two repeats, the seventh and eighth, at the carboxyl-terminus of Prp20. Carboxyl-terminal mutations in Prp20 differ from amino-terminal mutations in showing loss of dsDNA binding: their conditional lethal phenotype and the loss of dsDNA binding affinity are both suppressible by overproduction of Gsp1, a GTP-binding constituent of the Prp20 complex, homologous to the mammalian protein TC4/Ran. Although wild-type Prp20 does not bind to dsDNA on its own, two mutations in conserved residues were found that caused the isolated protein to bind dsDNA. These data imply that, in situ, the other components of the Prp20 complex regulate the conformation of Prp20 and thus its affinity for dsDNA. Gsp1 not only influences the dsDNA-binding ability of Prp20 but it also regulates other essential function(s) of the Prp20 complex. Overproduction of Gsp1 also suppresses the lethality of two conditional mutations in the penultimate carboxyl-terminal repeat of Prp20, even though these mutations do not eliminate the dsDNA binding activity of the Prp20 complex. Other site-directed mutants reveal that internal and carboxyl-terminal regions of Prp20 that lack homology to RCC1 are dispensable for dsDNA binding and growth.     

Analysis of homologous gene clusters in Caenorhabditis elegans reveals striking regional cluster domains

An algorithm for detecting local clusters of homologous genes was applied to the genome of Caenorhabditis elegans. Clusters of two or more homologous genes are abundant, totaling 1391 clusters containing 4607 genes, over one-fifth of all genes in C. elegans. Cluster genes are distributed unevenly in the genome, with the large majority located on autosomal chromosome arms, regions characterized by higher genetic recombination and more repeat sequences than autosomal centers and the X chromosome. Cluster genes are transcribed at much lower levels than average and very few have gross phenotypes as assayed by RNAi-mediated reduction of function. The molecular identity of cluster genes is unusual, with a preponderance of nematode-specific gene families that encode putative secreted and transmembrane proteins, and enrichment for genes implicated in xenobiotic detoxification and innate immunity. Gene clustering in Drosophila melanogaster is also substantial and the molecular identity of clustered genes follows a similar pattern. I hypothesize that autosomal chromosome arms in C. elegans undergo frequent local gene duplication and that these duplications support gene diversification and rapid evolution in response to environmental challenges. Although specific gene clusters have been documented in C. elegans, their abundance, genomic distribution, and unusual molecular identities were previously unrecognized.     

Genetic structure and gene flow in Gran Chaco populations of Argentina: evidence from Y-chromosome markers

The Gran Chaco region of central South America has been settled by humans for only the last 4,000-5,000 years. To investigate population structure and variation in this region's indigenous population, we scored males from tribes of the Argentinean part of the Gran Chaco (Pilagá, Wichí, and Toba, representing two major language groups, the Mataco and Guaycurú) for a number of Y-chromosome polymorphisms. The markers included eight microsatellites (DYS19, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, and DYS439) and the unique native American single nucleotide polymorphism, DYS199. Sixty males (77%) from the total sample carried the DYS199T chromosome, and these were the focus of the present analysis. Unlike most other native Americans, Gran Chaco males show a moderate level of diversity at the DYS19 locus but still less than that seen in non-native Americans. The FST value for Y-chromosome markers in Gran Chaco was 0.107, a value that is more than double that found for mtDNA haplogroups in the same tribes but is not particularly high compared with other Y-chromosome studies. Phylogenetic trees based on all eight microsatellites showed relatively poor correlation of the tribes with either geography or language, and this is possibly explained by their ecology. They are seasonal hunters living in small bands, and under such circumstances drift will be a powerful evolutionary force. An UPGMA tree based on five microsatellites (DYS19, DYS390, DYS391, DYS392, and DYS393), however, showed a more positive relationship, suggesting that DYS437, DYS438, and DYS439 may behave differently from the other microsatellites. No association was found between maternal and paternal lineage distributions. The time to the most recent common ancestor of the DYS199T chromosome is calculated to lie between 13,000 and 26,000 years. This range is consistent with estimates from other Y-chromosome studies as well as with estimates from mtDNA and the archeology of the colonization of South America. We conclude that the male lineages present in the contemporary Gran Chaco population reflect the level of diversity found in South America and that the region's male founders did not carry a restricted gene pool.     

Phylogeography of the human mitochondrial L1c haplogroup: genetic signatures of the prehistory of Central Africa

Interindividual variation of human mitochondrial DNA has been extensively studied over the last two decades, and its usefulness for reconstructing evolutionary relationships of extant populations has been proved. However, some mitochondrial lineages still need to be studied using a combination of larger and tailored datasets and increased level of resolution in order to shed light on their origin and on the processes underlying their present distribution. In this study, we analyze the phylogeny of the L1c haplogroup of human mitochondrial DNA using sequence data from hypervariable regions 1 and 2 obtained from 455 individuals (extracted from a total sampling of 2542 individuals) belonging to sub-Saharan African and African-American populations. We propose a substantial revision of L1c phylogeny, by introducing one new sub-haplogroup (L1c4), two new L1c1 clades (L1c1b and L1c1c), and by reassigning the previous L1c1a1 sequences to a clade which we termed L1c5. The new phylogeny encompasses distinct lineages with different evolutionary histories. In fact, based on population frequency, internal variation and mismatch distribution, we propose that L1c1b, L1c1c and L1c2 originated in Bantu ancestors, whereas L1c1a, L1c4 and L1c5 evolved among Western Pygmies. The population structure of L1c is not comparable to any known mitochondrial or, even, Y-chromosomal haplogroup, and challenges the current view that most of mtDNA variation in Pygmies might reflect admixture with Bantu or a persistence of plesiomorphic characters. In fact, the unique feature of the L1c is that it retains a signature of a phase common to the ancestors of the Bantu and Western Pygmies, while encompassing some specific sub-clades which can indicate their divergence. This allowed us to attempt a phylogenetically based assessment of the evolutionary relationships between the two groups. Taking into consideration estimates of the time to the most recent common ancestor of L1c and its clades together with archaeological and paleoclimatological evidence, we propose that the ancestors of Bantu and Western Pygmies separated between 60 and 30 kya.     

Y-chromosome variation in Altaian Kazakhs reveals a common paternal gene pool for Kazakhs and the influence of Mongolian expansions

Kazakh populations have traditionally lived as nomadic pastoralists that seasonally migrate across the steppe and surrounding mountain ranges in Kazakhstan and southern Siberia. To clarify their population history from a paternal perspective, we analyzed the non-recombining portion of the Y-chromosome from Kazakh populations living in southern Altai Republic, Russia, using a high-resolution analysis of 60 biallelic markers and 17 STRs. We noted distinct differences in the patterns of genetic variation between maternal and paternal genetic systems in the Altaian Kazakhs. While they possess a variety of East and West Eurasian mtDNA haplogroups, only three East Eurasian paternal haplogroups appear at significant frequencies (C3*, C3c and O3a3c*). In addition, the Y-STR data revealed low genetic diversity within these lineages. Analysis of the combined biallelic and STR data also demonstrated genetic differences among Kazakh populations from across Central Asia. The observed differences between Altaian Kazakhs and indigenous Kazakhs were not the result of admixture between Altaian Kazakhs and indigenous Altaians. Overall, the shared paternal ancestry of Kazakhs differentiates them from other Central Asian populations. In addition, all of them showed evidence of genetic influence by the 13(th) century CE Mongol Empire. Ultimately, the social and cultural traditions of the Kazakhs shaped their current pattern of genetic variation.     

Structure-function analysis of Escherichia coli DNA helicase I reveals non-overlapping transesterase and helicase domains

TraI (DNA helicase I) is an Escherichia coli F plasmid-encoded protein required for bacterial conjugative DNA transfer. The protein is a sequence-specific DNA transesterase that provides the site- and strand-specific nick required to initiate DNA strand transfer and a 5' to 3' DNA helicase that unwinds the F plasmid to provide the single-stranded DNA that is transferred from donor to recipient. Sequence comparisons with other transesterases and helicases suggest that these activities reside in the N- and C-terminal regions of TraI, respectively. Computer-assisted secondary structure probability analysis identified a potential interdomain region spanning residues 304-309. Proteins encoded by segments of traI, whose N or C terminus either flanked or coincided with this region, were purified and assessed for catalytic activity. Amino acids 1-306 contain the transesterase activity, whereas amino acids 309-1504 contain the helicase activity. The C-terminal 252 amino acids of the 1756-amino acid TraI protein are not required for either helicase or transesterase activity. Protein and nucleic acid sequence similarity searches indicate that the occurrence of both transesterase- and helicase-associated motifs in a conjugative DNA transfer initiator protein is rare. Only two examples (other than R100 plasmid TraI) were found: R388 plasmid TrwC and R46 plasmid (pKM101) TraH, belonging to the IncW and IncN groups of broad host range conjugative plasmids, respectively. The most significant structural difference between these proteins and TraI is that TraI contains an additional region of approximately 650 residues between the transesterase domain and the helicase-associated motifs. This region is required for helicase activity.     

Two distinct domains in Drosophila melanogaster telomeres

Telomeres are generally considered heterochromatic. On the basis of DNA composition, the telomeric region of Drosophila melanogaster contains two distinct subdomains: a subtelomeric region of repetitive DNA, termed TAS, and a terminal array of retrotransposons, which perform the elongation function instead of telomerase. We have identified several P-element insertions into this retrotransposon array and compared expression levels of transgenes with similar integrations into TAS and euchromatic regions. In contrast to insertions in TAS, which are silenced, reporter genes in the terminal HeT-A, TAHRE, or TART retroelements did not exhibit repressed expression in comparison with the same transgene construct in euchromatin. These data, in combination with cytological studies, provide evidence that the subtelomeric TAS region exhibits features resembling heterochromatin, while the terminal retrotransposon array exhibits euchromatic characteristics.     

Draft genome sequence of Xanthomonas fragariae reveals reductive evolution and distinct virulence-related gene content

Background

Xanthomonas fragariae (Xf) is a bacterial strawberry pathogen and an A2 quarantine organism on strawberry planting stock in the EU. It is taxonomically and metabolically distinct within the genus Xanthomonas, and known for its host specificity. As part of a broader pathogenicity study, the genome of a Belgian, virulent Xf strain (LMG 25863) was assembled to draft status and examined for its pathogenicity related gene content.

Results

The Xf draft genome (4.2 Mb) was considerably smaller than most known Xanthomonas genomes (~5 Mb). Only half of the genes coding for TonB-dependent transporters and cell-wall degrading enzymes that are typically present in other Xanthomonas genomes, were found in Xf. Other missing genes/regions with a possible impact on its plant-host interaction were: i) the three loci for xylan degradation and metabolism, ii) a locus coding for a ß-ketoadipate phenolics catabolism pathway, iii) xcs, one of two Type II Secretion System coding regions in Xanthomonas, and iv) the genes coding for the glyoxylate shunt pathway. Conversely, the Xf genome revealed a high content of externally derived DNA and several uncommon, possibly virulence-related features: a Type VI Secretion System, a second Type IV Secretion System and a distinct Type III Secretion System effector repertoire comprised of multiple rare effectors and several putative new ones.

Conclusions

The draft genome sequence of LMG 25863 confirms the distinct phylogenetic position of Xf within the genus Xanthomonas and reveals a patchwork of both lost and newly acquired genomic features. These features may help explain the specific, mostly endophytic association of Xf with the strawberry plant.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-829) contains supplementary material, which is available to authorized users.