False discovery rates of qpAdm-based screens for genetic admixture.

Although a broad range of methods exists for reconstructing population history from genome-wide single nucleotide polymorphism data, just a few methods gained popularity in archaeogenetics: principal component analysis (PCA); ADMIXTURE, an algorithm that models individuals as mixtures of multiple ancestral sources represented by actual or inferred populations; formal tests for admixture such as f3-statistics and D/f4-statistics; and qpAdm, a tool for fitting two-component and more complex admixture models to groups or individuals. Despite their popularity in archaeogenetics, which is explained by modest computational requirements and ability to analyze data of various types and qualities, protocols relying on qpAdm that screen numerous alternative models of varying complexity and find "fitting" models (often considering both estimated admixture proportions and p-values as a composite criterion of model fit) remain untested on complex simulated population histories in the form of admixture graphs of random topology. We analyzed genotype data extracted from such simulations and tested various types of high-throughput qpAdm protocols ("rotating" and "non-rotating", with or without temporal stratification of target groups and proxy ancestry sources, and with or without a "model competition" step). We caution that high-throughput qpAdm protocols may be inappropriate for exploratory analyses in poorly studied regions/periods since their false discovery rates varied between 12% and 68% depending on the details of the protocol and on the amount and quality of simulated data (i.e., >12% of fitting two-way admixture models imply gene flows that were not simulated). We demonstrate that for reducing false discovery rates of qpAdm protocols to nearly 0% it is advisable to use large SNP sets with low missing data rates, the rotating qpAdm protocol with a strictly enforced rule that target groups do not pre-date their proxy sources, and an unsupervised ADMIXTURE analysis as a way to verify feasible qpAdm models. Our study has a number of limitations: for instance, these recommendations depend on the assumption that the underlying genetic history is a complex admixture graph and not a stepping-stone model.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding.

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes.

BACKGROUND: Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role. RESULTS: We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication. CONCLUSIONS: Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

On the limits of fitting complex models of population history to f-statistics.

Our understanding of population history in deep time has been assisted by fitting admixture graphs (AGs) to data: models that specify the ordering of population splits and mixtures, which along with the amount of genetic drift and the proportions of mixture, is the only information needed to predict the patterns of allele frequency correlation among populations. The space of possible AGs relating populations is vast, and thus most published studies have identified fitting AGs through a manual process driven by prior hypotheses, leaving the majority of alternative models unexplored. Here, we develop a method for systematically searching the space of all AGs that can incorporate non-genetic information in the form of topology constraints. We implement this findGraphs tool within a software package, ADMIXTOOLS 2, which is a reimplementation of the ADMIXTOOLS software with new features and large performance gains. We apply this methodology to identify alternative models to AGs that played key roles in eight publications and find that in nearly all cases many alternative models fit nominally or significantly better than the published one. Our results suggest that strong claims about population history from AGs should only be made when all well-fitting and temporally plausible models share common topological features. Our re-evaluation of published data also provides insight into the population histories of humans, dogs, and horses, identifying features that are stable across the models we explored, as well as scenarios of populations relationships that differ in important ways from models that have been highlighted in the literature.

Water masses shape pico-nano eukaryotic communities of the Weddell Sea.

Polar oceans belong to the most productive and rapidly changing environments, yet our understanding of this fragile ecosystem remains limited. Here we present an analysis of a unique set of DNA metabarcoding samples from the western Weddell Sea sampled throughout the whole water column and across five water masses with different characteristics and different origin. We focus on factors affecting the distribution of planktonic pico-nano eukaryotes and observe an ecological succession of eukaryotic communities as the water masses move away from the surface and as oxygen becomes depleted with time. At the beginning of this succession, in the photic zone, algae, bacteriovores, and predators of small eukaryotes dominate the community, while another community develops as the water sinks deeper, mostly composed of parasitoids (syndinians), mesoplankton predators (radiolarians), and diplonemids. The strongly correlated distribution of syndinians and diplonemids along the depth and oxygen gradients suggests their close ecological link and moves us closer to understanding the biological role of the latter group in the ocean ecosystem.

Diplonemids - A Review on "New" Flagellates on the Oceanic Block.

Diplonemids are a group of flagellate protists, that belong to the phylum Euglenozoa alongside euglenids, symbiontids and kinetoplastids. They primarily inhabit marine environments, though are also found in freshwater lakes. Diplonemids have been considered as rare and unimportant eukaryotes for over a century, with only a handful of species described until recently. However, thanks to their unprecedented diversity and abundance in the world oceans, diplonemids now attract increased attention. Recent improvements in isolation and cultivation have enabled characterization of several new genera, warranting a re-examination of all available knowledge gathered so far. Here we summarize available data on diplonemids, focusing on the recent advances in the fields of diversity, ecology, genomics, metabolism, and endosymbionts. We illustrate the life stages of cultivated genera, and summarise all reported interspecies associations, which in turn suggest lifestyles of predation and parasitism. This review also includes the latest classification of diplonemids, with a taxonomic revision of the genus Diplonema. Ongoing efforts to sequence various diplonemids suggest the presence of large and complex genomes, which correlate with the metabolic versatility observed in the model species Paradiplonema papillatum. Finally, we highlight its successful transformation into one of few genetically tractable marine protists.

Indian genetic heritage in Southeast Asian populations.

The great ethnolinguistic diversity found today in mainland Southeast Asia (MSEA) reflects multiple migration waves of people in the past. Maritime trading between MSEA and India was established at the latest 300 BCE, and the formation of early states in Southeast Asia during the first millennium CE was strongly influenced by Indian culture, a cultural influence that is still prominent today. Several ancient Indian-influenced states were located in present-day Thailand, and various populations in the country are likely to be descendants of people from those states. To systematically explore Indian genetic heritage in MSEA populations, we generated genome-wide SNP data (using the Affymetrix Human Origins array) for 119 present-day individuals belonging to 10 ethnic groups from Thailand and co-analyzed them with published data using PCA, ADMIXTURE, and methods relying on f-statistics and on autosomal haplotypes. We found low levels of South Asian admixture in various MSEA populations for whom there is evidence of historical connections with the ancient Indian-influenced states but failed to find this genetic component in present-day hunter-gatherer groups and relatively isolated groups from the highlands of Northern Thailand. The results suggest that migration of Indian populations to MSEA may have been responsible for the spread of Indian culture in the region. Our results also support close genetic affinity between Kra-Dai-speaking (also known as Tai-Kadai) and Austronesian-speaking populations, which fits a linguistic hypothesis suggesting cladality of the two language families.

EukRef-excavates: seven curated SSU ribosomal RNA gene databases.

The small subunit ribosomal RNA (SSU rRNA) gene is a widely used molecular marker to study the diversity of life. Sequencing of SSU rRNA gene amplicons has become a standard approach for the investigation of the ecology and diversity of microbes. However, a well-curated database is necessary for correct classification of these data. While available for many groups of Bacteria and Archaea, such reference databases are absent for most eukaryotes. The primary goal of the EukRef project (eukref.org) is to close this gap and generate well-curated reference databases for major groups of eukaryotes, especially protists. Here we present a set of EukRef-curated databases for the excavate protists-a large assemblage that includes numerous taxa with divergent SSU rRNA gene sequences, which are prone to misclassification. We identified 6121 sequences, 625 of which were obtained from cultures, 3053 from cell isolations or enrichments and 2419 from environmental samples. We have corrected the classification for the majority of these curated sequences. The resulting publicly available databases will provide phylogenetically based standards for the improved identification of excavates in ecological and microbiome studies, as well as resources to classify new discoveries in excavate diversity.

Environmental determinants of the distribution of planktonic diplonemids and kinetoplastids in the oceans.

We analysed a widely used barcode, the V9 region of the 18S rRNA gene, to study the effect of environmental conditions on the distribution of two related heterotrophic protistan lineages in marine plankton, kinetoplastids and diplonemids. We relied on a major published dataset (Tara Oceans) where samples from the mesopelagic zone were available from just 32 of 123 locations, and both groups are most abundant in this zone. To close sampling gaps and obtain more information from the deeper ocean, we collected 57 new samples targeting especially the mesopelagic zone. We sampled in three geographic regions: the Arctic, two depth transects in the Adriatic Sea, and the anoxic Cariaco Basin. In agreement with previous studies, both protist groups are most abundant and diverse in the mesopelagic zone. In addition to that, we found that their abundance, richness, and community structure also depend on geography, oxygen concentration, salinity, temperature, and other environmental variables reflecting the abundance of algae and nutrients. Both groups studied here demonstrated similar patterns, although some differences were also observed. Kinetoplastids and diplonemids prefer tropical regions and nutrient-rich conditions and avoid high oxygen concentration, high salinity, and high density of algae.

Evolution of metabolic capabilities and molecular features of diplonemids, kinetoplastids, and euglenids.

BACKGROUND: The Euglenozoa are a protist group with an especially rich history of evolutionary diversity. They include diplonemids, representing arguably the most species-rich clade of marine planktonic eukaryotes; trypanosomatids, which are notorious parasites of medical and veterinary importance; and free-living euglenids. These different lifestyles, and particularly the transition from free-living to parasitic, likely require different metabolic capabilities. We carried out a comparative genomic analysis across euglenozoan diversity to see how changing repertoires of enzymes and structural features correspond to major changes in lifestyles. RESULTS: We find a gradual loss of genes encoding enzymes in the evolution of kinetoplastids, rather than a sudden decrease in metabolic capabilities corresponding to the origin of parasitism, while diplonemids and euglenids maintain more metabolic versatility. Distinctive characteristics of molecular machines such as kinetochores and the pre-replication complex that were previously considered specific to parasitic kinetoplastids were also identified in their free-living relatives. Therefore, we argue that they represent an ancestral rather than a derived state, as thought until the present. We also found evidence of ancient redundancy in systems such as NADPH-dependent thiol-redox. Only the genus Euglena possesses the combination of trypanothione-, glutathione-, and thioredoxin-based systems supposedly present in the euglenozoan common ancestor, while other representatives of the phylum have lost one or two of these systems. Lastly, we identified convergent losses of specific metabolic capabilities between free-living kinetoplastids and ciliates. Although this observation requires further examination, it suggests that certain eukaryotic lineages are predisposed to such convergent losses of key enzymes or whole pathways. CONCLUSIONS: The loss of metabolic capabilities might not be associated with the switch to parasitic lifestyle in kinetoplastids, and the presence of a highly divergent (or unconventional) kinetochore machinery might not be restricted to this protist group. The data derived from the transcriptomes of free-living early branching prokinetoplastids suggests that the pre-replication complex of Trypanosomatidae is a highly divergent version of the conventional machinery. Our findings shed light on trends in the evolution of metabolism in protists in general and open multiple avenues for future research.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America.

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.

Neobodonids are dominant kinetoplastids in the global ocean.

Kinetoplastid flagellates comprise basal mostly free-living bodonids and derived obligatory parasitic trypanosomatids, which belong to the best-studied protists. Due to their omnipresence in aquatic environments and soil, the bodonids are of ecological significance. Here, we present the first global survey of marine kinetoplastids and compare it with the strikingly different patterns of abundance and diversity in their sister clade, the diplonemids. Based on analysis of 18S rDNA V9 ribotypes obtained from 124 sites sampled during the Tara Oceans expedition, our results show generally low to moderate abundance and diversity of planktonic kinetoplastids. Although we have identified all major kinetoplastid lineages, 98% of kinetoplastid reads are represented by neobodonids, namely specimens of the Neobodo and Rhynchomonas genera, which make up 59% and 18% of all reads, respectively. Most kinetoplastids have small cell size (0.8-5 µm) and tend to be more abundant in the mesopelagic as compared to the euphotic zone. Some of the most abundant operational taxonomic units have distinct geographical distributions, and three novel putatively parasitic neobodonids were identified, along with their potential hosts.

Extreme Diversity of Diplonemid Eukaryotes in the Ocean.

The world's oceans represent by far the largest biome, with great importance for the global ecosystem [1-4]. The vast majority of ocean biomass and biodiversity is composed of microscopic plankton. Recent results from the Tara Oceans metabarcoding study revealed that a significant part of the plankton in the upper sunlit layer of the ocean is represented by an understudied group of heterotrophic excavate flagellates called diplonemids [5, 6]. We have analyzed the diversity and distribution patterns of diplonemid populations on the extended set of Tara Oceans V9 18S rDNA metabarcodes amplified from 850 size- fractionated plankton communities sampled across 123 globally distributed locations, for the first time also including samples from the mesopelagic zone, which spans the depth from about 200 to 1,000 meters. Diplonemids separate into four major clades, with the vast majority falling into the deep-sea pelagic diplonemid clade. Remarkably, diversity of this clade inferred from metabarcoding data surpasses even that of dinoflagellates, metazoans, and rhizarians, qualifying diplonemids as possibly the most diverse group of marine planktonic eukaryotes. Diplonemids display strong vertical separation between the photic and mesopelagic layers, with the majority of their relative abundance and diversity occurring in deeper waters. Globally, diplonemids display no apparent biogeographic structuring, with a few hyperabundant cosmopolitan operational taxonomic units (OTUs) dominating their communities. Our results suggest that the planktonic diplonemids are among the key heterotrophic players in the largest ecosystem of our biosphere, yet their roles in this ecosystem remain unknown.

Genome of Leptomonas pyrrhocoris: a high-quality reference for monoxenous trypanosomatids and new insights into evolution of Leishmania.

Many high-quality genomes are available for dixenous (two hosts) trypanosomatid species of the genera Trypanosoma, Leishmania, and Phytomonas, but only fragmentary information is available for monoxenous (single-host) trypanosomatids. In trypanosomatids, monoxeny is ancestral to dixeny, thus it is anticipated that the genome sequences of the key monoxenous parasites will be instrumental for both understanding the origin of parasitism and the evolution of dixeny. Here, we present a high-quality genome for Leptomonas pyrrhocoris, which is closely related to the dixenous genus Leishmania. The L. pyrrhocoris genome (30.4 Mbp in 60 scaffolds) encodes 10,148 genes. Using the L. pyrrhocoris genome, we pinpointed genes gained in Leishmania. Among those genes, 20 genes with unknown function had expression patterns in the Leishmania mexicana life cycle suggesting their involvement in virulence. By combining differential expression data for L. mexicana, L. major and Leptomonas seymouri, we have identified several additional proteins potentially involved in virulence, including SpoU methylase and U3 small nucleolar ribonucleoprotein IMP3. The population genetics of L. pyrrhocoris was also addressed by sequencing thirteen strains of different geographic origin, allowing the identification of 1,318 genes under positive selection. This set of genes was significantly enriched in components of the cytoskeleton and the flagellum.

Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry.

The Kets, an ethnic group in the Yenisei River basin, Russia, are considered the last nomadic hunter-gatherers of Siberia, and Ket language has no transparent affiliation with any language family. We investigated connections between the Kets and Siberian and North American populations, with emphasis on the Mal'ta and Paleo-Eskimo ancient genomes, using original data from 46 unrelated samples of Kets and 42 samples of their neighboring ethnic groups (Uralic-speaking Nganasans, Enets, and Selkups). We genotyped over 130,000 autosomal SNPs, identified mitochondrial and Y-chromosomal haplogroups, and performed high-coverage genome sequencing of two Ket individuals. We established that Nganasans, Kets, Selkups, and Yukaghirs form a cluster of populations most closely related to Paleo-Eskimos in Siberia (not considering indigenous populations of Chukotka and Kamchatka). Kets are closely related to modern Selkups and to some Bronze and Iron Age populations of the Altai region, with all these groups sharing a high degree of Mal'ta ancestry. Implications of these findings for the linguistic hypothesis uniting Ket and Na-Dene languages into a language macrofamily are discussed.

Diplonemids.

Lukes et al. introduce an enigmatic group of unicellular eukaryotes called the diplonemids, which according to recent surveys may be widespread in marine ecosystems.

Ocean plankton. Eukaryotic plankton diversity in the sunlit ocean.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

Lack of evidence for integration of Trypanosoma cruzi minicircle DNA in South American human genomes.

Horizontal gene transfer involving kinetoplast DNA minicircles between Trypanosoma cruzi and its mammalian hosts has recently been proposed as a usual consequence of infection (Hecht et al., 2010). However, we have found no sequences longer than 29 bp perfectly matching minicircles of T. cruzi in the unassembled reads from Colombian and Peruvian human populations provided by the 1,000 Genome project (129 individuals in total, coverage from 1.4× to 36.3×, read length from 42 to 101 bp). The weak sequence matches that were identified are shared with a Finnish population used as a control from a non-endemic area.

Specificity of genetic diversity in D1S80 revealed by SNP-VNTR haplotyping.

BACKGROUND: The allele frequency patterns of the D1S80 variable number tandem repeat (VNTR) locus have been shown to be multimodal in many different human populations. AIM: To explore the complex allele distribution of the D1S80 polymorphic locus in different populations comparing the derived single nucleotide polymorphism (SNP) rs16824398-D1S80 haplotype frequencies in samples of European (Russians), Asian (Yakuts) and sub-Saharan African origin. SUBJECTS AND METHODS: The D1S80 locus together with its 5'-flanking region including SNP rs16824398 was amplified using allele-specific polymerase chain reaction (PCR). RESULTS: Haplotype phase determination sub-divided the total D1S80 allele spectrum into two allele sets marked by the corresponding SNP rs16824398 alleles. In non-African samples, the most frequent D1S80 alleles had 24 and 18 repeats that were associated with different SNP backgrounds (T and G alleles, respectively). Both combinations also occurred in Africans, but these samples exhibited an expanded spectrum of VNTR alleles on both SNP backgrounds. CONCLUSIONS: The sub-division of the D1S80 allele spectrum shape on the linked SNP background is indicative of populations of the main human groups. The reported differences in D1S80 allele spectra between populations of different ethnic origins can be explained by the ratios of chromosomes with T and G alleles.

Haplotype frequencies at the DRD2 locus in populations of the East European Plain.

BACKGROUND: It was demonstrated previously that the three-locus RFLP haplotype, TaqI B-TaqI D-TaqI A (B-D-A), at the DRD2 locus constitutes a powerful genetic marker and probably reflects the most ancient dispersal of anatomically modern humans. RESULTS: We investigated TaqI B, BclI, MboI, TaqI D, and TaqI A RFLPs in 17 contemporary populations of the East European Plain and Siberia. Most of these populations belong to the Indo-European or Uralic language families. We identified three common haplotypes, which occurred in more than 90% of chromosomes investigated. The frequencies of the haplotypes differed according to linguistic and geographical affiliation. CONCLUSION: Populations in the northwestern (Byelorussians from Mjadel'), northern (Russians from Mezen' and Oshevensk), and eastern (Russians from Puchezh) parts of the East European Plain had relatively high frequencies of haplotype B2-D2-A2, which may reflect admixture with Uralic-speaking populations that inhabited all of these regions in the Early Middle Ages.