The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity.

The spike protein of SARS-CoV-2 has been undergoing mutations and is highly glycosylated. It is critically important to investigate the biological significance of these mutations. Here, we investigated 80 variants and 26 glycosylation site modifications for the infectivity and reactivity to a panel of neutralizing antibodies and sera from convalescent patients. D614G, along with several variants containing both D614G and another amino acid change, were significantly more infectious. Most variants with amino acid change at receptor binding domain were less infectious, but variants including A475V, L452R, V483A, and F490L became resistant to some neutralizing antibodies. Moreover, the majority of glycosylation deletions were less infectious, whereas deletion of both N331 and N343 glycosylation drastically reduced infectivity, revealing the importance of glycosylation for viral infectivity. Interestingly, N234Q was markedly resistant to neutralizing antibodies, whereas N165Q became more sensitive. These findings could be of value in the development of vaccine and therapeutic antibodies.

Phase Separation of Disease-Associated SHP2 Mutants Underlies MAPK Hyperactivation.

The non-receptor protein tyrosine phosphatase (PTP) SHP2, encoded by PTPN11, plays an essential role in RAS-mitogen-activated protein kinase (MAPK) signaling during normal development. It has been perplexing as to why both enzymatically activating and inactivating mutations in PTPN11 result in human developmental disorders with overlapping clinical manifestations. Here, we uncover a common liquid-liquid phase separation (LLPS) behavior shared by these disease-associated SHP2 mutants. SHP2 LLPS is mediated by the conserved well-folded PTP domain through multivalent electrostatic interactions and regulated by an intrinsic autoinhibitory mechanism through conformational changes. SHP2 allosteric inhibitors can attenuate LLPS of SHP2 mutants, which boosts SHP2 PTP activity. Moreover, disease-associated SHP2 mutants can recruit and activate wild-type (WT) SHP2 in LLPS to promote MAPK activation. These results not only suggest that LLPS serves as a gain-of-function mechanism involved in the pathogenesis of SHP2-associated human diseases but also provide evidence that PTP may be regulated by LLPS that can be therapeutically targeted.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits.

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

RNA Interactions Are Essential for CTCF-Mediated Genome Organization.

The function of the CCCTC-binding factor (CTCF) in the organization of the genome has become an important area of investigation, but the mechanisms by which CTCF dynamically contributes to genome organization are not clear. We previously discovered that CTCF binds to large numbers of endogenous RNAs, promoting its self-association. In this regard, we now report two independent features that disrupt CTCF association with chromatin: inhibition of transcription and disruption of CTCF-RNA interactions through mutations of 2 of its 11 zinc fingers that are not required for CTCF binding to its cognate DNA site: zinc finger 1 (ZF1) or zinc finger 10 (ZF10). These mutations alter gene expression profiles as CTCF mutants lose their ability to form chromatin loops and thus the ability to insulate chromatin domains and to mediate CTCF long-range genomic interactions. Our results point to the importance of CTCF-mediated RNA interactions as a structural component of genome organization.

Hsp70- and Hsp90-Mediated Regulation of the Conformation of p53 DNA Binding Domain and p53 Cancer Variants.

The activity of the tumor suppressor p53 has to be timed and balanced closely to prevent untimely induction of cell death. The stability of p53 depends on the ubiquitin ligase Mdm2 but also on Hsp70 and Hsp90 chaperones that interact with its DNA binding domain (DBD). Using hydrogen exchange mass spectrometry and biochemical methods, we analyzed conformational states of wild-type p53-DBD at physiological temperatures and conformational perturbations in three frequent p53 cancer mutants. We demonstrate that the Hsp70/Hdj1 system shifts the conformational equilibrium of p53 toward a flexible, more mutant-like, DNA binding inactive state by binding to the DNA binding loop. The analyzed cancer mutants are likewise destabilized by interaction with the Hsp70/Hdj1 system. In contrast, Hsp90 protects the DBD of p53 wild-type and mutant proteins from unfolding. We propose that the Hsp70 and Hsp90 chaperone systems assume complementary functions to optimally balance conformational plasticity with conformational stability.

Phosphorylation by the stress-activated MAPK Slt2 down-regulates the yeast TOR complex 2.

Saccharomyces cerevisiae target of rapamycin (TOR) complex 2 (TORC2) is an essential regulator of plasma membrane lipid and protein homeostasis. How TORC2 activity is modulated in response to changes in the status of the cell envelope is unclear. Here we document that TORC2 subunit Avo2 is a direct target of Slt2, the mitogen-activated protein kinase (MAPK) of the cell wall integrity pathway. Activation of Slt2 by overexpression of a constitutively active allele of an upstream Slt2 activator (Pkc1) or by auxin-induced degradation of a negative Slt2 regulator (Sln1) caused hyperphosphorylation of Avo2 at its MAPK phosphoacceptor sites in a Slt2-dependent manner and diminished TORC2-mediated phosphorylation of its major downstream effector, protein kinase Ypk1. Deletion of Avo2 or expression of a phosphomimetic Avo2 allele rendered cells sensitive to two stresses (myriocin treatment and elevated exogenous acetic acid) that the cell requires Ypk1 activation by TORC2 to survive. Thus, Avo2 is necessary for optimal TORC2 activity, and Slt2-mediated phosphorylation of Avo2 down-regulates TORC2 signaling. Compared with wild-type Avo2, phosphomimetic Avo2 shows significant displacement from the plasma membrane, suggesting that Slt2 inhibits TORC2 by promoting Avo2 dissociation. Our findings are the first demonstration that TORC2 function is regulated by MAPK-mediated phosphorylation.

Structural basis for specific DNA sequence recognition by the transcription factor NFIL3.

The CCAAT/enhancer-binding proteins (C/EBPs) constitute a family of pivotal transcription factors involved in tissue development, cellular function, proliferation, and differentiation. NFIL3, as one of them, plays an important role in regulating immune cell differentiation, circadian clock system, and neural regeneration, yet its specific DNA recognition mechanism remains enigmatic. In this study, we showed by the ITC binding experiments that NFIL3 prefers to bind to the TTACGTAA DNA motif. Our structural studies revealed that the α-helical NFIL3 bZIP domain dimerizes through its leucine zipper region, and binds to DNA via its basic region. The two basic regions of the NFIL3 bZIP dimer were pushed apart upon binding to DNA, facilitating the snug accommodation of the two basic regions within the major grooves of the DNA. Remarkably, our binding and structural data also revealed that both NFIL3 and C/EBPα/ß demonstrate a shared preference for the TTACGTAA sequence. Furthermore, our study revealed that disease-associated mutations within the NFIL3 bZIP domain result in either reduction or complete disruption of its DNA binding ability. These discoveries not only provide valuable insights into the DNA binding mechanisms of NFIL3 but also elucidate the causal role of NFIL3 mutations in disease pathogenesis.

A large sequenced mutant library - valuable reverse genetic resource that covers 98% of sorghum genes.

Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum, a drought and heat-tolerant C4 species, requires a vast, large-scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large-scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)-induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis-regulatory elements. A comparison of the distribution of mutations in the large-scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS-induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5' UTR, gene body, and 3'-UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the Ka /Ks value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB (https://www.depts.ttu.edu/igcast/sorbmutdb.php) and SorghumBase (https://sorghumbase.org/). This current large-scale sequence-indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

Combining Traditional Mutagenesis with New High-Throughput Sequencing and Genome Editing to Reveal Hidden Variation in Polyploid Wheat.

Induced mutations have been used to generate novel variation for breeding purposes since the early 1900s. However, the combination of this old technology with the new capabilities of high-throughput sequencing has resulted in powerful reverse genetic approaches in polyploid crops. Sequencing genomes or exomes of large mutant populations can generate extensive databases of mutations for most genes. These mutant collections, together with genome editing, are being used in polyploid species to combine mutations in all copies of a gene (homoeologs), and to expose phenotypic variation that was previously hidden by functional redundancy among homoeologs. This redundancy is more extensive in recently formed polyploids such as wheat, which can now benefit from the deployment of useful recessive mutations previously identified in its diploid relatives. Sequenced mutant populations and genome editing have changed the paradigm of what is possible in functional genetic analysis of wheat.

Cryo-EM structures of cancer-specific helical and kinase domain mutations of PI3Kα.

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that perform multiple and important cellular functions. The protein investigated here belongs to class IA of the PI3Ks; it is a dimer consisting of a catalytic subunit, p110α, and a regulatory subunit, p85α, and is referred to as PI3Kα. The catalytic subunit p110α is frequently mutated in cancer. The mutations induce a gain of function and constitute a driving force in cancer development. About 80% of these mutations lead to single-amino-acid substitutions in one of three sites of p110α: two in the helical domain of the protein (E542K and E545K) and one at the C-terminus of the kinase domain (H1047R). Here, we report the cryo-electron microscopy structures of these mutants in complex with the p110α-specific inhibitor BYL-719. The H1047R mutant rotates its sidechain to a new position and weakens the kα11 activation loop interaction, thereby reducing the inhibitory effect of p85α on p110α. E542K and E545K completely abolish the tight interaction between the helical domain of p110α and the N-terminal SH2 domain of p85α and lead to the disruption of all p85α binding and a dramatic increase in flexibility of the adaptor-binding domain (ABD) in p110α. Yet, the dimerization of PI3Kα is preserved through the ABD-p85α interaction. The local and global structural features induced by these mutations provide molecular insights into the activation of PI3Kα, deepen our understanding of the oncogenic mechanism of this important signaling molecule, and may facilitate the development of mutant-specific inhibitors.

A conserved mechanism affecting hydride shifting and deprotonation in the synthesis of hopane triterpenes as compositions of wax in oat.

Triterpenoids are biologically active metabolites synthesized from a common linear precursor catalyzed by 2,3-oxidosqualene cyclases (OSCs) to form diverse triterpenoid skeletons. OSCs corresponding to many discovered triterpene alcohols in nature have not been functionally and mechanistically characterized due to the diversity of chemical structures and complexity of the cyclization mechanism. We carried out a genome-wide investigation of OSCs from Avena strigosa and discovered two triterpene synthases, namely, AsHS1 and AsHS2, using a Nicotiana benthamiana expression system. These synthases produce hopenol B and hop-17(21)-en-3ß-ol, which are components of surface wax in oat panicles and sheathes, respectively. We demonstrated that substitutions of two to three amino acid residues in AsHS1 with corresponding residues from AsHS2 allowed it to be completely converted into a hop-17(21)-en-3ß-ol synthase. AsHS2 mutants with a substitution at site 410 could synthesize hopenol B alone or mixed with a side product isomotiol. The combined quantum mechanics and molecular mechanics calculation demonstrated that the side chain size of the residue at site 410 regulated the relative orientations between the hopyl C22 cation and Phe257, leading to a difference in deprotonation positions through providing or not providing cation­π interaction between the aromatic ring of F257 and the carbocation intermediate. A similar mechanism could be applied to a hopenol B synthase from a dicotyledonous plant Aquilegia. This study provided mechanistic insight into triterpenoid synthesis and discovered key amino acid residues acting on hydride transfer and a deprotonation site to differentiate between hopane-type scaffolds in diverse plant species.

Neutralization of Rubella Vaccine Virus and Immunodeficiency-Related Vaccine-Derived Rubella Viruses by Intravenous Immunoglobulins.

The association between granulomas and vaccine-derived rubella virus (VDRV) in people with primary immune deficiencies (PID) has raised concerns about the ability of immunoglobulin (IG) preparations to neutralize VDRVs. We investigated the capacity of IG to neutralize rubella vaccine virus and four VDRV strains. As expected, the rubella vaccine virus itself was potently neutralized by IG preparations; however, the VDRV isolates from patients after intra-host evolution, 2-6 times less so. Diagnosis of immune deficiencies before possible live-virus vaccination is thus of critical importance, while IG replacement therapy can be expected to provide protection from rubella virus infection.

The occurrence of granulomas associated with vaccine derived rubella viruses (VDRV) in people with primary immune deficiencies (PID) challenges immunoglobulin (IG) preparations regarding their rubella neutralizing ability. This study confirmed potent rubella virus neutralization capacity of IG preparations and thus suggests protection of IG-treated PID patients against rubella. The study also highlights the importance of early diagnosis and timely given IG to prevent possible systemic spread of VDRV persisting locally in granulomas.

Molecular insights into the binding interactions and energetics of the omicron spike variant with hACE2 and a neutralizing antibody.

The global spread of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) since 2019 has led to a continuous evolution of viral variants, with the latest concern being the Omicron (B.1.1.529) variant. In this study, classical molecular dynamics simulations were conducted to elucidate the biophysical aspects of the Omicron spike protein's receptor-binding domain (RBD) in its interaction with human angiotensin-converting enzyme 2 (hACE2) and a neutralizing antibody, comparing it to the wildtype (WT). To model the Omicron variant, 15 in silico mutations were introduced in the RBD region of WT (retrieved from PDB). The simulations of WT spike-hACE2 and Omicron spike-hACE2 complexes revealed comparable binding stability and dynamics. Notably, the Q493R mutation in the Omicron spike increased interactions with hACE2, particularly with ASP38 and ASP355. Additionally, mutations such as N417K, T478K, and Y505H contributed to enhanced structural stability in the Omicron variant. Conversely, when comparing WT with Omicron in complex with a neutralizing antibody, simulation results demonstrated poorer binding dynamics and stability for the Omicron variant. The E484K mutation significantly decreased binding interactions, resulting in an overall decrease in binding energy (∼-57 kcal/mol) compared to WT (∼-84 kcal/mol). This study provides valuable molecular insights into the heightened infectivity of the Omicron variant, shedding light on the specific mutations influencing its interactions with hACE2 and neutralizing antibodies.

The function of the phytoplasma effector SWP12 depends on the properties of two key amino acids.

Phytoplasmas are insect-borne bacterial pathogens capable of secreting effectors into host cells and interfering with host plant defense response processes. Previous studies have found that the Candidatus Phytoplasma tritici effector SWP12 binds to and destabilizes the wheat transcription factor TaWRKY74, increasing wheat susceptibility to phytoplasmas. Here, we used a Nicotiana benthamiana transient expression system to identify two key functional sites of SWP12 and screened a series of truncated mutants and amino acid substitution mutants to determine whether they inhibit Bax-induced cell death. Using a subcellular localization assay and online structure analysis websites, we found that structure rather than intracellular localization probably affects the function of SWP12. D33A and P85H are two inactive substitution mutants, neither of which interacts with TaWRKY74, and P85H does not inhibit Bax-induced cell death, suppress flg22-triggered reactive oxygen species (ROS) bursts, degrade TaWRKY74, or promote phytoplasma accumulation. D33A can weakly suppress Bax-induced cell death and flg22-triggered ROS bursts and degrade a portion of TaWRKY74 and weakly promote phytoplasma accumulation. S53L, CPP, and EPWB are three SWP12 homolog proteins from other phytoplasmas. Sequence analysis revealed that D33 was conserved in these proteins, and they exhibited the same polarity at P85. Transient expression in N. benthamiana showed that these proteins could inhibit Bax-induced cell death and suppress ROS bursts. Our findings clarified that P85 and D33 of SWP12 play critical and minor roles, respectively, in suppressing the plant defense response and that they play a preliminary role in determining the functions of homologous proteins.

The T1150A cancer mutant of the protein lysine dimethyltransferase NSD2 can introduce H3K36 trimethylation.

Protein lysine methyltransferases (PKMTs) play essential roles in gene expression regulation and cancer development. Somatic mutations in PKMTs are frequently observed in cancer cells. In biochemical experiments, we show here that the NSD1 mutations Y1971C, R2017Q, and R2017L observed mostly in solid cancers are catalytically inactive suggesting that NSD1 acts as a tumor suppressor gene in these tumors. In contrast, the frequently observed T1150A in NSD2 and its T2029A counterpart in NSD1, both observed in leukemia, are hyperactive and introduce up to three methyl groups in H3K36 in biochemical and cellular assays, while wildtype NSD2 and NSD1 only introduce up to two methyl groups. In Molecular Dynamics simulations, we determined key mechanistic and structural features controlling the product specificity of this class of enzymes. Simulations with NSD2 revealed that H3K36me3 formation is possible due to an enlarged active site pocket of T1150A and loss of direct contacts of T1150 to critical residues which regulate the product specificity of NSD2. Bioinformatic analyses of published data suggested that the generation of H3K36me3 by NSD2 T1150A could alter gene regulation by antagonizing H3K27me3 finally leading to the upregulation of oncogenes.

ZmXYL modulates auxin-induced maize growth.

Plant architecture, lodging resistance, and yield are closely associated with height. In this paper, we report the identification and characterization of two allelic EMS-induced mutants of Zea mays, xyl-1, and xyl-2 that display dwarf phenotypes. The mutated gene, ZmXYL, encodes an α-xylosidase which functions in releasing xylosyl residue from a ß-1,4-linked glucan chain. Total α-xylosidase activity in the two alleles is significantly decreased compared to wild-type plants. Loss-of-function mutants of ZmXYL resulted in a decreased xylose content, an increased XXXG content in xyloglucan (XyG), and a reduced auxin content. We show that auxin has an antagonistic effect with XXXG in promoting cell divisions within mesocotyl tissue. xyl-1 and xyl-2 were less sensitive to IAA compared to B73. Based on our study, a model is proposed that places XXXG, an oligosaccharide derived from XyG and the substrate of ZmXYL, as having a negative impact on auxin homeostasis resulting in the dwarf phenotypes of the xyl mutants. Our results provide a insight into the roles of oligosaccharides released from plant cell walls as signals in mediating plant growth and development.

Strigolactones repress nodule development and senescence in pea.

Strigolactones are a class of phytohormones that are involved in many different plant developmental processes, including the rhizobium-legume nodule symbiosis. Although both positive and negative effects of strigolactones on the number of nodules have been reported, the influence of strigolactones on nodule development is still unknown. Here, by means of the ramosus (rms) mutants of Pisum sativum (pea) cv Terese, we investigated the impact of strigolactone biosynthesis (rms1 and rms5) and signaling (rms3 and rms4) mutants on nodule growth. The rms mutants had more red, that is, functional, and larger nodules than the wild-type plants. Additionally, the increased nitrogen fixation and senescence zones with consequently reduced meristematic and infection zones indicated that the rms nodules developed faster than the wild-type nodules. An enhanced expression of the nodule zone-specific molecular markers for meristem activity and senescence supported the enlarged, fast maturing nodules. Interestingly, the master nodulation regulator, NODULE INCEPTION, NIN, was strongly induced in nodules of all rms mutants but not prior to inoculation. Determination of sugar levels with both bulk and spatial metabolomics in roots and nodules, respectively, hints at slightly increased malic acid levels early during nodule primordia formation and reduced sugar levels at later stages, possibly the consequence of an increased carbon usage of the enlarged nodules, contributing to the enhanced senescence. Taken together, these results suggest that strigolactones regulate the development of nodules, which is probably mediated through NIN, and available plant sugars.

The ALOG family members OsG1L1 and OsG1L2 regulate inflorescence branching in rice.

The architecture of the rice inflorescence is an important determinant of crop yield. The length of the inflorescence and the number of branches are among the key factors determining the number of spikelets, and thus grains, that a plant will develop. In particular, the timing of the identity transition from indeterminate branch meristem to determinate spikelet meristem governs the complexity of the inflorescence. In this context, the ALOG gene TAWAWA1 (TAW1) has been shown to delay the transition to determinate spikelet development in Oryza sativa (rice). Recently, by combining precise laser microdissection of inflorescence meristems with RNA-seq, we observed that two ALOG genes, OsG1-like 1 (OsG1L1) and OsG1L2, have expression profiles similar to that of TAW1. Here, we report that osg1l1 and osg1l2 loss-of-function CRISPR mutants have similar phenotypes to the phenotype of the previously published taw1 mutant, suggesting that these genes might act on related pathways during inflorescence development. Transcriptome analysis of the osg1l2 mutant suggested interactions of OsG1L2 with other known inflorescence architecture regulators and the data sets were used for the construction of a gene regulatory network (GRN), proposing interactions among genes potentially involved in controlling inflorescence development in rice. In this GRN, we selected the homeodomain-leucine zipper transcription factor encoding the gene OsHOX14 for further characterization. The spatiotemporal expression profiling and phenotypical analysis of CRISPR loss-of-function mutants of OsHOX14 suggests that the proposed GRN indeed serves as a valuable resource for the identification of new proteins involved in rice inflorescence development.

Stability of N-type inactivation and the coupling between N-type and C-type inactivation in the Aplysia Kv1 channel.

The voltage-dependent potassium channels (Kv channels) show several different types of inactivation. N-type inactivation is a fast inactivating mechanism, which is essentially an open pore blockade by the amino-terminal structure of the channel itself or the auxiliary subunit. There are several functionally discriminatable slow inactivation (C-type, P-type, U-type), the mechanism of which is supposed to include rearrangement of the pore region. In some Kv1 channels, the actual inactivation is brought about by coupling of N-type and C-type inactivation (N-C coupling). In the present study, we focused on the N-C coupling of the Aplysia Kv1 channel (AKv1). AKv1 shows a robust N-type inactivation, but its recovery is almost thoroughly from C-type inactivated state owing to the efficient N-C coupling. In the I8Q mutant of AKv1, we found that the inactivation as well as its recovery showed two kinetic components apparently correspond to N-type and C-type inactivation. Also, the cumulative inactivation which depends on N-type mechanism in AKv1 was hindered in I8Q, suggesting that N-type inactivation of I8Q is less stable. We also found that Zn 2 + specifically accelerates C-type inactivation of AKv1 and that H382 in the pore turret is involved in the Zn 2 + binding. Because the region around Ile 8 (I8) in AKv1 has been suggested to be involved in the pre-block binding of the amino-terminal structure, our results strengthen a hypothesis that the stability of the pre-block state is important for stable N-type inactivation as well as the N-C coupling in the Kv1 channel inactivation.

Genes required for survival and proliferation of Mycoplasma bovis in association with host cells.

Mycoplasma bovis is an important emerging pathogen of cattle and bison, but our understanding of the genetic basis of its interactions with its host is limited. The aim of this study was to identify genes of M. bovis required for interaction and survival in association with host cells. One hundred transposon-induced mutants of the type strain PG45 were assessed for their capacity to survive and proliferate in Madin-Darby bovine kidney cell cultures. The growth of 19 mutants was completely abrogated, and 47 mutants had a prolonged doubling time compared to the parent strain. All these mutants had a similar growth pattern to the parent strain PG45 in the axenic media. Thirteen genes previously classified as dispensable for the axenic growth of M. bovis were found to be essential for the growth of M. bovis in association with host cells. In most of the mutants with a growth-deficient phenotype, the transposon was inserted into a gene involved in transportation or metabolism. This included genes coding for ABC transporters, proteins related to carbohydrate, nucleotide and protein metabolism, and membrane proteins essential for attachment. It is likely that these genes are essential not only in vitro but also for the survival of M. bovis in infected animals. IMPORTANCE: Mycoplasma bovis causes chronic bronchopneumonia, mastitis, arthritis, keratoconjunctivitis, and reproductive tract disease in cattle around the globe and is an emerging pathogen in bison. Control of mycoplasma infections is difficult in the absence of appropriate antimicrobial treatment or effective vaccines. A comprehensive understanding of host-pathogen interactions and virulence factors is important to implement more effective control methods against M. bovis. Recent studies of other mycoplasmas with in vitro cell culture models have identified essential virulence genes of mycoplasmas. Our study has identified genes of M. bovis required for survival in association with host cells, which will pave the way to a better understanding of host-pathogen interactions and the role of specific genes in the pathogenesis of disease caused by M. bovis.