Riboswitches are RNAs that form complex, folded structures that selectively bind small molecules or ions. As with certain groups of protein enzymes and receptors, some riboswitch classes have evolved to change their ligand specificity. We developed a procedure to systematically analyze known riboswitch classes to find additional variants that have altered their ligand specificity. This approach uses multiple-sequence alignments, atomic-resolution structural information, and riboswitch gene associations. Among the discoveries are unique variants of the guanine riboswitch class that most tightly bind the nucleoside 2′-deoxyguanosine. In addition, we identified variants of the glycine riboswitch class that no longer recognize this amino acid, additional members of a rare flavin mononucleotide (FMN) variant class, and also variants of c-di-GMP-I and -II riboswitches that might recognize different bacterial signaling molecules. These findings further reveal the diverse molecular sensing capabilities of RNA, which highlights the potential for discovering a large number of additional natural riboswitch classes.iboswitches are structured noncoding RNA domains that regulate gene expression in response to the selective binding of small-molecule or ion ligands. The discovery of numerous classes of riboswitches has helped reveal how RNAs can form exquisitely precise ligand-binding pockets using only the four common RNA nucleotides (1-4). Furthermore, each discovery links the riboswitch ligand to the protein products of the genes under regulation. Recent riboswitch findings have exposed unique facets of biology, such as the widespread molecular mechanisms that confer fluoride (5) or guanidine (6) resistance, that maintain metal ion homeostasis (7-9), and that control important bacterial processes such as sporulation, biofilm formation, and chemotaxis (10-14). Thus, the identification of additional riboswitch classes promises to offer insights into otherwise hidden biological processes and their regulation.Riboswitch variants have been reported previously, wherein the ligand-binding "aptamer" domain has mutated to accommodate a different metabolite or signaling compound. The identification of such RNAs provides rare opportunities to study how small changes in RNA sequence can lead to major changes in smallmolecule ligand affinity. There have been seven examples, either experimentally validated or suspected, of ligand specificity changes reported to date. These include guanine aptamer variations present in riboswitches for adenine (15) and 2′-deoxyguanosine (2′-dG) (16), c-di-GMP-I aptamer variations that result in riboswitches that bind the recently discovered bacterial signaling molecule c-AMP-GMP (13, 14), and coenzyme B 12 aptamer changes (17, 18) that yield riboswitches selective for aquocobalamin (19). Three additional ligand specificity changes are suspected. Namely, some molybdenum cofactor riboswitches appear to exploit an altered aptamer structure to selectively recognize tungsten cofactor (20), certain flavin mononu...