Weed resistance to control measures, particularly herbicides, is a growing problem in agriculture. In the case of herbicides, resistance is sometimes connected to genetic changes that directly affect the target site of the herbicide. Other cases are less straightforward where resistance arises without such a clear-cut mechanism. Understanding the genetic and gene regulatory mechanisms that may lead to the rapid evolution of resistance in weedy species is critical to securing our food supply. To study this phenomenon, we exposed young Arabidopsis plants to sublethal levels of one of four weed management stressors, glyphosate herbicide, trifloxysulfuron herbicide, mechanical clipping, and shading. To evaluate responses to these stressors we collected data on gene expression and regulation via epigenetic modification (methylation) and small RNA (sRNA). For all of the treatments except shade, the stress was limited in duration, and the plants were allowed to recover until flowering, to identify changes that persist to reproduction. At flowering, DNA for methylation bisulfite sequencing, RNA, and sRNA were extracted from newly formed rosette leaf tissue. Analyzing the individual datasets revealed many differential responses when compared to the untreated control for gene expression, methylation, and sRNA expression. All three measures showed increases in differential abundance that were unique to each stressor, with very little overlap between stressors. Herbicide treatments tended to exhibit the largest number of significant differential responses, with glyphosate treatment most often associated with the greatest differences and contributing to overlap. To evaluate how large datasets from methylation, gene expression, and sRNA analyses could be connected and mined to link regulatory information with changes in gene expression, the information from each dataset and for each gene was united in a single large matrix and mined with classification algorithms. Although our models were able to differentiate patterns in a set of simulated data, the raw datasets were too noisy for the models to consistently identify differentially expressed genes. However, by focusing on responses at a local level, we identified several genes with differential expression, differential sRNA, and differential methylation. While further studies will be needed to determine whether these epigenetic changes truly influence gene expression at these sites, the changes detected at the treatment level could prime the plants for future incidents of stress, including herbicides. / Doctor of Philosophy / Growing resistance to herbicides, particularly glyphosate, is one of the many problems facing agriculture. The rapid rise of resistance across herbicide classes has caused some to wonder if there is a mechanism of adaptation that does not involve mutations. Epigenetics is the study of changes in the phenotype that cannot be attributed to changes in the genotype. Typically, studies revolve around two features of the chromosomes: cytosine methylation and histone modifications. The former can influence how proteins interact with DNA, and the latter can influence protein access to DNA. Both can affect each other in self-reinforcing loops. They can affect gene expression, and DNA methylation can be directed by small RNA (sRNA), which can also influence gene expression through other pathways. To study these processes and their role in abiotic stress response, we aimed to analyze sRNA, RNA, and DNA from Arabidopsis thaliana plants under stress. The stresses applied were sublethal doses of the herbicides, glyphosate and trifloxysulfuron, as well as mechanical clipping and shade to represent other weed management stressors. The focus of the project was to analyze these responses individually and together to find epigenetic responses to stresses routinely encountered by weeds. We tested RNA for gene expression changes under our stress conditions and identified many, including some pertaining to DNA methylation regulation. The herbicide treatments were associated with upregulated defense genes and downregulated growth genes. Shade treated plants had many downregulated defense and other stress response genes. We also detected differential methylation and sRNA responses when compared to the control plants. Changes to methylation and sRNA only accounted for about 20% of the variation in gene expression. While attempting to link the epigenetic process of methylation to gene expression, we connected all the data sets and developed computer programs to try to make correlations. While these methods worked on a simulated dataset, we did not detect broad patterns of changes to epigenetic pathways that correlated strongly with gene expression in our experiment's data. There are many factors that can influence gene expression that could create noise that would hinder the algorithms' abilities to detect differentially expressed genes. This does not, however, rule out the possibility of epigenetic influence on gene expression in local contexts. Through scoring the traits of individual genes, we found several that interest us for future studies.

Development of fungicide resistance in fungal and oomycete pathogens is a serious problem in grape production. Quinoxyfen is a fungicide widely used against grape powdery mildew (Erysiphe necator). In 2013, E. necator isolates with reduced quinoxyfen sensitivity (designated as quinoxyfen lab resistance or QLR) were detected in Virginia. Field trials were conducted in 2014, 2015, and 2016 at the affected vineyard to determine to what extent quinoxyfen might still contribute to disease control. Powdery mildew control by quinoxyfen was good, similar to, or only slightly less, than that provided by myclobutanil and boscalid in all three years. The frequency of QLR in vines not treated with quinoxyfen declined only slowly over the three years, from 65% to 46%. Information about the mode of action of quinoxyfen is limited; previous research suggests that quinoxyfen interferes with the signal transduction process. We profiled the transcriptomes of QLR and sensitive isolates in response to quinoxyfen treatment, providing support for this hypothesis. Additional transcriptional targets of quinoxyfen were revealed to be involved in the positive regulation of the MAPK signaling cascade, pathogenesis, and sporulation activity. Grape downy mildew (Plasmopara viticola), another important grape pathogen, is commonly controlled by phosphite fungicides. A field trial and laboratory bioassays were conducted to determine whether P. viticola isolates from vineyards with suspected control failures showed reduced sensitivity against phosphite fungicides. Prophyt applied at 14-day intervals under high disease pressure provided poor downy mildew control in the field. Next-generation sequencing technologies were utilized to identify 391,930 single nucleotide polymorphisms (SNPs) and generated a draft P. viticola genome assembly at ~130 megabase (Mb). Finally, field isolates of P. viticola collected from a Virginia vineyard with suspected mandipropamid control failure were bioassayed. The EC50 values of the isolates were >240 μg.ml-1 for mandipropamid, well above the field rate. The PvCesA3 gene of two resistant isolates was sequenced revealing that these isolates had a GGC-to-AGC substitution at codon 1105, the same mutation that has been found associated with CAA resistance elsewhere. / PHD

No / Whole genome sequencing (WGS) of Mycobacterium tuberculosis has rapidly progressed from a research tool to a clinical application for the diagnosis and management of tuberculosis and in public health surveillance. This development has been facilitated by drastic drops in cost, advances in technology and concerted efforts to translate sequencing data into actionable information. There is, however, a risk that, in the absence of a consensus and international standards, the widespread use of WGS technology may result in data and processes that lack harmonization, comparability and validation. In this Review, we outline the current landscape of WGS pipelines and applications, and set out best practices for M. tuberculosis WGS, including standards for bioinformatics pipelines, curated repositories of resistance-causing variants, phylogenetic analyses, quality control and standardized reporting. / European Research Council grant (INTERRUPTB; no. 311725), European Research Council grant (TB-ACCELERATE; no. 638553), Foundation for Innovative New Diagnostics, German Center for Infection Research (DZIF), Deutsche Forschungsgemeinschaft (German Research Foundation) under Germany’s Excellence Strategy (EXC 22167–390884018), FWO Odysseus G0F8316N, US National Institutes of Health BD2K K01 (MRF ES026835), Agence Nationale de la Recherche (ANR-16-CD35-0009)

Yes / Background: Mycobacterium tuberculosis rapid diagnostic tests (RDTs) are widely employed in routine laboratories and national surveys for detection of rifampicinresistant (RR)-TB. However, as next-generation sequencing technologies have become more commonplace in research and surveillance programs, RDTs are being increasingly complemented by whole genome sequencing (WGS). While comparison between RDTs is difficult, all RDT results can be derived from WGS data. This can facilitate continuous analysis of RR-TB burden regardless of the data generation technology employed. By converting WGS to RDT results, we enable comparison of data with different formats and sources particularly for low- and middle-income high TB-burden countries that employ different diagnostic algorithms for drug resistance surveys. This allows national TB control programs (NTPs) and epidemiologists to utilize all available data in the setting for improved RR-TB surveillance. Methods: We developed the Python-based MycTB Genome to Test (MTBGT) tool that transforms WGS-derived data into laboratory-validated results of the primary RDTs—Xpert MTB/RIF, XpertMTB/RIF Ultra, GenoType MDRTBplus v2.0, and GenoscholarNTM+MDRTB II. The tool was validated through RDT results of RR-TB strains with diverse resistance patterns and geographic origins and applied on routine-derived WGS data. Results: The MTBGT tool correctly transformed the single nucleotide polymorphism (SNP) data into the RDT results and generated tabulated frequencies of the RDT probes as well as rifampicin-susceptible cases. The tool supplemented the RDT probe reactions output with the RR-conferring mutation based on identified SNPs. The MTBGT tool facilitated continuous analysis of RR-TB and Xpert probe reactions from different platforms and collection periods in Rwanda. Conclusion: Overall, the MTBGT tool allows low- and middle-income countries to make sense of the increasingly generated WGS in light of the readily available RDT. / Erasmus Mundus Joint Doctorate Fellowship grant 2016- 1346.

Preimplantation genetic diagnosis (PGD) refers to the testing of embryos produced through in vitro fertilization (IVF) in order to identify those unaffected by a specific genetic disorder or chromosomal abnormality. In this study, different methodologies were examined and developed for performance of PGD. Investigation of various whole genome amplification (WGA) methods identified multiple displacement amplification as a reliable method for genotyping single cells. Furthermore, this technology was shown to be compatible with subsequent analysis using single nucleotide polymorphism (SNP) microarrays. Compared to conventional methods used in this study to perform single cell diagnosis (e.g. multiplex PCR), WGA techniques were found to be advantageous since they streamline the development of PGD protocols for couples at high risk of transmitting an inherited disorder and simultaneously offer the possibility of comprehensive chromosome screening (CCS). This study also aimed to develop a widely applicable protocol for accurate typing of the human leukocyte antigen (HLA) region with the purpose of identifying embryos that will be HLA-identical to an existing sibling affected by a disorder that requires haematopoietic stem cell transplantation. Additionally, a novel microarray platform was developed that, apart from accurate CCS, was capable of reliably determining the relative quantity of mitochondrial DNA in polar bodies removed from oocytes and single cells biopsied from embryos. Mitochondria are known to play an important role in oogenesis and preimplantation embryogenesis and their measurement may therefore be of clinical relevance. Moreover, real-time PCR was used for development of protocols for CCS, DNA fingerprinting of sperm samples and embryos and the relative quantitation of telomere length in embryos (since shortened telomeres might be associated with reduced viability). As well as considering the role of genetics in terms of oocyte and embryo viability assessment and the diagnosis of inherited genetic disorders, attention was given to a specific gene (Phospholipase C zeta) of relevance to male infertility. A novel mutation affecting the function of the resulting protein was discovered highlighting the growing importance of DNA sequence variants in the diagnosis and treatment of infertility.

Despite the global licensure of two live-attenuated rotavirus vaccines, Rotarix® and RotaTeq®, rotavirus remains the major cause of severe dehydrating diarrhoea in young mammals and the need for further development of additional rotavirus vaccines, especially vaccines effective against regional strains in developing country settings, is increasing. The design and formulation of new effective multivalent rotavirus vaccines is complicated by the wide rotavirus strain diversity. Novel rotavirus strains emerge periodically due to the propensity of rotaviruses to evolve using mechanisms such as point mutation, genome segment reassortment, genome segment recombination and interspecies transmission. Mutations occurring within the primer binding regions targeted by the current commonly employed sequence-dependent genotyping techniques lead to difficulties in genotyping novel mutant rotavirus strains. Therefore, use of sequence-independent techniques coupled with online rotavirus genotyping tools will help to understand the complete epidemiology of the circulating strains which, in turn, is vital for developing intervention measures such as vaccine and anti-viral therapies. In this study, sequence-independent cDNA synthesis that uses a single set of oligonucleotides that do not require prior sequence knowledge of the rotavirus strains, 454® pyrosequencing, and an online rotavirus genotyping tool, RotaC, were used to swiftly characterise the whole genome of rotaviruses. The robustness of this approach was demonstrated in characterising the complete genetic constellations and evolutionary origin of selected human rotavirus strains that emerged in the past two decades worldwide, human rotavirus strains frequently detected in Africa, and the whole genomes of some common strains frequently detected in bovine species. Most of the characterised strains emerged either through intra- or interspecies genome segment reassortment processes. The methods used in this study also allowed determination of the whole consensus genome sequence of multiple rotavirus variants present in a single stool sample and the elucidation of the evolutionary mechanisms that explained their origin. The 454® pyrosequence-generated data revealed evidence of intergenotype rotavirus genome segment recombination between the genome segments 6 (VP6), 8 (NSP2) and 10 (NSP4) of Wa-like and DS-1-like origin. The use of next generation sequencing technology combined with sequence-independent amplification of the rotavirus genomes allowed the determination of the consensus nucleotide sequence for each of the genome segments of the selected study strains directly from stool sample. The consensus nucleotide sequences of the genome segments encoding VP2, VP4, VP6 and VP7 of some of the study strains were codon optimised for insect cell expression and used to generate recombinant baculoviruses. The Bac-to-Bac baculovirus expression system was used to generate chimaeric rotavirus virus-like particles (RV-VLPs). These chimaeric RV-VLPs contained inner capsids (VP2 and VP6) derived from a South African RVA/Humanwt/ ZAF/GR10924/1999/G9P[6] strain, on to which outer capsid layer proteins composed of various combinations of VP4 and VP7 were assembled. The outer capsid proteins were derived from the dsRNA of G2, G8, G9 or G12 strains associated with either P[4], P[6] or P[8] genotypes that were directly extracted from human stool faecal specimens. The structures of these chimaeric RV-VLPs were morphologically evaluated using transmission electron microscopy (TEM). Based on the size and morphology of the particles, doublelayered (dRV-VLPs) and triple-layered RV-VLPs (tRV-VLPs) were produced. Recombinant rotavirus proteins readily assembled into dRV-VLPs, whereas approximately 10 – 30% of the assembled RV-VLPs from insect expressed recombinant VP2/6/7/4 were chimaeric tRVVLPs. These RV-VLPs will be evaluated in future animal studies as potential non-live rotavirus vaccine candidates. The novel approach of producing RV-VLPs introduced in this study, namely by using the consensus nucleotide sequence derived from dsRNA extracted directly from clinical specimens, should speed up vaccine research and development by bypassing the need to adapt the viruses to tissue culture and circumventing some other problems associated with cell culture adaptation as well. Thus, it is now possible to generate RV-VLPs for evaluation as non-live vaccine candidates for any human or animal field rotavirus strain. / Thesis (PhD (Biochemistry))--North-West University, Potchefstroom Campus, 2012

This work is part of a project focused on the study of the variability of human cytomegalovirus (HCMV) among clinical isolates with the aim to map the geographical distribution of HCMV genotypes, reveal the relationships between genotypes and the severity of HCMV-associated diseases, and identify regions in the HCMV genome with a potential for use as diagnostic and therapeutic targets. Attention was paid to the development of the methodology for the preparation of the material for next-generation sequencing (NGS) from HCMV clinical isolates and evaluation of the obtained sequencing data. Blood and urine samples collected from hematopoietic stem cell transplantat recipients and congenitally infected children were analyzed. Samples suitable for NGS were sequenced by the Illumina platform and sequences were created by de novo assembly followed by mapping assembly. Urine samples in comparison to blood samples had higher yield of material for NGS. Of the samples positive for HCMV DNA (7 of 50) after amplification in the cell cultures, only one sample had high purity of the viral DNA (98%) while six samples had purity of less than 7%. The sample containing 98% of the viral DNA was fully sequenced and the sequence was compared to the sequences of other clinical isolates from Belgium in 11 polymorphic...

Avec plus de 1000 génomes complets disponibles (la grande majorité venant de bactéries), les analyses comparatives de génomes deviennent indispensables pour leurs annotations fonctionnelles, ainsi que pour la compréhension de leur structure et leur évolution, et s'appliquent par exemple en phylogénomique ou au design des vaccins. L'une des approches de plus utilisées pour comparer des génomes est l'alignement de leurs séquences d'ADN, i.e. alignement de génomes complets, c'est à dire identifier les régions de similarité en s'affranchissant de toute annotation. Malgré des améliorations significatives durant les dernières années, des outils performants pour cette approche ainsi que des méthodes pour l'estimation de la qualité des résultats qu'elle produit, en particulier sur les génomes bactériens, restent encore à développer. Outre leurs grandes tailles qui rendent les solutions classiques basées sur la programmation dynamique inutilisables, l'alignement de génomes complets posent des difficultés supplémentaires dues à leur évolution particulière, comprenant: la divergence, qui estompe les similarités entre les séquences, le réordonnancent des portions génomiques (réarrangements), ou l'acquisition de matériel génétique extérieur, qui produit des régions non alignables entres les séquences, e.g. transfert horizontal des gènes, phages. En conséquence, les solutions pour l'alignement de génomes sont des heuristiques, dont la plus commune est appelée stratégie basée sur des ancres. Cette stratégie commence par identifier un ensemble initial de régions de similarité (phase 1). Ensuite une phase de chaînage sélectionne un sous-ensemble (non-chevauchantes et généralement colinéaires) de ces similarités de poids maximal, nommées ancres (phase 2). Les phases 1 et 2 sont appliquées de manière récursive sur les régions encore non-alignées (phase 3). La dernière phase consiste en l'application systématique des outils d'alignement classiques sur toutes les régions courtes qui n'ont pas encore été alignées. Cette thèse adresse plusieurs problèmes liés à l'alignement de génomes complets dont: l'évaluation de la qualité des résultats produits par les outils d'alignement et l'amélioration de la stratégie basée sur des ancres. Premièrement, nous avons créé un protocole pour évaluer la qualité des résultats d'alignement, contenant des mesures de calcul quantitatives et qualitatives, dont certaines basées sur des connaissances biologiques. Une analyse de la qualité des alignements produits par deux des principaux outils existants sur des paires de génomes bactériens intra-espèces révèle leurs limitations: des similarités non détectées et des portions d'alignement incorrectes. À partir de ces résultats, qui suggèrent un manque de sensibilité et spécificité, nous proposons un nouvel outil pour l'alignement deux à deux de génomes complets, YOC, qui implémente une version simplifiée de la stratégie basée sur des ancres, contenant seulement deux phases. Dans la phase 1, YOC améliore la sensibilité en utilisant comme ancres, pour la première fois dans cette stratégie, des similarités locales basées sur des graines espacées, capables de détecter des similarités plus longues dans des régions plus divergentes. Cette phase est suivie par une méthode de chainage adaptée aux similarités locales, un nouveau type de chaînage colinéaire, permettant des chevauchements proportionnels. Nous avons donné une formulation de ce nouveau problème et réalisé un premier algorithme. L'algorithme, qui adopte une approche de programmation dynamique basée sur le paradigme de la ``sweep-line'', donne une solution optimale, i.e. est exacte, et s'exécute en temps quadratique. Nous avons montré que cet algorithme, comparé au chainage colinéaire classique, améliore les résultats sur des génomes bactériens, tout en restant aussi efficace en pratique. / With more than 1000 complete genomes available (among which, the vast majority come from bacteria), comparative genomic analysis become essential for the functional annotation of genomes, the understanding of their structure and evolution and have applications in phylogenomics or vaccine design. One of the main approaches for comparing genomes is by aligning their DNA sequences, i.e. whole genome alignment (WGA), which means identifying the similarity regions without any prior annotation knowledge. Despite the significant improvements during the last years, reliable tools for WGA and methodology for estimating its quality, in particular for bacterial genomes, still need to be designed. Besides their extremely large lengths that make classical dynamic programming alignment methods unsuitable, aligning whole genomes involves several additional difficulties, due to the mechanisms through which genomes evolve: the divergence, which let sequence sim ilarity vanish over time, the reordering of genomic segments (rearrangements), or the acquisition of external genetic material generating regions that are unalignable between sequences, e.g. horizontal gene transfer, phages. Therefore, whole genome alignment tools implement heuristics, among which the most common is the anchor based strategy. It starts by detecting an initial set of similarity regions (phase 1), and, through a chaining phase (phase 2), selects a non-overlapping maximum-weighted, usually collinear, subset of those similarities, called anchors. Phases 1 and 2 are recursively applied on yet unaligned regions (phase 3). The last phase (phase 4) consists in systematically applying classical alignment tools to all short regions still left unaligned.This thesis addresses several problems related to whole genome alignment: the evaluation of the quality of results given by WGA tools and the improvement of the classical anchor based strategy. We first designed a protocol for evaluating the quality of alignment results, based on both computational and biological measures. An evaluation of the results given by two state of the art WGA tools on pairs of intra-species bacterial genomes revealed their shortcomings: the failure of detecting some of the similarities between sequences and the misalignment of some regions. Based on these results, which imply a lack in both sensitivity and specificity, we propose a novel, pairwise whole genome alignment tool, YOC, implementing a simplified two-phase version of the anchor strategy. In phase 1, YOC improves sensitivity by using as anchors, for the first time, local similarities based on spaced seeds that are capable of detecting larger similarity regions in divergent sequences. This ph ase is followed by a chaining method adapted to local similarities, a novel type of collinear chaining, allowing for proportional overlaps. We give a formulation for this novel problem and provide the first algorithm for it. The algorithm, implementing a dynamic programming approach based on the sweep-line paradigm, is exact and runs in quadratic time. We show that, compared to classical collinear chaining, chaining with overlaps improves on real bacterial data, while remaining almost as efficient in practice. Our novel tool, YOC, is evaluated together with other four WGA tools on a dataset composed of 694 pairs of intra-species bacterial genomes. The results show that YOC improves on divergent cases by detecting more distant similarities and by avoiding misaligned regions. In conclusion, YOC should be easier to apply automatically and systematically to incoming genomes, for it does not require a post-filtering step to detect misalignment and is less complex to calibrate.

The diploid amphibian Xenopus tropicalis represents a significant model organism for studies of early development, genes function and evolution. Such techniques as gynogenesis, injection of morpholino antisense oligonucleotide into fertilized eggs or transgenesis were established. In the recent ten years, many efforts have been made to complete the sequence information. X. tropicalis genome has been sequenced but the completion of its assembly only on the basis of sequence data has been impossible. Therefore, our first work was focused on one of approaches for a genome completing- genetic mapping. First of all, the genetic map of Xenopus tropicalis was established pursuant linkage and physical positions of markers. Since the map contained gaps, we developed a new method for genetic mapping based on the next generation sequencing of laser microdissected arm. Using Illumina next generation sequencing of fifteen copies of a short arm of chromosome 7, we obtained new insights into its genome by localizing previously unmapped genes and scaffolds as well as recognizing mislocalized portions of the genome assembly. This was the first time laser microdissection and sequencing of specific chromosomal regions has been used for the purpose of genome mapping. These data were also used in the evolution study of...

O objetivo deste trabalho foi investigar a presença de bactérias produtoras de β-lactamases adquiridas na microbiota Gram-negativa comensal de humanos e animais domésticos em uma comunidade remota na região da Floresta Amazônica. De março a julho de 2013 foram coletadas amostras de fezes de indivíduos atendidos e funcionários de um centro assistencial em saúde restrito a comunidades indígenas e de swab retal de animais de companhia da comunidade. Nas amostras de humanos foram detectados isolados de Escherichia coli, Klebsiella pneumoniae, Enterobacter aerogenes, Enterobacter kobei e Morganella morganii, carregando os genes blaTEM-1, blaCTX-M-15, blaCTX-M-14, blaCTX-M-8 e blaGES-5. Nas amostras de animais foram detectados apenas isolados de E. coli carregando os genes blaTEM-1, blaCTX-M-14, blaCTX-M-2 e blaCTX-M-8. Foi observada a relação clonal entre isolados de E. coli de origem humana e de origem animal. Estes resultados demonstram a disseminação de um problema endêmico em áreas urbanas para uma comunidade, em teoria, com baixa exposição a antibacterianos. / The aim of this study was to investigate the presence of acquired β-lactamase in the commensal Gram-negative microbiota of humans and domestic animals of a remote community in the Amazon Forest region. From March to July 2013 stool samples were collected from individuals attended in a health care center restricted to indigenous communities and from the local staff, and rectal swab samples were collected from companion animals in the community. In the human samples were detected Escherichia coli, Klebsiella pneumoniae, Enterobacter aerogenes, Enterobacter kobei and Morganella morganii isolates harboring blaTEM-1, blaCTX-M-15, blaCTX-M-14, blaCTX-M-8 e blaGES-5 genes. In the animal samples only E. coli strains harboring blaTEM-1, blaCTX-M-14, blaCTX-M-2 and blaCTX-M-8 were detected. The clonal relatedness between E. coli strains from human and animal samples was observed. These results demonstrate the dissemination of an urban endemic problem to a community, in theory, with low antimicrobial exposure.