Abstract

Shotgun metagenomic sequencing is a powerful approach to study microbiomes in an unbiased manner and of increasing relevance for identifying novel enzymatic functions. However, the potential of metagenomics to relate from microbiome composition to function has thus far been underutilized. Here, we introduce the Metagenomics Genome-Phenome Association (MetaGPA) study framework, which allows linking genetic information in metagenomes with a dedicated functional phenotype. We applied MetaGPA to identify enzymes associated with cytosine modifications in environmental samples. From the 2365 genes that met our significance criteria, we confirm known pathways for cytosine modifications and proposed novel cytosine-modifying mechanisms. Specifically, we characterized and identified a novel nucleic acid modifying enzyme, 5-hydroxymethylcytosine carbamoyltransferase, that catalyzes the formation of a previously unknown cytosine modification, 5-carbamoyloxymethylcytosine, in DNA and RNA. Our work introduces MetaGPA as a novel and versatile tool for advancing functional metagenomics.

Data availability

All raw and processed sequencing data generated in this study have been submitted to the NCBI Sequence Read Archive (SRA; https://www.ncbi.nlm.nih.gov/sra) under accession number PRJNA714147.

The following data sets were generated

Article and author information

Author details

  1. Weiwei Yang

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Weiwei Yang, The author is employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  2. Yu-Cheng Lin

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Yu-Cheng Lin, The author was an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  3. William Johnson

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    William Johnson, The author was an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  4. Nan Dai

    RNA Biology, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Nan Dai, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  5. Romualdas Vaisvila

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Romualdas Vaisvila, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  6. Peter Weigele

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Peter Weigele, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  7. Yan-Jiun Lee

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Yan-Jiun Lee, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  8. Ivan R Corrêa Jr

    RNA Biology, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Ivan R Corrêa, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3169-6878
  9. Ira Schildkraut

    Research department, New England Biolabs Inc, Ipswich, United States
    Competing interests
    Ira Schildkraut, The author is an employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
  10. Laurence Ettwiller

    Research department, New England Biolabs Inc, Ipswich, United States
    For correspondence
    laurence.ettwiller@gmail.com
    Competing interests
    Laurence Ettwiller, The author is employee of New England Biolabs Inc. a manufacturer of restriction enzymes and molecular reagents..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3957-6539

Funding

New England Biolabs (no data)

  • Weiwei Yang
  • Yu-Cheng Lin
  • William Johnson
  • Nan Dai
  • Romualdas Vaisvila
  • Peter Weigele
  • Yan-Jiun Lee
  • Ivan R Corrêa Jr
  • Ira Schildkraut
  • Laurence Ettwiller

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. María Mercedes Zambrano, CorpoGen, Colombia

Version history

  1. Preprint posted: March 23, 2021 (view preprint)
  2. Received: May 4, 2021
  3. Accepted: November 5, 2021
  4. Accepted Manuscript published: November 8, 2021 (version 1)
  5. Version of Record published: December 14, 2021 (version 2)

Copyright

© 2021, Yang et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,608
    views
  • 207
    downloads
  • 4
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Weiwei Yang
  2. Yu-Cheng Lin
  3. William Johnson
  4. Nan Dai
  5. Romualdas Vaisvila
  6. Peter Weigele
  7. Yan-Jiun Lee
  8. Ivan R Corrêa Jr
  9. Ira Schildkraut
  10. Laurence Ettwiller
(2021)
A Genome-Phenome Association study in native microbiomes identifies a mechanism for cytosine modification in DNA and RNA
eLife 10:e70021.
https://doi.org/10.7554/eLife.70021

Share this article

https://doi.org/10.7554/eLife.70021

Further reading

    1. Genetics and Genomics
    Gbolahan Bamgbose, Guillaume Bordet ... Alexei Tulin
    Research Article

    PARP-1 is central to transcriptional regulation under both normal and stress conditions, with the governing mechanisms yet to be fully understood. Our biochemical and ChIP-seq-based analyses showed that PARP-1 binds specifically to active histone marks, particularly H4K20me1. We found that H4K20me1 plays a critical role in facilitating PARP-1 binding and the regulation of PARP-1-dependent loci during both development and heat shock stress. Here, we report that the sole H4K20 mono-methylase, pr-set7, and parp-1 Drosophila mutants undergo developmental arrest. RNA-seq analysis showed an absolute correlation between PR-SET7- and PARP-1-dependent loci expression, confirming co-regulation during developmental phases. PARP-1 and PR-SET7 are both essential for activating hsp70 and other heat shock genes during heat stress, with a notable increase of H4K20me1 at their gene body. Mutating pr-set7 disrupts monomethylation of H4K20 along heat shock loci and abolish PARP-1 binding there. These data strongly suggest that H4 monomethylation is a key triggering point in PARP-1 dependent processes in chromatin.

    1. Cancer Biology
    2. Genetics and Genomics
    Ting Zhang, Alisa Ambrodji ... Steven M Offer
    Research Article

    Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.