# Construction of a complete set of Neisseria meningitidis mutants and its use for the phenotypic profiling of this human pathogen

A. Muir, I. Gurung, A. Cehovin, A. Bazin, D. Vallenet, V. Pelicic

► To cite this version:

A. Muir, I. Gurung, A. Cehovin, A. Bazin, D. Vallenet, et al. Construction of a complete set of Neisseria meningitidis mutants and its use for the phenotypic profiling of this human pathogen. Nature Communications, 2020, 11 (1), pp.5541. 10.1038/s41467-020-19347-y. hal-03512357

HAL Id: hal-03512357

https://hal.science/hal-03512357

Submitted on 13 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-020-19347-y

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Construction of a complete set of *Neisseria meningitidis* mutants and its use for the phenotypic profiling of this human pathogen

Alastair Muir<fnref n="1" />, Ishwori Gurung<fnref n="1" />, Ana Cehovin<fnref n="1" />, Adelme Bazin ID ², David Vallenet ID ² & Vladimir Pelicic ID <fnref n="1" />✉

The bacterium *Neisseria meningitidis* causes life-threatening meningitis and sepsis. Here, we construct a complete collection of defined mutants in protein-coding genes of this organism, identifying all genes that are essential under laboratory conditions. The collection, named NeMeSys 2.0, consists of individual mutants in 1584 non-essential genes. We identify 391 essential genes, which are associated with basic functions such as expression and preservation of genome information, cell membrane structure and function, and metabolism. We use this collection to shed light on the functions of diverse genes, including a gene encoding a member of a previously unrecognised class of histidinol-phosphatases; a set of 20 genes required for type IV pili function; and several conditionally essential genes encoding anti-toxins and/or immunity proteins. We expect that NeMeSys 2.0 will facilitate the phenotypic profiling of a major human bacterial pathogen.

<footnote n="1">¹MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, UK. ²LABGeM, Génomique Métabolique, CEA, Genoscope, Institut François Jacob, Université d'Evry, Université Paris-Saclay, CNRS, Evry, France. ✉email: v.pelicic@imperial.ac.uk</footnote>
Low-cost massively parallel sequencing methods have transformed biology, with the availability of complete genome sequences increasing exponentially. In bacteria, genome sequences can be determined for hundreds of isolates in parallel in just a matter of days. For bacterial pathogens, genome sequences are often available for thousands of clinical isolates and represent a key resource for epidemiology and comparative genomics¹. This explosion has also led to the identification of millions of new genes, with the important insight that many are of unknown function. Even in a model species such as *Escherichia coli* K-12, perhaps the most thoroughly studied biological entity, only 60% of 4623 protein-coding genes have experimentally verified functions², while 35% have either no predicted function or a predicted function that is yet to be verified². In less well-studied species, this partition is dramatically skewed towards the latter class, with an overwhelming majority of genes having predicted but not experimentally verified function or no predicted function at all. This indicates that there are still major gaps in our understanding of basic bacterial physiology, including well-studied metabolic pathways³. This has fundamental consequences for systems and synthetic biology, by hindering holistic understanding of the life-supporting functions and processes allowing bacteria to grow and replicate. It also has practical consequences by slowing down the identification of new means for controlling bacterial pathogens, which are needed more than ever in this era of antimicrobial resistance. Therefore, developing methods for systematically elucidating gene function is a critical biological endeavour.

Molecular genetics, which involves the creation and phenotypic analysis of large collections of mutants, is the most direct path to determining gene function on a genome scale by providing a crucial link between phenotype and genotype. In many bacteria, genome-wide mutant collections in protein-coding genes have been constructed using either targeted mutagenesis or random transposon (Tn) mutagenesis. Because of its time and cost-effectiveness, Tn mutagenesis is more widely used⁴⁻⁶, especially since it has been coupled with massively parallel sequencing of Tn insertion sites (Tn-Seq)⁷, which allows simultaneous monitoring of the fitness in pools of hundreds of thousands of mutants⁸,⁹. However, Tn-Seq is not without drawbacks. (1) Not all intragenic Tn insertions lead to a loss of function. (2) Some properties cannot be assessed in a pooled assay format. (3) Specific mutants cannot be recovered from pools for further individual study. (4) Libraries of Tn mutant are never complete, no matter how large. In contrast, complete collections of mutants constructed by targeted mutagenesis do not suffer from these limitations, but their construction is laborious and expensive. For this reason, such collections are available for only a handful of bacterial species: the model Gram-negative *E. coli* K-12 (Gammaproteobacteria)¹⁰, the model Gram-positive *Bacillus subtilis* (Firmicute)¹¹, *Acinetobacter baylyi* (Gammaproteobacteria) studied for its metabolic properties and possible biotechnological applications¹², and *Streptococcus sanguinis* (Firmicute) an opportunistic pathogen causing infective endocarditis¹³. The utility of such resources is best exemplified by a large-scale study using the Keio collection of ~4000 mutants in *E. coli* K-12, which identified mutants presenting altered colony size upon growth in more than 300 different conditions¹⁴. More than 10,000 phenotypes were identified for ~2000 mutants, including many in genes of unknown function¹⁴. Another major advantage of mutant collections constructed by targeted mutagenesis is that they also reveal complete lists of essential genes, those that cannot be mutated because their loss is lethal. These genes that encode proteins essential for cellular life, many of which are expected to be common to all living organisms¹⁵, have attracted particular interest because they are expected to shed light on the origin of life. Additionally, essential

genes specific to bacteria are also interesting because they are prime targets for the design of new drugs to tackle antimicrobial resistance.

Since the few species in which complete collections of mutants have been constructed represent only a tiny fraction of bacterial diversity, additional complete collections of mutants in diverse species are needed to tackle more effectively the challenge of genes of unknown function. Human bacterial pathogens are particularly attractive candidates for this venture because of the potential practical implications for human health. Among these, *Neisseria meningitidis*—a Gram-negative Betaproteobacteria causing life-threatening meningitis and sepsis—stands out as an ideal candidate for several reasons. (1) It has a small genome with ~2000 protein-coding genes, and yet displays a robust metabolism allowing rapid growth (doubling time of 40 min), even on minimal medium. (2) It has been the subject of intensive investigation for decades and there are several thousand publicly available meningococcal genome sequences. (3) It is naturally competent, which makes it a model for molecular genetics. Previously, we have used these advantages to design a modular toolbox for *Neisseria meningitidis* systematic analysis of gene function, named NeMeSys¹⁶. The central module of NeMeSys was an ordered library of ~4500 Tn mutants in which the Tn insertion sites were sequence-defined¹⁶,¹⁷, revealing insertions in more than 900 genes¹⁶.

In the present study, we have (1) used the original library of Tn mutants as a starting point to generate a complete collection of defined mutants in *N. meningitidis*—named NeMeSys 2.0— comprising one mutant in each non-essential protein-coding gene, (2) defined and analysed the essential meningococcal genome and (3) illustrated the potential of NeMeSys 2.0 for the phenotypic profiling of the meningococcal genome by identifying the function of multiple genes of unknown function.

## Results and discussion

**Construction of the complete NeMeSys 2.0 collection of meningococcal mutants.** The present study gave new impetus to a systematic re-annotation of the genome of *N. meningitidis* 8013 published in 2009¹⁶. We took into account new information in the public databases, new publications, an RNA-Seq analysis performed in 8013¹⁸ and findings in the present study. In brief, we dropped one previously annotated gene, updated the gene product annotations of 521 genes, added 69 non-coding RNAs and changed 176 gene start sites. The genome of *N. meningitidis* 8013 contains 2060 protein-coding genes (labelled NMV_.). In parallel, we re-sequenced the genome of the 8013 wild-type (WT) strain by Illumina whole-genome sequencing (WGS), which confirmed the accuracy of the original genome sequence determined by Sanger sequencing more than 11 years ago¹⁶, with an error rate of ~1 per 134,000 bases. We identified 17 differences with the published sequence, which are likely to be genuine since they were also found in the few mutants that were verified by WGS (see below). Of these differences, 11 correspond to single-nucleotide polymorphisms, five are single-nucleotide indels and one is the deletion of a GT dinucleotide within a repeated GT tract in NMV_0677, which therefore corresponds to a phase variation event¹⁹.

In keeping with previous studies describing the construction of complete collections of mutants¹⁰⁻¹³, 85 genes (4.1%) were not considered to be meaningful targets for mutagenesis (Fig. 1) because they encode transposases for repeated insertion sequences, or they correspond to short remnants of truncated genes and/or non-expressed cassettes (Supplementary Data 1). We, therefore, set out to assemble a complete collection of mutants in the remaining 1975 genes (95.9%) following a two-step
<figure><img src="image_4.png" /><figcaption>**Fig. 1 Flowchart of the construction of the NeMeSys 2.0 complete collection of mutants in *N. meningitidis* 8013.** As for similar efforts in other bacteria¹⁰,¹²,¹³, we first selected protein-coding genes to be targeted by systematic mutagenesis, excluding 85 genes (4.1%, highlighted in black in the first pie chart) because they encode transposases of repeated insertion sequences, or correspond to short remnants of truncated genes or cassettes (Supplementary Data 1). We then followed a two-step mutagenesis approach explained in the text and in Supplementary Fig. 1. In brief, we first selected a subset of sequence-verified Tn mutants from a previously constructed arrayed library¹⁶,¹⁷. Mutations were re-transformed in strain 8013 and PCR-verified. We thus selected 801 Tn mutants with a disrupting transposon in the corresponding target genes (highlighted in black in the second and third pie charts). Next, we systematically mutagenised the remaining 1174 target genes using a validated no-cloning mutagenesis method relying on sPCR²⁰. For each successful transformation, two colonies were isolated and PCR-verified. To minimise false-positive identification of essential genes, each transformation that yielded no transformants was repeated at least three times. In total, we could construct an additional 783 mutants (highlighted in grey in the third pie chart), generating an ordered library of defined mutants in 1584 meningococcal genes (Supplementary Data 3). This effort also identified 391 candidate essential genes, which could not be disrupted, encoding proteins required for *N. meningitidis* growth on rich medium (Supplementary Data 4).</figcaption></figure>

approach described below. In brief, we first selected a subset of suitable Tn mutants from the NeMeSys library, and then systematically mutagenised the remaining target genes by allelic exchange (Supplementary Fig. 1).

First, since we previously constructed an arrayed library of Tn mutants in strain 8013 in which Tn insertion sites were sequence-defined¹⁶,¹⁷, we selected a subset of suitable mutants from this library. Specifically, mutations with a Tn insertion closer to the centre of the genes (excluding the first and last 15% that are often not gene-inactivating) (Supplementary Fig. 1) were first re-transformed into strain 8013, which is naturally competent. After extraction of the genomic DNA of the transformants, we PCR-verified that they contained a Tn inserted in the expected target gene using suitable pairs of flanking primers (Supplementary Data 2). We thus selected a subset of 801 Tn mutants (40.6% of target genes), harbouring Tn insertions expected to be gene-inactivating (Fig. 1) (Supplementary Data 3). As a control, Illumina WGS of one mutant, in *nlaB*, confirmed that it contained the expected mutation, a Tn inserted after a TA dinucleotide at position 2,157,750 in the genome.

Next, we set out to systematically mutagenise the remaining 1174 target genes (59.4%) by allelic exchange. For this, we used a previously validated no-cloning mutagenesis method²⁰, based on splicing PCR (sPCR), to construct non-polar mutants in which a central portion of the target gene (between the first and last 30%) would be deleted and replaced by the same kanamycin (Km) resistance cassette present in the above Tn (Supplementary Fig. 1). In brief, three PCR products were first amplified, corresponding to the Km cassette and regions upstream and downstream target genes. For failed reactions, we undertook as many rounds of primer design as needed, until all PCRs were successful. The three PCR products were then combined and spliced together. The final sPCR product was transformed directly into *N. meningitidis* 8013. Typically, for non-essential genes, we obtained hundreds of colonies resistant to kanamycin (Km⁴). For

each successful transformation, two Km⁴ colonies were isolated and PCR-verified to contain the expected mutation. To minimise false-positive identification of essential genes, each transformation that yielded no transformants was repeated at least three times with different sPCR products. Only when all transformations failed to yield transformants, was the target gene deemed to be essential. In summary (Fig. 1), out of the 1174 genes that were targeted, we could disrupt 783 (Supplementary Data 3), while 391 could not be mutated and are thus essential for growth on rich medium (Supplementary Data 4). As above, WGS of a few mutants constructed by sPCR, in *Int*, *tatB* and *secB*, confirmed that each contained the expected mutation.

Taken together, using the above two-step approach, we constructed a comprehensive set of 1584 defined mutants in strain 8013, named the NeMeSys 2.0 collection. In addition, we also identified 391 candidate essential genes (19% of all protein-coding genes), which are required for growth of the meningococcus on rich medium. For each gene that was successfully mutated, one mutant was stored at -80 °C in glycerol, while the corresponding genomic DNA was stored at -20 °C. This allows for easy distribution of mutants to the community and/or re-transformation of the corresponding mutations in 8013, other meningococcal strains, or even genetically closely related species like the gonococcus, which share most of their genome with the meningococcus.

**Analysis of the meningococcal essential genome.** Considering their fundamental and practical importance, we first focused on the 391 genes essential for meningococcal life. Since corresponding proteins are expected to be broadly conserved, we determined the partition of the 391 essential genes between conserved and variable genomes in meningococci. To do this, we first used the recently described PPanGGOLiN software²¹—an expectation-maximisation algorithm based on multivariate Bernoulli mixture model coupled with a Markov random field—to
compute the pangenome of *N. meningitidis*, based on complete genome sequences of 108 meningococcal isolates available in RefSeq (Supplementary Data 5). We thereby classified genes in three categories²¹, persistent (gene families present in almost all genomes), shell (present at intermediate frequencies), or cloud (present at low frequency). We then determined how the 2060 protein-coding genes of *N. meningitidis* 8013 partitioned between these three classes (Supplementary Data 5). Only 1664 genes (80.8%) belong to the persistent genome, while 396 genes (19.2%) correspond to gene families present at intermediate/low frequency in meningococci since they are part of the shell (241 genes) and cloud genomes (155 genes) (Fig. 2a). Critically, besides highlighting the well-known genomic plasticity of this naturally competent species, this analysis confirmed our original prediction by revealing that essential genes are overwhelmingly conserved in meningococcal isolates. Indeed, 382 essential genes (97.4%) are part of the persistent meningococcal genome (Fig. 2b and Supplementary Data 5). Of the remaining nine essential genes, six are part of the shell genome, while three belong to the cloud genome (Supplementary Data 5).

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Partition of essential meningococcal genes into persistent, shell and cloud genomes. To perform this analysis, we used the PPanGGOLiN²¹ method as explained in the text. **a** Partition of the 2060 genes in the genome of *N. meningitidis* 8013: persistent (not highlighted), shell (highlighted in grey), and cloud (highlighted in black). The corresponding datasets are listed in Supplementary Data 5. **b** Partition of the subset of 391 essential genes identified in this study (same colour code than in **a**). The corresponding datasets are listed in Supplementary Data 5.</figcaption></figure>

Furthermore, many meningococcal essential genes (30.2%) are conserved and essential in other bacteria where systematic mutagenesis efforts have been performed (Fig. 3a), including in phylogenetically distant species such as *S. sanguinis* (Supplementary Data 6)¹⁰,¹²,¹³. When conservation was assessed with Proteobacteria only (*E. coli* and *A. baylyi*), which are more closely related to *N. meningitidis*, half of the 391 essential meningococcal genes (195 in total) are conserved and essential in these three Gram-negative species (Fig. 3a). Interestingly, a similar proportion of meningococcal essential genes (175 in total) was conserved in JCVI-syn3.0 (Fig. 3b) (Supplementary Data 7), a synthetic *Mycoplasma mycoides* bacterium engineered with a minimal genome²² smaller than that of any naturally occurring autonomously replicating cell. Taken together, these observations are consistent with the notion that a sizeable portion of the genome essential for cellular life is widely conserved in bacteria.

Next, to have a more global understanding of the biological functions that are essential in the meningococcus, we analysed the distribution of the 391 essential genes in functional categories (Supplementary Data 8). Genes were classified in specific functional categories using the bioinformatic tools embedded in the MicroScope platform²³, which hosts NeMeSys 2.0. In particular, we used predictions from MultiFun²⁴, MetaCyc²⁵, eggNOG²⁶, COG²⁷, FIGfam²⁸ and/or InterProScan²⁹. As can be seen in Table 1, the essential meningococcal genes could be distributed in a surprisingly limited number of pathways. In keeping with expectations, many essential genes are involved in key cellular processes such as (1) transcription (11 genes encoding subunits of the RNA polymerase, a series of factors modulating transcription, and two transcriptional regulators), (2) RNA modification/degradation (19 genes) and (3) protein biosynthesis (120 genes) (Table 1). The large class of genes involved in protein biosynthesis encode (1) all 24 proteins involved in tRNA charging, (2) virtually all ribosomal proteins (51/56), (3) five proteins involved in ribosome biogenesis/ maturation, (4) 10 factors modulating translation, (5) two enzymes introducing post-translational modifications in proteins, (6) seven factors facilitating protein folding and (7) 21 factors involved in protein export (signal peptidase, Lol system involved in lipoprotein export, general secretion Sec system, Tat translocation system and TAM translocation and assembly module specific for autotransporters) (Table 1). Furthermore, many other essential genes encode proteins involved in additional key processes such as (1) genome replication and maintenance

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Comparison of *N. meningitidis* essential genome to that of other bacteria. **a** Edwards-Venn diagram displaying overlaps between essential genes in bacteria in which complete libraries of mutants have been constructed: *N. meningitidis* (black line), *E. coli* (green line)¹⁰, *A. baylyi* (red line)¹² and *S. sanguinis* (blue line)¹³. To perform this analysis, we queried the DEG database⁶⁵ of essential genes using our set of essential genes. The corresponding datasets are listed in Supplementary Data 6. **b** Comparison of *N. meningitidis* essential genome (black line) to JCVI-syn3.0 (red line), a synthetic *M. mycoides* designed with a minimal genome²². The corresponding datasets are listed in Supplementary Data 7.</figcaption></figure>
<table><caption>**Table 1 Essential meningococcal genes listed by functional category/sub-category.**</caption><thead><tr><th>Functional category/sub-category</th><th>Number of genes</th><th>% of total</th></tr></thead><tbody><tr><td>Gene/protein expression</td><td>150</td><td>38.4</td></tr><tr><td>Ribosomal protein</td><td>51</td><td>13</td></tr><tr><td>Aminoacyl-tRNA synthesis</td><td>24</td><td>6.1</td></tr><tr><td>Protein export</td><td>21</td><td>5.4</td></tr><tr><td>RNA modification/ degradation</td><td>19</td><td>4.9</td></tr><tr><td>DNA transcription</td><td>11</td><td>2.8</td></tr><tr><td>Translation factors</td><td>10</td><td>2.6</td></tr><tr><td>Protein folding</td><td>7</td><td>1.8</td></tr><tr><td>Ribosome biogenesis/ maturation</td><td>5</td><td>1.3</td></tr><tr><td>Protein post-translational modification</td><td>2</td><td>0.5</td></tr><tr><td>Genome/cell replication</td><td>33</td><td>8.4</td></tr><tr><td>Chromosome replication/ maintenance</td><td>20</td><td>5.1</td></tr><tr><td>Cell division</td><td>13</td><td>3.3</td></tr><tr><td>Cell membrane/wall</td><td>54</td><td>13.8</td></tr><tr><td>biogenesis</td><td></td><td></td></tr><tr><td>Peptidoglycan biosynthesis/ recycling</td><td>17</td><td>4.3</td></tr><tr><td>LOS biosynthesis</td><td>14</td><td>3.6</td></tr><tr><td>Fatty acid biosynthesis</td><td>12</td><td>3.1</td></tr><tr><td>Phospholipids biosynthesis</td><td>8</td><td>2</td></tr><tr><td>UDP-GlcNAc biosynthesis</td><td>3</td><td>0.8</td></tr><tr><td>Cytosolic metabolism</td><td>120</td><td>30.7</td></tr><tr><td>Vitamins biosynthesis</td><td>24</td><td>6.1</td></tr><tr><td>Nucleotide biosynthesis</td><td>15</td><td>3.8</td></tr><tr><td>Aerobic respiration/ cytochrome c</td><td>12</td><td>3.1</td></tr><tr><td>Amino acid metabolism</td><td>9</td><td>2.3</td></tr><tr><td>Isoprenoid biosynthesis</td><td>9</td><td>2.3</td></tr><tr><td>Electron transport</td><td>9</td><td>2.3</td></tr><tr><td>Heme biosynthesis</td><td>9</td><td>2.3</td></tr><tr><td>Sugar metabolism</td><td>9</td><td>2.3</td></tr><tr><td>Fe-S cluster biosynthesis</td><td>6</td><td>1.5</td></tr><tr><td>CoA/ac-CoA biosynthesis</td><td>6</td><td>1.5</td></tr><tr><td>Ubiquinol biosynthesis</td><td>4</td><td>1</td></tr><tr><td>NAD/NADP biosynthesis</td><td>3</td><td>0.8</td></tr><tr><td>Lipoate biosynthesis</td><td>2</td><td>0.5</td></tr><tr><td>SAM biosynthesis</td><td>1</td><td>0.3</td></tr><tr><td>Phosphate metabolism</td><td>1</td><td>0.3</td></tr><tr><td>Sulfur metabolism</td><td>1</td><td>0.3</td></tr><tr><td>Unassigned</td><td>34</td><td>8.7</td></tr><tr><td>Total</td><td>391</td><td>100</td></tr></tbody></table>

For this classification, we used predictions from MultiFun²⁴, MetaCyc²⁵, eggNOG²⁶, COG²⁷, FIGfam²⁸ and/or InterProScan²⁹. The corresponding datasets are listed in Supplementary Data 8.

(20 proteins including subunits of the DNA polymerase, topoisomerases, gyrase, ligase etc.), (2) cell division (13 proteins), (3) peptidoglycan cell wall biogenesis (17 proteins) and (4) membrane biogenesis (37 proteins) (Table 1). The latter class includes all 12 proteins required for fatty acid biosynthesis, eight proteins required for the biosynthesis of phospholipids (phosphatidylethanolamine (PE) and phosphatidylglycerol (PG)), and 17 proteins involved in the biosynthesis and export of the lipopolysaccharide (LOS), which constitutes the external leaflet of the outer membrane. Finally, most of the remaining genes of known function (Table 1) are involved in energy generation, production of key metabolic intermediates (such as dihydroxyacetone phosphate (DHAP), or 5-phosphoribosyl diphosphate (PRPP)) and a variety of metabolic pathways leading to the biosynthesis of

(1) vitamins (B1, B2, B6 and B9), (2) nucleotides, (3) amino acids (notably meso-diaminopimelate (M-DAP) a component of the peptidoglycan), (4) isoprenoid, (5) heme, (6) CoA/ac-CoA, (7) NAD/NADP, (8) iron-sulfur (Fe-S) clusters, (9) ubiquinol, (10) lipoate and (11) S-adenosyl-methionine (SAM). Many of these compounds are cofactors and/or coenzymes known to be essential for the activity of many bacterial enzymes.

When essential genes were further integrated into networks and metabolic pathways—using primarily MetaCyc²⁵—the above picture became dramatically simpler. As many as 91.3% of the meningococcal essential genes partitioned into just four major functional categories (Supplementary Fig. 2). Namely, (1) gene/protein expression, (2) genome/cell replication, (3) membrane/cell wall biogenesis and (4) cytosolic metabolism. Strikingly, for most multi-step enzymatic pathways, most, and often all, the corresponding genes were identified as essential (Fig. 4), although they are most often scattered throughout the genome. We view this as important quality control of our library. Only 34 essential genes (8.7%) could not be clearly assigned to these four categories. Critically, this functional partition is coherent with the one previously determined for the minimal synthetic bacterium JCVI-syn3.0²². As a result, we were able to generate a concise overview of the meningococcal essential genome in the context of the cell (Fig. 4), in which many of the above pathways are detailed and often linked by critical metabolic intermediates. This overview provides a useful blueprint for systems and synthetic biology and a global understanding of meningococcal biology.

**Essential genes that are part of the variable genome are conditionally essential.** A puzzling finding was that some essential genes, which could not be assigned to the above four functional groups, are part of shell/cloud genomes present at intermediate/low frequency in the meningococcus (Supplementary Data 5). Interestingly, many of the corresponding genes are located in regions of genome plasticity (RGP)—often thought to have been acquired by horizontal gene transfer—as confirmed by the recently described predictive method panRGP³⁰ that identifies RGP using pangenome data generated by PPanGGOLiN²¹. A panRGP analysis identified 32 RGP in the genome of strain 8013, encompassing a total of 348,885 bp (15.3% of the genome) (Supplementary Data 9). Interestingly, one of these genes offered a plausible explanation to the apparent paradox of essential genes that are not part of the persistent meningococcal genome. NMV_2289 is predicted to encode a DNA-methyltransferase, based on the presence of a DNA methylase N-4/N-6 domain (InterPro IPR002941). Since the neighbouring NMV_2288 is predicted to encode a restriction enzyme, the probable role of the NMV_2289 methyltransferase is to protect meningococcal DNA against degradation by this restriction enzyme. This would explain the lethal phenotype of the ΔNMV_2289 mutant. A closer examination of the other essential genes that are part of shell/cloud genomes and often found in RGP suggests that a similar scenario might be often applicable. We noticed that NMV_1478, which encodes a protein of unknown function in RGP_2, is likely to be co-transcribed with NMV_1479 that is predicted to encode the toxin of a toxin-antitoxin (TA) system³¹ (Fig. 5a). Therefore, since a toxin and its cognate antitoxin are often encoded by closely linked genes, we hypothesised that NMV_1478 might be the antitoxin for the neighbouring NMV_1479 toxin. This would explain why NMV_1478 is essential, i.e. in its absence, the NMV_1479 toxin would kill the cell. If this is true, we reasoned that it should be possible to delete NMV_1478 together with NMV_1479, which was attempted. As predicted, in contrast to NMV_1478 which could not be deleted individually despite
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Concise cellular overview of the essential meningococcal genome. Essential genes were integrated into networks and metabolic pathways using primarily MetaCyc²⁵. The four basic functional groups are highlighted using the same colour code as in Supplementary Fig. 2, i.e. orange (gene/protein expression), green (genome/cell replication), blue (cell membrane/wall biogenesis) and yellow (cytosolic metabolism). The 34 essential genes that could not be clearly assigned to one of these four categories are not represented on the figure. Genes involved in the different reactions are indicated by their name or NMV\_label, in red when essential, in black when dispensable. \*Genes involved in more than one pathway. Key compounds/proteins are abbreviated as follows. ACP acyl carrier protein, AIR aminoimidazole ribotide, CDP-DAG CDP-diacylglycerol, CMP-Kdo CMP-ketodeoxyoctonate, DHAP dihydroxyacetone phosphate, DPP dimethylallyl diphosphate, DXP 1-deoxyxylulose-5P, Fe-S iron-sulfur, FMN flavin mononucleotide, G3P glyceraldehyde-3P, Glu-6P glucose-6P, IMP inosine monophosphate, IPP isopentenyl diphosphate, LOS lipo-oligosaccharide, M-DAP meso-diaminopimelate, NaMN nicotinate d-ribonucleotide, OPP all-trans-octaprenyl diphosphate, PE phosphatidylethanolamine, PG phosphatidylglycerol, PHBA p-hydroxybenzoate, PRPP 5-phosphoribosyl diphosphate, Rib-5P ribulose-5P, SAM S-adenosyl-methionine, UDP-GlcNAc UDP-N-acetyl-glucosamine, UPP di-trans-poly-cis-undecaprenyl diphosphate. The corresponding datasets are listed in Supplementary Data 8.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Essential genes in RGP—NMV_1479 and NMV_0559—are conditionally essential. **a** Gene organisation of the putative antitoxin NMV_1478 (highlighted in red) with its neighbouring NMV_1479 toxin. Results of the mutagenesis (represented by scissors) are shown. Viable mutant (represented by a grey diplococcus); lethal phenotype (represented by a skull). **b** Gene organisation of the *tps* RGP (RGP_0) to which NMV_0559 (highlighted in red) belongs (Supplementary Data 9). Genes on the + strand are in white, genes on the—strand are in black. Results of the mutagenesis (scissors) are shown. Viable mutant (grey diplococcus); lethal phenotype (skull). In contrast to NMV_0559, each of the other target genes in the *tps* RGP could be mutated individually (not shown for readability).</figcaption></figure>

repeated attempts, a double deletion mutant ΔNMV_1478/1479 could be readily obtained (Fig. 5a). This confirms that the NMV_1478 gene is essential only when the toxin product of NMV_1479 is also present. Such a scenario was also tested for NMV_0559, which is found in a large RGP of 86.8 kbp (RGP_0),

encompassing genes NMV_0527 to NMV_0618/0619 (Fig. 5b). We named this RGP *tps* because it contains multiple *tpsA* genes predicted to encode large haemagglutinin/haemolysin-related proteins of two-partner secretion systems³². Since NMV_0559 is the only essential gene in the *tps* RGP, we tested whether it might
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 6** NMV_1317 encodes a novel histidinol-phosphatase. **a** Histidine biosynthesis pathway in the meningococcus, with genes highlighted in red. PRPP 5-phosphoribosyl diphosphate, PRFAR phosphoribulosylformimino-AICAR-P, IGP erythro-imidazole-glycerol-P. **b** Growth on M9 minimal medium, with or without added histidine (His). The plates also contained 0.5 mM IPTG for inducing expression of the complementing genes. WT, strain 8013; Δ1317, ΔNMV_1317 mutant; Δ1317::1317, ΔNMV_1317 complemented with NMV_1317; Δ1317::hisBEC, ΔNMV_1317 cross-complemented with hisBEC from E. coli, which encodes the unrelated histidinol-phosphatase present in this species; Δ1718, ΔNMV_1718 mutant; Δ1718::hisBEC, control showing that ΔNMV_1718 (hisH) cannot be cross-complemented by hisBEC. Source data are provided as a Source Data file.</figcaption></figure>

be conditionally essential by attempting to delete most of RGP_0. As predicted, while NMV_0559 could not be deleted on its own despite repeated attempts, an 80 kbp Δtps deletion in RGP_0 encompassing NMV_0559 could be readily obtained (Fig. 5b). This suggests that NMV_0559 is essential only in the presence of another gene in RGP_0, which remains to be identified. We predict that several other essential genes of unknown function that are not part of the persistent meningococcal genome may be similarly conditionally essential, including (1) NMV_1305 and NMV_1310 part of a putative prophage (RGP_1), (2) NMV_1541 the putative antitoxin for neighbouring NMV_1539 toxin, (3) NMV_1544 putative antitoxin for neighbouring NMV_1543 toxin, (4) NMV_1761 putative immunity protein against neighbouring MafB toxins (RGP_4), which represent a family of secreted toxins in pathogenic *Neisseria* species³³, (5) NMV_1918 part of a second tps RGP (RGP_7), (6) NMV_2010 the putative immunity protein against neighbouring bacteriocins (RGP_3) and (7) NMV_2333 and NMV_2335 putative immunity proteins in a second *maf* RGP (RGP_11).

**Filling holes in metabolic pathways: identification of a novel and widespread histidinol-phosphatase.** Another major utility of the NeMeSys 2.0 collection of mutants is to determine the function of genes of unknown function by reverse genetics. We, therefore, wished to illustrate this aspect. We noticed that the 10-step metabolic pathway leading to histidine biosynthesis presents an apparent “hole” (Fig. 6a). The “missing” enzyme in 8013 corresponds to histidinol-phosphatase (EC 3.1.3.15), which catalyses the dephosphorylation of histidinol-P into histidinol (Fig. 6a). The meningococcal genome does not have homologues of known histidinol-phosphatases³⁴, a feature it shares with many other bacterial species according to MetaCyc²⁵. First, we excluded the possibility that the meningococcus might be auxotrophic for histidine by showing that 8013 can grow on M9 minimal medium without added histidine (Fig. 6b). Next, we sought to identify the unknown histidinol-phosphatase using NeMeSys 2.0, by specifically testing mutants in genes encoding putative phosphatases of unknown function, in search of a mutant that would not grow on M9 without added histidine. Using this approach, we identified

ΔNMV_1317 as a histidine auxotroph, growing on M9 plates only when histidine was added (Fig. 6b). The product of this gene was annotated as a putative hydrolase of unknown function, belonging to the IB subfamily of the haloacid dehalogenase superfamily of aspartate-nucleophile hydrolases (IPR006385). The corresponding sequences include a variety of phosphatases none of which is annotated as a histidinol-phosphatase. Complementation of ΔNMV_1317 mutant with NMV_1317 restored growth on M9 without added histidine (Fig. 6b), confirming that the auxotrophic phenotype was due to the mutation in NMV_1317. To confirm that NMV_1317 encodes the missing enzyme in histidine biosynthesis, we performed a cross-species complementation assay with the *hisB* gene from *E. coli* DH5α (*hisBEC*), which encodes a histidinol-phosphatase unrelated to NMV_1317. Complementation of ΔNMV_1317 with *hisBEC* restored growth on M9 without added histidine (Fig. 6b). Complementation was gene-specific, since *hisBEC* could not complement the growth deficiency of a different histidine auxotrophic mutant in NMV_1778 (*hisH*) (Fig. 6b). Together, these results show that NMV_1317 defines a previously unrecognised class of histidinol-phosphatases, allowing us to fill a hole in the histidine biosynthesis pathway in the meningococcus. This finding, which illustrates the utility of NeMeSys 2.0 for annotating genes of unknown function by reverse genetics, has implications for other species lacking known histidinol-phosphatases. Indeed, NMV_1317 homologues are widespread in Betaproteobacteria and Gammaproteobacteria in which no histidinol-phosphatase has been identified according to MetaCyc²⁵. We expect that NMV_1317 homologues will encode the elusive histidinol-phosphatase in a variety of Burkholderiales and Pseudomonadales, including many pathogenic species of *Bordetella*, *Pseudomonas* and *Burkholderia* (there is some evidence for this in *Burkholderia phytofirmans*³⁵).

**Identification of a comprehensive set of genes involved in type IV pilus biology in N. meningitidis.** Another utility of NeMeSys 2.0 we wished to illustrate is the possibility to perform whole-genome phenotypic screens to identify all the genes responsible for a phenotype of interest. We chose to focus on type IV pili
<table><caption>**Table 2 Mutants in the complete NeMeSys 2.0 library of meningococcal mutants affected for two functions mediated by T4P: aggregation and twitching motility.**</caption><thead><tr><th>Mutated gene</th><th>Product</th><th>Aggregation</th><th>Twitching</th></tr></thead><tbody><tr><td colspan="4">Genes previously known to be involved in T4P biology in <i>N. meningitidis</i></td></tr><tr><td>pilE</td><td>Major pilin PilE</td><td>—</td><td>N/A</td></tr><tr><td>pilT</td><td>Type IV pilus retraction ATPase PilT</td><td>+ (Irregular)</td><td>—</td></tr><tr><td>pilF</td><td>Type IV pilus extension ATPase PilF</td><td>—</td><td>N/A</td></tr><tr><td>pilD</td><td>Leader peptidase/N-methyltransferase PilD</td><td>—</td><td>N/A</td></tr><tr><td>pilG</td><td>Type IV pilus biogenesis protein PilG</td><td>—</td><td>N/A</td></tr><tr><td>pilW</td><td>Type IV pilus biogenesis lipoprotein PilW</td><td>—</td><td>N/A</td></tr><tr><td>pilX</td><td>Minor pilin PilX</td><td>—</td><td>N/A</td></tr><tr><td>pilK</td><td>Type IV pilus biogenesis protein PilK</td><td>—</td><td>N/A</td></tr><tr><td>pilJ</td><td>Type IV pilus biogenesis protein PilJ</td><td>—</td><td>N/A</td></tr><tr><td>pilI</td><td>Type IV pilus biogenesis protein PilI</td><td>—</td><td>N/A</td></tr><tr><td>pilH</td><td>Type IV pilus biogenesis protein PilH</td><td>—</td><td>N/A</td></tr><tr><td>pilZ</td><td>PilZ protein</td><td>—</td><td>N/A</td></tr><tr><td>pilM</td><td>Type IV pilus biogenesis protein PilM</td><td>—</td><td>N/A</td></tr><tr><td>pilN</td><td>Type IV pilus biogenesis protein PilN</td><td>—</td><td>N/A</td></tr><tr><td>pilO</td><td>Type IV pilus biogenesis protein PilO</td><td>—</td><td>N/A</td></tr><tr><td>pilP</td><td>Type IV pilus biogenesis lipoprotein PilP</td><td>—</td><td>N/A</td></tr><tr><td>pilQ</td><td>Type IV pilus secretin PilQ</td><td>—</td><td>N/A</td></tr><tr><td colspan="4">Genes not previously known to be involved in T4P biology in <i>N. meningitidis</i></td></tr><tr><td>tsaP</td><td>Secretin-associated protein TsaP</td><td>—</td><td>N/A</td></tr><tr><td>NMV_1205</td><td>Conserved hypothetical periplasmic protein</td><td>—</td><td>N/A</td></tr><tr><td>NMV_2228</td><td>Putative carbonic anhydrase</td><td>+/-</td><td>N/A</td></tr></tbody></table>

The other 1569 mutants can form aggregates and exhibit twitching motility. N/A not assayable because twitching motility could only be assessed in mutants forming aggregates. The corresponding datasets are listed in Supplementary Data 10.

(T4P), which are pivotal virulence factors in the meningococcus³⁶, for two reasons. First, T4P, are ubiquitous in prokar-yotes, which has made them a hot topic for research for the past 40 years³⁷. Second, although many aspects of T4P biology remain incompletely understood, most genes composing the multi-protein machinery involved in the assembly of these filaments and/or their multiple functions have been identified, including in the meningococcus that is one of the mainstream T4P models³⁶. Hence, we could readily benchmark the results of our screen against previous mutational analyses, especially for the recovery rate of expected mutants. The two T4P-linked phenotypes we decided to study are the formation of bacterial aggregates and twitching motility. Both phenotypes can be simultaneously assessed—allowing thus a dual screen—by observing the mutants growing in liquid medium by phase-contrast microscopy³⁸. We thus scored the 1589 mutants for the presence of round aggregates, and for the continuous and vigorous jerky movements of cells within these aggregates, which corresponds to twitching motility³⁸ (Supplementary Data 10). This analysis revealed that 20 mutants (1.2%) present phenotypic defects in these two T4P-linked phenotypes (Table 2). Significantly, we identified 100% of the expected single mutants in 17 *pil* genes known to affect these phenotypes in the meningococcus³⁸⁻⁴⁰. This is another quality control of the NeMeSys 2.0 collection of mutants, which confirms excellent correlation between phenotype and genotype, and shows that when a robust screening method is used all the genes involved in a given phenotype can be identified. Furthermore, we could readily identify mutants such as Δ*pilT* able to form (irregular) aggregates, but showing no twitching motility³⁸. PilT is known to encode the motor powering pilus retraction, which is directly responsible for twitching motility⁴¹. No other mutant in the library exhibited a similar phenotype. Strikingly, we also identified mutants in three genes not previously associated with T4P biology in *N. meningitidis* (Table 2), which was unexpected considering that T4P has been studied for decades in this species. TsaP has been shown in the closely related pathogen Neisseria gonorrhoeae to interact with the secretin PilQ⁴², which forms a

pore in the outer membrane through which T4P translocate onto the cell surface. TsaP has a poorly understood role in T4P biology in *N. gonorrhoeae* and *Myxococcus xanthus*⁴², but has apparently no role in *Pseudomonas aeruginosa*⁴³. In addition, we identified NMV_1205 and NMV_2228, which to the best of our knowledge have never been previously linked to T4P biology. NMV_1205 is predicted to encode a periplasmic protein of unknown function and is found mainly in Neisseriales, where it is widespread. NMV_2228 is predicted to encode a carbonic anhydrase, an enzyme catalysing the reversible hydration of carbon dioxide⁴⁴.

Next, we analysed the three new mutants in detail using an approach previously validated in the meningococcus with many *pil* genes³⁸,⁴⁵. We first tested whether piliation was affected in the mutants by purifying T4P using a procedure in which filaments sheared by vortexing are precipitated using ammonium sulfate³⁸. Pilus preparations, obtained from equivalent numbers of cells, were separated by SDS-PAGE and either stained using Coomassie blue, or tested by immunoblotting using an anti-PilE antibody. This revealed the major pilin PilE as a 17 kDa species (Fig. 7a) and confirmed that all three mutants were piliated. However, when compared to the WT, pilus yields differed between the mutants. While piliation in ΔNMV_2228 was apparently normal, it was decreased in ΔtsaP and ΔNMV_1205. The decrease was dramatic in ΔNMV_1205, where filaments could be detected only by immunoblotting (Fig. 7a). Quantification of piliation using a whole-cell ELISA procedure³⁸ confirmed these findings (Fig. 7b) and showed that piliation levels were 40, 31 and 154% of WT, in ΔtsaP, ΔNMV_1205 and ΔNMV_2228, respectively. Therefore, none of the corresponding proteins is required for T4P biogenesis.

We then further characterised the new mutants for T4P-linked functions. Since it is known that in some meningococcal mutants piliation and/or aggregation defects can be restored when filament retraction is abolished by a concurrent mutation in *pilT*³⁸,⁴⁵, we tested if that might be the case for the new mutants. First, we showed that aggregation, which is abolished in ΔtsaP and ΔNMV_1205 and dramatically affected in ΔNMV_2228,
<figure><img src="image_11.png" /><figcaption>**Fig. 7 Assaying piliation in the mutants in genes not previously associated with T4P biology in *N. meningitidis*.** The WT strain and a non-piliated Δ*pilD* mutant were included as positive and negative controls, respectively. Δ*pisaP*, Δ*pisaP* mutant; Δ1205, Δ*NMV*_1205 mutant; Δ2228, Δ*NMV*_2228 mutant. **a** T4P purified using a shearing/precipitation method were separated by SDS-PAGE and either stained with Coomassie blue (upper panel) or analysed by immunoblotting using an antibody against the major pilin PilE (lower panel)⁵⁹. Samples were prepared from equivalent numbers of cells and identical volumes were loaded in each lane. MW molecular weight marker lane, with values in kDa. Source data are provided as a Source Data file. **b** T4P were quantified by whole-cell ELISA using a monoclonal antibody specific for the filaments of strain 8013⁵⁸. Equivalent numbers of cells, based on OD₆₀₀ readings, were applied to the wells of microtiter plates. Results are expressed in % piliation (ratio to WT) and are the average ± standard deviations from five independent experiments. Source data are provided as a Source Data file. For statistical analysis, one-way ANOVA followed by Dunnett's multiple comparison tests were performed (****P < 0.0001).</figcaption></figure>

could be restored when mutants were complemented with the corresponding WT alleles (Fig. 8a). This confirmed that the phenotypic defects in these mutants were indeed due to the mutations in the above three genes. Similarly, some aggregation was restored in Δ*pisaP*/Δ*pilT* and Δ*NMV* _1205/Δ*pilT* (Δ*NMV* _2228/Δ*pilT* could not be constructed), which harboured a concurrent mutation in *pilT* that encodes the retraction motor PilT. However, aggregates were morphologically distinct from the round aggregates formed by the WT or the irregular aggregates formed by Δ*pilT*. These findings suggest that *isaP* and *NMV* _1205 participate in the formation of aggregates by counterbalancing PilT-mediated pilus retraction. Lastly, we checked whether the mutants were affected for another T4P-linked phenotype, competence for DNA uptake, which makes the meningococcus naturally transformable. We quantified competence in the three mutants as described⁴⁶ by transformation to rifampicin resistance (Fig. 8b). We found that Δ*pisaP* was almost as transformable as the WT, Δ*NMV* _1205 showed a 26-fold decrease in transformation, while competence was completely abolished in Δ*NMV* _2228 (explaining why the Δ*NMV* _2228/Δ*pilT* double mutant could not be constructed). These findings show that *isaP*, *NMV* _1205 and *NMV* _2228 contribute to the fine-tuning of multiple functions mediated by T4P in *N. meningitidis*.

Taken together, these findings are of significance in two ways. The identification of three new genes playing a role in T4P biology (how exactly remains to be understood), will contribute to a better understanding of these filaments. This has broad implications because T4P and T4P-related filamentous nanomachines are ubiquitous in Bacteria and Archaea⁴⁷. In addition, the identification of new genes contributing to phenotypes that have been extensively studied in several species for the past 30 years is unambiguous evidence of the potential of NeMeSys 2.0 to lead to a global phenotypic profiling of the meningococcal genome.

**Concluding remarks.** Here, we describe the construction of a complete collection of defined mutants in *N. meningitidis*, which has been fully integrated in our modular NeMeSys toolbox¹⁶, accessible online in *MicroScope*²³. Furthermore, we illustrate

NeMeSys 2.0 utility for tackling the challenge of genes of unknown function, following a variety of approaches. Although they were limited at this stage to a few different biological properties, these experiments provide significant information on meningococcal biology, with possible implications for many other species. They also clearly demonstrate the potential of NeMeSys 2.0 to make a significant contribution to ongoing efforts directed towards a comprehensive understanding of a bacterial cell. Such potential is further amplified by several useful properties of *N. meningitidis* such as (1) its small genome with limited functional redundancy, which in species with larger genomes can obscure the link between phenotype and genotype, (2) its hardy nature with a robust metabolism allowing it to grow on minimal medium, (3) its taxonomy making it the first Betaproteobacteria in which a complete library of mutants is available and (4) the fact that it is a major human pathogen, which allows addressing virulence properties absent in non-pathogenic model species. This latter point makes the meningococcus a prime model for the identification of new means for controlling bacterial pathogens, which would have practical implications for human health.

The other achievement in the present study is the identification of the essential meningococcal genome, comprising 391 genes that could not be disrupted. The finding that more than 90% of these genes are associated with just four basic biological functions, helped us generate a coherent concise cellular overview of the meningococcal essential metabolism. We surmise that most, if not all, of the remaining essential genes of unknown function, which are conserved in other bacterial species, will be involved in either expression of genome information, preservation of genome information, cell membrane structure/function, or cytosolic metabolism (in particular the production of key metabolic intermediates, vitamins and crucial cofactors/coenzymes). Our cellular overview provides a useful starting point for systems biology, which together with previous genome-scale metabolic network models⁴⁸,⁴⁹, could help us progress towards a global understanding of meningococcal biology. This, in turn, would have widespread implications, because many of the meningococcal genes, especially the essential ones, are widely conserved.
a

<figure><img src="image_12.png" /><figcaption>**Fig. 8** Functional analysis of the mutants in genes not previously associated with T4P biology in *N. meningitidis*. The WT strain and a non-piliated *pilD* mutant were included as positive and negative controls, respectively. **a** Aggregation in liquid culture as assessed by phase-contrast microscopy³⁸. Δ*tsaP*::*tsaP*, Δ*tsaP* complemented with *tsaP*; Δ*tsaP*/Δ*pilT*, double mutant in *tsaP* and *pilT*; Δ1025::1205, ΔNMV_1205 complemented with NMV_1205; Δ1205/Δ*pilT*, double mutant in NMV_1205 and *pilT*; Δ2228::2228, ΔNMV_2228 complemented with NMV_2228. Scale bar, 20 µm. Source data are provided as a Source Data file. **b** Quantification of the competence for DNA transformation. Equivalent numbers of recipient cells were transformed using a rpoB PCR product containing a point mutation leading to rifampicin resistance. Results are expressed as transformation frequencies and are the average ± standard deviations from four independent experiments. Source data are provided as a Source Data file. For statistical analysis, one-way ANOVA followed by Dunnett's multiple comparison tests were performed (****P < 0.0001).</figcaption></figure>

b

<figure><img src="image_13.png" /></figure>

Methods

**Strains and growth conditions.** All the *N. meningitidis* strains that were generated and used in this study were derived from a highly adhesive variant—sometimes called clone 12 or 2C43—of the clinical isolate 8013⁵⁰. This serogroup C strain, which belongs to the sequence type 177 and clonal complex ST-18, has been previously sequenced¹⁶. Meningococci were routinely grown on plates with GC agar Base (GCB) (Difco) containing 5 g/l agar (all chemicals were from Sigma unless stated otherwise) and Kellogg's supplements (4 g/l glucose, 0.59 µM thiamine hydrochloride, 12.37 µM Fe(NO₃)₃·9H₂O, 68.4 µM l-glutamine). Plates were incubated overnight (O/N) at 37 °C in a moist atmosphere containing 5% CO₂. Alternatively, we used agar plates with M9 minimal medium (4 g/l glucose, 4.78 mM Na₂HPO₄, 2.2 mM KH₂PO₄, 1.87 mM NH₄Cl, 8.56 mM NaCl, 2 mM MgSO₄, 0.1 mM CaCl₂). When required, plates contained 100 µg/ml kanamycin, 3 µg/ml erythromycin, 5 µg/ml rifampicin, 25 µg/ml l-histidine and/or 0.5 mM isopropyl-β-D-thiogalactopyranoside (IPTG) (Merck). Strains were stored at -80 °C in 10% glycerol in liquid GC (15 g/l protease peptone No. 3, 23 mM K₂HPO₄, 7.34 mM KH₂PO₄, 85.6 mM NaCl). *E. coli* TOP10 and DH5α were grown at 37 °C in liquid or solid lysogeny broth (LB), which contained 100 µg/ml spectinomycin or 50 µg/ml kanamycin, when appropriate.

**Construction of strains.** Genomic DNA from *N. meningitidis* and *E. coli* strains were prepared using the Wizard Genomic DNA Purification kit (Promega) following the manufacturer's instructions. Plasmid DNA from *E. coli* strains was purified using the QIAprep Spin Miniprep Kit (Qiagen) following the manufacturer's instructions. Bacteria were transformed as follows. *N. meningitidis* is naturally competent for transformation. A loopful of bacteria grown on GCB plates was resuspended in liquid GC containing 5 mM MgCl₂ (GC transfo), 200 µl was aliquoted in the well of a 24-well plate, and DNA was added. After incubation for 30 min at 37 °C on an orbital shaker, 0.8 ml GC transfo was added to the wells and the plates were further incubated for 3 h at 37 °C, without shaking. Transformants were selected on GCB plates containing suitable antibiotics. For the transformation of *E. coli*, ultra-competent cells were prepared as described elsewhere⁵¹ and transformed by a standard heat shock procedure⁵². Transformants were selected on LB plates containing the suitable antibiotic.

The construction of the NeMeSys 2.0 library of meningococcal mutants followed a two-step procedure explained in the Results section (Supplementary Fig. 1). First, we selected a subset of potentially suitable mutants from an archived library of Tn mutants with sequence-defined Tn insertion sites, which was described previously¹⁶,¹⁷. The corresponding genomic DNAs were PCR-verified
queries (by BLAST or keyword searches) and download datasets/annotations in a variety of formats (including EMBL and GenBank).

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

**Data availability**

Whole-genome sequencing data that support the findings in this study have been deposited in the European Nucleotide Archive under accession number PRJEB39197. All the datasets generated during this study are either included in this paper and its Supplementary Information files, or available in MicroScope (http://mage.genoscope.cns.fr/microscope/mage/viewer.php?O_id=99). The RefSeq database (http://www.ncbi.nlm.nih.gov/refseq/) was used to access annotations of publicly available genomes of *N. meningitidis*. Lists of essential genes in other bacteria were obtained from the DEG repository (http://tubic.tju.edu.cn/deg/). Source data are provided with this paper.

Received: 27 July 2020; Accepted: 5 October 2020;
Published online: 02 November 2020

**References**

1. Jolley, K. A. & Maiden, M. C. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics **11**, 595 (2010).

2. Ghatak, S., King, Z. A., Sastry, A. & Palsson, B. O. The y-ome defines the 35% of *Escherichia coli* genes that lack experimental evidence of function. *Nucleic Acids Res.* **47**, 2446–2454 (2019).

3. Galperin, M. Y. & Koonin, E. V. From complete genome sequence to 'complete' understanding? *Trends Biotechnol.* **28**, 398–406 (2010).

4. Akerley, B. J. et al. A genome-scale analysis for identification of genes required for growth or survival of *Haemophilus influenzae*. *Proc. Natl Acad. Sci. USA* **99**, 966–971 (2002).

5. Liberati, N. T. et al. An ordered, nonredundant library of *Pseudomonas aeruginosa* strain PA14 transposon insertion mutants. *Proc. Natl Acad. Sci. USA* **103**, 2833–2838 (2006).

6. Deutschbauer, A. et al. Evidence-based annotation of gene function in *Shewanella oneidensis* MR-1 using genome-wide fitness profiling across 121 conditions. *PLoS Genet.* **7**, e1002385 (2011).

7. van Oprijnen, T., Lazinski, D. W. & Camilli, A. Genome-wide fitness and genetic interactions determined by Tn-seq, a high-throughput massively parallel sequencing method for microorganisms. *Curr. Protoc. Mol. Biol.* **106**, 7.16.11–7.16.24 (2014).

8. Langridge, G. C. et al. Simultaneous assay of every *Salmonella* Typhi gene using one million transposon mutants. *Genome Res.* **19**, 2308–2316 (2009).

9. Christen, B. et al. The essential genome of a bacterium. *Mol. Syst. Biol.* **7**, 528 (2011).

10. Baba, T. et al. Construction of *Escherichia coli* K-12 in-frame, single-gene knockout mutants: the Keio collection. *Mol. Syst. Biol.* **2**, 2006.0008 (2006).

11. Kobayashi, K. et al. Essential *Bacillus subtilis* genes. *Proc. Natl Acad. Sci. USA* **100**, 4678–4683 (2003).

12. de Berardinis, V. et al. A complete collection of single-gene deletion mutants of *Acinetobacter baylyi* ADP1. *Mol. Syst. Biol.* **4**, 174 (2008).

13. Xu, P. et al. Genome-wide essential gene identification in *Streptococcus sanguinis*. *Sci. Rep.* **1**, 125 (2011).

14. Nichols, R. J. et al. Phenotypic landscape of a bacterial cell. *Cell* **144**, 143–156 (2011).

15. Mushegian, A. R. & Koonin, E. V. A minimal gene set for cellular life derived by comparison of complete bacterial genomes. *Proc. Natl Acad. Sci. USA* **93**, 10268–10273 (1996).

16. Rusniok, C. et al. NeMeSys: a resource for narrowing the gap between sequence and function in the human pathogen *Neisseria meningitidis*. *Genome Biol.* **10**, R110 (2009).

17. Geoffroy, M. C., Floquet, S., Métais, A., Nassif, X. & Pelicic, V. Large-scale analysis of the meningococcus genome by gene disruption: resistance to complement-mediated lysis. *Genome Res.* **13**, 391–398 (2003).

18. Heidrich, N. et al. The primary transcriptome of *Neisseria meningitidis* and its interaction with the RNA chaperone Hfq. *Nucleic Acids Res.* **45**, 6147–6167 (2017).

19. Martin, P. et al. Experimentally revised repertoire of putative contingency loci in *Neisseria meningitidis* strain MC58: evidence for a novel mechanism of phase variation. *Mol. Microbiol.* **50**, 245–257 (2003).

20. Georgiadou, M., Castagnini, M., Karimova, G., Ladant, D. & Pelicic, V. Large-scale study of the interactions between proteins involved in type IV pilus biology in *Neisseria meningitidis*: characterization of a subcomplex involved in pilus assembly. *Mol. Microbiol.* **84**, 857–873 (2012).

21. Gautreau, G. et al. PpanGGOLiN: depicting microbial diversity via a partitioned pangenome graph. *PLoS Comput. Biol.* **16**, e1007732 (2020).

22. Hutchison, C. A. et al. Design and synthesis of a minimal bacterial genome. *Science* **351**, aa6253 (2016).

23. Vallenet, D. et al. MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis. *Nucleic Acids Res.* **48**, D579–D589 (2020).

24. Serres, M. H. & Riley, M. MultiFun, a multifunctional classification scheme for *Escherichia coli* K-12 gene products. *Microb. Comp. Genomics* **5**, 205–222 (2000).

25. Caspi, R. et al. The MetaCyc database of metabolic pathways and enzymes - a 2019 update. *Nucleic Acids Res.* **48**, D445–D453 (2020).

26. Huerta-Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. *Nucleic Acids Res.* **47**, D309–D314 (2019).

27. Galperin, M. Y., Makarova, K. S., Wolf, Y. I. & Koonin, E. V. Expanded microbial genome coverage and improved protein family annotation in the COG database. *Nucleic Acids Res.* **43**, D261–D269 (2015).

28. Meyer, F., Overbeek, R. & Rodriguez, A. FIGfams: yet another set of protein families. *Nucleic Acids Res.* **37**, 6643–6654 (2009).

29. Jones, P. et al. InterProScan 5: genome-scale protein function classification. *Bioinformatics* **30**, 1236–1240 (2014).

30. Bazin, A., Gautreau, G., Médigue, C., Vallenet, D. & Calteau, A. panRGP: a pangenome-based method to predict genomic islands and explore their diversity. *Bioinformatics* https://doi.org/10.1093/bioinformatics/btaa792 (2020).

31. Yamaguchi, Y., Park, J. H. & Inouye, M. Toxin-antitoxin systems in bacteria and archaea. *Annu. Rev. Genet.* **45**, 61–79 (2011).

32. van Ulsen, P., Rutten, L., Feller, M., Tommassen, J. & van der Ende, A. Two-partner secretion systems of *Neisseria meningitidis* associated with invasive clonal complexes. *Infect. Immun.* **76**, 4649–4658 (2008).

33. Jamet, A. et al. A new family of secreted toxins in pathogenic *Neisseria* species. *PLoS Pathog.* **11**, e1004592 (2015).

34. Kulis-Horn, R. K., Ruckert, C., Kalinowski, J. & Persicke, M. Sequence-based identification of inositol monophosphatase-like histidinol-phosphate phosphatases (HisN) in *Corynebacterium glutamicum*, *Actinobacteria*, and beyond. *BMC Microbiol.* **17**, 161 (2017).

35. Price, M. N. et al. Mutant phenotypes for thousands of bacterial genes of unknown function. *Nature* **557**, 503–509 (2018).

36. Pelicic, V. Type IV pili: *e pluribus unum? Mol. Microbiol.* **68**, 827–837 (2008).

37. Berry, J. L. & Pelicic, V. Exceptionally widespread nano-machines composed of type IV pilins: the prokaryotic Swiss Army knives. *FEMS Microbiol. Rev.* **39**, 134–154 (2015).

38. Brown, D., Helaine, S., Carbonnelle, E. & Pelicic, V. Systematic functional analysis reveals that a set of 7 genes is involved in fine tuning of the multiple functions mediated by type IV pili in *Neisseria meningitidis*. *Infect. Immun.* **78**, 3053–3063 (2010).

39. Carbonnelle, E., Helaine, S., Prouvensier, L., Nassif, X. & Pelicic, V. Type IV pilus biogenesis in *Neisseria meningitidis*: PilW is involved in a step occurring after pilus assembly, essential for fiber stability and function. *Mol. Microbiol.* **55**, 54–64 (2005).

40. Helaine, S. et al. PilX, a pilus-associated protein essential for bacterial aggregation, is a key to pilus-facilitated attachment of *Neisseria meningitidis* to human cells. *Mol. Microbiol.* **55**, 65–77 (2005).

41. Merz, A. J., So, M. & Sheetz, M. P. Pilus retraction powers bacterial twitching motility. *Nature* **407**, 98–102 (2000).

42. Siewering, K. et al. Peptidoglycan-binding protein TsAP functions in surface assembly of type IV pili. *Proc. Natl Acad. Sci. USA* **111**, E953–E961 (2014).

43. Koo, J., Lamers, R. P., Rubinstein, J. L., Burrows, L. L. & Howell, P. L. Structure of the *Pseudomonas aeruginosa* type IVa pilus secretin at 7.4 Å. *Structure* **24**, 1778–1787 (2016).

44. Smith, K. S. & Ferry, J. G. Prokaryotic carbonic anhydrases. *FEMS Microbiol. Rev.* **24**, 335–366 (2000).

45. Carbonnelle, E., Helaine, S., Nassif, X. & Pelicic, V. A systematic genetic analysis in *Neisseria meningitidis* defines the Pil proteins required for assembly, functionality, stabilization and export of type IV pili. *Mol. Microbiol.* **61**, 1510–1522 (2006).

46. Berry, J. L., Cehovin, A., McDowell, M. A., Lea, S. M. & Pelicic, V. Functional analysis of the interdependence between DNA uptake sequence and its cognate ComP receptor during natural transformation in *Neisseria* species. *PLoS Genet.* **9**, e1004014 (2013).

47. Denise, R., Abby, S. S. & Rocha, E. P. C. Diversification of the type IV filament superfamily into machines for adhesion, protein secretion, DNA uptake, and motility. *PLoS Biol.* **17**, e3000390 (2019).

48. Baart, G. J. et al. Modeling *Neisseria meningitidis* metabolism: from genome to metabolic fluxes. *Genome Biol.* **8**, R136 (2007).

49. Mendum, T. A., Newcombe, J., Mannan, A. A., Kierzek, A. M. & McFadden, J. Interrogation of global mutagenesis data with a genome scale model of
*Neisseria meningitidis* to assess gene fitness *in vitro* and in sera. *Genome Biol.* **12**, R127 (2011).

50. Nassif, X. et al. Antigenic variation of pilin regulates adhesion of *Neisseria meningitidis* to human epithelial cells. *Mol. Microbiol.* **8**, 719–725 (1993).

51. Inoue, H., Nojima, H. & Okayama, H. High efficiency transformation of *Escherichia coli* with plasmids. *Gene* **96**, 23–28 (1990).

52. Sambrook, J. & Russell, D. W. *Molecular Cloning: A Laboratory Manual* (Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York, 2001).

53. Untergasser, A. et al. Primer3 – new capabilities and interfaces. *Nucleic Acids Res.* **40**, e115 (2012).

54. Pelicic, V., Morelle, S., Lampe, D. & Nassif, X. Mutagenesis of *Neisseria meningitidis* by *in vitro* transposition of *Himar1 mariner*. *J. Bacteriol.* **182**, 5391–5398 (2000).

55. Metzgar, D. et al. *Acinetobacter* sp. ADP1: an ideal model organism for genetic analysis and genome engineering. *Nucleic Acids Res.* **32**, 5780–5790 (2004).

56. Pujol, C., Eugène, E., Marceau, M. & Nassif, X. The meningococcal PilT protein is required for induction of intimate attachment to epithelial cells following pilus-mediated adhesion. *Proc. Natl Acad. Sci. USA* **96**, 4017–4022 (1999).

57. Mehr, I. J., Long, C. D., Serkin, C. D. & Seifert, H. S. A homologue of the recombination-dependent growth gene, *rdgC*, is involved in gonococcal pilin antigenic variation. *Genetics* **154**, 523–532 (2000).

58. Pujol, C., Eugène, E., de Saint Martin, L. & Nassif, X. Interaction of *Neisseria meningitidis* with a polarized monolayer of epithelial cells. *Infect. Immun.* **65**, 4836–4842 (1997).

59. Morand, P. C. et al. Type IV pilus retraction in pathogenic *Neisseria* is regulated by the PilC proteins. *EMBO J.* **23**, 2009–2017 (2004).

60. Cehovin, A. et al. Specific DNA recognition mediated by a type IV pilin. *Proc. Natl Acad. Sci. USA* **110**, 3065–3070 (2013).

61. Carter, P. E., Abadi, F. J., Yakubu, D. E. & Pennington, T. H. Molecular characterization of rifampin-resistant *Neisseria meningitidis*. *Antimicrob. Agents Chemother.* **38**, 1256–1261 (1994).

62. Carver, T., Harris, S. R., Berriman, M., Parkhill, J. & McQuillan, J. A. Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. *Bioinformatics* **28**, 464–469 (2012).

63. Marck, C. 'DNA Strider': a 'C' program for the fast analysis of DNA and protein sequences on the Apple Macintosh family of computers. *Nucleic Acids Res.* **16**, 1829–1836 (1988).

64. Haft, D. H. et al. RefSeq: an update on prokaryotic genome annotation and curation. *Nucleic Acids Res.* **46**, D851–D860 (2018).

65. Luo, H., Lin, Y., Gao, F., Zhang, C. T. & Zhang, R. DEG 10, an update of the database of essential genes that includes both protein-coding genes and noncoding genomic elements. *Nucleic Acids Res.* **42**, D574–D580 (2014).

## Acknowledgements

This work was supported by funding from the Wellcome Trust (092290/Z/10/Z) to V.P. The France Génomique and French Bioinformatics Institute national infrastructures—

funded as part of Investissement d'Avenir programme managed by the Agence Nationale pour la Recherche (contracts ANR-10-INBS-09 and ANR-11-INBS-0013)—are acknowledged for support of the MicroScope annotation platform. We thank Angelika Gründling (Imperial College London) and Christoph Tang (University of Oxford) for critical reading of the manuscript.

## Author contributions

V.P. was responsible for conception and supervision of the work, interpretation of data and writing of the manuscript. A.M., I.G., A.C. and V.P. performed the experimental studies. A.B. and D.V. performed computational studies.

## Competing interests

The authors declare no competing interests.

## Additional information

*Supplementary information* is available for this paper at https://doi.org/10.1038/s41467-020-19347-y.

**Correspondence** and requests for materials should be addressed to V.P.

**Peer review information** *Nature Communications* thanks Adam Deutschbauer and Hank Seifert for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_14.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020