Sialotranscriptomics of the brown ear ticks, Rhipicephalus appendiculatus Neumann, 1901 and R. Zambeziensis Walker, Norval and Corwin, 1981, vectors of Corridor disease

Text in EnglishCorridor disease is an economically important tick-borne disease of cattle in southern Africa. The disease is caused by Theileria parva and transmitted by the vectors, Rhipicephalus appendiculatus and R. zambeziensis. There is currently no vaccine to protect cattle against T. parva that is permitted in South Africa. To develop recombinant anti-tick vaccines against Corridor disease, comprehensive databases of genes expressed in the tick’s salivary glands are required. Therefore, in Chapters 2 and 3, mRNA from the salivary glands of R. appendiculatus and R. zambeziensis was sequenced and assembled using next generation sequencing technologies. Respectively, 12 761 and 13 584 non-redundant protein sequences were predicted from the sialotranscriptomes of R. appendiculatus and R. zambeziensis and uploaded to public sequence domains. This greatly expanded the number of sequences available for the two vectors, which will be invaluable resources for the selection of vaccine candidates in future. Further, in Chapter 3, differential gene expression analysis in R. zambeziensis revealed dynamic expression of secretory protein transcripts during feeding, suggestive of stringent transcriptional regulation of these proteins. Knowledge of these intricate expression profiles will further assist vaccine development in future. In Chapter 4, comparative sialotranscriptomic analyses were performed between R. appendiculatus and R. zambeziensis. The ticks have previously shown varying vector competence for T. parva and this chapter presents the search for correlates of this variance. Phylogenetic analyses were performed using these and other publically available tick transcriptomes, which indicated that R. appendiculatus and R. zambeziensis are closely related but distinct species. However, significant expression differences were observed between the two ticks, specifically of genes involved in tick immunity or pathogen transmission, signifying potential bioinformatic signatures of vector competence. Furthermore, nearly four thousand putative long non-coding RNAs (lncRNAs) were predicted in each of the two ticks. A large number of these showed differential expression and suggested a potential transcriptional regulatory function of lncRNA in tick blood feeding. LncRNAs are completely unexplored in ticks. Finally, in Chapter 5, concluding remarks are given on the potential impact the R. appendiculatus and R. zambeziensis sialotranscriptomes may have on future vaccine developments and some future research endeavours are discussed.Life and Consumer SciencesPh. D. (Life Sciences

Allelic diversity in the CAD2 and LIM1 lignin biosynthetic genes of Eucalyptus grandis Hill ex Maiden and E. smithii R.T. Baker

Lignin is a highly abundant aromatic biopolymer deposited during the final stages of secondary cell wall formation in plants and it constitutes a substantial proportion of the dry weight of woody plant stems. Lignin contributes structural support to xylem cell walls and hydrophobisity to water-conducting vessels and forms a defence mechanism against pathogen invasion. Although being an essential part of normal plant cell development, lignin content and composition are targets for tree improvement, because residual lignin in paper pulp has negative effects on paper quality and lignin therefore has to be removed using treatments that are expensive and often detrimental to the environment. At present, little is known about the amount of allelic diversity in lignin biosynthetic genes and whether such diversity may be associated with variation in lignin content and composition. However, the identification of alleles associated with desirable lignin phenotypes is dependent on a detailed understanding of the molecular evolution and population genetics of these genes. This M.Sc. study was aimed at analysing nucleotide and allelic diversity in two lignin biosynthetic genes of Eucalyptus trees. Additionally, the study aimed to develop single nucleotide polymorphism (SNP) markers that could be used to assay allelic diversity for these genes in populations of two target species, E. grandis and E. smithii. Orthologues of the tobacco LIM-domain1 (NtLIM1) transcription factor gene involved in the regulation of lignin biosynthesis were isolated from E. grandis and E. smithii. Approximately 3 kb of genomic sequence including the promoter and full-length gene regions were isolated for the two orthologues, respectively labeled EgrLIM1 and EsLIM1. The predicted amino acid sequences of EgrLIM1 and EsLIM1 were 99.4% identical to each other and indicated that LIM1 is a small protein of only 188 residues in eucalypt trees and has a predicted molecular weight of 21.0 kDa. Quantitative, real-time RT-PCR analysis confirmed the expression of LIM1 in wood-forming tissues undergoing lignification. Ten putative cis-regulatory elements were observed in the promoter regions of EgrLIM1 and EsLIM1including a GA-dinucleotide microsatellite that appears to be specific to LIM1 promoters of Eucalyptus tree species. The full-length LIM1 gene sequences could subsequently be used in the assessment of nucleotide and allelic diversity, together with the full-length CAD2 sequences that were already available in the public domain. The level of nucleotide and allelic diversity and the distribution and decay of linkage disequilibrium (LD) were surveyed in 5’ and 3’ derived gene fragments of CAD2 and LIM1 obtained from 20 E. grandis and 20 E. smithii individuals. Each gene displayed a unique genetic diversity profile, but for the most part, nucleotide diversity (π) was estimated at approximately 0.0010 except for the E. grandis LIM1 gene where π lower than 0.0040 was observed. Generally, except for the high amounts of LD observed in the CAD2 gene of E. grandis (> 2.5 kb), LD decayed within 500 bp. A large number (13 to 45) of SNP sites (defined as single nucleotide changes with minor allele frequencies of at least 0.10 in each species) were observed in each gene of each species. The high SNP density (ranging from one per 45 to one per 155 bp) observed in the two genes facilitated the efficient development of SNP markers to be used in future aspects of LD mapping, association genetics and marker-assisted breeding. The allele sequences obtained for the CAD2 and LIM1 genes were used as templates for the development of SNP marker panels (a series of six or seven SNP markers analysed together) for the analysis (tagging) of SNP haplotype diversity in species-wide reference populations (100 E. grandis and 137E. smithii individuals) of the two species. Each tag SNP was assayed using a single base extension assay and capillary gel electrophoresis. High polymorphism information content (average PIC of 0.836) was observed for the SNP marker panels. Four SNPs in the CAD2 and two in the LIM1 genes were found to be polymorphic in E. grandis and E. smithii (i.e. trans-specific SNPs), suggesting a possible ancestral origin for these polymorphisms. Assessment of candidate gene variation in the genomes of forest trees is of importance to ultimately be able to predict the amount and structure of nucleotide diversity available for the future design of SNP assays at the whole-genome level. Such assays will be useful to study differentiation among tree species and populations, to associate nucleotide polymorphisms with desirable phenotypes and to increase the efficiency of tree improvement approaches.Dissertation (MSc (Genetics))--University of Pretoria, 2009.Geneticsunrestricte

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.The Economic Competitive Support Programme (30/01/V010) and the National Research Foundation (NRF) Incentive Funding (IFR2011032400016) for Rated Researchers (NRF-Mans). MdC was supported by an NRF/Department of Science and Technology—Professional Development Program (NRF/DST-PDP) studentship.http://www.frontiersin.org/Cellular_and_Infection_Microbiologyam2017Veterinary Tropical Disease

Next-generation sequencing as means to retrieve tick systematic markers, with the focus on Nuttalliella namaqua (Ixodoidea : Nuttalliellidae)

Nuclear ribosomal RNA (18S and 28S rRNA) and mitochondrial genomes are commonly used in tick systematics. The ability to retrieve these markers using next-generation sequencing was investigated using the tick Nuttalliella namaqua. Issues related to nuclear markers may be resolved using this approach, notably, the monotypic status of N. namaqua and its basal relationship to other tick families. Four different Illumina datasets (∼55 million, 100 bp reads each) were generated from a single tick specimen and assembled to give 350k-390k contigs. A genome size of ∼1 Gbp was estimated with low levels of repetitive elements. Contigs (>1000 bp, >50-fold coverage) present in most assemblies (n=69), included host-derived 18S and 28S rRNA, tick and host-derived transposable elements, full-length tick 18S and 28 rRNA, the mitochondrial genome in single contig assemblies and the histone cassette. Coverage for the nuclear rRNA genes was above 1000-fold confirming previous sequencing errors in the 18S rRNA gene, thereby maintaining the monotypic status of this tick. Nuclear markers for the soft tick Argas africolumbae were also retrieved from next-generation data. Phylogenetic analysis of a concatenated 18S-28S rRNA dataset supported the grouping of N. namaqua at the base of the tick tree and the two main tick families in separate clades. This study confirmed the monotypic status of N. namaqua and its basal relationship to other tick families. Next-generation sequencing of genomic material to retrieve high quality nuclear and mitochondrial systematic markers for ticks is viable and may resolve issues around conventional sequencing errors when comparing closely related tick species.The Joy Liebenberg Trust (21/19/JT02) allocated to BM, incentive funding for rated researchers grant from the National Research Foundation of South Africa (NRF-Mans 76499), and the ARC Tick Vaccine project (30/01/V010).http://www.elsevier.com/locate/ttbdishb2017Veterinary Tropical Disease

The Mitochondrial genomes of Nuttalliella namaqua (Ixodoidea: Nuttalliellidae) and Argas africolumbae (Ixodoidae: Argasidae) : estimation of divergence dates for the major tick lineages and reconstruction of ancestral blood-feeding characters

Ixodida are composed of hard (Ixodidae), soft (Argasidae) and the monotypic Nuttalliellidae (Nuttalliella namaqua) tick families. Nuclear 18S rRNA analysis suggested that N. namaqua was the closest extant relative to the last common ancestral tick lineage. The mitochondrial genomes of N. namaqua and Argas africolumbae were determined using next generation sequencing and de novo assembly to investigate this further. The latter was included since previous estimates on the divergence times of argasids lacked data for this major genus. Mitochondrial gene order for both was identical to that of the Argasidae and Prostriata. Bayesian analysis of the COI, Cytb, ND1, ND2 and ND4 genes confirmed the monophyly of ticks, the basal position of N. namaqua to the other tick families and the accepted systematic relationships of the other tick genera. Molecular clock estimates were derived for the divergence of the major tick lineages and supported previous estimates on the origins of ticks in the Carboniferous. N. namaqua larvae fed successfully on lizards and mice in a prolonged manner similar to many argasids and all ixodids. Excess blood meal-derived water was secreted via the salivary glands, similar to ixodids. We propose that this prolonged larval feeding style eventually gave rise to the long feeding periods that typify the single larval, nymphal and adult stages of ixodid ticks and the associated secretion of water via the salivary glands. Ancestral reconstruction of characters involved in blood-feeding indicates that most of the characteristics unique to either hard or soft tick families were present in the ancestral tick lineage.The Joy Liebenberg Trust (21/19/JT02) allocated to BM and a South African National Research Foundation grant allocated to AL (NRF-Spain).http://www.plosone.orgam2013ab201

Diversity and cis-element architecture of the promoter regions of cellulose synthase genes in Eucalyptus

Lignocellulosic biomass from fast-growing plantation trees is composed of carbohydrate-rich materials deposited into plant cell walls in a coordinated manner during wood formation. The diversity and evolution of the transcriptional networks regulating this process have not been studied extensively.We investigated patterns of species-level nucleotide diversity in the promoters of cellulose synthase (CesA) genes from different Eucalyptus tree species and assessed the possible roles of DNA sequence polymorphism in the gain or loss of cis-elements harboured within the promoters. Promoter regions of three primary and three secondary cell wall-associated CesA genes were isolated from 13 Eucalyptus species and were analysed for nucleotide and cis-element diversity. Species-level nucleotide diversity (π) ranged from 0.014 to 0.068, and different CesA promoters exhibited distinct patterns of sequence conservation. A set of 22 putative cis-elements were mapped to the CesA promoters using in silico methods. Forty-two percent of the mapped cis-element occurrences contained singleton polymorphisms which resulted in either gain or loss of a ciselement in a particular Eucalyptus species. The promoters of Eucalyptus CesA genes contained regions that are highly conserved at the species (Eucalyptus) and genus (with Arabidopsis and Populus) level, suggesting the presence of regulatory modules imposing functional constraint on such regions. Nucleotide polymorphisms in the CesA promoters more frequently created new cis-element occurrences than disrupted existing cis-element occurrences, a process which may be important for the maintenance and evolution of cellulose gene regulation in plants.Mondi and Sappi, through the Forest Molecular Genetics (FMG) Programme, the Technology and Human Resources for Industry Programme (THRIP)and the National Research Foundation of South Africa (NRF).http://link.springer.com/journal/11295hb201

Ancestral reconstruction of tick lineages

Ancestral reconstruction in its fullest sense aims to describe the complete evolutionary history of a lineage. This depends on accurate phylogenies and an understanding of the key characters of each parental lineage. An attempt is made to delineate our current knowledge with regard to the ancestral reconstruction of the tick (Ixodida) lineage. Tick characters may be assigned to Core of Life, Lineages of Life or Edges of Life phenomena depending on how far back these characters may be assigned in the evolutionary Tree of Life. These include housekeeping genes, sub-cellular systems, heme processing (Core of Life), development, moulting, appendages, nervous and organ systems, homeostasis, respiration (Lineages of Life), specific adaptations to a blood-feeding lifestyle, including the complexities of salivary gland secretions and tick–host interactions (Edges of Life). The phylogenetic relationships of lineages, their origins and importance in ancestral reconstruction are discussed. Uncertainties with respect to systematic relationships, ancestral reconstruction and the challenges faced in comparative transcriptomics (next-generation sequencing approaches) are highlighted. While almost 150 years of information regarding tick biology have been assembled, progress in recent years indicates that we are in the infancy of understanding tick evolution. Even so, broad reconstructions can be made with relation to biological features associated with various lineages. Conservation of characters shared with sister and parent lineages are evident, but appreciable differences are present in the tick lineage indicating modification with descent, as expected for Darwinian evolutionary theory. Many of these differences can be related to the hematophagous lifestyle of ticks.This project was funded by the Joy Liebenberg Trust (21/19/JT02)allocated to BM, incentive funding for rated researchers grant from the National Research Foundation of South Africa (NRF-Mans76499), and the ARC Tick Vaccine project (30/01/V010).http://www.elsevier.com/locate/ttbdis2017-06-30hb2016Veterinary Tropical Disease

De novo assembly and annotation of the salivary gland transcriptome of Rhipicephalus appendiculatus male and female ticks during blood feeding

Tick secretory proteins modulate haemostasis, inflammation and immune responses of the host and are attractive recombinant anti-tick vaccine candidates. Yet, many of the proteins have not been characterised due to the limited sequence availability for ticks and other arthropods for homology-based annotation. To address this limitation, we sequenced the salivary glands of the economically important adult male and female Rhipicephalus appendiculatus ticks during feeding. The quality-filtered Illumina sequencing reads were de novo assembled to generate a R. appendiculatus sialotranscriptome of 21,410 transcripts. A non-redundant set of 12,761 R. appendiculatus proteins was predicted from the transcripts, including 2134 putative secretory and 8237 putative housekeeping proteins. Secretory proteins accounted for most of the expression in the salivary gland transcriptome (63%). Of the secretory protein class, the Glycine-rich superfamily contributed 66% and the Lipocalin family 12% of the transcriptome expression. Differential expression analysis identified 1758 female and 2346 male up-regulated transcripts, suggesting varying blood-feeding mechanisms employed between female and male ticks. The sialotranscriptome assembled in this work, greatly improves on the sequence information available for R. appendiculatus and is a valuable resource for potential future vaccine candidate selection.The Economic Competitive Support Programme (30/01/V010) and Incentive Funding for Rated Researchers (NRF-Mans).http://www.elsevier.com/locate/ttbdis2017-06-30hb2016Veterinary Tropical Disease

The genetic legacy of the expansion of Bantu-speaking peoples in Africa.

The expansion of people speaking Bantu languages is the most dramatic demographic event in Late Holocene Africa and fundamentally reshaped the linguistic, cultural and biological landscape of the continent1-7. With a comprehensive genomic dataset, including newly generated data of modern-day and ancient DNA from previously unsampled regions in Africa, we contribute insights into this expansion that started 6,000-4,000 years ago in western Africa. We genotyped 1,763 participants, including 1,526 Bantu speakers from 147 populations across 14 African countries, and generated whole-genome sequences from 12 Late Iron Age individuals8. We show that genetic diversity amongst Bantu-speaking populations declines with distance from western Africa, with current-day Zambia and the Democratic Republic of Congo as possible crossroads of interaction. Using spatially explicit methods9 and correlating genetic, linguistic and geographical data, we provide cross-disciplinary support for a serial-founder migration model. We further show that Bantu speakers received significant gene flow from local groups in regions they expanded into. Our genetic dataset provides an exhaustive modern-day African comparative dataset for ancient DNA studies10 and will be important to a wide range of disciplines from science and humanities, as well as to the medical sector studying human genetic variation and health in African and African-descendant populations