Search

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., Competing Interests: Declaration of interests The authors have no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Malaria-associated pathogenesis such as parasite invasion, egress, host cell remodelling and antigenic variation requires concerted action by many proteins, but the molecular regulation is poorly understood. Here we have characterized an essential Plasmodium-specific Apicomplexan AP2 transcription factor in Plasmodium falciparum (PfAP2-P; pathogenesis) during the blood-stage development with two peaks of expression. An inducible knockout of gene function showed that PfAP2-P is essential for trophozoite development, and critical for var gene regulation, merozoite development and parasite egress. Chromatin immunoprecipitation sequencing data collected at timepoints matching the two peaks of pfap2-p expression demonstrate PfAP2-P binding to promoters of genes controlling trophozoite development, host cell remodelling, antigenic variation and pathogenicity. Single-cell RNA sequencing and fluorescence-activated cell sorting revealed de-repression of most var genes in Δpfap2-p parasites. Δpfap2-p parasites also overexpress early gametocyte marker genes, indicating a regulatory role in sexual stage conversion. We conclude that PfAP2-P is an essential upstream transcriptional regulator at two distinct stages of the intra-erythrocytic development cycle., (© 2023. The Author(s).)

Background: The excessive inflammatory responses provoked by SARS-CoV-2 infection are critical factors affecting the severity and mortality of COVID-19. Previous work found that two adjacent co-occurring mutations R203K and G204R (KR) on the nucleocapsid (N) protein correlate with increased disease severity in COVID-19 patients. However, links with the host immune response remain unclear., Methods: Here, we grouped nasopharyngeal swab samples of COVID-19 patients into two cohorts based on the presence and absence of SARS-CoV-2 nucleocapsid KR mutations. We performed nasopharyngeal transcriptome analysis of age, gender, and ethnicity-matched COVID-19 patients infected with either SARS-CoV-2 with KR mutations in the N protein (KR patients n = 39) or with the wild-type N protein (RG patients n = 39) and compared to healthy controls (n = 34). The impact of KR mutation on immune response was further characterized experimentally by transcriptomic and proteomic profiling of virus-like-particle (VLP) incubated cells., Results: We observed markedly elevated expression of proinflammatory cytokines, chemokines, and interferon-stimulated (ISGs) genes in the KR patients compared to RG patients. Using nasopharyngeal transcriptome data, we found significantly higher levels of neutrophils and neutrophil-to-lymphocyte (NLR) ratio in KR patients than in the RG patients. Furthermore, transcriptomic and proteomic profiling of VLP incubated cells confirmed a similar hyper-inflammatory response mediated by the KR variant., Conclusions: Our data demonstrate an unforeseen connection between nucleocapsid KR mutations and augmented inflammatory immune response in severe COVID-19 patients. These findings provide insights into how mutations in SARS-CoV-2 modulate host immune output and pathogenesis and may contribute to more efficient therapeutics and vaccine development., (© 2023. The Author(s).)

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.

The outcome of transplant recipients is variable depending on the study population, vaccination status and COVID-19 variants. Our aim was to study the impact of Omicron subvariants on the mortality of transplant recipients. We reviewed the results of SARS-CoV-2 whole genome sequence of random isolates collected from 29 December 2021 until 17 May 2022 in King Faisal Specialist Hospital and Research center, Jeddah (KFSHRC-J), Saudi Arabia performed as hospital genomic surveillance program for COVID-19 variants. We included 25 transplant patients infected with confirmed Omicron variants.17 (68%) and 8 (32%) patients had Omicron BA.1 and BA.2, respectively. 12 (68%) patients had renal transplants. Only 36% of patients received three doses of COVID-19 vaccines. 23 (92%) patients required hospitalization. 20 (80%) patients survived and 6 (25%) required intensive care unit (ICU) admission. Among ICU patients, 66.7% were more than 50 years, 50% had two to three comorbidities and 5 out of 6 (83%) died. The mortality of transplant patients infected with Omicron variants in our cohort was higher than other centers as a limited number of patients received booster vaccines. Optimizing booster vaccination is the most efficient method to improve the mortality of COVID-19 in transplant recipients recognizing the inefficacy of monoclonal antibodies in the presence of SARS-CoV-2 emerging variants. We did not show a difference in mortality in transplant patients infected with Omicron BA.1 and BA.2 knowing the limitation of our sample size., (© 2023. The Author(s).)

A fungal metabolite, FR235222, specifically inhibits a histone deacetylase of the apicomplexan parasite Toxoplasma gondii and TgHDAC3 has emerged as a key factor regulating developmental stage transition in this species. Here, we exploited FR235222 to ask if changes in histone acetylation regulate developmental stage transition of Theileria annulata, another apicomplexan species. We found that FR235222 treatment of T. annulata-infected transformed leukocytes induced a proliferation arrest. The blockade in proliferation was due to drug-induced conversion of intracellular schizonts to merozoites that lack the ability to maintain host leukocyte cell division. Induction of merogony by FR235222 leads to an increase in expression of merozoite-marker (rhoptry) proteins. RNA-seq of FR235222-treated T. annulata-infected B cells identified deregulated expression of 468 parasite genes including a number encoding parasite ApiAP2 transcription factors. Thus, similar to T. gondii, FR235222 inhibits T. annulata HDAC (TaHDAC1) activity and places parasite histone acetylation as a major regulatory event of the transition from schizonts to merozoites., (© 2022. The Author(s).)

Monitoring SARS-CoV-2 spread and evolution through genome sequencing is essential in handling the COVID-19 pandemic. Here, we sequenced 892 SARS-CoV-2 genomes collected from patients in Saudi Arabia from March to August 2020. We show that two consecutive mutations (R203K/G204R) in the nucleocapsid (N) protein are associated with higher viral loads in COVID-19 patients. Our comparative biochemical analysis reveals that the mutant N protein displays enhanced viral RNA binding and differential interaction with key host proteins. We found increased interaction of GSK3A kinase simultaneously with hyper-phosphorylation of the adjacent serine site (S206) in the mutant N protein. Furthermore, the host cell transcriptome analysis suggests that the mutant N protein produces dysregulated interferon response genes. Here, we provide crucial information in linking the R203K/G204R mutations in the N protein to modulations of host-virus interactions and underline the potential of the nucleocapsid protein as a drug target during infection., (© 2022. The Author(s).)

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., (© 2021. The Author(s).)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., Competing Interests: Declaration of competing interest None declared., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Apicomplexan parasites cause major human disease and food insecurity. They owe their considerable success to highly specialized cell compartments and structures. These adaptations drive their recognition, nondestructive penetration, and elaborate reengineering of the host's cells to promote their growth, dissemination, and the countering of host defenses. The evolution of unique apicomplexan cellular compartments is concomitant with vast proteomic novelty. Consequently, half of apicomplexan proteins are unique and uncharacterized. Here, we determine the steady-state subcellular location of thousands of proteins simultaneously within the globally prevalent apicomplexan parasite Toxoplasma gondii. This provides unprecedented comprehensive molecular definition of these unicellular eukaryotes and their specialized compartments, and these data reveal the spatial organizations of protein expression and function, adaptation to hosts, and the underlying evolutionary trajectories of these pathogens., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Cell cycle transitions are generally triggered by variation in the activity of cyclin-dependent kinases (CDKs) bound to cyclins. Malaria-causing parasites have a life cycle with unique cell-division cycles, and a repertoire of divergent CDKs and cyclins of poorly understood function and interdependency. We show that Plasmodium berghei CDK-related kinase 5 (CRK5), is a critical regulator of atypical mitosis in the gametogony and is required for mosquito transmission. It phosphorylates canonical CDK motifs of components in the pre-replicative complex and is essential for DNA replication. During a replicative cycle, CRK5 stably interacts with a single Plasmodium -specific cyclin (SOC2), although we obtained no evidence of SOC2 cycling by transcription, translation or degradation. Our results provide evidence that during Plasmodium male gametogony, this divergent cyclin/CDK pair fills the functional space of other eukaryotic cell-cycle kinases controlling DNA replication., Competing Interests: AB, MZ, ER, CP, LB, TM, AS, NK, PA, RP, DB, SV, AH, AP, DF, AH, RT, MB No competing interests declared, (© 2020, Balestra et al.)

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., Competing Interests: The authors have declared that no competing interests exist.

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., (© The Author(s) 2019. Published by Oxford University Press for the Infectious Diseases Society of America.)

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions., (© 2018 John Wiley & Sons Ltd.)

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., Competing Interests: The authors have declared that no competing interests exist.

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.

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.

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.

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., Competing Interests: The authors declare no conflict of interest.

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.

Circular DNAs are extra-chromosomal fragments that become circularized by genomic recombination events. We have recently shown that yeast LTR elements generate circular DNAs through recombination events between their flanking long terminal repeats (LTRs). Similarly, circular DNAs can be generated by recombination between LTRs residing at different genomic loci, in which case the circular DNA will contain the intervening sequence. In yeast, this can result in gene copy number variations when circles contain genes and origins of replication. Here, I speculate on the potential and implications of circular DNAs generated through recombination between human transposable elements.

Outbreaks of zoonotic diseases in humans and livestock are not uncommon, and an important component in containment of such emerging viral diseases is rapid and reliable diagnostics. Such methods are often PCR-based and hence require the availability of sequence data from the pathogen. Rattus norvegicus (R. norvegicus) is a known reservoir for important zoonotic pathogens. Transmission may be direct via contact with the animal, for example, through exposure to its faecal matter, or indirectly mediated by arthropod vectors. Here we investigated the viral content in rat faecal matter (n=29) collected from two continents by analyzing 2.2 billion next-generation sequencing reads derived from both DNA and RNA. Among other virus families, we found sequences from members of the Picornaviridae to be abundant in the microbiome of all the samples. Here we describe the diversity of the picornavirus-like contigs including near-full-length genomes closely related to the Boone cardiovirus and Theiler's encephalomyelitis virus. From this study, we conclude that picornaviruses within R. norvegicus are more diverse than previously recognized. The virome of R. norvegicus should be investigated further to assess the full potential for zoonotic virus transmission.

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., (Copyright © 2016 Mollerup et al.)

Virus discovery from high throughput sequencing data often follows a bottom-up approach where taxonomic annotation takes place prior to association to disease. Albeit effective in some cases, the approach fails to detect novel pathogens and remote variants not present in reference databases. We have developed a species independent pipeline that utilises sequence clustering for the identification of nucleotide sequences that co-occur across multiple sequencing data instances. We applied the workflow to 686 sequencing libraries from 252 cancer samples of different cancer and tissue types, 32 non-template controls, and 24 test samples. Recurrent sequences were statistically associated to biological, methodological or technical features with the aim to identify novel pathogens or plausible contaminants that may associate to a particular kit or method. We provide examples of identified inhabitants of the healthy tissue flora as well as experimental contaminants. Unmapped sequences that co-occur with high statistical significance potentially represent the unknown sequence space where novel pathogens can be identified.

Extrachromosomal circular DNA (eccDNA) derived from chromosomal Ty retrotransposons in yeast can be generated in multiple ways. Ty eccDNA can arise from the circularization of extrachromosomal linear DNA during the transpositional life cycle of retrotransposons, or from circularization of genomic Ty DNA. Circularization may happen through nonhomologous end-joining (NHEJ) of long terminal repeats (LTRs) flanking Ty elements, by Ty autointegration, or by LTR-LTR recombination. By performing an in-depth investigation of sequence reads stemming from Ty eccDNAs obtained from populations of Saccharomyces cerevisiae S288c, we find that eccDNAs predominantly correspond to full-length Ty1 elements. Analyses of sequence junctions reveal no signs of NHEJ or autointegration events. We detect recombination junctions that are consistent with yeast Ty eccDNAs being generated through recombination events within the genome. This opens the possibility that retrotransposable elements could move around in the genome without an RNA intermediate directly through DNA circularization., (Copyright © 2016 by the Genetics Society of America.)

From Illumina sequencing of DNA from brain and liver tissue from the lion, Panthera leo, and tumor samples from the pike-perch, Sander lucioperca, we obtained two assembled sequence contigs with similarity to known retroviruses. Phylogenetic analyses suggest that the pike-perch retrovirus belongs to the epsilonretroviruses, and the lion retrovirus to the gammaretroviruses. To determine if these novel retroviral sequences originate from an endogenous retrovirus or from a recently integrated exogenous retrovirus, we assessed the genetic diversity of the parental sequences from which the short Illumina reads are derived. First, we showed by simulations that we can robustly infer the level of genetic diversity from short sequence reads. Second, we find that the measures of nucleotide diversity inferred from our retroviral sequences significantly exceed the level observed from Human Immunodeficiency Virus infections, prompting us to conclude that the novel retroviruses are both of endogenous origin. Through further simulations, we rule out the possibility that the observed elevated levels of nucleotide diversity are the result of co-infection with two closely related exogenous retroviruses.

Although nearly one fifth of all human cancers have an infectious aetiology, the causes for the majority of cancers remain unexplained. Despite the enormous data output from high-throughput shotgun sequencing, viral DNA in a clinical sample typically constitutes a proportion of host DNA that is too small to be detected. Sequence variation among virus genomes complicates application of sequence-specific, and highly sensitive, PCR methods. Therefore, we aimed to develop and characterize a method that permits sensitive detection of sequences despite considerable variation. We demonstrate that our low-stringency in-solution hybridization method enables detection of <100 viral copies. Furthermore, distantly related proviral sequences may be enriched by orders of magnitude, enabling discovery of hitherto unknown viral sequences by high-throughput sequencing. The sensitivity was sufficient to detect retroviral sequences in clinical samples. We used this method to conduct an investigation for novel retrovirus in samples from three cancer types. In accordance with recent studies our investigation revealed no retroviral infections in human B-cell lymphoma cells, cutaneous T-cell lymphoma or colorectal cancer biopsies. Nonetheless, our generally applicable method makes sensitive detection possible and permits sequencing of distantly related sequences from complex material.

The eukaryotic phylum Apicomplexa encompasses thousands of obligate intracellular parasites of humans and animals with immense socio-economic and health impacts. We sequenced nuclear genomes of Chromera velia and Vitrella brassicaformis, free-living non-parasitic photosynthetic algae closely related to apicomplexans. Proteins from key metabolic pathways and from the endomembrane trafficking systems associated with a free-living lifestyle have been progressively and non-randomly lost during adaptation to parasitism. The free-living ancestor contained a broad repertoire of genes many of which were repurposed for parasitic processes, such as extracellular proteins, components of a motility apparatus, and DNA- and RNA-binding protein families. Based on transcriptome analyses across 36 environmental conditions, Chromera orthologs of apicomplexan invasion-related motility genes were co-regulated with genes encoding the flagellar apparatus, supporting the functional contribution of flagella to the evolution of invasion machinery. This study provides insights into how obligate parasites with diverse life strategies arose from a once free-living phototrophic marine alga.

Mycoplasma capricolum subsp. capricolum is a well-known pathogen of small ruminants. A recent human case of septicemia involving this agent raised the question of its potential pathogenicity to humans. We present the first draft genome sequence of a human Mycoplasma capricolum subsp. capricolum isolate., (Copyright © 2015 Seersholm et al.)

Toxoplasma gondii is an important protozoan parasite that infects all warm-blooded animals and causes opportunistic infections in immuno-compromised humans. Its closest relative, Neospora caninum, is an important veterinary pathogen that causes spontaneous abortion in livestock. Comparative genomics of these two closely related coccidians has been of particular interest to identify genes that contribute to varied host cell specificity and disease. Here, we describe a manual evaluation of these genomes based on strand-specific RNA sequencing and shotgun proteomics from the invasive tachyzoite stages of these two parasites. We have corrected predicted structures of over one third of the previously annotated gene models and have annotated untranslated regions (UTRs) in over half of the predicted protein-coding genes. We observe distinctly long UTRs in both the organisms, almost four times longer than other model eukaryotes. We have also identified a putative set of cis-natural antisense transcripts (cis-NATs) and long intergenic non-coding RNAs (lincRNAs). We have significantly improved the annotation quality in these genomes that would serve as a manually curated dataset for Toxoplasma and Neospora research communities.

Viral infections cause many different diseases stemming both from well-characterized viral pathogens but also from emerging viruses, and the search for novel viruses continues to be of great importance. High-throughput sequencing is an important technology for this purpose. However, viral nucleic acids often constitute a minute proportion of the total genetic material in a sample from infected tissue. Techniques to enrich viral targets in high-throughput sequencing have been reported, but the sensitivity of such methods is not well established. This study compares different library preparation techniques targeting both DNA and RNA with and without virion enrichment. By optimizing the selection of intact virus particles, both by physical and enzymatic approaches, we assessed the effectiveness of the specific enrichment of viral sequences as compared to non-enriched sample preparations by selectively looking for and counting read sequences obtained from shotgun sequencing. Using shotgun sequencing of total DNA or RNA, viral targets were detected at concentrations corresponding to the predicted level, providing a foundation for estimating the effectiveness of virion enrichment. Virion enrichment typically produced a 1000-fold increase in the proportion of DNA virus sequences. For RNA virions the gain was less pronounced with a maximum 13-fold increase. This enrichment varied between the different sample concentrations, with no clear trend. Despite that less sequencing was required to identify target sequences, it was not evident from our data that a lower detection level was achieved by virion enrichment compared to shotgun sequencing.

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits., (Copyright © 2014, American Association for the Advancement of Science.)

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., (© 2014 Reid et al.; Published by Cold Spring Harbor Laboratory Press.)

Transposable elements (TEs) are ubiquitous in eukaryotic genomes. Barbara McClintock's famous notion of TEs acting as controlling elements modifying the genetic response of an organism upon exposure to stressful environments has since been solidly supported in a series of model organisms. This requires the TE activity response to possess an element of specificity and be targeted toward certain parts of the genome. We propose that a similar TE response is present in human cells, and that this stress response may drive the onset of human cancers. As such, TE-driven cancers may be viewed as an evolutionary by-product of organisms' abilities to genetically adapt to environmental stress.

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.

Molecular detection of viruses has been aided by high-throughput sequencing, permitting the genomic characterization of emerging strains. In this study, we comprehensively screened 500 respiratory secretions from children with upper and/or lower respiratory tract infections for viral pathogens. The viruses detected are described, including a divergent human parainfluenza virus type 4 from GS FLX pyrosequencing of 92 specimens. Complete full-genome characterization of the virus followed, using Single Molecule, Real-Time (SMRT) sequencing. Subsequent "primer walking" combined with Sanger sequencing validated the RS platform's utility in viral sequencing from complex clinical samples. Comparative genomics reveals the divergent strain clusters with the only completely sequenced HPIV4a subtype. However, it also exhibits various structural features present in one of the HPIV4b reference strains, opening questions regarding their lifecycle and evolutionary relationships among these viruses. Clinical data from patients infected with the strain, as well as viral prevalence estimates using real-time PCR, is also described.

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.

In 2010, a chimpanzee died at Copenhagen Zoo following an outbreak of respiratory disease among chimpanzees in the zoo. Identification of coxsackie B3 virus, a common human pathogen, as the causative agent, and its severe manifestation, raise questions about pathogenicity and transmissibility among humans and other primates.

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.

Background: Piwi-associated RNAs (piRNAs) bind transcripts from retrotransposable elements (RTE) in mouse germline cells and seemingly act as guides for genomic methylation, thereby repressing the activity of RTEs. It is currently unknown if and how Piwi proteins distinguish RTE transcripts from other cellular RNAs. During germline development, the main target of piRNAs switch between different types of RTEs. Using the piRNA targeting of RTEs as an indicator of RTE activity, and considering the entire population of genomic RTE loci along with their age and location, this study aims at further elucidating the dynamics of RTE activity during mouse germline development., Results: Due to the inherent sequence redundancy between RTE loci, assigning piRNA targeting to specific loci is problematic. This limits the analysis, although certain features of piRNA targeting of RTE loci are apparent. As expected, young RTEs display a much higher level of piRNA targeting than old RTEs. Further, irrespective of age, RTE loci near protein-coding coding genes are targeted to a greater extent than RTE loci far from genes. During development, a shift in piRNA targeting is observed, with a clear increase in the relative piRNA targeting of RTEs residing within boundaries of protein-coding gene transcripts., Conclusions: Reanalyzing published piRNA sequences and taking into account the features of individual RTE loci provide novel insight into the activity of RTEs during development. The obtained results are consistent with some degree of proportionality between what transcripts become substrates for Piwi protein complexes and the level by which the transcripts are present in the cell. A transition from active transcription of RTEs to passive co-transcription of RTE sequences residing within protein-coding transcripts appears to take place in postnatal development. Hence, the previously reported increase in piRNA targeting of SINEs in postnatal testis development does not necessitate widespread active transcription of SINEs, but may simply be explained by the prevalence of SINEs residing in introns.

Background: Retrotransposons are transposable elements that proliferate within eukaryotic genomes through a process involving reverse transcription. The numbers of retrotransposons within genomes and differences between closely related species may yield insight into the evolutionary history of the elements. Less is known about the ongoing dynamics of retrotransposons, as analysis of genome sequences will only reveal insertions of retrotransposons that are fixed--or near fixation--in the population or strain from which genetic material has been extracted for sequencing. One pre-requisite for retrotransposition is transcription of the elements. Given their intrinsic sequence redundancy, transcriptome-level analyses of transposable elements are scarce. We have used recently published transcriptome data from the fission yeast Schizosaccharomyces pombe to assess the ability to detect and describe transcriptional activity from Long Terminal Repeat (LTR) retrotransposons. LTR retrotransposons are normally flanked by two LTR sequences. However, the majority of LTR sequences in S. pombe exist as solitary LTRs, i.e. as single terminal repeat sequences not flanking a retrotransposon. Transcriptional activity was analysed for both full-length LTR retrotransposons and solitary LTRs., Results: Two independent sets of transcriptome data reveal the presence of full-length, polyadenylated transcripts from LTR retrotransposons in S. pombe during growth phase in rich medium. The redundancy of retrotransposon sequences makes it difficult to assess which elements are transcriptionally active, but data strongly indicates that only a subset of the LTR retrotransposons contribute significantly to the detected transcription. A considerable level of reverse strand transcription is also detected. Equal levels of transcriptional activity are observed from both strands of solitary LTR sequences. Transcriptome data collected during meiosis suggests that transcription of solitary LTRs is correlated with the transcription of nearby protein-coding genes., Conclusions: Presumably, the host organism negatively regulates proliferation of LTR retrotransposons. The finding of considerable transcriptional activity of retrotransposons suggests that part of this regulation is likely to take place at a post-transcriptional level. Alternatively, the transcriptional activity may signify a hitherto unrecognized activity level of retrotransposon proliferation. Our findings underline the usefulness of transcriptome data in elucidating dynamics in retrotransposon transcription.

Retrotransposons constitute a significant fraction of mammalian genomes. Considering the finding of widespread transcriptional activity across entire genomes, it is not surprising that retrotransposons contribute to the collective RNA pool. However, the transcriptional output from retrotransposons does not merely represent spurious transcription. We review examples of functional RNAs transcribed from retrotransposons, and address the collection of non-protein coding RNAs derived from transposable element sequences, including numerous human microRNAs and the neuronal BC RNAs. Finally, we review the emerging understanding of how retrotransposons themselves are regulated by small RNAs.

We undertook a genome-wide search for novel noncoding RNAs (ncRNA) in the malaria parasite Plasmodium falciparum. We used the RNAz program to predict structures in the noncoding regions of the P. falciparum 3D7 genome that were conserved with at least one of seven other Plasmodium spp. genome sequences. By using Northern blot analysis for 76 high-scoring predictions and microarray analysis for the majority of candidates, we have verified the expression of 33 novel ncRNA transcripts including four members of a ncRNA family in the asexual blood stage. These transcripts represent novel structured ncRNAs in P. falciparum and are not represented in any RNA databases. We provide supporting evidence for purifying selection acting on the experimentally verified ncRNAs by comparing the nucleotide substitutions in the predicted ncRNA candidate structures in P. falciparum with the closely related chimp malaria parasite P. reichenowi. The high confirmation rate within a single parasite life cycle stage suggests that many more of the predictions may be expressed in other stages of the organism's life cycle.