{'fraction_non_alphanumeric': 0.03741280913126189, 'fraction_numerical': 0.01077996195307546, 'mean_word_length': 4.725952813067151, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Null models for gene enrichment in plasmids

The recent article by Che et al. (1) presents an analysis of a dataset of 14,029 plasmids. I enjoyed the paper but have some concerns about its analysis. While it is clear that more conjugative plasmids carry antimicrobial resistance (AMR) genes, it is unclear whether this is simply because conjugative plasmids tend to carry more genes in general. Specifically, Che et al. aim to determine whether the AMR gene pool is “homogeneously distributed” or “constrained” by plasmid categories. To do this, they perform a Fisher’s exact test of the proportion of AMR gene-bearing plasmids between plasmid categories. In my view, this test is insufficient by itself for investigating “enrichment”—a term which suggests differential processes of horizontal gene transfer (HGT) between categories. To highlight the problem with an analogy, consider two types of dice trial: one where a die is rolled 10 times and one where it is rolled only twice. Trials of the first type will have more 6s. It is wrong to conclude, from this observation alone, that trials with 6s are “enriched” in the first type; one should also control for the number of rolls. [Che et al.’s (1) subsequent analysis of enrichment of AMR gene types in sets of genes within plasmid categories is not subject to this problem.] For the case in hand, assuming independence of association in Fisher’s exact test is equivalent to assuming that plasmids of different sizes have equal proportions of AMR gene-bearing plasmids. But there are better null models of homogeneous distribution by HGT. For example, one could assume AMR genes are Poisson distributed, occurring at a mean rate λ per base pair (bp) irrespective of plasmid category. Assuming 13,899 AMR genes distributed across the total length of all plasmids (~1.414 Gbp) gives λ = 9.8× 10−6 per bp. A simulation shows Fisher’s exact test gives a highly significant result for conjugative plasmids (Table 1), which might lead to the conclusion that “AMR gene-bearing plasmids were significantly enriched in this category” (1). However, we know the AMR genes have been homogeneously distributed; a greater proportion of conjugative plasmids are AMR gene bearing simply due to their greater length. Table 1 shows the Poisson model is poor (probably a zero-inflated Poissonmodel would be a better choice). But we can still use it as a guide for thinking about odds ratios (ORs). Since the model produces an OR of around three, that suggests (e.g.) OR > 3.5 as a more appropriate threshold thanOR> 1. In fact, we observeOR> 5, so there is a proportional difference for conjugative plasmids that remains unexplained. By this rough rule of thumb, Che et al.’s (1) original conclusion stands. Nevertheless, I raise this general point because I think language should match analysis. For example, elsewhere, Che et al. (1) describe detection of the same AMR gene on two plasmids as “when an AMR gene was transferred between two plasmids” without further analysis—I think this must be a mistake. But this is not just semantics. When investigating enrichment, considering null models carefully is important.