Engaging large numbers of undergraduates in authentic scientific discovery is desirable but difficult to achieve. We have developed a general model in which faculty and teaching assistants from diverse academic institutions are trained to teach a research course for first-year undergraduate students focused on bacteriophage discovery and genomics. The course is situated within a broader scientific context aimed at understanding viral diversity, such that faculty and students are collaborators with established researchers in the field. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) course has been widely implemented and has been taken by over 4,800 students at 73 institutions. We show here that this alliance-sourced model not only substantially advances the field of phage genomics but also stimulates students’ interest in science, positively influences academic achievement, and enhances persistence in science, technology, engineering, and mathematics (STEM) disciplines. Broad application of this model by integrating other research areas with large numbers of early-career undergraduate students has the potential to be transformative in science education and research training.

Temperate phages are common and prophages are abundant residents of sequenced bacterial genomes. Mycobacteriophages are viruses infecting mycobacterial hosts including Mycobacterium tuberculosis and Mycobacterium smegmatis, encompass substantial genetic diversity, and are commonly temperate. Characterization of ten Cluster N temperate mycobacteriophages reveals at least five distinct prophage-expressed viral defense systems that interfere with infection of lytic and temperate phages that are either closely-related (homotypic defense) or unrelated (heterotypic defense). Target specificity is unpredictable, ranging from a single target phage to one-third of those tested. The defense systems include a single-subunit restriction system, a heterotypic exclusion system, and a predicted (p)ppGpp synthetase, which blocks lytic phage growth, promotes bacterial survival, and enables efficient lysogeny. The predicted (p)ppGpp synthetase coded by the Phrann prophage defends against phage Tweety infection, but Tweety codes for a tetrapeptide repeat protein, gp54, that acts as a highly effective counter-defense system. Prophage-mediated viral defense offers an efficient mechanism for bacterial success in host-virus dynamics, and counter-defense promotes phage co-evolution.

While course-based research in genomics can generate both knowledge gains and a greater appreciation for how science is done, a significant investment of course time is required to enable students to show gains commensurate to a summer research experience. Nonetheless, this is a very cost-effective way to reach larger numbers of students.

There have been numerous calls to engage students in science as science is done. A survey of 90-plus faculty members explores barriers and incentives when developing a research-based genomics course. The results indicate that a central core supporting a national experiment can help overcome local obstacles.

The Muller F element (4.2 Mb, ~80 protein-coding genes) is an unusual autosome of Drosophila melanogaster; it is mostly heterochromatic with a low recombination rate. To investigate how these properties impact the evolution of repeats and genes, we manually improved the sequence and annotated the genes on the D. erecta, D. mojavensis, and D. grimshawi F elements and euchromatic domains from the Muller D element. We find that F elements have greater transposon density (25–50%) than euchromatic reference regions (3–11%). Among the F elements, D. grimshawi has the lowest transposon density (particularly DINE-1: 2% vs. 11–27%). F element genes have larger coding spans, more coding exons, larger introns, and lower codon bias. Comparison of the Effective Number of Codons with the Codon Adaptation Index shows that, in contrast to the other species, codon bias in D. grimshawi F element genes can be attributed primarily to selection instead of mutational biases, suggesting that density and types of transposons affect the degree of local heterochromatin formation. F element genes have lower estimated DNA melting temperatures than D element genes, potentially facilitating transcription through heterochromatin. Most F element genes (~90%) have remained on that element, but the F element has smaller syntenic blocks than genome averages (3.4–3.6 vs. 8.4–8.8 genes per block), indicating greater rates of inversion despite lower rates of recombination. Overall, the F element has maintained characteristics that are distinct from other autosomes in the Drosophila lineage, illuminating the constraints imposed by a heterochromatic milieu.