The saponin bomb: a nucleolar-localized ß-glucosidase hydrolyzes triterpene saponins in Medicago truncatula.

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a ß-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned ß-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Siliplant1 protein precipitates silica in sorghum silica cells.

Silicon is absorbed by plant roots as silicic acid. The acid moves with the transpiration stream to the shoot, and mineralizes as silica. In grasses, leaf epidermal cells called silica cells deposit silica in most of their volume using an unknown biological factor. Using bioinformatics tools, we identified a previously uncharacterized protein in Sorghum bicolor, which we named Siliplant1 (Slp1). Slp1 is a basic protein with seven repeat units rich in proline, lysine, and glutamic acid. We found Slp1 RNA in sorghum immature leaf and immature inflorescence. In leaves, transcription was highest just before the active silicification zone (ASZ). There, Slp1 was localized specifically to developing silica cells, packed inside vesicles and scattered throughout the cytoplasm or near the cell boundary. These vesicles fused with the membrane, releasing their content in the apoplastic space. A short peptide that is repeated five times in Slp1 precipitated silica in vitro at a biologically relevant silicic acid concentration. Transient overexpression of Slp1 in sorghum resulted in ectopic silica deposition in all leaf epidermal cell types. Our results show that Slp1 precipitates silica in sorghum silica cells.

T-DNA-genome junctions form early after infection and are influenced by the chromatin state of the host genome.

Agrobacterium tumefaciens mediated T-DNA integration is a common tool for plant genome manipulation. However, there is controversy regarding whether T-DNA integration is biased towards genes or randomly distributed throughout the genome. In order to address this question, we performed high-throughput mapping of T-DNA-genome junctions obtained in the absence of selection at several time points after infection. T-DNA-genome junctions were detected as early as 6 hours post-infection. T-DNA distribution was apparently uniform throughout the chromosomes, yet local biases toward AT-rich motifs and T-DNA border sequence micro-homology were detected. Analysis of the epigenetic landscape of previously isolated sites of T-DNA integration in Kanamycin-selected transgenic plants showed an association with extremely low methylation and nucleosome occupancy. Conversely, non-selected junctions from this study showed no correlation with methylation and had chromatin marks, such as high nucleosome occupancy and high H3K27me3, that correspond to three-dimensional-interacting heterochromatin islands embedded within euchromatin. Such structures may play a role in capturing and silencing invading T-DNA.

The 'TranSeq' 3'-end sequencing method for high-throughput transcriptomics and gene space refinement in plant genomes.

High-throughput RNA sequencing has proven invaluable not only to explore gene expression but also for both gene prediction and genome annotation. However, RNA sequencing, carried out on tens or even hundreds of samples, requires easy and cost-effective sample preparation methods using minute RNA amounts. Here, we present TranSeq, a high-throughput 3'-end sequencing procedure that requires 10- to 20-fold fewer sequence reads than the current transcriptomics procedures. TranSeq significantly reduces costs and allows a greater increase in size of sample sets analyzed in a single experiment. Moreover, in comparison with other 3'-end sequencing methods reported to date, we demonstrate here the reliability and immediate applicability of TranSeq and show that it not only provides accurate transcriptome profiles but also produces precise expression measurements of specific gene family members possessing high sequence similarity. This is difficult to achieve in standard RNA-seq methods, in which sequence reads cover the entire transcript. Furthermore, mapping TranSeq reads to the reference tomato genome facilitated the annotation of new transcripts improving >45% of the existing gene models. Hence, we anticipate that using TranSeq will boost large-scale transcriptome assays and increase the spatial and temporal resolution of gene expression data, in both model and non-model plant species. Moreover, as already performed for tomato (ITAG3.0; www.solgenomics.net), we strongly advocate its integration into current and future genome annotations.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping.

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Network Modeling Unravels Mechanisms of Crosstalk between Ethylene and Salicylate Signaling in Potato.

To develop novel crop breeding strategies, it is crucial to understand the mechanisms underlying the interaction between plants and their pathogens. Network modeling represents a powerful tool that can unravel properties of complex biological systems. In this study, we aimed to use network modeling to better understand immune signaling in potato (Solanum tuberosum). For this, we first built on a reliable Arabidopsis (Arabidopsis thaliana) immune signaling model, extending it with the information from diverse publicly available resources. Next, we translated the resulting prior knowledge network (20,012 nodes and 70,091 connections) to potato and superimposed it with an ensemble network inferred from time-resolved transcriptomics data for potato. We used different network modeling approaches to generate specific hypotheses of potato immune signaling mechanisms. An interesting finding was the identification of a string of molecular events illuminating the ethylene pathway modulation of the salicylic acid pathway through Nonexpressor of PR Genes1 gene expression. Functional validations confirmed this modulation, thus supporting the potential of our integrative network modeling approach for unraveling molecular mechanisms in complex systems. In addition, this approach can ultimately result in improved breeding strategies for potato and other sensitive crops.

MYB107 and MYB9 Homologs Regulate Suberin Deposition in Angiosperms.

Suberin, a polymer composed of both aliphatic and aromatic domains, is deposited as a rough matrix upon plant surface damage and during normal growth in the root endodermis, bark, specialized organs (e.g., potato [Solanum tuberosum] tubers), and seed coats. To identify genes associated with the developmental control of suberin deposition, we investigated the chemical composition and transcriptomes of suberized tomato (Solanum lycopersicum) and russet apple (Malus x domestica) fruit surfaces. Consequently, a gene expression signature for suberin polymer assembly was revealed that is highly conserved in angiosperms. Seed permeability assays of knockout mutants corresponding to signature genes revealed regulatory proteins (i.e., AtMYB9 and AtMYB107) required for suberin assembly in the Arabidopsis thaliana seed coat. Seeds of myb107 and myb9 Arabidopsis mutants displayed a significant reduction in suberin monomers and altered levels of other seed coat-associated metabolites. They also exhibited increased permeability, and lower germination capacities under osmotic and salt stress. AtMYB9 and AtMYB107 appear to synchronize the transcriptional induction of aliphatic and aromatic monomer biosynthesis and transport and suberin polymerization in the seed outer integument layer. Collectively, our findings establish a regulatory system controlling developmentally deposited suberin, which likely differs from the one of stress-induced polymer assembly recognized to date.

Characterization of a New Pink-Fruited Tomato Mutant Results in the Identification of a Null Allele of the SlMYB12 Transcription Factor.

The identification and characterization of new tomato (Solanum lycopersicum) mutants affected in fruit pigmentation and nutritional content can provide valuable insights into the underlying biology, as well as a source of new alleles for breeding programs. To date, all characterized pink-pigmented tomato fruit mutants appear to result from low SlMYB12 transcript levels in the fruit skin. Two new mutant lines displaying a pink fruit phenotype (pf1 and pf2) were characterized in this study. In the pf mutants, SlMYB12 transcripts accumulated to wild-type levels but exhibited the same truncation, which resulted in the absence of the essential MYB activation domain coding region. Allelism and complementation tests revealed that both pf mutants were allelic to the y locus and showed the same recessive null allele in homozygosis: Δy A set of molecular and metabolic effects, reminiscent of those observed in the Arabidopsis (Arabidopsis thaliana) myb11 myb12 myb111 triple mutant, were found in the tomato Δy mutants. To our knowledge, these have not been described previously, and our data support the idea of their being null mutants, in contrast to previously described transcriptional hypomorphic pink fruit lines. We detected a reduction in the expression of several flavonol glycosides and some associated glycosyl transferases. Transcriptome analysis further revealed that the effects of the pf mutations extended beyond the flavonoid pathway into the interface between primary and secondary metabolism. Finally, screening for Myb-binding sites in the candidate gene promoter sequences revealed that 141 of the 152 co-down-regulated genes may be direct targets of SlMYB12 regulation.

Rewiring Host Lipid Metabolism by Large Viruses Determines the Fate of Emiliania huxleyi, a Bloom-Forming Alga in the Ocean.

Marine viruses are major ecological and evolutionary drivers of microbial food webs regulating the fate of carbon in the ocean. We combined transcriptomic and metabolomic analyses to explore the cellular pathways mediating the interaction between the bloom-forming coccolithophore Emiliania huxleyi and its specific coccolithoviruses (E. huxleyi virus [EhV]). We show that EhV induces profound transcriptome remodeling targeted toward fatty acid synthesis to support viral assembly. A metabolic shift toward production of viral-derived sphingolipids was detected during infection and coincided with downregulation of host de novo sphingolipid genes and induction of the viral-encoded homologous pathway. The depletion of host-specific sterols during lytic infection and their detection in purified virions revealed their novel role in viral life cycle. We identify an essential function of the mevalonate-isoprenoid branch of sterol biosynthesis during infection and propose its downregulation as an antiviral mechanism. We demonstrate how viral replication depends on the hijacking of host lipid metabolism during the chemical "arms race" in the ocean.

Mapping the diatom redox-sensitive proteome provides insight into response to nitrogen stress in the marine environment.

Diatoms are ubiquitous marine photosynthetic eukaryotes responsible for approximately 20% of global photosynthesis. Little is known about the redox-based mechanisms that mediate diatom sensing and acclimation to environmental stress. Here we used a quantitative mass spectrometry-based approach to elucidate the redox-sensitive signaling network (redoxome) mediating the response of diatoms to oxidative stress. We quantified the degree of oxidation of 3,845 cysteines in the Phaeodactylum tricornutum proteome and identified approximately 300 redox-sensitive proteins. Intriguingly, we found redox-sensitive thiols in numerous enzymes composing the nitrogen assimilation pathway and the recently discovered diatom urea cycle. In agreement with this finding, the flux from nitrate into glutamine and glutamate, measured by the incorporation of (15)N, was strongly inhibited under oxidative stress conditions. Furthermore, by targeting the redox-sensitive GFP sensor to various subcellular localizations, we mapped organelle-specific oxidation patterns in response to variations in nitrogen quota and quality. We propose that redox regulation of nitrogen metabolism allows rapid metabolic plasticity to ensure cellular homeostasis, and thus is essential for the ecological success of diatoms in the marine ecosystem.

Combined correlation-based network and mQTL analyses efficiently identified loci for branched-chain amino acid, serine to threonine, and proline metabolism in tomato seeds.

Correlation-based network analysis (CNA) of the metabolic profiles of seeds of a tomato introgression line mapping population revealed a clique of proteinogenic amino acids: Gly, Ile, Pro, Ser, Thr, and Val. Correlations between profiles of these amino acids exhibited a statistically significant average correlation coefficient of 0.84 as compared with an average correlation coefficient of 0.39 over the 16 119 other metabolite cliques containing six metabolites. In silico removal of cliques was used to quantify their importance in determining seminal network properties, highlighting the strong effects of the amino acid clique. Quantitative trait locus analysis revealed co-localization for the six amino acids on chromosome 2, 4 and 10. Sequence analysis identified a unique set of 10 genes on chromosome 2 only, which were associated with amino acid metabolism and specifically the metabolism of Ser-Gly and their conversion into branched-chain amino acids. Metabolite profiling of a set of sublines, with introgressions on chromosome 2, identified a significant change in the abundance of the six amino acids in comparison with M82. Expression analysis of candidate genes affecting Ser metabolism matched the observation from the metabolite data, suggesting a coordinated behavior of the level of these amino acids at the genetic level. Analysis of transcription factor binding sites in the promoter regions of the identified genes suggested combinatorial response to light and the circadian clock.

The MORPH algorithm: ranking candidate genes for membership in Arabidopsis and tomato pathways.

Closing gaps in our current knowledge about biological pathways is a fundamental challenge. The development of novel computational methods along with high-throughput experimental data carries the promise to help in the challenge. We present an algorithm called MORPH (for module-guided ranking of candidate pathway genes) for revealing unknown genes in biological pathways. The method receives as input a set of known genes from the target pathway, a collection of expression profiles, and interaction and metabolic networks. Using machine learning techniques, MORPH selects the best combination of data and analysis method and outputs a ranking of candidate genes predicted to belong to the target pathway. We tested MORPH on 230 known pathways in Arabidopsis thaliana and 93 known pathways in tomato (Solanum lycopersicum) and obtained high-quality cross-validation results. In the photosynthesis light reactions, homogalacturonan biosynthesis, and chlorophyll biosynthetic pathways of Arabidopsis, genes ranked highly by MORPH were recently verified to be associated with these pathways. MORPH candidates ranked for the carotenoid pathway from Arabidopsis and tomato are derived from pathways that compete for common precursors or from pathways that are coregulated with or regulate the carotenoid biosynthetic pathway.

The MORPH-R web server and software tool for predicting missing genes in biological pathways.

A biological pathway is the set of molecular entities involved in a given biological process and the interrelations among them. Even though biological pathways have been studied extensively, discovering missing genes in pathways remains a fundamental challenge. Here, we present an easy-to-use tool that allows users to run MORPH (MOdule-guided Ranking of candidate PatHway genes), an algorithm for revealing missing genes in biological pathways, and demonstrate its capabilities. MORPH supports the analysis in tomato, Arabidopsis and the two new species: rice and the newly sequenced potato genome. The new tool, called MORPH-R, is available both as a web server (at http://bioinformatics.psb.ugent.be/webtools/morph/) and as standalone software that can be used locally. In the standalone version, the user can apply the tool to new organisms using any proprietary and public data sources.

Lycopene cyclase paralog CruP protects against reactive oxygen species in oxygenic photosynthetic organisms.

In photosynthetic organisms, carotenoids serve essential roles in photosynthesis and photoprotection. A previous report designated CruP as a secondary lycopene cyclase involved in carotenoid biosynthesis [Maresca J, et al. (2007) Proc Natl Acad Sci USA 104:11784-11789]. However, we found that cruP KO or cruP overexpression plants do not exhibit correspondingly reduced or increased production of cyclized carotenoids, which would be expected if CruP was a lycopene cyclase. Instead, we show that CruP aids in preventing accumulation of reactive oxygen species (ROS), thereby reducing accumulation of ß-carotene-5,6-epoxide, a ROS-catalyzed autoxidation product, and inhibiting accumulation of anthocyanins, which are known chemical indicators of ROS. Plants with a nonfunctional cruP accumulate substantially higher levels of ROS and ß-carotene-5,6-epoxide in green tissues. Plants overexpressing cruP show reduced levels of ROS, ß-carotene-5,6-epoxide, and anthocyanins. The observed up-regulation of cruP transcripts under photoinhibitory and lipid peroxidation-inducing conditions, such as high light stress, cold stress, anoxia, and low levels of CO(2), fits with a role for CruP in mitigating the effects of ROS. Phylogenetic distribution of CruP in prokaryotes showed that the gene is only present in cyanobacteria that live in habitats characterized by large variation in temperature and inorganic carbon availability. Therefore, CruP represents a unique target for developing resilient plants and algae needed to supply food and biofuels in the face of global climate change.

Evaluation and integration of functional annotation pipelines for newly sequenced organisms: the potato genome as a test case.

BACKGROUND: For most organisms, even if their genome sequence is available, little functional information about individual genes or proteins exists. Several annotation pipelines have been developed for functional analysis based on sequence, 'omics', and literature data. However, researchers encounter little guidance on how well they perform. Here, we used the recently sequenced potato genome as a case study. The potato genome was selected since its genome is newly sequenced and it is a non-model plant even if there is relatively ample information on individual potato genes, and multiple gene expression profiles are available. RESULTS: We show that the automatic gene annotations of potato have low accuracy when compared to a "gold standard" based on experimentally validated potato genes. Furthermore, we evaluate six state-of-the-art annotation pipelines and show that their predictions are markedly dissimilar (Jaccard similarity coefficient of 0.27 between pipelines on average). To overcome this discrepancy, we introduce a simple GO structure-based algorithm that reconciles the predictions of the different pipelines. We show that the integrated annotation covers more genes, increases by over 50% the number of highly co-expressed GO processes, and obtains much higher agreement with the gold standard. CONCLUSIONS: We find that different annotation pipelines produce different results, and show how to integrate them into a unified annotation that is of higher quality than each single pipeline. We offer an improved functional annotation of both PGSC and ITAG potato gene models, as well as tools that can be applied to additional pipelines and improve annotation in other organisms. This will greatly aid future functional analysis of '-omics' datasets from potato and other organisms with newly sequenced genomes. The new potato annotations are available with this paper.

Single nucleotide variation catalog from clinical isolates mapped on tertiary and quaternary structures of ESX-1-related proteins reveals critical regions as putative Mtb therapeutic targets.

Proteins encoded by the ESX-1 genes of interest are essential for full virulence in all Mycobacterium tuberculosis complex (Mtbc) lineages, the pathogens causing the highest mortality worldwide. Identifying critical regions in these ESX-1-related proteins could provide preventive or therapeutic targets for Mtb infection, the game changer needed for tuberculosis control. We analyzed a compendium of whole genome sequences of clinical Mtb isolates from all lineages from >32,000 patients and identified single nucleotide polymorphisms. When mutations corresponding to all non-synonymous single nucleotide polymorphisms were mapped on structural models of the ESX-1 proteins, fully conserved regions emerged. Some could be assigned to known quaternary structures, whereas others could be predicted to be involved in yet-to-be-discovered interactions. Some mutants had clonally expanded (found in >1% of the isolates); these mutants were mostly located at the surface of globular domains, remote from known intra- and inter-molecular protein-protein interactions. Fully conserved intrinsically disordered regions of proteins were found, suggesting that these regions are crucial for the pathogenicity of the Mtbc. Altogether, our findings highlight fully conserved regions of proteins as attractive vaccine antigens and drug targets to control Mtb virulence. Extending this approach to the whole Mtb genome as well as other microorganisms will enhance vaccine development for various pathogens. IMPORTANCE: We mapped all non-synonymous single nucleotide polymorphisms onto each of the experimental and predicted ESX-1 proteins' structural models and inspected their placement. Varying sizes of conserved regions were found. Next, we analyzed predicted intrinsically disordered regions within our set of proteins, finding two putative long stretches that are fully conserved, and discussed their potential essential role in immunological recognition. Combined, our findings highlight new targets for interfering with Mycobacterium tuberculosis complex virulence.

Volatile organic compound detection of Buruli ulcer disease: Headspace analysis of Mycobacterium ulcerans and used gauzes of Buruli-compatible ulcers.

Diagnosing Buruli ulcer (BU) is complicated by limited access to the sensitive IS2404 qPCR. Experienced clinicians report a distinct odour of Buruli ulcers. We explored the potential of headspace analysis by thermal desorption-gas chromatography-mass spectrometry to detect volatile organic compounds (VOCs) from Mycobacterium ulcerans both in vitro and clinically. This study was conducted in two phases: a discovery and validation phase. During the discovery phase, VOCs that enable identification of M. ulcerans cultures were determined. During the validation phase, these VOCs were evaluated in clinical samples for which we used gauzes from patients with skin ulcerations in the Democratic Republic of Congo. Seven M. ulcerans headspace samples were compared with four from sterile growth medium and laboratory environmental air. The univariate analysis resulted in the selection of 24 retained VOC fragments and a perfect differentiation between cultures and controls. Sixteen of 24 fragments were identified, resulting in eleven unique compounds, mainly alkanes. Methylcyclohexane was the best performing compound. Based on these 24 fragments, headspace samples originating from gauzes of 50 open skin lesions (12 qPCR positive and 38 negative) were analysed and an AUC of 0.740 (95%-CI 0.583-0.897) was obtained. As this is an experimental study, future research has to confirm whether the identified compounds can serve as novel biomarkers.

Transcriptome and metabolome profiling identify factors potentially involved in pro-vitamin A accumulation in cassava landraces.

Cassava (Manihot esculenta Crantz) is a predominant food security crop in several developing countries. Its storage roots, rich in carbohydrate, are deficient in essential micronutrients, including provitamin A carotenoids. Increasing carotenoid content in cassava storage roots is important to reduce the incidence of vitamin A deficiency, a public health problem in sub-Saharan Africa. However, cassava improvement advances slowly, mainly due to limited information on the molecular factors influencing ß-carotene accumulation in cassava. To address this problem, we performed comparative transcriptomic and untargeted metabolic analyses of roots and leaves of eleven African cassava landraces ranging from white to deep yellow colour, to uncover regulators of carotenoid biosynthesis and accumulation with conserved function in yellow cassava roots. Sequence analysis confirmed the presence of a mutation, known to influence ß-carotene content, in PSY transcripts of deep yellow but not of pale yellow genotypes. We identified genes and metabolites with expression and accumulation levels significantly associated with ß-carotene content. Particularly an increased activity of the abscisic acid catabolism pathway together with a reduced amount of L-carnitine, may be related to the carotenoid pathway flux, higher in yellow than in white storage roots. In fact, NCED_3.1 was specifically expressed at a lower level in all yellow genotypes suggesting that it could be a potential target for increasing carotenoid accumulation in cassava. These results expand the knowledge on metabolite compositions and molecular mechanisms influencing carotenoid biosynthesis and accumulation in cassava and provide novel information for biotechnological applications and genetic improvement of cassava with high nutritional values.

Multidrug-resistant tuberculosis control in Rwanda overcomes a successful clone that causes most disease over a quarter century.

SUMMARY BACKGROUND: Multidrug-resistant (MDR) tuberculosis (TB) poses an important challenge in TB management and control. Rifampicin resistance (RR) is a solid surrogate marker of MDR-TB. We investigated the RR-TB clustering rates, bacterial population dynamics to infer transmission dynamics, and the impact of changes to patient management on these dynamics over 27 years in Rwanda. METHODS: We analysed whole genome sequences of a longitudinal collection of nationwide RR-TB isolates. The collection covered three important periods: before programmatic management of MDR-TB (PMDT; 1991-2005), the early PMDT phase (2006-2013), in which rifampicin drug-susceptibility testing (DST) was offered to retreatment patients only, and the consolidated phase (2014-2018), in which all bacteriologically confirmed TB patients had rifampicin DST done mostly via Xpert MTB/RIF assay. We constructed clusters based on a 5 SNP cut-off and resistance conferring SNPs. We used Bayesian modelling for dating and population size estimations, TransPhylo to estimate the number of secondary cases infected by each patient, and multivariable logistic regression to assess predictors of being infected by the dominant clone. RESULTS: Of 308 baseline RR-TB isolates considered for transmission analysis, the clustering analysis grouped 259 (84.1%) isolates into 13 clusters. Within these clusters, a single dominant clone was discovered containing 213 isolates (82.2% of clustered and 69.1% of all RR-TB), which we named the "Rwanda Rifampicin-Resistant clone" (R3clone). R3clone isolates belonged to Ugandan sub-lineage 4.6.1.2 and its rifampicin and isoniazid resistance were conferred by the Ser450Leu mutation in rpoB and Ser315Thr in katG genes, respectively. All R3clone isolates had Pro481Thr, a putative compensatory mutation in the rpoC gene that likely restored its fitness. The R3clone was estimated to first arise in 1987 and its population size increased exponentially through the 1990s', reaching maximum size (∼84%) in early 2000 s', with a declining trend since 2014. Indeed, the highest proportion of R3clone (129/157; 82·2%, 95%CI: 75·3-87·8%) occurred between 2000 and 13, declining to 64·4% (95%CI: 55·1-73·0%) from 2014 onward. We showed that patients with R3clone detected after an unsuccessful category 2 treatment were more likely to generate secondary cases than patients with R3clone detected after an unsuccessful category 1 treatment regimen. CONCLUSIONS: RR-TB in Rwanda is largely transmitted. Xpert MTB/RIF assay as first diagnostic test avoids unnecessary rounds of rifampicin-based TB treatment, thus preventing ongoing transmission of the dominant R3clone. As PMDT was intensified and all TB patients accessed rifampicin-resistance testing, the nationwide R3clone burden declined. To our knowledge, our findings provide the first evidence supporting the impact of universal DST on the transmission of RR-TB.