Whole genome sequencing of Theileria parva using target capture.

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.

The 4-Celled Tetrabaena socialis Nuclear Genome Reveals the Essential Components for Genetic Control of Cell Number at the Origin of Multicellularity in the Volvocine Lineage.

Multicellularity is the premier example of a major evolutionary transition in individuality and was a foundational event in the evolution of macroscopic biodiversity. The volvocine chlorophyte lineage is well suited for studying this process. Extant members span unicellular, simple colonial, and obligate multicellular taxa with germ-soma differentiation. Here, we report the nuclear genome sequence of one of the most morphologically simple organisms in this lineage-the 4-celled colonial Tetrabaena socialis and compare this to the three other complete volvocine nuclear genomes. Using conservative estimates of gene family expansions a minimal set of expanded gene families was identified that associate with the origin of multicellularity. These families are rich in genes related to developmental processes. A subset of these families is lineage specific, which suggests that at a genomic level the evolution of multicellularity also includes lineage-specific molecular developments. Multiple points of evidence associate modifications to the ubiquitin proteasomal pathway (UPP) with the beginning of coloniality. Genes undergoing positive or accelerating selection in the multicellular volvocines were found to be enriched in components of the UPP and gene families gained at the origin of multicellularity include components of the UPP. A defining feature of colonial/multicellular life cycles is the genetic control of cell number. The genomic data presented here, which includes diversification of cell cycle genes and modifications to the UPP, align the genetic components with the evolution of this trait.

Molecular characterization of Morogoro maize-associated virus, a nucleorhabdovirus detected in maize (Zea mays) in Tanzania.

RNAtag-seq of maize samples collected in Tanzania revealed the presence of a previously undescribed nucleorhabdovirus, tentatively named "Morogoro maize-associated virus" (MMaV), in three samples. The MMaV genome is 12,185-12,187 nucleotides long and shares a 69-70% nucleotide sequence identity with taro vein chlorosis virus. Annotation of the genomes showed a typical nucleorhabdovirus gene organization. PCR was unable to detect the same virus in the remaining 35 samples collected in the region.

Diversity and distribution of Maize-associated totivirus strains from Tanzania.

Typically associated with fungal species, members of the viral family Totiviridae have recently been shown to be associated with plants, including important crop species, such as Carica papaya (papaya) and Zea mays (maize). Maize-associated totivirus (MATV) was first described in China and more recently in Ecuador, where it has been found to co-occur with other viruses known to elicit maize lethal necrosis disease (MLND). In a survey for maize-associated viruses, 35 samples were selected for Illumina HiSeq sequencing, from the Tanzanian maize producing regions of Mara, Arusha, Manyara, Kilimanjaro, Morogoro and Pwani. Libraries were prepared using an RNA-tag-seq methodology. Taxonomic classification of the resulting datasets showed that 6 of the 35 samples from the regions of Arusha, Kilimanjaro, Morogoro and Mara, contained reads that were assigned to MATV reference sequences. This was confirmed with PCR and Sanger sequencing. Read assembly of the six MATV-associated datasets yielded partial MATV genomes, two of which were selected for further characterization, using RACE. This yielded two full-length MATV genomes, one of which is divergent from other available MATV genomes.

Evolution of cytokinesis-related protein localization during the emergence of multicellularity in volvocine green algae.

BACKGROUND: The volvocine lineage, containing unicellular Chlamydomonas reinhardtii and differentiated multicellular Volvox carteri, is a powerful model for comparative studies aiming at understanding emergence of multicellularity. Tetrabaena socialis is the simplest multicellular volvocine alga and belongs to the family Tetrabaenaceae that is sister to more complex multicellular volvocine families, Goniaceae and Volvocaceae. Thus, T. socialis is a key species to elucidate the initial steps in the evolution of multicellularity. In the asexual life cycle of C. reinhardtii and multicellular volvocine species, reproductive cells form daughter cells/colonies by multiple fission. In embryogenesis of the multicellular species, daughter protoplasts are connected to one another by cytoplasmic bridges formed by incomplete cytokinesis during multiple fission. These bridges are important for arranging the daughter protoplasts in appropriate positions such that species-specific integrated multicellular individuals are shaped. Detailed comparative studies of cytokinesis between unicellular and simple multicellular volvocine species will help to elucidate the emergence of multicellularity from the unicellular ancestor. However, the cytokinesis-related genes between closely related unicellular and multicellular species have not been subjected to a comparative analysis. RESULTS: Here we focused on dynamin-related protein 1 (DRP1), which is known for its role in cytokinesis in land plants. Immunofluorescence microscopy using an antibody against T. socialis DRP1 revealed that volvocine DRP1 was localized to division planes during cytokinesis in unicellular C. reinhardtii and two simple multicellular volvocine species T. socialis and Gonium pectorale. DRP1 signals were mainly observed in the newly formed division planes of unicellular C. reinhardtii during multiple fission, whereas in multicellular T. socialis and G. pectorale, DRP1 signals were observed in all division planes during embryogenesis. CONCLUSIONS: These results indicate that the molecular mechanisms of cytokinesis may be different in unicellular and multicellular volvocine algae. The localization of DRP1 during multiple fission might have been modified in the common ancestor of multicellular volvocine algae. This modification may have been essential for the re-orientation of cells and shaping colonies during the emergence of multicellularity in this lineage.

Organelle genome complexity scales positively with organism size in volvocine green algae.

It has been argued that for certain lineages, noncoding DNA expansion is a consequence of the increased random genetic drift associated with long-term escalations in organism size. But a lack of data has prevented the investigation of this hypothesis in most plastid-bearing protists. Here, using newly sequenced mitochondrial and plastid genomes, we explore the relationship between organelle DNA noncoding content and organism size within volvocine green algae. By looking at unicellular, colonial, and differentiated multicellular algae, we show that organelle DNA complexity scales positively with species size and cell number across the volvocine lineage. Moreover, silent-site genetic diversity data suggest that the volvocine species with the largest cell numbers and most bloated organelle genomes have the smallest effective population sizes. Together, these findings support the view that nonadaptive processes, like random genetic drift, promote the expansion of noncoding regions in organelle genomes.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam.

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

De novo assembled salivary gland transcriptome and expression pattern analyses for Rhipicephalus evertsi evertsi Neuman, 1897 male and female ticks.

Ticks secrete proteins in their saliva that change over the course of feeding to modulate the host inflammation, immune responses, haemostasis or may cause paralysis. RNA next generation sequencing technologies can reveal the complex dynamics of tick salivary glands as generated from various tick life stages and/or males and females. The current study represents 15,115 Illumina sequenced contigs of the salivary gland transcriptome from male and female Rhipicephalus evertsi evertsi ticks of early, mid and late feeding stages from 1320 separate assemblies using three short read assemblers. The housekeeping functional class contributed to the majority of the composition of the transcriptome (80%) but with lower expression (51%), while the secretory protein functional class represented only 14% of the transcriptome but 46% of the total coverage. Six percent had an unknown status contributing 3% of the overall expression in the salivary glands. Platelet aggregation inhibitors, blood clotting inhibitors and immune-modulators orthologous to the ancestral tick lineages were confirmed in the transcriptome and their differential expression during feeding in both genders observed. This transcriptome contributes data of importance to salivary gland biology and blood feeding physiology of non-model organisms.

Transcriptomic Analysis of Salivary Glands of Ornithodoros brasiliensis Aragão, 1923, the Agent of a Neotropical Tick-Toxicosis Syndrome in Humans.

Tick salivary glands produce and secrete a variety of compounds that modulate host responses and ensure a successful blood meal. Despite great progress made in the identification of ticks salivary compounds in recent years, there is still a paucity of information concerning salivary molecules of Neotropical argasid ticks. Among this group of ticks, considering the number of human cases of parasitism, including severe syndromes and hospitalization, Ornithodoros brasiliensis can be considered one of the major Neotropical argasid species with impact in public health. Here, we describe the transcriptome analysis of O. brasiliensis salivary glands (ObSG). The transcriptome yielded ~14,957 putative contigs. A total of 368 contigs were attributed to secreted proteins (SP), which represent approximately 2.5% of transcripts but ~53% expression coverage transcripts per million. Lipocalins are the major protein family among the most expressed SP, accounting for ~16% of the secretory transcripts and 51% of secretory protein abundance. The most expressed transcript is an ortholog of TSGP4 (tick salivary gland protein 4), a lipocalin first identified in Ornithodoros kalahariensis that functions as a leukotriene C4 scavenger. A total of 55 lipocalin transcripts were identified in ObSG. Other transcripts potentially involved in tick-host interaction included as: basic/acid tail secretory proteins (second most abundant expressed group), serine protease inhibitors (including Kunitz inhibitors), 5' nucleotidases (tick apyrases), phospholipase A2, 7 disulfide bond domain, cystatins, and tick antimicrobial peptides. Another abundant group of proteins in ObSG is metalloproteases. Analysis of these major protein groups suggests that several duplication events after speciation were responsible for the abundance of redundant compounds in tick salivary glands. A full mitochondrial genome could be assembled from the transcriptome data and confirmed the close genetic identity of the tick strain sampled in the current study, to a tick strain previously implicated in tick toxicoses. This study provides novel information on the molecular composition of ObSG, a Brazilian endemic tick associated with several human cases of parasitism. These results could be helpful in the understanding of clinical findings observed in bitten patients, and also, could provide more information on the evolution of Neotropical argasids.

Nuclear (18S-28S rRNA) and mitochondrial genome markers of Carios (Carios) vespertilionis (Argasidae) support Carios Latreille, 1796 as a lineage embedded in the Ornithodorinae: re-classification of the Carios sensu Klompen and Oliver (1993) clade into its respective subgenera.

Argasid systematics remains controversial with widespread adherence to the Hoogstraal (1985) classification scheme, even though it does not reflect evolutionary relationships and results in paraphyly for the main genera of soft ticks (Argasidae), namely Argas and Ornithodoros. The alternative classification scheme, proposed by Klompen and Oliver (1993), has problems of its own: most notably paraphyly of the subgenus Pavlovskyella and the controversial grouping together of the subgenera Alectorobius, Antricola, Carios, Chiropterargas, Nothoaspis, Parantricola, Reticulinasus and Subparmatus into the genus Carios. Recent phylogenetic analyses of 18S/28S rRNA sequences and mitochondrial genomes agree with the scheme of Klompen and Oliver (1993), with regard to the paraphyly of Pavlovskyella, placement of Alveonasus, Ogadenus, Proknekalia and Secretargas in the Argasinae and placement of Carios and Chiropterargas in the Ornithodorinae (Mans et al., 2019). The Carios clade and its constituent subgenera remain controversial, since the phylogenetic position of its type species Carios (Carios) vespertilionis Latreille, 1796 (formerly Argas vespertilionis) has not been determined with confidence. The current study aimed to resolve Carios sensu lato Klompen and Oliver, 1993, and Carios sensu stricto Hoogstraal, 1985, by determining and analysing phylogenetic nuclear and mitochondrial markers for C. (C.) vespertilionis. Both the nuclear and mitochondrial markers support placement of Carios s.s. within the subfamily Ornithodorinae, but to the exclusion of the clade that includes the 6 other subgenera that are part of Carios s.l. Klompen and Oliver (1993), namely Alectorobius, Antricola, Nothoaspis, Parantricola, Reticulinasus and Subparmatus. These 6 subgenera form a monophyletic clade that might be placed as new subgenera within the genus Alectorobius, or elevated to genera. Given the substantial differences in biology among these subgenera, we propose that these 6 subgenera be elevated to genera. Thus, we propose to modify the classification scheme of Mans et al. (2019) so that the subfamily Argasinae now has six genera, Alveonasus, Argas (subgenera Argas and Persicargas), Navis, Ogadenus, Proknekalia and Secretargas, and the subfamily Ornithodorinae has nine genera, Alectorobius, Antricola (subgenera Antricola and Parantricola), Carios, Chiropterargas, Nothoaspis, Ornithodoros (subgenera Microargas, Ornamentum, Ornithodoros, Pavlovskyella and Theriodoros), Otobius, Reticulinasus and Subparmatus (genera indicated in bold).

Next generation sequencing and RNA-seq characterization of adipose tissue in the Nile crocodile (Crocodylus niloticus) in South Africa: Possible mechanism(s) of pathogenesis and pathophysiology of pansteatitis.

BACKGROUND: Concerted efforts to identify the pathogenesis and mechanism(s) involved in pansteatitis, (a generalized inflammation of the adipose tissue), that was attributed to the recent crocodile die off in the Olifants River and Loskop Dam in Kruger National Park, Mpumalanga, South Africa have been in the forefront of research in recent time. As part of the efforts, molecular characterization of healthy and pansteatitis adipose tissue was carried out by RNA sequencing (RNA-Seq) using Next Generation Sequencing (NGS) and de novo assembly of the adipose transcriptome, followed by differential gene expression analysis. METHODOLOGY: Healthy adipose tissue consisting of fifty samples was collected from the subcutaneous, visceral, intermuscular adipose tissues and the abdominal fat body of ten 4 years old juvenile crocodiles from a local crocodile farm in Pretoria, South Africa. Ten pansteatitis samples were collected from visceral and intermuscular adipose tissues of five crocodiles that were dying of pansteatitis. RESULTS: Forty-two thousand, two hundred and one (42,201) transcripts were assembled, out of which 37, 835 had previously been characterized. The de novo assembled transcriptome had an N50 (average sequence) of 436 bp, percentage GC content of 43.92, which compared well with previously assembled transcripts in the saltwater crocodile. Seventy genes were differentially expressed and upregulated in pansteatitis. These included genes coding for extracellular matrix (ECM) signaling ligands, inflammatory cytokines and tumour necrosis factor alpha (TNFα) receptors, fatty acid synthase and fatty acid binding proteins, peroxisome proliferator-activated receptor gamma (PPARγ), nuclear factor and apoptosis signaling ligands, and mitogen activated protein kinase enzymes among others. Majority (88.6%) of the upregulated genes were found to be involved in hypoxia inducible pathways for activation of NFkß and inflammation, apoptosis, Toll-like receptor pathway and PPARγ. Bicaudal homologous 2 Drosophila gene (BICD2) associated with spinal and lower extremity muscle atrophy was also upregulated in pansteatitis while Sphingosine -1-phosphate phosphatase 2 (SGPP2) involved in Sphingosine -1- phosphate metabolism was downregulated. Futhermore, Doublesex-mab-related transcription factor 1 (DMRT1) responsible for sex gonad development and germ cell differentiation was also downregulated. CONCLUSION: Thus, from the present study, based on differentially expressed genes in pansteatitis, affected Nile crocodiles might have died partly due to their inability to utilize stored triglycerides as a result of inflammation induced insulin resistance, leading to starvation in the midst of plenty. Affected animals may have also suffered muscular atrophy of the lower extremities and poor fertility.

Argasid and ixodid systematics: Implications for soft tick evolution and systematics, with a new argasid species list.

The systematics of the genera and subgenera within the soft tick family Argasidae is not adequately resolved. Different classification schemes, reflecting diverse schools of scientific thought that elevated or downgraded groups to genera or subgenera, have been proposed. In the most recent classification scheme, Argas and Ornithodoros are paraphyletic and the placement of various subgenera remains uncertain because molecular data are lacking. Thus, reclassification of the Argasidae is required. This will enable an understanding of soft tick systematics within an evolutionary context. This study addressed that knowledge gap using mitochondrial genome and nuclear (18S and 28S ribosomal RNA) sequence data for representatives of the subgenera Alectorobius, Argas, Chiropterargas, Ogadenus, Ornamentum, Ornithodoros, Navis (subgen. nov.), Pavlovskyella, Persicargas, Proknekalia, Reticulinasus and Secretargas, from the Afrotropical, Nearctic and Palearctic regions. Hard tick species (Ixodidae) and a new representative of Nuttalliella namaqua (Nuttalliellidae), were also sequenced with a total of 83 whole mitochondrial genomes, 18S rRNA and 28S rRNA genes generated. The study confirmed the utility of next-generation sequencing to retrieve systematic markers. Paraphyly of Argas and Ornithodoros was resolved by systematic analysis and a new species list is proposed. This corresponds broadly with the morphological cladistic analysis of Klompen and Oliver (1993). Estimation of divergence times using molecular dating allowed dissection of phylogeographic patterns for argasid evolution. The discovery of cryptic species in the subgenera Chiropterargas, Ogadenus and Ornithodoros, suggests that cryptic speciation is common within the Argasidae. Cryptic speciation has implications for past biological studies of soft ticks. These are discussed in particular for the Ornithodoros (Ornithodoros) moubata and Ornithodoros (Ornithodoros) savignyi groups.

Draft Genomes of Two Artocarpus Plants, Jackfruit (A. heterophyllus) and Breadfruit (A. altilis).

Two of the most economically important plants in the Artocarpus genus are jackfruit (A. heterophyllus Lam.) and breadfruit (A. altilis (Parkinson) Fosberg). Both species are long-lived trees that have been cultivated for thousands of years in their native regions. Today they are grown throughout tropical to subtropical areas as an important source of starch and other valuable nutrients. There are hundreds of breadfruit varieties that are native to Oceania, of which the most commonly distributed types are seedless triploids. Jackfruit is likely native to the Western Ghats of India and produces one of the largest tree-borne fruit structures (reaching up to 45 kg). To-date, there is limited genomic information for these two economically important species. Here, we generated 273 Gb and 227 Gb of raw data from jackfruit and breadfruit, respectively. The high-quality reads from jackfruit were assembled into 162,440 scaffolds totaling 982 Mb with 35,858 genes. Similarly, the breadfruit reads were assembled into 180,971 scaffolds totaling 833 Mb with 34,010 genes. A total of 2822 and 2034 expanded gene families were found in jackfruit and breadfruit, respectively, enriched in pathways including starch and sucrose metabolism, photosynthesis, and others. The copy number of several starch synthesis-related genes were found to be increased in jackfruit and breadfruit compared to closely-related species, and the tissue-specific expression might imply their sugar-rich and starch-rich characteristics. Overall, the publication of high-quality genomes for jackfruit and breadfruit provides information about their specific composition and the underlying genes involved in sugar and starch metabolism.

A draft genome sequence of the rose black spot fungus Diplocarpon rosae reveals a high degree of genome duplication.

BACKGROUND: Black spot is one of the most severe and damaging diseases of garden roses. We present the draft genome sequence of its causative agent Diplocarpon rosae as a working tool to generate molecular markers and to analyze functional and structural characteristics of this fungus. RESULTS: The isolate DortE4 was sequenced with 191x coverage of different read types which were assembled into 2457 scaffolds. By evidence supported genome annotation with the MAKER pipeline 14,004 gene models were predicted and transcriptomic data indicated that 88.5% of them are expressed during the early stages of infection. Analyses of k-mer distributions resulted in unexpectedly large genome size estimations between 72.5 and 91.4 Mb, which cannot be attributed to its repeat structure and content of transposable elements alone, factors explaining such differences in other fungal genomes. In contrast, different lines of evidences demonstrate that a huge proportion (approximately 80%) of genes are duplicated, which might indicate a whole genome duplication event. By PCR-RFLP analysis of six paralogous gene pairs of BUSCO orthologs, which are expected to be single copy genes, we could show experimentally that the duplication is not due to technical error and that not all isolates tested possess all of the paralogs. CONCLUSIONS: The presented genome sequence is still a fragmented draft but contains almost the complete gene space. Therefore, it provides a useful working tool to study the interaction of D. rosae with the host and the influence of a genome duplication outside of the model yeast in the background of a phytopathogen.

Gene Duplication and Protein Evolution in Tick-Host Interactions.

Ticks modulate their hosts' defense responses by secreting a biopharmacopiea of hundreds to thousands of proteins and bioactive chemicals into the feeding site (tick-host interface). These molecules and their functions evolved over millions of years as ticks adapted to blood-feeding, tick lineages diverged, and host-shifts occurred. The evolution of new proteins with new functions is mainly dependent on gene duplication events. Central questions around this are the rates of gene duplication, when they occurred and how new functions evolve after gene duplication. The current review investigates these questions in the light of tick biology and considers the possibilities of ancient genome duplication, lineage specific expansion events, and the role that positive selection played in the evolution of tick protein function. It contrasts current views in tick biology regarding adaptive evolution with the more general view that neutral evolution may account for the majority of biological innovations observed in ticks.

Draft Whole-Genome Sequence and Annotation of Xenorhabdus griffiniae Strain BMMCB Associated with the South African Entomopathogenic Nematode Steinernema khoisanae Strain BMMCB.

Xenorhabdus griffiniae strain BMMCB (LDNM00000000) belongs to the family Enterobacteriaceae and was isolated from the South African entomopathogenic nematode Steinernema khoisanae strain BMMCB (GenBank accession no. KT027382). Here, we report the draft whole-genome sequence of X. griffinae strain BMMCB with a genome size of 4,183,779 bp and 44.7% G+C content. The NCBI Prokaryotic Automatic Annotation Pipeline (PGAAP) revealed 3,970 genes.

Draft Whole-Genome Sequence of Serratia sp. Strain TEL, Associated with Oscheius sp. TEL-2014 (Nematoda: Rhabditidae) Isolated from a Grassland in South Africa.

Here, we report on the draft genome sequence of Serratia sp. strain TEL, associated with Oscheius sp. TEL-2014 (Nematoda: Rhabditidae, KM492926) isolated from a grassland in Suikerbosrand Nature Reserve near Johannesburg in South Africa. Serratia sp. strain TEL has a genome size of 5,000,541 bp with 4,647 genes and a G+C content of 59.1%.

Draft Whole-Genome Sequence and Annotation of the Entomopathogenic Bacterium Xenorhabdus khoisanae Strain MCB.

We report here the draft genome sequence of Xenorhabdus khoisanae strain MCB, a Gram-negative bacterium and symbiont of a Steinernema entomopathogenic nematode. The genome assembly consists of 266 contigs covering 4.68 Mb. Genome annotation revealed 3,869 protein-coding sequences, with a G+C content of 43.5%.

Draft Whole-Genome Sequence of Xenorhabdus sp. Strain GDc328, Isolated from the Indigenous South African Nematode Host Steinernema khoisanae.

Here, we describe the draft genome sequence of Xenorhabdus sp. GDc328, an endosymbiont of the native South African entomopathogenic nematode host, Steinernema khoisanae. The total genome size of the bacteria is 4.09 Mb. The genome comprises a total of 3,608 genes with a molecular G+C content of 44.64%.

Draft Genome Sequence and Assembly of Photorhabdus heterorhabditis Strain VMG, a Bacterial Symbiont Associated with the Entomopathogenic Nematode Heterorhabditis zealandica.

Here, we report the draft genome sequence of Photorhabdus heterorhabditis strain VMG, a symbiont of the entomopathogenic nematode Heterorhabditis zealandica in South Africa. The draft genome sequence is 4,878,919 bp long and contains 4,023 protein-coding genes. The genome assembly contains 262 contigs with a G+C content of 42.22%.