Reconstruction of Microbial Haplotypes by Integration of Statistical and Physical Linkage in Scaffolding.

DNA sequencing technologies provide unprecedented opportunities to analyze within-host evolution of microorganism populations. Often, within-host populations are analyzed via pooled sequencing of the population, which contains multiple individuals or "haplotypes." However, current next-generation sequencing instruments, in conjunction with single-molecule barcoded linked-reads, cannot distinguish long haplotypes directly. Computational reconstruction of haplotypes from pooled sequencing has been attempted in virology, bacterial genomics, metagenomics, and human genetics, using algorithms based on either cross-host genetic sharing or within-host genomic reads. Here, we describe PoolHapX, a flexible computational approach that integrates information from both genetic sharing and genomic sequencing. We demonstrated that PoolHapX outperforms state-of-the-art tools tailored to specific organismal systems, and is robust to within-host evolution. Importantly, together with barcoded linked-reads, PoolHapX can infer whole-chromosome-scale haplotypes from 50 pools each containing 12 different haplotypes. By analyzing real data, we uncovered dynamic variations in the evolutionary processes of within-patient HIV populations previously unobserved in single position-based analysis.

Correction: Highly diverged novel subunit composition of apicomplexan F-type ATP synthase identified from Toxoplasma gondii.

[This corrects the article DOI: 10.1371/journal.pbio.2006128.].

The genome of the zoonotic malaria parasite Plasmodium simium reveals adaptations to host switching.

BACKGROUND: Plasmodium simium, a malaria parasite of non-human primates (NHP), was recently shown to cause zoonotic infections in humans in Brazil. We sequenced the P. simium genome to investigate its evolutionary history and to identify any genetic adaptions that may underlie the ability of this parasite to switch between host species. RESULTS: Phylogenetic analyses based on whole genome sequences of P. simium from humans and NHPs reveals that P. simium is monophyletic within the broader diversity of South American Plasmodium vivax, suggesting P. simium first infected NHPs as a result of a host switch of P. vivax from humans. The P. simium isolates show the closest relationship to Mexican P. vivax isolates. Analysis of erythrocyte invasion genes reveals differences between P. vivax and P. simium, including large deletions in the Duffy-binding protein 1 (DBP1) and reticulocyte-binding protein 2a genes of P. simium. Analysis of P. simium isolated from NHPs and humans revealed a deletion of 38 amino acids in DBP1 present in all human-derived isolates, whereas NHP isolates were multi-allelic. CONCLUSIONS: Analysis of the P. simium genome confirmed a close phylogenetic relationship between P. simium and P. vivax, and suggests a very recent American origin for P. simium. The presence of the DBP1 deletion in all human-derived isolates tested suggests that this deletion, in combination with other genetic changes in P. simium, may facilitate the invasion of human red blood cells and may explain, at least in part, the basis of the recent zoonotic infections.

Host-directed editing of the SARS-CoV-2 genome.

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.

Plasmodium kinesin-8X associates with mitotic spindles and is essential for oocyst development during parasite proliferation and transmission.

Kinesin-8 proteins are microtubule motors that are often involved in regulation of mitotic spindle length and chromosome alignment. They move towards the plus ends of spindle microtubules and regulate the dynamics of these ends due, at least in some species, to their microtubule depolymerization activity. Plasmodium spp. exhibit an atypical endomitotic cell division in which chromosome condensation and spindle dynamics in the different proliferative stages are not well understood. Genome-wide shared orthology analysis of Plasmodium spp. revealed the presence of two kinesin-8 motor proteins, kinesin-8X and kinesin-8B. Here we studied the biochemical properties of kinesin-8X and its role in parasite proliferation. In vitro, kinesin-8X has motility and depolymerization activities like other kinesin-8 motors. To understand the role of Plasmodium kinesin-8X in cell division, we used fluorescence-tagging and live cell imaging to define its location, and gene targeting to analyse its function, during all proliferative stages of the rodent malaria parasite P. berghei life cycle. The results revealed a spatio-temporal involvement of kinesin-8X in spindle dynamics and an association with both mitotic and meiotic spindles and the putative microtubule organising centre (MTOC). Deletion of the kinesin-8X gene revealed a defect in oocyst development, confirmed by ultrastructural studies, suggesting that this protein is required for oocyst development and sporogony. Transcriptome analysis of Δkinesin-8X gametocytes revealed modulated expression of genes involved mainly in microtubule-based processes, chromosome organisation and the regulation of gene expression, supporting a role for kinesin-8X in cell division. Kinesin-8X is thus required for parasite proliferation within the mosquito and for transmission to the vertebrate host.

Highly diverged novel subunit composition of apicomplexan F-type ATP synthase identified from Toxoplasma gondii.

The mitochondrial F-type ATP synthase, a multisubunit nanomotor, is critical for maintaining cellular ATP levels. In T. gondii and other apicomplexan parasites, many subunit components necessary for proper assembly and functioning of this enzyme appear to be missing. Here, we report the identification of 20 novel subunits of T. gondii F-type ATP synthase from mass spectrometry analysis of partially purified monomeric (approximately 600 kDa) and dimeric (>1 MDa) forms of the enzyme. Despite extreme sequence diversification, key FO subunits a, b, and d can be identified from conserved structural features. Orthologs for these proteins are restricted to apicomplexan, chromerid, and dinoflagellate species. Interestingly, their absence in ciliates indicates a major diversion, with respect to subunit composition of this enzyme, within the alveolate clade. Discovery of these highly diversified novel components of the apicomplexan F-type ATP synthase complex could facilitate the development of novel antiparasitic agents. Structural and functional characterization of this unusual enzyme complex will advance our fundamental understanding of energy metabolism in apicomplexan species.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

High-Throughput Sequencing-Based Investigation of Viruses in Human Cancers by Multienrichment Approach.

BACKGROUND: Viruses and other infectious agents cause more than 15% of human cancer cases. High-throughput sequencing-based studies of virus-cancer associations have mainly focused on cancer transcriptome data. METHODS: In this study, we applied a diverse selection of presequencing enrichment methods targeting all major viral groups, to characterize the viruses present in 197 samples from 18 sample types of cancerous origin. Using high-throughput sequencing, we generated 710 datasets constituting 57 billion sequencing reads. RESULTS: Detailed in silico investigation of the viral content, including exclusion of viral artefacts, from de novo assembled contigs and individual sequencing reads yielded a map of the viruses detected. Our data reveal a virome dominated by papillomaviruses, anelloviruses, herpesviruses, and parvoviruses. More than half of the included samples contained 1 or more viruses; however, no link between specific viruses and cancer types were found. CONCLUSIONS: Our study sheds light on viral presence in cancers and provides highly relevant virome data for future reference.

Substitutions of short heterologous DNA segments of intragenomic or extragenomic origins produce clustered genomic polymorphisms.

In a screen for unexplained mutation events we identified a previously unrecognized mechanism generating clustered DNA polymorphisms such as microindels and cumulative SNPs. The mechanism, short-patch double illegitimate recombination (SPDIR), facilitates short single-stranded DNA molecules to invade and replace genomic DNA through two joint illegitimate recombination events. SPDIR is controlled by key components of the cellular genome maintenance machinery in the gram-negative bacterium Acinetobacter baylyi. The source DNA is primarily intragenomic but can also be acquired through horizontal gene transfer. The DNA replacements are nonreciprocal and locus independent. Bioinformatic approaches reveal occurrence of SPDIR events in the gram-positive human pathogen Streptococcus pneumoniae and in the human genome.

Correction to: Cancer associated fibroblasts (CAFs) are activated in cutaneous basal cell carcinoma and in the peritumoural skin.

After publication of the original article [1] it was identified that order of the author list had been presented incorrectly. The author Robert Gniadecki's surname was also incorrect in the original article.

Cutavirus in Cutaneous Malignant Melanoma.

A novel human protoparvovirus related to human bufavirus and preliminarily named cutavirus has been discovered. We detected cutavirus in a sample of cutaneous malignant melanoma by using viral enrichment and high-throughput sequencing. The role of cutaviruses in cutaneous cancers remains to be investigated.

Genomic analysis of the causative agents of coccidiosis in domestic chickens.

Global production of chickens has trebled in the past two decades and they are now the most important source of dietary animal protein worldwide. Chickens are subject to many infectious diseases that reduce their performance and productivity. Coccidiosis, caused by apicomplexan protozoa of the genus Eimeria, is one of the most important poultry diseases. Understanding the biology of Eimeria parasites underpins development of new drugs and vaccines needed to improve global food security. We have produced annotated genome sequences of all seven species of Eimeria that infect domestic chickens, which reveal the full extent of previously described repeat-rich and repeat-poor regions and show that these parasites possess the most repeat-rich proteomes ever described. Furthermore, while no other apicomplexan has been found to possess retrotransposons, Eimeria is home to a family of chromoviruses. Analysis of Eimeria genes involved in basic biology and host-parasite interaction highlights adaptations to a relatively simple developmental life cycle and a complex array of co-expressed surface proteins involved in host cell binding.

Cancer associated fibroblasts (CAFs) are activated in cutaneous basal cell carcinoma and in the peritumoural skin.

BACKGROUND: Cutaneous basal cell carcinoma (BCC) is the commonest cancer worldwide. BCC is locally invasive and the surrounding stromal microenvironment is pivotal for tumourigenesis. Cancer associated fibroblasts (CAFs) in the microenvironment are essential for tumour growth in a variety of neoplasms but their role in BCC is poorly understood. METHODS: Material included facial BCC and control skin from the peritumoural area and from the buttocks. With next-generation sequencing (NGS) we compared mRNA expression between BCC and peritumoural skin. qRT-PCR, immunohistochemical and immunofluorescent staining were performed to validate the NGS results and to investigate CAF-related cyto-and chemokines. RESULTS: NGS revealed upregulation of 65 genes in BCC coding for extracellular matrix components pointing at CAF-related matrix remodeling. qRT-PCR showed increased mRNA expression of CAF markers FAP-α, PDGFR-ß and prolyl-4-hydroxylase in BCC. Peritumoural skin (but not buttock skin) also exhibited high expression of PDGFR-ß and prolyl-4-hydroxylase but not FAP-α. We found a similar pattern for the CAF-associated chemokines CCL17, CCL18, CCL22, CCL25, CXCL12 and IL6 with high expression in BCC and peritumoural skin but absence in buttock skin. Immunofluorescence revealed correlation between FAP-α and PDGFR-ß and CXCL12 and CCL17. CONCLUSION: Matrix remodeling is the most prominent molecular feature of BCC. CAFs are present within BCC stroma and associated with increased expression of chemokines involved in tumour progression and immunosuppression (CXCL12, CCL17). Fibroblasts from chronically sun-exposed skin near tumours show gene expression patterns resembling that of CAFs, indicating that stromal fibroblasts in cancer-free surgical BCC margins exhibit a tumour promoting phenotype.

Potential movement of transposable elements through DNA circularization.

The generation of circular DNAs is a relatively unrecognized type of genomic structural variation, but recent findings point to a possible role of circular DNAs in the movement of transposable elements. Circularization of genomic DNA is observed across eukaryotic species, in a range of different cell types, and from all parts of the genome. A recent study on circular DNAs in yeast found that transposable element sequence residing in circular structures mostly corresponded to full-length transposable elements. Transposable elements are mobile genetic elements scattered across eukaryotic genomes. Different classes of transposable elements move either through a copy-and-paste or a cut-and-paste. As circular DNA structures may recombine with the genome and re-integrate into a novel genomic locus, transposable elements could move through circularization. In yeast, the predominant type of transposable element is a so-called LTR (long terminal repeats) retrotransposable element that moves through a copy-and-paste mechanism. The observed circularization of this element means it potentially could move through a cut-and-paste mechanism as well. Although further experimental evidence is needed to establish the extent to which movement of transposable elements through DNA circularization takes place, such movement is likely to have a functional impact on the genomic context.

Propionibacterium acnes: Disease-Causing Agent or Common Contaminant? Detection in Diverse Patient Samples by Next-Generation Sequencing.

Propionibacterium acnesis the most abundant bacterium on human skin, particularly in sebaceous areas.P. acnesis suggested to be an opportunistic pathogen involved in the development of diverse medical conditions but is also a proven contaminant of human clinical samples and surgical wounds. Its significance as a pathogen is consequently a matter of debate. In the present study, we investigated the presence ofP. acnesDNA in 250 next-generation sequencing data sets generated from 180 samples of 20 different sample types, mostly of cancerous origin. The samples were subjected to either microbial enrichment, involving nuclease treatment to reduce the amount of host nucleic acids, or shotgun sequencing. We detected high proportions ofP. acnesDNA in enriched samples, particularly skin tissue-derived and other tissue samples, with the levels being higher in enriched samples than in shotgun-sequenced samples.P. acnesreads were detected in most samples analyzed, though the proportions in most shotgun-sequenced samples were low. Our results show thatP. acnescan be detected in practically all sample types when molecular methods, such as next-generation sequencing, are employed. The possibility of contamination from the patient or other sources, including laboratory reagents or environment, should therefore always be considered carefully whenP. acnesis detected in clinical samples. We advocate that detection ofP. acnesalways be accompanied by experiments validating the association between this bacterium and any clinical condition.

Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: Coccidia differing in host range and transmission strategy.

Toxoplasma gondii is a zoonotic protozoan parasite which infects nearly one third of the human population and is found in an extraordinary range of vertebrate hosts. Its epidemiology depends heavily on horizontal transmission, especially between rodents and its definitive host, the cat. Neospora caninum is a recently discovered close relative of Toxoplasma, whose definitive host is the dog. Both species are tissue-dwelling Coccidia and members of the phylum Apicomplexa; they share many common features, but Neospora neither infects humans nor shares the same wide host range as Toxoplasma, rather it shows a striking preference for highly efficient vertical transmission in cattle. These species therefore provide a remarkable opportunity to investigate mechanisms of host restriction, transmission strategies, virulence and zoonotic potential. We sequenced the genome of N. caninum and transcriptomes of the invasive stage of both species, undertaking an extensive comparative genomics and transcriptomics analysis. We estimate that these organisms diverged from their common ancestor around 28 million years ago and find that both genomes and gene expression are remarkably conserved. However, in N. caninum we identified an unexpected expansion of surface antigen gene families and the divergence of secreted virulence factors, including rhoptry kinases. Specifically we show that the rhoptry kinase ROP18 is pseudogenised in N. caninum and that, as a possible consequence, Neospora is unable to phosphorylate host immunity-related GTPases, as Toxoplasma does. This defense strategy is thought to be key to virulence in Toxoplasma. We conclude that the ecological niches occupied by these species are influenced by a relatively small number of gene products which operate at the host-parasite interface and that the dominance of vertical transmission in N. caninum may be associated with the evolution of reduced virulence in this species.

Evolutionary analysis identifies a Golgi pathway and correlates lineage-specific factors with endomembrane organelle emergence in apicomplexans.

The organelle paralogy hypothesis (OPH) aims to explain the evolution of non-endosymbiotically derived organelles. It predicts that lineage-specific pathways or organelles should result when identity-encoding membrane-trafficking components duplicate and co-evolve. Here, we investigate the presence of such lineage-specific membrane-trafficking machinery paralogs in Apicomplexa, a globally important parasitic lineage. We are able to identify 18 paralogs of known membrane-trafficking machinery, in several cases co-incident with the presence of new endomembrane organelles in apicomplexans or their parent lineage, the Alveolata. Moreover, focused analysis of the apicomplexan Arf-like small GTPases (i.e., ArlX3) revealed a specific post-Golgi trafficking pathway. This pathway appears involved in delivery of proteins to micronemes and rhoptries, with knockdown demonstrating reduced invasion capacity. Overall, our data have identified an unforeseen post-Golgi trafficking pathway in apicomplexans and are consistent with the OPH mechanism acting to produce endomembrane pathways or organelles at various evolutionary stages across the alveolate lineage.

Statistical guidelines for detecting past population shifts using ancient DNA.

Populations carry a genetic signal of their demographic past, providing an opportunity for investigating the processes that shaped their evolution. Our ability to infer population histories can be enhanced by including ancient DNA data. Using serial-coalescent simulations and a range of both quantitative and temporal sampling schemes, we test the power of ancient mitochondrial sequences and nuclear single-nucleotide polymorphisms (SNPs) to detect past population bottlenecks. Within our simulated framework, mitochondrial sequences have only limited power to detect subtle bottlenecks and/or fast post-bottleneck recoveries. In contrast, nuclear SNPs can detect bottlenecks followed by rapid recovery, although bottlenecks involving reduction of less than half the population are generally detected with low power unless extensive genetic information from ancient individuals is available. Our results provide useful guidelines for scaling sampling schemes and for optimizing our ability to infer past population dynamics. In addition, our results suggest that many ancient DNA studies may face power issues in detecting moderate demographic collapses and/or highly dynamic demographic shifts when based solely on mitochondrial information.

miR-34c-3p Regulates Protein Kinase A Activity Independent of cAMP by Dicing prkar2b Transcripts in Theileria annulata-Infected Leukocytes.

MicroRNAs (miRNAs) are small noncoding RNAs that can play critical roles in regulating various cellular processes, including during many parasitic infections. Here, we report a regulatory role for miR-34c-3p in cAMP-independent regulation of host cell protein kinase A (PKA) activity in Theileria annulata-infected bovine leukocytes. We identified prkar2b (cAMP-dependent protein kinase A type II-beta regulatory subunit) as a novel miR-34c-3p target gene and demonstrate how infection-induced upregulation of miR-34c-3p repressed PRKAR2B expression to increase PKA activity. As a result, the disseminating tumorlike phenotype of T. annulata-transformed macrophages is enhanced. Finally, we extend our observations to Plasmodium falciparum-parasitized red blood cells, where infection-induced augmentation in miR-34c-3p levels led to a drop in the amount of prkar2b mRNA and increased PKA activity. Collectively, our findings represent a novel cAMP-independent way of regulating host cell PKA activity in infections by Theileria and Plasmodium parasites. IMPORTANCE Small microRNA levels are altered in many diseases, including those caused by parasites. Here, we describe how infection by two important animal and human parasites, Theileria annulata and Plasmodium falciparum, induce changes in infected host cell miR-34c-3p levels to regulate host cell PKA kinase activity by targeting mammalian prkar2b. Infection-induced changes in miR-34c-3p levels provide a novel epigenetic mechanism for regulating host cell PKA activity independent of fluxes in cAMP to both aggravate tumor dissemination and improve parasite fitness.