1. Beres, S. B. et al. Transcriptome remodeling contributes to epidemic disease caused by the human pathogen Streptococcus pyogenes. mBio 7, e00403-16 (2016).

2. Chewapreecha, C. et al. Comprehensive identification of single nucleotide polymorphisms associated with beta-lactam resistance within pneumococcal mosaic genes. PLoS Genet. 10, e1004547 (2014).

3. Fernandez-Romero, N. et al. Uncoupling between core genome and virulome in extraintestinal pathogenic Escherichia coli. Can. J. Microbiol. 61, 647–652 (2015).

4. Long, S. W. et al. Population genomic analysis of 1,777 extended-spectrum beta-lactamase-producing Klebsiella pneumoniae isolates, Houston, Texas: unexpected abundance of clonal group 307. mBio 8, e00489-17 (2017).

5. Mukherjee, S. et al. 1,003 reference genomes of bacterial and archaeal isolates expand coverage of the tree of life. Nat. Biotechnol. 35, 676–683 (2017).

6. Nasser, W. et al. Evolutionary pathway to increased virulence and epidemic group A Streptococcus disease derived from 3,615 genome sequences. Proc. Natl Acad. Sci. USA 111, E1768–E1776 (2014).

7. Bruchmann, S. et al. Deep transcriptome profiling of clinical Klebsiella pneumoniae isolates reveals strain and sequence type-specific adaptation. Environ. Microbiol. 17, 4690–4710 (2015).

8. Dotsch, A. et al. The Pseudomonas aeruginosa transcriptional landscape is shaped by environmental heterogeneity and genetic variation. mBio 6, e00749 (2015).

9. Sharma-Kuinkel, B. K. et al. Potential influence of Staphylococcus aureus clonal complex 30 genotype and transcriptome on hematogenous infections. Open Forum Infect. Dis. 2, ofv093 (2015).

10. Felek, S., Tsang, T. M. & Krukonis, E. S. Three Yersinia pestis adhesins facilitate Yop delivery to eukaryotic cells and contribute to plague virulence. Infect. Immun. 78, 4134–4150 (2010).

11. Swearingen, M. C., Porwollik, S., Desai, P. T., McClelland, M. & Ahmer, B. M. Virulence of 32 Salmonella strains in mice. PLoS One 7, e36043 (2012).

12. Schreiber, H. L. T. et al. Bacterial virulence phenotypes of Escherichia coli and host susceptibility determine risk for urinary tract infections.Sci. Transl. Med. 9, eaaf1283 (2017).

13. Carapetis, J. R., Steer, A. C., Mulholland, E. K. & Weber, M. The global burden of group A streptococcal diseases. Lancet Infect. Dis. 5, 685–694 (2005).

14. Carapetis, J. R. et al. Acute rheumatic fever and rheumatic heart disease. Nat. Rev. Dis. Primers 2, 15084 (2016).

15. Zhu, L. et al. A molecular trigger for intercontinental epidemics of group A Streptococcus. J. Clin. Invest. 125, 3545–3559 (2015).

16. Zhu, L., Olsen, R. J., Nasser, W., de la Riva Morales, I. & Musser, J. M. Trading capsule for increased cytotoxin production: contribution to virulence of a newly emerged clade of emm89 Streptococcus pyogenes. mBio 6, e01378-15 (2015).

17. Colman, G., Tanna, A., Efstratiou, A. & Gaworzewska, E. T. The serotypes of Streptococcus pyogenes present in Britain during 1980–1990 and their association with disease. J. Med. Microbiol. 39, 165–178 (1993).

18. Gherardi, G., Vitali, L. A. & Creti, R. Prevalent emm types among invasive GAS in Europe and North America since year 2000. Front. Public Health 6, 59 (2018).

19. Smit, P. W. et al. Epidemiology and emm types of invasive group A streptococcal infections in Finland, 2008–2013. Eur. J. Clin. Microbiol. Infect. Dis. 34, 2131–2136 (2015).

20. Ikebe, T. et al. Increased prevalence of group A Streptococcus isolates in streptococcal toxic shock syndrome cases in Japan from 2010 to 2012. Epidemiol. Infect. 143, 864–872 (2015).

21. Naseer, U., Steinbakk, M., Blystad, H. & Caugant, D. A. Epidemiology of invasive group A streptococcal infections in Norway 2010–2014: a retrospective cohort study.Eur. J. Clin. Microbiol. Infect. Dis. 35, 1639–1648 (2016).

22. Nelson, G. E. et al. Epidemiology of invasive group A streptococcal infections in the United States, 2005–2012. Clin. Infect. Dis. 63, 478–486 (2016).

23. Plainvert, C. et al. Invasive group A streptococcal infections in adults, France (2006–2010). Clin. Microbiol. Infect. 18, 702–710 (2012).

24. Al-Shahib, A. et al. Emergence of a novel lineage containing a prophage in emm/M3 group A Streptococcus associated with upsurge in invasive disease in the UK. Microb. Genom. 2, e000059 (2016).

25. Davies, M. R. et al. Emergence of scarlet fever Streptococcus pyogenes emm12 clones in Hong Kong is associated with toxin acquisition and multidrug resistance. Nat. Genet. 47, 84–87 (2015).

26. Fittipaldi, N. et al. Full-genome dissection of an epidemic of severe invasive disease caused by a hypervirulent, recently emerged clone of group A Streptococcus. Am. J. Pathol. 180, 1522–1534 (2012).

27. Hamilton, S. M., Stevens, D. L. & Bryant, A. E. Pregnancy-related group a streptococcal infections: temporal relationships between bacterial acquisition, infection onset, clinical findings, and outcome. Clin. Infect. Dis. 57, 870–876 (2013).

28. Johnson, D. R., Stevens, D. L. & Kaplan, E. L. Epidemiologic analysis of group A streptococcal serotypes associated with severe systemic infections, rheumatic fever, or uncomplicated pharyngitis. J. Infect. Dis. 166, 374–382 (1992).

29. Shea, P. R. et al. Group A Streptococcus emm gene types in pharyngeal isolates, Ontario, Canada, 2002–2010. Emerg. Infect. Dis. 17, 2010–2017 (2011).

30. Smoot, J. C. et al. Genome sequence and comparative microarray analysis of serotype M18 group A Streptococcus strains associated with acute rheumatic fever outbreaks. Proc. Natl Acad. Sci. USA 99, 4668–4673 (2002).

31. Ben Zakour, N. L., Venturini, C., Beatson, S. A. & Walker, M. J. Analysis of a Streptococcus pyogenes puerperal sepsis cluster by use of whole-genome sequencing. J. Clin. Microbiol. 50, 2224–2228 (2012).

32. Chuang, I., Van Beneden, C., Beall, B. & Schuchat, A. Population-based surveillance for postpartum invasive group A Streptococcus infections, 1995–2000. Clin. Infect. Dis. 35, 665–670 (2002).

33. Gaworzewska, E. & Colman, G. Changes in the pattern of infection caused by Streptococcus pyogenes. Epidemiol. Infect. 100, 257–269 (1988).

34. Raymond, J., Schlegel, L., Garnier, F. & Bouvet, A. Molecular characterization of Streptococcus pyogenes isolates to investigate an outbreak of puerperal sepsis. Infect. Control Hosp. Epidemiol. 26, 455–461 (2005).

35. Sims, D., Sudbery, I., Ilott, N. E., Heger, A. & Ponting, C. P. Sequencing depth and coverage: key considerations in genomic analyses. Nat. Rev. Genet. 15, 121–132 (2014).

36. Bricker, A. L., Carey, V. J. & Wessels, M. R. Role of NADase in virulence in experimental invasive group A streptococcal infection. Infect. Immun. 73, 6562–6566 (2005).

37. Bricker, A. L., Cywes, C., Ashbaugh, C. D. & Wessels, M. R. NAD+-glycohydrolase acts as an intracellular toxin to enhance the extracellular survival of group A streptococci. Mol. Microbiol. 44, 257–269 (2002).

38. Sumby, P. et al. Evolutionary origin and emergence of a highly successful clone of serotype M1 group A Streptococcus involved multiple horizontal gene transfer events. J. Infect. Dis. 192, 771–782 (2005).

39. Zhu, L. et al. Contribution of secreted NADase and streptolysin O to the pathogenesis of epidemic serotype M1 Streptococcus pyogenes infections. Am. J. Pathol. 187, 605–613 (2017).

40. Meehl, M. A., Pinkner, J. S., Anderson, P. J., Hultgren, S. J. & Caparon, M. G. A novel endogenous inhibitor of the secreted streptococcal NAD-glycohydrolase. PLoS Pathog. 1, e35 (2005).

41. Tatsuno, I. et al. Characterization of the NAD-glycohydrolase in streptococcal strains. Microbiology 153, 4253–4260 (2007).

42. Shimomura, Y. et al. Complete genome sequencing and analysis of a Lancefield group G Streptococcus dysgalactiae subsp. equisimilis strain causing streptococcal toxic shock syndrome (STSS). BMC Genomics 12, 17 (2011).

43. Carroll, R. K. et al. Naturally occurring single amino acid replacements in a regulatory protein alter streptococcal gene expression and virulence in mice. J. Clin. Invest. 121, 1956–1968 (2011).

44. Graham, M. R. et al. Virulence control in group A Streptococcus by a two-component gene regulatory system: global expression profiling and in vivo infection modeling. Proc. Natl Acad. Sci. USA 99, 13855–13860 (2002).

45. Ribardo, D. A. & McIver, K. S. Defining the Mga regulon: comparative transcriptome analysis reveals both direct and indirect regulation by Mga in the group A Streptococcus. Mol. Microbiol. 62, 491–508 (2006).

46. Ramalinga, A., Danger, J. L., Makthal, N., Kumaraswami, M. & Sumby, P. Multimerization of the virulence-enhancing group A Streptococcus transcription factor RivR is required for regulatory activity.J. Bacteriol. 199, e00452-16 (2017).

47. Trevino, J., Liu, Z., Cao, T. N., Ramirez-Pena, E. & Sumby, P. RivR is a negative regulator of virulence factor expression in group A Streptococcus. Infect. Immun. 81, 364–372 (2013).

48. Nyberg, P., Rasmussen, M. & Bjorck, L. α 2 -Macroglobulin-proteinase complexes protect Streptococcus pyogenes from killing by the antimicrobial peptide LL-37. J. Biol. Chem. 279, 52820–52823 (2004).

49. Rasmussen, M., Muller, H. P. & Bjorck, L. Protein GRAB of Streptococcus pyogenes regulates proteolysis at the bacterial surface by binding α 2 -macroglobulin. J. Biol. Chem. 274, 15336–15344 (1999).

50. Toppel, A. W., Rasmussen, M., Rohde, M., Medina, E. & Chhatwal, G. S. Contribution of protein G-related α 2 -macroglobulin-binding protein to bacterial virulence in a mouse skin model of group A streptococcal infection. J. Infect. Dis. 187, 1694–1703 (2003).

51. Haas, B. J., Chin, M., Nusbaum, C., Birren, B. W. & Livny, J. How deep is deep enough for RNA-Seq profiling of bacterial transcriptomes? BMC Genomics 13, 734 (2012).

52. Shishkin, A. A. et al. Simultaneous generation of many RNA-Seq libraries in a single reaction. Nat. Methods 12, 323–325 (2015).

53. Engleberg, N. C., Heath, A., Miller, A., Rivera, C. & DiRita, V. J. Spontaneous mutations in the CsrRS two-component regulatory system of Streptococcus pyogenes result in enhanced virulence in a murine model of skin and soft tissue infection. J. Infect. Dis. 183, 1043–1054 (2001).

54. Li, J. et al. Neutrophils select hypervirulent CovRS mutants of M1T1 group A Streptococcus during subcutaneous infection of mice. Infect. Immun. 82, 1579–1590 (2014).

55. Mayfield, J. A. et al. Mutations in the control of virulence sensor gene from Streptococcus pyogenes after infection in mice lead to clonal bacterial variants with altered gene regulatory activity and virulence. PLoS One 9, e100698 (2014).

56. Sumby, P., Whitney, A. R., Graviss, E. A., DeLeo, F. R. & Musser, J. M. Genome-wide analysis of group A streptococci reveals a mutation that modulates global phenotype and disease specificity. PLoS Pathog. 2, e5 (2006).

57. Tatsuno, I., Okada, R., Zhang, Y., Isaka, M. & Hasegawa, T. Partial loss of CovS function in Streptococcus pyogenes causes severe invasive disease. BMC Res. Notes 6, 126 (2013).

58. Trevino, J. et al. CovS simultaneously activates and inhibits the CovR-mediated repression of distinct subsets of group A Streptococcus virulence factor-encoding genes. Infect. Immun. 77, 3141–3149 (2009).

59. Breiman, L. Random Forests. Mach. Learn. 45, 5–32 (2001).

60. Stalhammar-Carlemalm, M., Areschoug, T., Larsson, C. & Lindahl, G. The R28 protein of Streptococcus pyogenes is related to several group B streptococcal surface proteins, confers protective immunity and promotes binding to human epithelial cells. Mol. Microbiol. 33, 208–219 (1999).

61. Stalhammar-Carlemalm, M., Stenberg, L. & Lindahl, G. Protein rib: a novel group B streptococcal cell surface protein that confers protective immunity and is expressed by most strains causing invasive infections. J. Exp. Med. 177, 1593–1603 (1993).

62. Beres, S. B. & Musser, J. M. Contribution of exogenous genetic elements to the group A Streptococcus metagenome. PLoS One 2, e800 (2007).

63. Green, N. M. et al. Genome sequence of a serotype M28 strain of group A Streptococcus: potential new insights into puerperal sepsis and bacterial disease specificity. J. Infect. Dis. 192, 760–770 (2005).

64. Coll, F. et al. Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis. Nat. Genet. 50, 307–316 (2018).

65. Earle, S. G. et al. Identifying lineage effects when controlling for population structure improves power in bacterial association studies. Nat. Microbiol. 1, 16041 (2016).

66. Gibson, G., Powell, J. E. & Marigorta, U. M. Expression quantitative trait locus analysis for translational medicine. Genome Med. 7, 60 (2015).

67. Nicolae, D. L. et al. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet. 6, e1000888 (2010).

68. Olsen, R. J. & Musser, J. M. Molecular pathogenesis of necrotizing fasciitis. Annu. Rev. Pathol. 5, 1–31 (2010).

69. Rodriguez-Ortega, M. J. et al. Characterization and identification of vaccine candidate proteins through analysis of the group A Streptococcus surface proteome. Nat. Biotechnol. 24, 191–197 (2006).

70. Zhu, L. et al. Intergenic variable-number tandem-repeat polymorphism upstream of rocA alters toxin production and enhances virulence in Streptococcus pyogenes. Infect. Immun. 84, 2086–2093 (2016).

71. Hammarlof, D. L. et al. Role of a single noncoding nucleotide in the evolution of an epidemic African clade of Salmonella.Proc. Natl Acad. Sci. USA 115, E2614–E2623 (2018).

72. Blount, Z. D., Barrick, J. E., Davidson, C. J. & Lenski, R. E. Genomic analysis of a key innovation in an experimental Escherichia coli population. Nature 489, 513–518 (2012).

73. Zaunbrecher, M. A., Sikes, R. D. Jr, Metchock, B., Shinnick, T. M. & Posey, J. E. Overexpression of the chromosomally encoded aminoglycoside acetyltransferase eis confers kanamycin resistance in Mycobacterium tuberculosis. Proc. Natl Acad. Sci. USA 106, 20004–20009 (2009).

74. Puopolo, K. M. & Madoff, L. C. Upstream short sequence repeats regulate expression of the alpha C protein of group B Streptococcus. Mol. Microbiol. 50, 977–991 (2003).

75. Stalhammar-Carlemalm, M., Areschoug, T., Larsson, C. & Lindahl, G. Cross-protection between group A and group B streptococci due to cross-reacting surface proteins. J. Infect. Dis. 182, 142–149 (2000).

76. Weckel, A. et al. The N-terminal domain of the R28 protein promotes emm28 group A Streptococcus adhesion to host cells via direct binding to three integrins. J. Biol. Chem. 293, 16006–16018 (2018).

77. Valdes, K. M. et al. The fruRBA operon is necessary for group A streptococcal growth in fructose and for resistance to neutrophil killing during growth in whole human blood. Infect. Immun. 84, 1016–1031 (2016).

78. Jeukens, J. et al. Genomics of antibiotic-resistance prediction in Pseudomonas aeruginosa.Ann. NY Acad. Sci. 1435, 5–17 (2017).

79. Nguyen, M. et al. Developing an in silico minimum inhibitory concentration panel test for Klebsiella pneumoniae. Sci. Rep. 8, 421 (2018).

80. Pesesky, M. W. et al. Evaluation of machine learning and rules-based approaches for predicting antimicrobial resistance profiles in Gram-negative bacilli from whole genome sequence data. Front. Microbiol. 7, 1887 (2016).

81. Rishishwar, L., Petit, R. A. 3rd, Kraft, C. S. & Jordan, I. K. Genome sequence-based discriminator for vancomycin-intermediate Staphylococcus aureus. J. Bacteriol. 196, 940–948 (2014).

82. Li, Y. et al. Validation of beta-lactam minimum inhibitory concentration predictions for pneumococcal isolates with newly encountered penicillin binding protein (PBP) sequences. BMC Genomics 18, 621 (2017).

83. Li, Y. et al. Penicillin-binding protein transpeptidase signatures for tracking and predicting beta-lactam resistance levels in Streptococcus pneumoniae. mBio 7, e00756-16 (2016).

84. Hao, K. et al. Lung eQTLs to help reveal the molecular underpinnings of asthma. PLoS Genet. 8, e1003029 (2012).

85. Naranbhai, V. et al. Genomic modulators of gene expression in human neutrophils. Nat. Commun. 6, 7545 (2015).

86. Ongen, H. et al. Estimating the causal tissues for complex traits and diseases. Nat. Genet. 49, 1676–1683 (2017).

87. Tung, J., Zhou, X., Alberts, S. C., Stephens, M. & Gilad, Y. The genetic architecture of gene expression levels in wild baboons. eLife https://doi.org/10.7554/eLife.04729.001 (2015).

88. Albert, F. W., Treusch, S., Shockley, A. H., Bloom, J. S. & Kruglyak, L. Genetics of single-cell protein abundance variation in large yeast populations. Nature 506, 494–497 (2014).

89. Parker, C. C. et al. Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice. Nat. Genet. 48, 919–926 (2016).

90. Francesconi, M. & Lehner, B. The effects of genetic variation on gene expression dynamics during development. Nature 505, 208–211 (2014).

91. Beres, S. B. et al. Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics. Proc. Natl Acad. Sci. USA 107, 4371–4376 (2010).

92. Olsen, R. J. et al. The majority of 9,729 group A Streptococcus strains causing disease secrete SpeB cysteine protease: pathogenesis implications. Infect. Immun. 83, 4750–4758 (2015).

93. Beres, S. B. et al. Genome sequence analysis of emm89 Streptococcus pyogenes strains causing infections in Scotland, 2010–2016. J. Med. Microbiol. 66, 1765–1773 (2017).

94. Liu, Y., Schroder, J. & Schmidt, B. Musket: a multistage k-mer spectrum-based error corrector for Illumina sequence data. Bioinformatics 29, 308–315 (2013).

95. Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9, e112963 (2014).

96. Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).

97. Inouye, M. et al. SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 6, 90 (2014).

98. Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43, e15 (2015).

99. Cheng, L., Connor, T. R., Siren, J., Aanensen, D. M. & Corander, J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. Mol. Biol. Evol. 30, 1224–1228 (2013).

100. Huson, D. H. SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14, 68–73 (1998).

101. Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 13, e1005595 (2017).

102. Long, S. W., Kachroo, P., Musser, J. M. & Olsen, R. J. Whole-genome sequencing of a human clinical isolate of emm28 Streptococcus pyogenes causing necrotizing fasciitis acquired contemporaneously with Hurricane Harvey.Genome Announc. 5, e01269-17 (2017).

103. Lees, J. A. et al. Sequence element enrichment analysis to determine the genetic basis of bacterial phenotypes. Nat. Commun. 7, 12797 (2016).

104. Bishop, C. Pattern Recognition and Machine Learning (Springer, New York, 2006).

105. Eraso, J. M. et al. Genomic landscape of intrahost variation in group A Streptococcus: repeated and abundant mutational inactivation of the fabT gene encoding a regulator of fatty acid synthesis. Infect. Immun. 84, 3268–3281 (2016).

106. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

107. Magoc, T., Wood, D. & Salzberg, S. L. EDGE-pro: estimated degree of gene expression in prokaryotic genomes. Evol. Bioinform. Online 9, 127–136 (2013).

108. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome. Biol. 15, 550 (2014).

109. Hoffman, G. E. & Schadt, E. E. variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics 17, 483 (2016).

110. Tarazona, S. et al. Data quality aware analysis of differential expression in RNA-Seq with NOISeq R/Bioc package. Nucleic Acids Res. 43, e140 (2015).

111. Ondov, B. D. et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome. Biol. 17, 132 (2016).

112. Shabalin, A. A. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353–1358 (2012).