Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

植物中异源转座因子标签的研究

把已知的异源转座因子转化植物，通过转座因子标签法定向诱导突变体，己成为克隆基因的有效途径之一．本研究在充分了解该领域研究进展的同时，把玉米转座因子Ac导入单倍体烟草，尝试用单倍体细胞的遗传体系进行转座因子标签研究;并把Ac转化普通烟草双二倍体，分析Ac在转基因烟草及其后代中的行为，探索突变体筛选的有效途径，另外．构建了一系列用于单子叶植物的农杆菌转化．Ds．，Ac转座因子标签及Cre／lox位点特异性重墨的双元载体系统．上述体系为禾谷类植物的遗传转化．基因突变与克隆，以及基因组图谱的分析提供了研究基础．

Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis.

Cytosine methylation is important for transposon silencing and epigenetic regulation of endogenous genes, although the extent to which this DNA modification functions to regulate the genome is still unknown. Here we report the first comprehensive DNA methylation map of an entire genome, at 35 base pair resolution, using the flowering plant Arabidopsis thaliana as a model. We find that pericentromeric heterochromatin, repetitive sequences, and regions producing small interfering RNAs are heavily methylated. Unexpectedly, over one-third of expressed genes contain methylation within transcribed regions, whereas only approximately 5% of genes show methylation within promoter regions. Interestingly, genes methylated in transcribed regions are highly expressed and constitutively active, whereas promoter-methylated genes show a greater degree of tissue-specific expression. Whole-genome tiling-array transcriptional profiling of DNA methyltransferase null mutants identified hundreds of genes and intergenic noncoding RNAs with altered expression levels, many of which may be epigenetically controlled by DNA methylation.

A miniature insertion element transposable in Microcystis sp FACHB 854

A 186-bp sequence with imperfect terminal inverted repeats and target direct repeats but without any transposase-encoding capacity was found to be transposable in an isolate derived from Microcystis sp. FACHB 854. This miniature insertion element, designated as ISM854-1, and with its homologues present at least 10 copies in the genome of Microcystis FACHB 854, is inserted into the 8-bp long and AT-rich target sequences, but none or few in other Microcystis strains. A variant of ISM854-1, denoted ISM854-1A, has perfect inverted repeat sequences and may transpose in pairs in a structure like a composite transposon. This is the first report of non-autonomous transposition of a mini-IS in a cyanobacterium.

Excision of Sleeping Beauty transposons: parameters and applications to gene therapy

A major problem in gene therapy is the determination of the rates at which gene transfer has occurred. Our work has focused on applications of the Sleeping Beauty (SB) transposon system as a non-viral vector for gene therapy. Excision of a transposon from a donor molecule and its integration into a cellular chromosome are catalyzed by SB transposase. In this study, we used a plasmid-based excision assay to study the excision step of transposition. We used the excision assay to evaluate the importance of various sequences that border the sites of excision inside and outside the transposon in order to determine the most active sequences for transposition from a donor plasmid. These findings together with our previous results in transposase binding to the terminal repeats suggest that the sequences in the transposon-junction of SB are involved in steps subsequent to DNA binding but before excision, and that they may have a role in transposase-transposon interaction. We found that SB transposons leave characteristically different footprints at excision sites in different cell types, suggesting that alternative repair machineries operate in concert with transposition. Most importantly, we found that the rates of excision correlate with the rates of transposition. We used this finding to assess transposition in livers of mice that were injected with the SB transposon and transposase. The excision assay appears to be a relatively quick and easy method to optimize protocols for delivery of genes in SB transposons to mammalian chromosomes in living animals. Copyright (C) 2004 John Wiley Sons, Ltd.

alr0117, a two-component histidine kinase gene, is involved in heterocyst development in Anabaena sp PCC 7120

Anabaena sp. PCC; 7120 was mutagenized by transposon Tn5-1087b, generating a mutant whose heterocysts lack the envelope polysaccharide layer. The transposon was located between nucleotides 342 and 343 of alr0117, a 918 bp gene encoding a histidine kinase for a two-component regulatory system. Complementation of the mutant with a DNA fragment containing alr0117 and targeted inactivation of the gene confirmed that alr0117 is involved in heterocyst development. RT-PCR showed that alr0117 was constitutively expressed in the presence or absence of a combined-nitrogen source. hepA and patB, the two genes turned on during wild-type heterocyst development, were no longer activated in an alr0117-null mutant. The two-component signal transduction system involving alr0117 may control the formation of the envelope polysaccharide layer and certain late events essential to the function of heterocysts.

Construction and analysis of a high-CO2-requiring mutant of the cyanobacterium Anabaena sp strain PCC7120

A mutant of Anabaena sp. strain PCC7120 requiring high CO2 was generated using Tn5 mutagenesis. This is the first data for a filamentous cyanobacterium. The mutant was capable of growing at 5% CO2, but incapable of growing at air levels of CO2. Southern hybridization analysis indicated that the Anabaena genome was inserted by the transposon at one site. The apparent photosynthetic affinity of the mutant to external dissolved inorganic carbon (DIC) was about 300 times lower that of the wild type (WT), and the medium alkalization rate as well as the carboxysomal carbonic anhydrase activity of the mutant was also lower than those of the WT. When the mutant was transferred from the culture medium bubbled with 5% CO2 to higher DIC (8.4% CO2) or 1% CO2, it showed similar responses to the WT. However, aberrant carboxysomes were found in the mutant cells through ultrastructural analysis, indicating it was most probably the wrong organization of the carboxysomes that eventually led to the inefficient operation of carboxysomal carbonic anhydrase and the subsequent defectiveness of the mutant in utilizing DIC.

转座子Minos介导的果蝇转基因平台构建及水稻高表达启动子的克隆 任娟 指

本实验室果蝇研究工作，主要集中在黑腹果蝇的新基因起源的研究。新基因起源的分子机制主要包括：外显子重排、基因复制、基因逆转座、移动元件介导、基因水平转移、基因从头起源、基因的断裂融合。为了阐述这些新基因的产生和它们所带来的物种适应性，我们对这些新近起源的基因进行了功能研究。但是，仅仅限于新基因所在物种的功能研究并不能完全解释新基因产生的进化原因，我们需要了解它是否能够给没有该基因的果蝇物种带来一定的适应性。例如一些生殖相关新基因，如果我们将它们转入没有该基因的果蝇，那是否能够给该果蝇带来生殖能力的提高？无论结果如何，这都为我们研究新基因的起源提供一个重要线索。由此，黑腹果蝇以外的其它果蝇物种中实现转基因成为该研究的重要技术环节。但是，实验室目前的转基因系统仅限于P转座子介导的黑腹果蝇转基因系统，因而我们需要建立一种新的转基因平台。而转座子Minos打破物种范围的转基因特性，以及它的转座特点为我们提供了选择。转座子Minos是从果蝇D. hydei中克隆出来长约1.8kb的Ⅱ型转座子，Tc1家族转座元件成员。Minos的转座机制与大部分转座子一样，在宿主基因组里面实行着剪切和粘贴的运作机制。Minos在转座时，偏向插入TA位点并且主要集中于内含子区域，这样可以减少对插入位置基因的影响。此外，Minos在黑腹果蝇中的转座效率约30%，并且拥有一套成熟的选择标记。因此，Minos成为我们解决非黑腹果蝇转基因技术难题的首选。 在本文的工作中，我们采用由希腊Savakis教授（希腊分子生物学与生物技术研究所）提供的Minos转基因系统，完成果蝇的转基因实验。在这套转基因系统中，非自主的转座子Minos和转座酶基因被克隆到了不同载体当中。其中Minos转座子序列中插入了由3xP3眼睛特异表达的启动子介导表达的eGFP报告基因，而转座酶基因则由热激蛋白hsp70启动子调控表达。实验过程中，我们在果蝇D. melanogaster 和D. yakuba的胚胎中分别同时显微注射入含有转座子和转座酶本实验室果蝇研究工作，主要集中在黑腹果蝇的新基因起源的研究。新基因起源的分子机制主要包括：外显子重排、基因复制、基因逆转座、移动元件介导、基因水平转移、基因从头起源、基因的断裂融合。为了阐述这些新基因的产生和它们所带来的物种适应性，我们对这些新近起源的基因进行了功能研究。但是，仅仅限于新基因所在物种的功能研究并不能完全解释新基因产生的进化原因，我们需要了解它是否能够给没有该基因的果蝇物种带来一定的适应性。例如一些生殖相关新基因，如果我们将它们转入没有该基因的果蝇，那是否能够给该果蝇带来生殖能力的提高？无论结果如何，这都为我们研究新基因的起源提供一个重要线索。由此，黑腹果蝇以外的其它果蝇物种中实现转基因成为该研究的重要技术环节。但是，实验室目前的转基因系统仅限于P转座子介导的黑腹果蝇转基因系统，因而我们需要建立一种新的转基因平台。而转座子Minos打破物种范围的转基因特性，以及它的转座特点为我们提供了选择。转座子Minos是从果蝇D. hydei中克隆出来长约1.8kb的Ⅱ型转座子，Tc1家族转座元件成员。Minos的转座机制与大部分转座子一样，在宿主基因组里面实行着剪切和粘贴的运作机制。Minos在转座时，偏向插入TA位点并且主要集中于内含子区域，这样可以减少对插入位置基因的影响。此外，Minos在黑腹果蝇中的转座效率约30%，并且拥有一套成熟的选择标记。因此，Minos成为我们解决非黑腹果蝇转基因技术难题的首选。 在本文的工作中，我们采用由希腊Savakis教授（希腊分子生物学与生物技术研究所）提供的Minos转基因系统，完成果蝇的转基因实验。在这套转基因系统中，非自主的转座子Minos和转座酶基因被克隆到了不同载体当中。其中Minos转座子序列中插入了由3xP3眼睛特异表达的启动子介导表达的eGFP报告基因，而转座酶基因则由热激蛋白hsp70启动子调控表达。实验过程中，我们在果蝇D. melanogaster 和D. yakuba的胚胎中分别同时显微注射入含有转座子和转座酶所在的质粒。转座酶在37度条件诱导下进行表达，协助Minos完成转座过程。在转基因果蝇的阳性筛选中，我们利用眼睛特异表达的绿色荧光蛋作为选择标记。并且，我们通过PCR实验进一步验证了转基因果蝇的真实性。本研究中，我们对转基因实验条件进行了初步优化。我们通过对黑腹果蝇白眼突变品系W1118和D. yakuba注射后胚胎进行保湿，对D. yakuba注射胚胎进行非退壳处理。在改进条件下W1118和D. yakuba的存活率分别为10%和3%左右。通过筛选转基因阳性果蝇，我们得出Minos在W1118和D. yakuba中的转座效率分别在32%和20%左右。我们的实验结果再一次证实了Minos在果蝇D. melanogaster中可行性。同时，该工作也初步完成了在果蝇D. yakuba 中的第一次Minos介导的转基因实验，为新基因的跨物种功能研究奠定了实验基础。在未来的工作计划中，我们将采用Minos转基因系统，把实验室目前研究的黑腹果蝇新基因导入其它物种果蝇进行功能研究。 水稻是一种重要的世界粮食作物，世界上过半的人口以水稻为主食。水稻相对别的粮食作物来讲具有较小的基因组，并且拥有较好的基因组注释，是一种理想的单子叶模式生物。植物转基因技术的发展推动着水稻功能基因组学的研究，目前水稻的转基因技术主要依赖于土壤细菌农杆菌（Agrobacterium tumefaciens）T-DNA介导的外源基因染色体插入。在自然状态下，农杆菌的T-DNA位于Ti致瘤质粒当中。它包括了一些转座元件和一些帮助T-DNA转座的毒性蛋白基因和调节基因。由于Ti质粒上的T-DNA太长，并且没有太多的酶切位点，因此自然状态的T-DNA不适合进行转基因实验。为了方便T-DNA的实际应用，研究人员创立了双载体转基因系统。T-DNA转座区被分离到出Ti载体，并且装载到另外一个适合实验操作的质粒当中，而毒性蛋白表达基因等则保留在Ti质粒上。因此，在进行T-DNA介导的转基因实验时，需要同时存在T-DNA载体和Ti质粒。 本文以“水稻注释计划数据库RAP-DB”的表达数据为参考，选择了60个高表达基因的启动子区域进行克隆。通过对T-DNA载体pCAMBIA1301 进行改造，去掉其原来的35S启动子，将预测的基因启动子克隆到该载体中并与报告基 摘要 因GUS 基因融合。通过分子克隆实验，我们得到了45个高表达基因的启动子载体。最终，为了测试这45个启动子的启动效率，我们会将它们转化到水稻愈伤组织中通过启动子融合的GUS基于表达情况来判断我们启动子的启动效率。

Characterization of amphioxus GDF8/11 gene, an archetype of vertebrate MSTN and GDF11

MSTN, also known as growth and differentiation factor 8 (GDF8), and GDF11 are members of the transforming growth factor-beta (TGF-beta) subfamily. They have been thought to be derived from one ancestral gene. In the present study, we report the isolation and characterization of an invertebrate GDF8/11 homolog from the amphioxus (Branchiostoma belcheri tsingtauense). The amphioxus GDF8/11 gene consists of five exons flanked by four introns, which have two more exons and introns than that of other species. In intron III, a possible transposable element was identified. This suggested that this intron might be derived from transposon. The amphioxus GDF8/11 cDNA encodes a polypeptide of 419 amino acid residues. Phologenetic analysis shows that the GDF8/11 is at the base of vertebrate MSTNs and GDF11s. This result might prove that the GDF8/11 derived from one ancestral gene and the amphioxus GDF8/11 may be the common ancestral gene, and also the gene duplication event generating MSTN and GDF11 occurred before the divergence of vertebrates and after or at the divergence of amphioxus from vertebrates. Reverse transcriptase polymerase chain reaction results showed that the GDF8/11 gene was expressed in new fertilized cell, early gastrulation, and knife-shaped embryo, which was different from that in mammals. It suggested that the GDF8/11 gene might possess additional functions other than regulating muscle growth in amphioxus.

Molecular analysis of the anaerobic rumen fungus Orpinomyces - insights into an AT-rich genome

Matthew J. Nicholson, Michael K. Theodorou and Jayne L. Brookman. (2005). Molecular analysis of the anaerobic rumen fungus Orpinomyces - insights into an AT-rich genome. Microbiology, 151 (1), 121-133. Sponsorship: BBSRC RAE2008

Investigation of global regulators influencing styrene metabolism and bioplastic synthesis in Pseudomonas putida CA-3

The genetics and biochemistry involved in the biodegradation of styrene and the production of polyhydroxyalkanoates in Pseudomonas putida CA-3 have been well characterised to date. Knowledge of the role played by global regulators in controlling these pathways currently represents a critical knowledge gap in this area. Here we report on our efforts to identify such regulators using mini-Tn5 transposon mutagenesis of the P. putida CA-3 genome. The library generated was subjected to phenotypic screening to identify mutants exhibiting a reduced sensitivity to the effects of carbon catabolite repression of aromatic pathway activity. Our efforts identified a clpX disrupted mutant which exhibited wild-type levels of growth on styrene but significantly reduced growth on phenylacetic acid. RT-PCR analysis of key PACoA catabolon genes necessary for phenylacetic acid metabolism, and SDS-PAGE protein profile analyses suggest that no direct alteration of PACoA pathway transcriptional or translational activity was involved. The influence of global regulators affecting the accumulation of PHAs in P. putida CA-3 was also studied. Phenotypic screening of the mini-Tn5 library revealed a gacS sensor kinase gene disruption resulting in the loss of PHA accumulation capacity in P. putida CA-3. Subsequent SDS-PAGE protein analyses of the wild type and gacS mutant strains identified post-transcriptional control of phaC1 synthase as a key point of control of PHA synthesis in P. putida CA-3. Disruption of the gacS gene in another PHA accumulating organism, P. putida S12, also demonstrated a reduction of PHA accumulation capacity. PHA accumulation was observed to be disrupted in the CA-3 gacS mutant under phosphorus limited growth conditions. Over-expression studies in both wild type CA-3 and gacS mutant demonstrated that rsmY over-expression in gacS disrupted P. putida CA-3 is insufficient to restore PHA accumulation in the cell however in wild type cells, over-expression of rsmY results in an altered PHA monomer compositions.

Investigation of growth and stress tolerance characteristics of Cronobacter spp.

Cronobacter spp. are opportunistic pathogens which can be isolated from a wide variety of foods and environments. They are Gram negative, motile, non-spore forming, peritrichous rods of the Enterobacteriaceae family. This food-borne pathogen is associated with the ingestion of contaminated infant milk formula (IMF), causing necrotizing enterocolitis, sepsis and meningitis in neonatal infants. The work presented in this thesis involved the investigation and characterisation of a bank of Cronobacter strains for their ability to tolerate physiologically relevant stress conditions that are commonly encountered in the gastrointestinal tract. While all strains were able to endure the suboptimal conditions tested, noteworthy variations were observed between strains. A collection of these strains were Lux-tagged to determine if their growth could be tracked in IMF by measuring bioluminescence. The resulting strains could be easily and reproducibly monitored in real time by measuring light emission. Following this a transposon mutagenesis library was created in one of the Lux-tagged strains of Cronobacter sakazakii. This library was screened for mutants with affected growth in milk. The majority of mutants identified were associated with amino acid metabolism. The final section of this thesis identified genes involved in the tolerance of C. sakazakii to the milk derived antimicrobial peptide, Lactoferricin B (Lfcin B). This was achieved by creating a transposon mutagenesis library in C. sakazakii and screening for mutants with increased susceptibility to Lfcin B. Overall this thesis demonstrates the variation between Cronobacter strains. It also identifies genes required for growth of the bacteria in milk, as well as genes needed for antimicrobial peptide tolerance.

Environmental and genetic determinants of colony morphology in yeast.

Nutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Overexpression of the rhamnose catabolism regulatory protein, RhaR: a novel mechanism for metronidazole resistance in Bacteroides thetaiotaomicron

Objectives: The aim of the investigation was to use in vitro transposon mutagenesis to generate metronidazole resistance in the obligately anaerobic pathogenic bacterium Bacteroides thetaiotaomicron, and to identify the genes involved to enable investigation of potential mechanisms for the generation of metronidazole resistance.
Methods: The genes affected by the transposon insertion were identified by plasmid rescue and sequencing. Expression levels of the relevant genes were determined by semi-quantitative RNA hybridization and catabolic activity by lactate dehydrogenase/pyruvate oxidoreductase assays.
Results: A metronidazole-resistant mutant was isolated and the transposon insertion site was identified in an intergenic region between the rhaO and rhaR genes of the gene cluster involved in the uptake and catabolism of rhamnose. Metronidazole resistance was observed during growth in defined medium containing either rhamnose or glucose. The metronidazole-resistant mutant showed improved growth in the presence of rhamnose as compared with the wild-type parent. There was increased transcription of all genes of the rhamnose gene cluster in the presence of rhamnose and glucose, likely due to the transposon providing an additional promoter for the rhaR gene, encoding the positive transcriptional regulator of the rhamnose operon. The B. thetaiotaomicron metronidazole resistance phenotype was recreated by overexpressing the rhaR gene in the B. thetaiotaomicron wild-type parent. Both the metronidazole-resistant transposon mutant and RhaR overexpression strains displayed a phenotype of higher lactate dehydrogenase and lower pyruvate oxidoreductase activity in comparison with the parent strain during growth in rhamnose.
Conclusions: These data indicate that overexpression of the rhaR gene generates metronidazole resistance in B. thetaiotaomicron