Roving methyltransferases generate a mosaic epigenetic landscape and influence evolution in Bacteroides fragilis group

Three types of DNA methyl modifications have been detected in bacterial genomes, and mechanistic studies have demonstrated roles for DNA methylation in physiological functions ranging from phage defense to transcriptional control of virulence and host-pathogen interactions. Despite the ubiquity of methyltransferases and the immense variety of possible methylation patterns, epigenomic diversity remains unexplored for most bacterial species. Members of the Bacteroides fragilis group (BFG) reside in the human gastrointestinal tract as key players in symbiotic communities but also can establish anaerobic infections that are increasingly multi-drug resistant. In this work, we utilize long-read sequencing technologies to perform pangenomic (n = 383) and panepigenomic (n = 268) analysis of clinical BFG isolates cultured from infections seen at the NIH Clinical Center over four decades. Our analysis reveals that single BFG species harbor hundreds of DNA methylation motifs, with most individual motif combinations occurring uniquely in single isolates, implying immense unsampled methylation diversity within BFG epigenomes. Mining of BFG genomes identified more than 6000 methyltransferase genes, approximately 1000 of which were associated with intact prophages. Network analysis revealed substantial gene flow among disparate phage genomes, implying a role for genetic exchange between BFG phages as one of the ultimate sources driving BFG epigenome diversity.

The Host Adapted Fungal Pathogens of Pneumocystis Genus Utilize Genic Regional Centromeres

Centromeres are genomic regions that coordinate accurate chromosomal segregation during mitosis and meiosis. Yet, despite their essential function, centromeres evolve rapidly across eukaryotes. Centromeres are often the sites of chromosomal breaks which contribute to genome shuffling and promote speciation by inhibiting gene flow. How centromeres form in strongly host-adapted fungal pathogens has yet to be investigated. Here, we characterized the centromere structures in closely related species of mammalian-specific pathogens of the fungal phylum of Ascomycota. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of Schizosaccharomyces pombe. Using organisms from a short-term in vitro culture or infected animal models and ChIP-seq, we identified centromeres in three Pneumocystis species that diverged ~100 million years ago. Each species has a unique short regional centromere (< 10kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. CENP-C, a scaffold protein that links the inner centromere to the kinetochore appears dispensable in one species, suggesting a kinetochore rewiring. Despite the loss of DNA methyltransferases, 5-methylcytosine DNA methylation occurs in these species, though not related to centromere function. These features suggest an epigenetic specification of centromere function.

Hypermutator strains of Pseudomonas aeruginosa reveal novel pathways of resistance to combinations of cephalosporin antibiotics and beta-lactamase inhibitors

Hypermutation due to DNA mismatch repair (MMR) deficiencies can accelerate the development of antibiotic resistance in Pseudomonas aeruginosa. Whether hypermutators generate resistance through predominantly similar molecular mechanisms to wild-type (WT) strains is not fully understood. Here, we show that MMR-deficient P. aeruginosa can evolve resistance to important broad-spectrum cephalosporin/beta-lactamase inhibitor combination antibiotics through novel mechanisms not commonly observed in WT lineages. Using whole-genome sequencing (WGS) and transcriptional profiling of isolates that underwent in vitro adaptation to ceftazidime/avibactam (CZA), we characterized the detailed sequence of mutational and transcriptional changes underlying the development of resistance. Surprisingly, MMR-deficient lineages rapidly developed high-level resistance (>256 μg/mL) largely without corresponding fixed mutations or transcriptional changes in well-established resistance genes. Further investigation revealed that these isolates had paradoxically generated an early inactivating mutation in the mexB gene of the MexAB-OprM efflux pump, a primary mediator of CZA resistance in P. aeruginosa, potentially driving an evolutionary search for alternative resistance mechanisms. In addition to alterations in a number of genes not known to be associated with resistance, 2 mutations were observed in the operon encoding the RND efflux pump MexVW. These mutations resulted in a 4- to 6-fold increase in resistance to ceftazidime, CZA, cefepime, and ceftolozane-tazobactam when engineered into a WT strain, demonstrating a potentially important and previously unappreciated mechanism of resistance to these antibiotics in P. aeruginosa. Our results suggest that MMR-deficient isolates may rapidly evolve novel resistance mechanisms, sometimes with complex dynamics that reflect gene inactivation that occurs with hypermutation. The apparent ease with which hypermutators may switch to alternative resistance mechanisms for which antibiotics have not been developed may carry important clinical implications.

Disseminated Mycoplasma orale infection in patients with Phosphoinositide-3-Kinase Regulatory Subunit 1 Mutations

Mycoplasma orale is a rare cause of invasive infection in immunodeficient hosts. Phosphatidylinositol 3-Kinase, Regulatory Subunit 1 (PI3KR1) mutations predispose patients to sinopulmonary infections, alongside bronchiectasis autoimmunity and lymphoproliferation. We report the two cases of PI3KR1 deficiency with invasive M. orale and effective treatment options.

Genomic insights into the host specific adaptation of the Pneumocystis genus

Pneumocystis jirovecii, the fungal agent of human Pneumocystis pneumonia, is closely related to macaque Pneumocystis. Little is known about other Pneumocystis species in distantly related mammals, none of which are capable of establishing infection in humans. The molecular basis of host specificity in Pneumocystis remains unknown as experiments are limited due to an inability to culture any species in vitro. To explore Pneumocystis evolutionary adaptations, we have sequenced the genomes of species infecting macaques, rabbits, dogs and rats and compared them to available genomes of species infecting humans, mice and rats. Complete whole genome sequence data enables analysis and robust phylogeny, identification of important genetic features of the host adaptation, and estimation of speciation timing relative to the rise of their mammalian hosts. Our data reveals insights into the evolution of P. jirovecii, the sole member of the genus able to infect humans.

Cissé, Ma et al. utilize genomic data from Pneumocystis species infecting macaques, rabbit, dogs and rats to investigate the molecular basis of host specificity in Pneumocystis. Their analyses provide insight to the specific adaptations enabling the infection of humans by P. jirovecii.

ACKR1 Alleles at 5.6 kb in a Well-Characterized Renewable US Food and Drug Administration (FDA) Reference Panel for Standardization of Blood Group Genotyping

The glycoprotein encoded by the ACKR1 gene expresses the Duffy blood group antigens and is a receptor for malaria parasites. We recently described 18 long-range ACKR1 alleles in an autochthonous population of a malaria endemic region. Extending this work, we sequenced the gene in a 53-sample repository established by the US Food and Drug Administration (FDA) as reference reagents for blood group genotyping. The FDA samples have been characterized for 19 genes; however, long-range haplotype information for these genes, including ACKR1, was lacking. We used a hybrid approach, novel for this type of gene, to characterize ACKR1 by combining two next-generation sequencing technologies, the short-read massively parallel sequencing and the long-read nanopore sequencing. The expedient integration of data from both next-generation sequencing systems were necessary and sufficient to allow determination of all 25 long-range ACKR1 alleles found in the 53 samples accurately. All 25 alleles identified in our current FDA cohort were novel and, unexpectedly, none had been observed among the 18 alleles in our previous study. The alleles will be useful for validation, calibration, and proficiency testing of red cell genotyping. The lack of any overlap between the ACKR1 alleles in the two studies documents differences in mutation rate and recombination frequency among populations. The exact haplotype and their interethnic or interpopulation dissimilarities can influence disease susceptibility and therapy.

232. Genomic Evidence for Dissemination of Mycobacterium marinum in an HIV Patient with Multifocal Cutaneous Disease

Background

Hematogenous dissemination has been proposed to explain multifocal cutaneous disease caused by Mycobacterium marinum in certain patients. Treatment duration for disseminated disease is often months longer than for skin infection alone. However, distinguishing multiple independent inoculation events from dissemination has relied primarily on clinical judgement. Additionally, whether temperature-sensitive non-tuberculous mycobacteria such as M. marinum are indeed capable of invading the vascular space at core body temperature is debated. Here we used whole-genome sequencing (WGS) of serial isolates from a single patient with multifocal cutaneous M. marinum infection to distinguish dissemination of a clonal strain from multiple inoculation events.

Methods

A 35-year-old male with HIV (CD4 of 66 cells/µL) presented with a two-month history of a non-healing M. marinum wound on his left elbow (isolate MM0). This was followed a month later after initiation of antiretroviral therapy by a second M. marinum lesion on the right heel (MM1) without history of repeat inoculation, and increased swelling and erythema of the wound on the left arm (MM2) consistent with paradoxical immune reconstitution inflammatory syndrome. A PacBio genome was generated for MM0 and short read Illumina genomes were generated for MM1 and MM2.

Results

All isolates were found to be closely related, with MM1 and MM2 distinguished from MM0 by one and five single-nucleotide variants (SNVs), respectively. Given the substantial genetic heterogeneity among environmental M. marinum strains, such close relatedness of these isolates suggests common origin, and provides strong evidence for dissemination of a clonal strain in this patient. The SNVs included a frameshift mutation in the purT gene, which encodes a formate-dependent phosphoribosylglycinamide formyltransferase involved in de novo purine synthesis, and missense mutations in atsA and the DNA methylase hsdM. All isolates grew at 35°C, compared with the optimal growth temperature of 30°C typically observed for M. marinum, suggesting thermotolerance permissive for dissemination.

Conclusion

These results demonstrate the potential role of WGS for providing supportive evidence of disseminated infection with M. marinum.

Disclosures

All authors: No reported disclosures.

Dynamic Emergence of Mismatch Repair Deficiency Facilitates Rapid Evolution of Ceftazidime-Avibactam Resistance in Pseudomonas aeruginosa Acute Infection

Strains of Pseudomonas aeruginosa with deficiencies in DNA mismatch repair have been studied in the context of chronic infection, where elevated mutational rates ("hypermutation") may facilitate the acquisition of antimicrobial resistance. Whether P. aeruginosa hypermutation can also play an adaptive role in the more dynamic context of acute infection remains unclear. In this work, we demonstrate that evolved mismatch repair deficiencies may be exploited by P. aeruginosa to facilitate rapid acquisition of antimicrobial resistance in acute infection, and we directly document rapid clonal succession by such a hypermutating lineage in a patient. Whole-genome sequencing (WGS) was performed on nine serially cultured blood and respiratory isolates from a patient in whom ceftazidime-avibactam (CZA) resistance emerged in vivo over the course of days. The CZA-resistant clone was differentiated by 14 mutations, including a gain-of-function G183D substitution in the PDC-5 chromosomal AmpC cephalosporinase conferring CZA resistance. This lineage also contained a substitution (R656H) at a conserved position in the ATPase domain of the MutS mismatch repair (MMR) protein, and elevated mutational rates were confirmed by mutational accumulation experiments with WGS of evolved lineages in conjunction with rifampin resistance assays. To test whether MMR-deficient hypermutation could facilitate rapid acquisition of CZA resistance, in vitro adaptive evolution experiments were performed with a mutS-deficient strain. These experiments demonstrated rapid hypermutation-facilitated acquisition of CZA resistance compared with the isogenic wild-type strain. Our results suggest a possibly underappreciated role for evolved MMR deficiency in facilitating rapid adaptive evolution of P. aeruginosa in the context of acute infection.IMPORTANCE Antimicrobial resistance in bacteria represents one of the most consequential problems in modern medicine, and its emergence and spread threaten to compromise central advances in the treatment of infectious diseases. Ceftazidime-avibactam (CZA) belongs to a new class of broad-spectrum beta-lactam/beta-lactamase inhibitor combinations designed to treat infections caused by multidrug-resistant bacteria. Understanding the emergence of resistance to this important new drug class is of critical importance. In this work, we demonstrate that evolved mismatch repair deficiency in P. aeruginosa, an important pathogen responsible for significant morbidity and mortality among hospitalized patients, may facilitate rapid acquisition of resistance to CZA in the context of acute infection. These findings are relevant for both diagnosis and treatment of antimicrobial resistance emerging in acute infection in the hypermutator background and additionally have implications for the emergence of more virulent phenotypes.

Protracted course of disseminated adenovirus disease with necrotizing granulomas in the liver

A 52- year-old male with chronic lymphocytic leukemia was hospitalized with disseminated adenovirus disease. More than a month following recovery, hepatic necrotizing granulomas secondary to adenovirus were found. This case illustrates the protracted course that adenovirus disease may take and emphasizes an unusual presentation with hepatic necrotizing granulomas.

Severe BCG-osis Misdiagnosed as Multidrug-Resistant Tuberculosis in an IL-12Rβ1-Deficient Peruvian Girl

PURPOSE

Mendelian suceptibility to mycobacterial disease (MSMD) is a rare primary immunodeficiency predisposing to severe disease caused by mycobacteria and other intracellular pathogens. Delay in diagnosis can have an impact on the patient's prognosis.

METHODS

We evaluated the IFN-γ circuit by studying IFN-γ production after mycobacterial challenge as well as IL-12Rβ1 expression and STAT4 phosphorylation in response to IL-12p70 stimulation in whole blood of a 6-year-old Peruvian girl with disseminated recurrent mycobacterial infection diagnosed as multidrug-resistant tuberculosis. Genetic studies with Sanger sequencing were used to identify the causative mutation. Microbiological studies based on PCR reactions were used to diagnose the specific mycobacterial species.

RESULTS

We identified a homozygous mutation in the IL12RB1 gene (p. Arg211*) causing abolished expression of IL-12Rβ1 and IL-12 response. MSMD diagnosis led to a microbiological reevaluation of the patient, revealing a BCG vaccine-related infection instead of tuberculosis. Treatment was then adjusted, with good response.

CONCLUSIONS

We report the first Peruvian patient with IL-12Rβ1 deficiency. Specific mycobacterial species diagnosis within Mycobacterium tuberculosis complex is still challenging in countries with limited access to PCR-based microbiological diagnostic techniques. Awareness of MSMD warning signs and accurate microbiological diagnosis of mycobacterial infections are of the utmost importance for optimal diagnosis and management of affected patients.

Rapid Nanopore Sequencing of Plasmids and Resistance Gene Detection in Clinical Isolates

Recent advances in nanopore sequencing technology have led to a substantial increase in throughput and sequence quality. Together, these improvements may permit real-time benchtop genomic sequencing and antimicrobial resistance gene detection in clinical isolates. In this study, we evaluated workflows and turnaround times for a benchtop long-read sequencing approach in the clinical microbiology laboratory using the Oxford Nanopore Technologies MinION sequencer. We performed genomic and plasmid sequencing of three clinical isolates with both MinION and Illumina MiSeq, using different library preparation methods (2D and rapid 1D) with the goal of antimicrobial resistance gene detection. We specifically evaluated the advantages of using plasmid DNA for sequencing and the value of supplementing MinION sequences with MiSeq reads for increasing assembly accuracy. Resequencing of three plasmids in a reference Klebsiella pneumoniae isolate demonstrated ∼99% accuracy of draft MinION-only assembly and >99.9% accuracy of assembly polished with MiSeq reads. Plasmid DNA sequencing of previously uncharacterized clinical extended-spectrum β-lactamase (ESBL)-producing Escherichia coli and K. pneumoniae isolates using MinION allowed successful identification of antimicrobial resistance genes in the draft assembly corresponding to all classes of observed plasmid-based phenotypic resistance. Importantly, use of plasmid DNA enabled lower depth sequencing, and assemblies sufficient for full antimicrobial resistance gene annotation were obtained with as few as 2,000 to 5,000 reads, which could be acquired in 20 min of sequencing. With a MinION-only workflow that balances accuracy against turnaround time, full annotation of plasmid resistance gene content could be obtained in under 6 h from a subcultured isolate, less time than traditional phenotypic susceptibility testing.

Integrated transcriptome analysis of mouse spermatogenesis

BACKGROUND

Differentiation of primordial germ cells into mature spermatozoa proceeds through multiple stages, one of the most important of which is meiosis. Meiotic recombination is in turn a key part of meiosis. To achieve the highly specialized and diverse functions necessary for the successful completion of meiosis and the generation of spermatozoa thousands of genes are coordinately regulated through spermatogenesis. A complete and unbiased characterization of the transcriptome dynamics of spermatogenesis is, however, still lacking.

RESULTS

In order to characterize gene expression during spermatogenesis we sequenced eight mRNA samples from testes of juvenile mice from 6 to 38 days post partum. Using gene expression clustering we defined over 1,000 novel meiotically-expressed genes. We also developed a computational de-convolution approach and used it to estimate cell type-specific gene expression in pre-meiotic, meiotic and post-meiotic cells. In addition, we detected 13,000 novel alternative splicing events around 40% of which preserve an open reading frame, and found experimental support for 159 computational gene predictions. A comparison of RNA polymerase II (Pol II) ChIP-Seq signals with RNA-Seq coverage shows that gene expression correlates well with Pol II signals, both at promoters and along the gene body. However, we observe numerous instances of non-canonical promoter usage, as well as intergenic Pol II peaks that potentially delineate unannotated promoters, enhancers or small RNA clusters.

CONCLUSIONS

Here we provide a comprehensive analysis of gene expression throughout mouse meiosis and spermatogenesis. Importantly, we find over a thousand of novel meiotic genes and over 5,000 novel potentially coding isoforms. These data should be a valuable resource for future studies of meiosis and spermatogenesis in mammals.

Sensitive mapping of recombination hotspots using sequencing-based detection of ssDNA

Meiotic DNA double-stranded breaks (DSBs) initiate genetic recombination in discrete areas of the genome called recombination hotspots. DSBs can be directly mapped using chromatin immunoprecipitation followed by sequencing (ChIP-seq). Nevertheless, the genome-wide mapping of recombination hotspots in mammals is still a challenge due to the low frequency of recombination, high heterogeneity of the germ cell population, and the relatively low efficiency of ChIP. To overcome these limitations we have developed a novel method—single-stranded DNA (ssDNA) sequencing (SSDS)—that specifically detects protein-bound single-stranded DNA at DSB ends. SSDS comprises a computational framework for the specific detection of ssDNA-derived reads in a sequencing library and a new library preparation procedure for the enrichment of fragments originating from ssDNA. The use of our technique reduces the nonspecific double-stranded DNA (dsDNA) background >10-fold. Our method can be extended to other systems where the identification of ssDNA or DSBs is desired.

Genetic crossovers are predicted accurately by the computed human recombination map

Hotspots of meiotic recombination can change rapidly over time. This instability and the reported high level of inter-individual variation in meiotic recombination puts in question the accuracy of the calculated hotspot map, which is based on the summation of past genetic crossovers. To estimate the accuracy of the computed recombination rate map, we have mapped genetic crossovers to a median resolution of 70 Kb in 10 CEPH pedigrees. We then compared the positions of crossovers with the hotspots computed from HapMap data and performed extensive computer simulations to compare the observed distributions of crossovers with the distributions expected from the calculated recombination rate maps. Here we show that a population-averaged hotspot map computed from linkage disequilibrium data predicts well present-day genetic crossovers. We find that computed hotspot maps accurately estimate both the strength and the position of meiotic hotspots. An in-depth examination of not-predicted crossovers shows that they are preferentially located in regions where hotspots are found in other populations. In summary, we find that by combining several computed population-specific maps we can capture the variation in individual hotspots to generate a hotspot map that can predict almost all present-day genetic crossovers.

In eukaryotes genetic crossovers are responsible for generating genetic diversity and ensuring the proper segregation of chromosomes. Genetic crossovers are tightly clustered in hotspots. Although the existence of hotspots in humans is clearly proven, mechanisms of their formation and the regulation of meiotic recombination in general remain poorly understood. An additional complication in studies of meiotic recombination is the fact that the direct experimental mapping of human hotspots on a genome-wide scale is not feasible with current methods. The best available indirect methods compute the position of hotspots from patterns of historic associations between genetic markers in population samples. In this study we determined the positions of genetic crossovers in ten pedigrees of European origin and then compared the positions of crossovers with the hotspots computed from HapMap data. Importantly, we find that the population-averaged computed map is in close agreement with the observed distribution of genetic crossovers. We also find that cryptic hotspots that are not easily detected in the computed European map can be more effectively identified if other populations are included in the analysis. Our analysis shows that high-resolution recombination profiles are highly similar between distantly related populations and that by including computed hotspots from several populations we can predict nearly all crossovers.

Molecular Features and Functional Constraints in the Evolution of the Mammalian X Chromosome

Recent advances in genomic sequencing of multiple organisms have fostered significant advances in our understanding of the evolution of the sex chromosomes. The integration of this newly available sequence information with functional data has facilitated a considerable refinement of our conceptual framework of the forces driving this evolution. Here we address multiple functional constraints that were encountered in the evolution of the X chromosome and the impact that this evolutionary history has had on its modern behavior.

Variation in patterns of human meiotic recombination

In the last 30 years it has become evident that patterns of meiotic recombination can be highly variable among individuals. The evidence comes from both low and high resolution analyses of hotspots of recombination in human and other species. In addition, a comparison of the recombination profiles in closely related species such as human and chimpanzee reveals essentially no correlation in the position of hotspots. Although the variation in hotspots of meiotic recombination is clearly documented, the mechanisms responsible for such variation are far from being understood. Here we will review the available evidence of natural variation in meiotic recombination and will discuss potential implications of this variation on the functional mechanisms of crossover formation and control.

Gene expression profiles of Spo11-/- mouse testes with spermatocytes arrested in meiotic prophase I

Spo11, a meiosis-specific protein, introduces double-strand breaks on chromosomal DNA and initiates meiotic recombination in a wide variety of organisms. Mouse null Spo11 spermatocytes fail to synapse chromosomes and progress beyond the zygotene stage of meiosis. We analyzed gene expression profiles in Spo11(-/ -)adult and juvenile wild-type testis to describe genes expressed before and after the meiotic arrest resulting from the knocking out of Spo11. These genes were characterized using the Gene Ontology data base. To focus on genes involved in meiosis, we performed comparative gene expression analysis of Spo11(-/ -)and wild-type testes from 15-day mice, when spermatocytes have just entered pachytene. We found that the knockout of Spo11 causes dramatic changes in the level of expression of genes that participate in meiotic recombination (Hop2, Brca2, Mnd1, FancG) and in the meiotic checkpoint (cyclin B2, Cks2), but does not affect genes encoding protein components of the synaptonemal complex. Finally, we discovered unknown genes that are affected by the disruption of the Spo11 gene and therefore may be specifically involved in meiosis and spermatogenesis.

X for intersection: retrotransposition both on and off the X chromosome is more frequent

As the heteromorphic sex chromosomes evolved from a pair of autosomes, the sex chromosomes became increasingly different in gene content and structure from each other and from the autosomes. Although recently there has been progress in documenting and understanding these differences, the molecular mechanisms that have fashioned some of these changes remain unclear. A new study addresses the differential distribution of retroposed genes in human and mouse genomes. Surprisingly, chromosome X is a major source and a preferred target for retrotransposition.

The mouse X chromosome is enriched for sex-biased genes not subject to selection by meiotic sex chromosome inactivation

Sex chromosomes are subject to sex-specific selective evolutionary forces. One model predicts that genes with sex-biased expression should be enriched on the X chromosome. In agreement with Rice's hypothesis, spermatogonial genes are over-represented on the X chromosome of mice and sex- and reproduction-related genes are over-represented on the human X chromosome. Male-biased genes are under-represented on the X chromosome in worms and flies, however. Here we show that mouse spermatogenesis genes are relatively under-represented on the X chromosome and female-biased genes are enriched on it. We used Spo11(-/-) mice blocked in spermatogenesis early in meiosis to evaluate the temporal pattern of gene expression in sperm development. Genes expressed before the Spo11 block are enriched on the X chromosome, whereas those expressed later in spermatogenesis are depleted. Inactivation of the X chromosome in male meiosis may be a universal driving force for X-chromosome demasculinization.

Crystal structure of Escherichia coli protein ybgI, a toroidal structure with a dinuclear metal site

BACKGROUND

The protein encoded by the gene ybgI was chosen as a target for a structural genomics project emphasizing the relation of protein structure to function.

RESULTS

The structure of the ybgI protein is a toroid composed of six polypeptide chains forming a trimer of dimers. Each polypeptide chain binds two metal ions on the inside of the toroid.

CONCLUSION

The toroidal structure is comparable to that of some proteins that are involved in DNA metabolism. The di-nuclear metal site could imply that the specific function of this protein is as a hydrolase-oxidase enzyme.

Characterization of the DNA damage-inducible helicase DinG from Escherichia coli

The dinG promoter was first isolated in a genetic screen scoring for damage-inducible loci in Escherichia coli (Lewis, L. K., Jenkins, M. E., and Mount, D. W. (1992) J. Bacteriol. 174, 3377-3385). Sequence analysis suggests that the dinG gene encodes a putative helicase related to a group of eukaryotic helicases that includes mammalian XPD (Koonin, E. V. (1993) Nucleic Acids Res. 21, 1497), an enzyme involved in transcription-coupled nucleotide excision repair and basal transcription. We have characterized the dinG gene product from E. coli using genetic and biochemical approaches. Deletion of dinG has no severe phenotype, indicating that it is non-essential for cell viability. Both dinG deletion and over-expression of the DinG protein from a multicopy plasmid result in a slight reduction of UV resistance. DinG, purified as a fusion protein from E. coli cells, behaves as a monomer in solution, as judged from gel filtration experiments. DinG is an ATP-hydrolyzing enzyme; single-stranded (ss) DNA stimulates the ATPase activity 15-fold. Kinetic data yield a Hill coefficient of 1, consistent with one ATP-hydrolyzing site per DinG molecule. DinG possesses a DNA helicase activity; it translocates along ssDNA in a 5' --> 3' direction, as revealed in experiments with substrates containing non-natural 5'-5' and 3'-3' linkages. The ATP-dependent DNA helicase activity of DinG requires divalent cations (Mg2+, Ca2+, and Mn2+) but is not observed in the presence of Zn2+. The DinG helicase does not discriminate between ribonucleotide and deoxyribonucleotide triphosphates, and it unwinds duplex DNA with similar efficiency in the presence of ATP or dATP. We discuss the possible involvement of the DinG helicase in DNA replication and repair processes.

Over 1000 genes are involved in the DNA damage response of Escherichia coli

Changes in gene expression after treatment of Escherichia coli cultures with mitomycin C were assessed using gene array technology. Unexpectedly, a large number of genes (nearly 30% of all genes) displayed significant changes in their expression level. Analysis and classification of expression profiles of the corresponding genes allowed us to assign this large number of genes into one or two dozen small clusters of genes with similar expression profiles. This assignment allowed us to describe systematically the changes in the level of gene expression in response to DNA damage. Among the damage-induced genes, more than 100 are novel. From those genes involved in DNA metabolism that have not previously been shown to be induced by DNA damage, the mutS gene involved in mismatch repair is especially noteworthy. In addition to the SOS response, we observed the induction of other stress response pathways, such as those of oxidative stress and osmotic protection. Among the genes that are downregulated in response to DNA damage are numerous protein biosynthesis genes. Analysis of the gene expression data highlighted the essential involvement of sigma(s)-regulated genes and the general stress response network in the response to DNA damage.

Differences in HERV-K LTR insertions in orthologous loci of humans and great apes

The classification of the long terminal repeats (LTRs) of the human endogenous retrovirus HERV-K (HML-2) family was refined according to diagnostic differences between the LTR sequences. The mutation rate was estimated to be approximately equal for LTRs belonging to different families and branches of human endogenous retroviruses (HERVs). An average mutation rate value was calculated based on differences between LTRs of the same HERV and was found to be 0.13% per million years (Myr). Using this value, the ages of different LTR groups belonging to the LTR HML-2 subfamily were found to vary from 3 to 50Myr. Orthologous potential LTR-containing loci from different primate species were PCR amplified using primers corresponding to the genomic sequences flanking LTR integration sites. This allowed us to calculate the phylogenetic times of LTR integrations in primate lineages in the course of the evolution and to demonstrate that they are in good agreement with the LTR ages calculated from the mutation rates. Human-specific integrations for some very young LTRs were demonstrated. The possibility of LTRs and HERVs involvement in the evolution of primates is discussed.

A human endogenous retrovirus-like (HERV) LTR formed more than 10 million years ago due to an insertion of HERV-H LTR into the 5' LTR of HERV-K is situated on human chromosomes 10, 19 and Y

A chimeric long terminal repeat (LTR) containing the whole LTR of a human endogenous retrovirus-like element of the H family (HERV-H) inserted downstream of the core enhancer region of the 5' LTR of a HERV-K retroelement was detected and sequenced in the human 19p12 locus, known to be enriched with genes encoding zinc finger proteins. Similar chimeras were also detected in human chromosomes 10 and Y in human-hamster hybrid cells containing individual human chromosomes. This finding was interpreted as evidence of transpositions of the chimera in the genome. PCR analyses detected the chimera in the genomes of chimpanzee and gorilla, but not in that of orangutan. These data demonstrate that the chimera appeared in the primate germ cells more than 10 million years ago, before divergence of the human/chimpanzee and the gorilla lineages. The combination of the two LTRs forms a new regulatory system that can be involved in nearby gene expression.

High resolution mapping DNAs by R-loop atomic force microscopy

R-loops formed by short RNA transcripts have been imaged by atomic force microscopy (AFM) at a constant force in the height mode. The technique was applied to mapping the human endogenous retrovirus K10 family (HERV-K10) long terminal repeats (LTR) within individual plasmids and cosmids. RNA probes specific for the U3 (384 nt) and U5 (375 nt) LTR regions separated by a span of 200 bp were used for R-loop formation with LTRs located within plasmid (3.8 kb) or cosmid ( approximately 40 kb) DNAs. R-loops stabilized by glyoxal treatment and adsorbed onto the mica surface in the presence of magnesium ions looked like looped out segments of RNA:DNA hybrids. The total yield of R-loops was usually approximately 95%. The RNA:DNA hybrids were found to be 12-15% shorter than the corresponding DNA:DNA duplex. The two regions of the LTR could be easily discerned in the AFM images as clearly separated loops. R-loop positions determined on cosmids by AFM were accurate to approximately 0.5% of the cosmid length. This technique might be easily adapted for mapping various sequences such as gene exons or regulatory regions and for detecting insertions, deletions and rearrangements that cause human genetic diseases.

[Long terminal repeat of the human endogenous retrovirus HERV-K in the intron of the ZNF91 gene]

The relative positions of the ZNF91 gene (exon-intron organization) and the HERV-K human endogenous virus long terminal repeat (LTR), which was earlier found to be located in the ZNF91 gene locus, on the EcoRI restriction metric map of the human chromosome 19 were determined with a high resolution. The direction of the ZNF91 gene transcription relative to the chromosome 19 telomeres was determined. The HERV-K LTR was localized to the ZNF91 gene intron 19. The role of retroviral sequences in the evolution of the ZNF91 gene family is discussed.

Long terminal repeats of human endogenous retrovirus K family (HERV-K) specifically bind host cell nuclear proteins

Solitary long terminal repeats (LTRs) of the human endogenous retroviruses, scattered in several thousand copies throughout the human genome, are potentially capable of affecting the expression of closely located genes. To assess their regulatory potential, the LTR sequences of one of the most abundant HERV families (HERV-K) were screened for the presence of binding sites for the host cell nuclear factors using mobility shift and UV-crosslinking assays. It was shown that the LTR sequences of two subfamilies harbor a specific binding site for a complex consisting of at least three proteins, ERF1, ERF2 and ERF3 of 98, 91 and 88 kDa apparent molecular mass, respectively. This binding site is located in the 5' region of the LTR U3 element. The preservation of the specific protein binding site in different HERV-K LTR sequences suggests their possible role in regulation of nearby located genes.