Relaxed selection is a precursor to the evolution of phenotypic plasticity

Phenotypic plasticity allows organisms to produce alternative phenotypes under different conditions and represents one of the most important ways by which organisms adaptively respond to the environment. However, the relationship between phenotypic plasticity and molecular evolution remains poorly understood. We addressed this issue by investigating the evolution of genes associated with phenotypically plastic castes, sexes, and developmental stages of the fire ant Solenopsis invicta. We first determined if genes associated with phenotypic plasticity in S. invicta evolved at a rapid rate, as predicted under theoretical models. We found that genes differentially expressed between S. invicta castes, sexes, and developmental stages all exhibited elevated rates of evolution compared with ubiquitously expressed genes. We next investigated the evolutionary history of genes associated with the production of castes. Surprisingly, we found that orthologs of caste-biased genes in S. invicta and the social bee Apis mellifera evolved rapidly in lineages without castes. Thus, in contrast to some theoretical predictions, our results suggest that rapid rates of molecular evolution may not arise primarily as a consequence of phenotypic plasticity. Instead, genes evolving under relaxed purifying selection may more readily adopt new forms of biased expression during the evolution of alternate phenotypes. These results suggest that relaxed selective constraint on protein-coding genes is an important and underappreciated element in the evolutionary origin of phenotypic plasticity.     

Intercontinental differences in resource use reveal the importance of mutualisms in fire ant invasions

Mutualisms play key roles in the functioning of ecosystems. However, reciprocally beneficial interactions that involve introduced species also can enhance invasion success and in doing so compromise ecosystem integrity. For example, the growth and competitive ability of introduced plant species can increase when fungal or microbial associates provide limiting nutrients. Mutualisms also may aid animal invasions, but how such systems may promote invasion success has received relatively little attention. Here we examine how access to food-for-protection mutualisms involving the red imported fire ant (Solenopsis invicta) aids the success of this prominent invader. Intense interspecific competition in its native Argentina constrained the ability of S. invicta to benefit from honeydew-producing Hemiptera (and other accessible sources of carbohydrates), whereas S. invicta dominated these resources in its introduced range in the United States. Consistent with this strong pattern, nitrogen isotopic data revealed that fire ants from populations in the United States occupy a lower trophic position than fire ants from Argentina. Laboratory and field experiments demonstrated that honeydew elevated colony growth, a crucial determinant of competitive performance, even when insect prey were not limiting. Carbohydrates, obtained largely through mutualistic partnerships with other organisms, thus represent critical resources that may aid the success of this widespread invasive species. These results illustrate the potential for mutualistic interactions to play a fundamental role in the establishment and spread of animal invasions.     

Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile)

Ants are some of the most abundant and familiar animals on Earth, and they play vital roles in most terrestrial ecosystems. Although all ants are eusocial, and display a variety of complex and fascinating behaviors, few genomic resources exist for them. Here, we report the draft genome sequence of a particularly widespread and well-studied species, the invasive Argentine ant (Linepithema humile), which was accomplished using a combination of 454 (Roche) and Illumina sequencing and community-based funding rather than federal grant support. Manual annotation of >1,000 genes from a variety of different gene families and functional classes reveals unique features of the Argentine ant's biology, as well as similarities to Apis mellifera and Nasonia vitripennis. Distinctive features of the Argentine ant genome include remarkable expansions of gustatory (116 genes) and odorant receptors (367 genes), an abundance of cytochrome P450 genes (>110), lineage-specific expansions of yellow/major royal jelly proteins and desaturases, and complete CpG DNA methylation and RNAi toolkits. The Argentine ant genome contains fewer immune genes than Drosophila and Tribolium, which may reflect the prominent role played by behavioral and chemical suppression of pathogens. Analysis of the ratio of observed to expected CpG nucleotides for genes in the reproductive development and apoptosis pathways suggests higher levels of methylation than in the genome overall. The resources provided by this genome sequence will offer an abundance of tools for researchers seeking to illuminate the fascinating biology of this emerging model organism.     

Unconscious Woman in Shock and Covered with Ants Pulled from an Abandoned Automobile

BackgroundA middle-aged woman was taken from an abandoned automobile unconscious and covered with ants in Tucson, Arizona. When hospitalized in July 2018, she had an extensive papular-pustular skin eruption on her abdomen and thigh and disseminated intravascular coagulation. She was stung innumerable times by native golden fire ants (Solenopsis aurea) while sleeping in the vehicle. The large amount of venom injected by stings into this individual may have triggered dissemnated intravascualar coagulation because the venom contains powerful hemolytic factors.MethodsThe patient history is presented and ants were captured and identified.ResultsClinical findings of fire ant stings are presented and the importance of recognizing the distinctive skin lesions that occur is emphasized. Stings of the red imported fire ant, Solenopsis invicta, and the black imported fire ant, Solenopsis richteri, cause skin lesions recognized by physicians and victims alike in the southern and southeastern United States. Native fire ant stings are documented much less often. However, there is significant cross-reactivity among the venoms of Solenopsis species.ConclusionIt is important for clinicians to recognize the characteristic skin lesions of fire ant envenomation as fire ant populations are expanding and they sting millions of people each year.     

Laccase detoxification mediates the nutritional alliance between leaf-cutting ants and fungus-garden symbionts

Leaf-cutting ants combine large-scale herbivory with fungus farming to sustain advanced societies. Their stratified colonies are major evolutionary achievements and serious agricultural pests, but the crucial adaptations that allowed this mutualism to become the prime herbivorous component of neotropical ecosystems has remained elusive. Here we show how coevolutionary adaptation of a specific enzyme in the fungal symbiont has helped leaf-cutting ants overcome plant defensive phenolic compounds. We identify nine putative laccase-coding genes in the fungal genome of Leucocoprinus gongylophorus cultivated by the leaf-cutting ant Acromyrmex echinatior. One of these laccases (LgLcc1) is highly expressed in the specialized hyphal tips (gongylidia) that the ants preferentially eat, and we confirm that these ingested laccase molecules pass through the ant guts and remain active when defecated on the leaf pulp that the ants add to their gardens. This accurate deposition ensures that laccase activity is highest where new leaf material enters the fungus garden, but where fungal mycelium is too sparse to produce extracellular enzymes in sufficient quantities to detoxify phenolic compounds. Phylogenetic analysis of LgLcc1 ortholog sequences from symbiotic and free-living fungi revealed significant positive selection in the ancestral lineage that gave rise to the gongylidia-producing symbionts of leaf-cutting ants and their non–leaf-cutting ant sister group. Our results are consistent with fungal preadaptation and subsequent modification of a particular laccase enzyme for the detoxification of secondary plant compounds during the transition to active herbivory in the ancestor of leaf-cutting ants between 8 and 12 Mya.     

Distant Mimivirus relative with a larger genome highlights the fundamental features of Megaviridae

Mimivirus, a DNA virus infecting acanthamoeba, was for a long time the largest known virus both in terms of particle size and gene content. Its genome encodes 979 proteins, including the first four aminoacyl tRNA synthetases (ArgRS, CysRS, MetRS, and TyrRS) ever found outside of cellular organisms. The discovery that Mimivirus encoded trademark cellular functions prompted a wealth of theoretical studies revisiting the concept of virus and associated large DNA viruses with the emergence of early eukaryotes. However, the evolutionary significance of these unique features remained impossible to assess in absence of a Mimivirus relative exhibiting a suitable evolutionary divergence. Here, we present Megavirus chilensis, a giant virus isolated off the coast of Chile, but capable of replicating in fresh water acanthamoeba. Its 1,259,197-bp genome is the largest viral genome fully sequenced so far. It encodes 1,120 putative proteins, of which 258 (23%) have no Mimivirus homologs. The 594 Megavirus/Mimivirus orthologs share an average of 50% of identical residues. Despite this divergence, Megavirus retained all of the genomic features characteristic of Mimivirus, including its cellular-like genes. Moreover, Megavirus exhibits three additional aminoacyl-tRNA synthetase genes (IleRS, TrpRS, and AsnRS) adding strong support to the previous suggestion that the Mimivirus/Megavirus lineage evolved from an ancestral cellular genome by reductive evolution. The main differences in gene content between Mimivirus and Megavirus genomes are due to (i) lineages specific gains or losses of genes, (ii) lineage specific gene family expansion or deletion, and (iii) the insertion/migration of mobile elements (intron, intein).     

Identification of a pathogenicity island required for Salmonella survival in host cells.

We have identified a region unique to the Salmonella typhimurium chromosome that is essential for virulence in mice. This region harbors at least three genes: two (spiA and spiB) encode products that are similar to proteins found in type III secretion systems, and a third (spiR) encodes a putative regulator. A strain with a mutation in spiA was unable to survive within macrophages but displayed wild-type levels of epithelial cell invasion. The culture supernatants of the spi mutants lacked a modified form of flagellin, which was present in the supernatant of the wild-type strain. This suggests that the Spi secretory apparatus exports a protease, or a protein that can alter the activity of a secreted protease. The "pathogenicity island" harboring the spi genes may encode the virulence determinants that set Salmonella apart from other enteric pathogens.     

Searching for the autoimmune thyroid disease susceptibility genes: from gene mapping to gene function

The autoimmune thyroid diseases (AITD) are complex diseases that are caused by an interaction between susceptibility genes and environmental triggers. Genetic susceptibility, in combination with external factors (e.g., dietary iodine), is believed to initiate the autoimmune response to thyroid antigens. Abundant epidemiological data, including family and twin studies, point to a strong genetic influence on the development of AITD. Various techniques have been used to identify the genes contributing to the etiology of AITD, including candidate gene analysis and whole genome screening. These studies have enabled the identification of several loci (genetic regions) that are linked with AITD, and in some of these loci putative AITD susceptibility genes have been identified. Some of these genes/loci are unique to Graves' disease (GD) and Hashimoto's thyroiditis (HT), and some are common to both diseases, indicating that there is a shared genetic susceptibility to GD and HT. The putative GD and HT susceptibility genes include both immune modifying genes (e.g., human leukocyte antigen, cytotoxic T lymphocyte antigen-4) and thyroid-specific genes (e.g., TSH receptor, thyroglobulin). Most likely these loci interact, and their interactions may influence disease phenotype and severity. It is hoped that in the near future additional AITD susceptibility genes will be identified and the mechanisms by which they induce AITD will be unraveled.     

Chromosome segregation in Archaea mediated by a hybrid DNA partition machine

Eukarya and, more recently, some bacteria have been shown to rely on a cytoskeleton-based apparatus to drive chromosome segregation. In contrast, the factors and mechanisms underpinning this fundamental process are underexplored in archaea, the third domain of life. Here we establish that the archaeon Sulfolobus solfataricus harbors a hybrid segrosome consisting of two interacting proteins, SegA and SegB, that play a key role in genome segregation in this organism. SegA is an ortholog of bacterial, Walker-type ParA proteins, whereas SegB is an archaea-specific factor lacking sequence identity to either eukaryotic or bacterial proteins, but sharing homology with a cluster of uncharacterized factors conserved in both crenarchaea and euryarchaea, the two major archaeal sub-phyla. We show that SegA is an ATPase that polymerizes in vitro and that SegB is a site-specific DNA-binding protein contacting palindromic sequences located upstream of the segAB cassette. SegB interacts with SegA in the presence of nucleotides and dramatically affects its polymerization dynamics. Our data demonstrate that SegB strongly stimulates SegA polymerization, possibly by promoting SegA nucleation and accelerating polymer growth. Increased expression levels of segAB resulted in severe growth and chromosome segregation defects, including formation of anucleate cells, compact nucleoids confined to one half of the cell compartment and fragmented nucleoids. The overall picture emerging from our findings indicates that the SegAB complex fulfills a crucial function in chromosome segregation and is the prototype of a DNA partition machine widespread across archaea.     

Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle

As an obligatory parasite of humans, the body louse (Pediculus humanus humanus) is an important vector for human diseases, including epidemic typhus, relapsing fever, and trench fever. Here, we present genome sequences of the body louse and its primary bacterial endosymbiont Candidatus Riesia pediculicola. The body louse has the smallest known insect genome, spanning 108 Mb. Despite its status as an obligate parasite, it retains a remarkably complete basal insect repertoire of 10,773 protein-coding genes and 57 microRNAs. Representing hemimetabolous insects, the genome of the body louse thus provides a reference for studies of holometabolous insects. Compared with other insect genomes, the body louse genome contains significantly fewer genes associated with environmental sensing and response, including odorant and gustatory receptors and detoxifying enzymes. The unique architecture of the 18 minicircular mitochondrial chromosomes of the body louse may be linked to the loss of the gene encoding the mitochondrial single-stranded DNA binding protein. The genome of the obligatory louse endosymbiont Candidatus Riesia pediculicola encodes less than 600 genes on a short, linear chromosome and a circular plasmid. The plasmid harbors a unique arrangement of genes required for the synthesis of pantothenate, an essential vitamin deficient in the louse diet. The human body louse, its primary endosymbiont, and the bacterial pathogens that it vectors all possess genomes reduced in size compared with their free-living close relatives. Thus, the body louse genome project offers unique information and tools to use in advancing understanding of coevolution among vectors, symbionts, and pathogens.     

Identification of neuropeptide-like protein gene families in Caenorhabditiselegans and other species.

Neuropeptides play critical roles in synaptic signaling in all nervous systems. Unlike classical neurotransmitters, peptidergic neurotransmitters are encoded as preproproteins that are posttranslationally processed to yield bioactive neuropeptides. To identify novel peptidergic neurotransmitters, the Caenorhabditis elegans genome was searched for predicted proteins with the structural hallmarks of neuropeptide preproproteins. Thirty-two C. elegans neuropeptide-like protein (nlp) genes were identified. The nlp genes define at least 11 families of putative neuropeptides with unique motifs; similar expressed sequence tags were identified in other invertebrate species for all 11 families. Six of these families are defined by putative bioactive motifs (FAFA, GGxYamide, MRxamide, LQFamide, LxDxamide, and GGARAF); the remaining five families are related to allatostatin, myomodulin, buccalin/drosulfakinin, orcokinin, and APGWamide neuropeptides (MGL/Famide, FRPamide, MSFamide, GFxGF, and YGGWamide families, respectively). Most C. elegans nlp gene expression is in neurons. The C. elegans nlp genes and similar genes encoding putative neuropeptides in other species are likely to play diverse roles in nervous system function.     

Identification and chromosomal distribution of 5S rRNA genes in Neurospora crassa.

The 5S rRNA genes of Neurospora crassa, unlike those of most organisms, are not tandemly arranged, and they are found imbedded in a variety of unique sequences. The 5S rRNA regions of most of the genes are of one type, alpha; however, several other "isotypes" (beta, gamma, delta, zeta, and eta) are also found. We asked whether Neurospora 5S rRNA genes are dispersed on a chromosomal scale and whether genes of different isotypes are spatially segregated. We identified, by DNA sequencing, 5S rRNA genes in 22 5S DNA clones, and we mapped these genes by conventional crosses by using restriction fragment length polymorphisms in their flanking sequences as genetic markers. The results show that the 5S rRNA genes are distributed on at least six of the seven chromosomes. Their location does not appear to be completely random. Some of them are closely linked. One of the chromosomes carries a disproportionate number of 5S rRNA genes of the most common structural type, alpha; another chromosome carries three of the four mapped beta 5S rRNA genes. None of the 5S rRNA genes studied maps close to the nucleolus organizer, the site of the genes that code for the three larger rRNAs.     

New target region of allelic loss in hepatocellular carcinomas within a 1-cM interval on chromosome 6q23

BACKGROUND/AIMS: Frequent allelic losses on the long arm of chromosome 16 in several types of human cancers have suggested that 16q harbors one or more genes that are important for suppressing tumorigenesis in the tissues in question. METHODS: To identify the locations of putative tumor suppressor genes involved in hepatocellular carcinoma, we examined 96 primary hepatocellular carcinomas for their patterns of allelic loss at 18 microsatellite marker loci distributed along this chromosome arm. RESULTS: Allelic loss at one or more loci was observed in 48 (50%) of these tumors. The highest frequency of loss of heterozygosity (42%) was observed with marker D6S311 on chromosome 6q23. Through detailed deletion mapping of tumors having partial or interstitial deletions, we identified two commonly deleted regions at 6q23 and at 6q26-27. CONCLUSIONS: The common region at 6q23 lay within a 1-cM interval, flanked by D6S977 and D6S311. The previously documented deletion region that includes the M6P/IGF2R locus was confined to a 20-cM region at band 6q26-27 in our panel of tumors. The location we defined at 6q23 for a putative suppressor of hepatocellular carcinoma has not been reported before.     

The genome of Syntrophus aciditrophicus: life at the thermodynamic limit of microbial growth

Biochemically, the syntrophic bacteria constitute the missing link in our understanding of anaerobic flow of carbon in the biosphere. The completed genome sequence of Syntrophus aciditrophicus SB, a model fatty acid- and aromatic acid-degrading syntrophic bacterium, provides a glimpse of the composition and architecture of the electron transfer and energy-transducing systems needed to exist on marginal energy economies of a syntrophic lifestyle. The genome contains 3,179,300 base pairs and 3,169 genes where 1,618 genes were assigned putative functions. Metabolic reconstruction of the gene inventory revealed that most biosynthetic pathways of a typical Gram-negative microbe were present. A distinctive feature of syntrophic metabolism is the need for reverse electron transport; the presence of a unique Rnf-type ion-translocating electron transfer complex, menaquinone, and membrane-bound Fe-S proteins with associated heterodisulfide reductase domains suggests mechanisms to accomplish this task. Previously undescribed approaches to degrade fatty and aromatic acids, including multiple AMP-forming CoA ligases and acyl-CoA synthetases seem to be present as ways to form and dissipate ion gradients by using a sodium-based energy strategy. Thus, S. aciditrophicus, although nutritionally self-sufficient, seems to be a syntrophic specialist with limited fermentative and respiratory metabolism. Genomic analysis confirms the S. aciditrophicus metabolic and regulatory commitment to a nonconventional mode of life compared with our prevailing understanding of microbiology.     

Giant virus with a remarkable complement of genes infects marine zooplankton

As major consumers of heterotrophic bacteria and phytoplankton, microzooplankton are a critical link in aquatic foodwebs. Here, we show that a major marine microflagellate grazer is infected by a giant virus, Cafeteria roenbergensis virus (CroV), which has the largest genome of any described marine virus (≈730 kb of double-stranded DNA). The central 618-kb coding part of this AT-rich genome contains 544 predicted protein-coding genes; putative early and late promoter motifs have been detected and assigned to 191 and 72 of them, respectively, and at least 274 genes were expressed during infection. The diverse coding potential of CroV includes predicted translation factors, DNA repair enzymes such as DNA mismatch repair protein MutS and two photolyases, multiple ubiquitin pathway components, four intein elements, and 22 tRNAs. Many genes including isoleucyl-tRNA synthetase, eIF-2γ, and an Elp3-like histone acetyltransferase are usually not found in viruses. We also discovered a 38-kb genomic region of putative bacterial origin, which encodes several predicted carbohydrate metabolizing enzymes, including an entire pathway for the biosynthesis of 3-deoxy-d-manno-octulosonate, a key component of the outer membrane in Gram-negative bacteria. Phylogenetic analysis indicates that CroV is a nucleocytoplasmic large DNA virus, with Acanthamoeba polyphaga mimivirus as its closest relative, although less than one-third of the genes of CroV have homologs in Mimivirus. CroV is a highly complex marine virus and the only virus studied in genetic detail that infects one of the major groups of predators in the oceans.     

The genome of Clostridium kluyveri, a strict anaerobe with unique metabolic features

Clostridium kluyveri is unique among the clostridia; it grows anaerobically on ethanol and acetate as sole energy sources. Fermentation products are butyrate, caproate, and H₂. We report here the genome sequence of C. kluyveri, which revealed new insights into the metabolic capabilities of this well studied organism. A membrane-bound energy-converting NADH:ferredoxin oxidoreductase (RnfCDGEAB) and a cytoplasmic butyryl-CoA dehydrogenase complex (Bcd/EtfAB) coupling the reduction of crotonyl-CoA to butyryl-CoA with the reduction of ferredoxin represent a new energy-conserving module in anaerobes. The genes for NAD-dependent ethanol dehydrogenase and NAD(P)-dependent acetaldehyde dehydrogenase are located next to genes for microcompartment proteins, suggesting that the two enzymes, which are isolated together in a macromolecular complex, form a carboxysome-like structure. Unique for a strict anaerobe, C. kluyveri harbors three sets of genes predicted to encode for polyketide/nonribosomal peptide synthetase hybrides and one set for a nonribosomal peptide synthetase. The latter is predicted to catalyze the synthesis of a new siderophore, which is formed under iron-deficient growth conditions.     

Burkholderia xenovorans LB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility

Burkholderia xenovorans LB400 (LB400), a well studied, effective polychlorinated biphenyl-degrader, has one of the two largest known bacterial genomes and is the first nonpathogenic Burkholderia isolate sequenced. From an evolutionary perspective, we find significant differences in functional specialization between the three replicons of LB400, as well as a more relaxed selective pressure for genes located on the two smaller vs. the largest replicon. High genomic plasticity, diversity, and specialization within the Burkholderia genus are exemplified by the conservation of only 44% of the genes between LB400 and Burkholderia cepacia complex strain 383. Even among four B. xenovorans strains, genome size varies from 7.4 to 9.73 Mbp. The latter is largely explained by our findings that >20% of the LB400 sequence was recently acquired by means of lateral gene transfer. Although a range of genetic factors associated with in vivo survival and intercellular interactions are present, these genetic factors are likely related to niche breadth rather than determinants of pathogenicity. The presence of at least eleven "central aromatic" and twenty "peripheral aromatic" pathways in LB400, among the highest in any sequenced bacterial genome, supports this hypothesis. Finally, in addition to the experimentally observed redundancy in benzoate degradation and formaldehyde oxidation pathways, the fact that 17.6% of proteins have a better LB400 paralog than an ortholog in a different genome highlights the importance of gene duplication and repeated acquirement, which, coupled with their divergence, raises questions regarding the role of paralogs and potential functional redundancies in large-genome microbes.     

Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach

Escherichia coli is a model laboratory bacterium, a species that is widely distributed in the environment, as well as a mutualist and pathogen in its human hosts. As such, E. coli represents an attractive organism to study how environment impacts microbial genome structure and function. Uropathogenic E. coli (UPEC) must adapt to life in several microbial communities in the human body, and has a complex life cycle in the bladder when it causes acute or recurrent urinary tract infection (UTI). Several studies designed to identify virulence factors have focused on genes that are uniquely represented in UPEC strains, whereas the role of genes that are common to all E. coli has received much less attention. Here we describe the complete 5,065,741-bp genome sequence of a UPEC strain recovered from a patient with an acute bladder infection and compare it with six other finished E. coli genome sequences. We searched 3,470 ortholog sets for genes that are under positive selection only in UPEC strains. Our maximum likelihood-based analysis yielded 29 genes involved in various aspects of cell surface structure, DNA metabolism, nutrient acquisition, and UTI. These results were validated by resequencing a subset of the 29 genes in a panel of 50 urinary, periurethral, and rectal E. coli isolates from patients with UTI. These studies outline a computational approach that may be broadly applicable for studying strain-specific adaptation and pathogenesis in other bacteria.     

Temporal profile of replication of human chromosomes

Chromosomes in human cancer cells are expected to initiate replication from predictably localized origins, firing reproducibly at discrete times in S phase. Replication products obtained from HeLa cells at different stages of S phase were hybridized to cDNA and genome tiling oligonucleotide microarrays to determine the temporal profile of replication of human chromosomes on a genome-wide scale. About 1,000 genes and chromosomal segments were identified as sites containing efficient origins that fire reproducibly. Early replication was correlated with high gene density. An acute transition of gene density from early to late replicating areas suggests that discrete chromatin states dictate early versus late replication. Surprisingly, at least 60% of the interrogated chromosomal segments replicate equally in all quarters of S phase, suggesting that large stretches of chromosomes are replicated by inefficient, variably located and asynchronous origins and forks, producing a pan-S phase pattern of replication. Thus, at least for aneuploid cancer cells, a typical discrete time of replication in S phase is not seen for large segments of the chromosomes.     

Genome sequence of Streptococcus mutans UA159, a cariogenic dental pathogen

Streptococcus mutans is the leading cause of dental caries (tooth decay) worldwide and is considered to be the most cariogenic of all of the oral streptococci. The genome of S. mutans UA159, a serotype c strain, has been completely sequenced and is composed of 2,030,936 base pairs. It contains 1,963 ORFs, 63% of which have been assigned putative functions. The genome analysis provides further insight into how S. mutans has adapted to surviving the oral environment through resource acquisition, defense against host factors, and use of gene products that maintain its niche against microbial competitors. S. mutans metabolizes a wide variety of carbohydrates via nonoxidative pathways, and all of these pathways have been identified, along with the associated transport systems whose genes account for almost 15% of the genome. Virulence genes associated with extracellular adherent glucan production, adhesins, acid tolerance, proteases, and putative hemolysins have been identified. Strain UA159 is naturally competent and contains all of the genes essential for competence and quorum sensing. Mobile genetic elements in the form of IS elements and transposons are prominent in the genome and include a previously uncharacterized conjugative transposon and a composite transposon containing genes for the synthesis of antibiotics of the gramicidin/bacitracin family; however, no bacteriophage genomes are present.