Flavonols improve thermotolerance in tomato pollen during germination and tube elongation by maintaining ROS homeostasis.

Elevated temperatures impair pollen performance and reproductive success, resulting in lower crop yields. The Solanum lycopersicum anthocyanin reduced ( are ) mutant has a FLAVANONE 3 HYDROXYLASE ( F3H ) gene mutation resulting in impaired synthesis of flavonol antioxidants. The are mutant has reduced pollen performance and seed set relative to the VF36 parental line, which is accentuated at elevated temperatures. Transformation of are with the wild-type F3H gene, or chemical complementation with flavonols, prevented temperature-dependent ROS accumulation in pollen and reversed are's reduced viability, germination, and tube elongation to VF36 levels. VF36 transformed with an F3H overexpression construct prevented temperature driven ROS increases and impaired pollen performance, revealing thermotolerance results from elevated flavonol synthesis. Although stigmas of are had reduced flavonols and elevated ROS, the growth of are pollen tubes were similarly impaired in both are and VF36 pistils. RNA-Seq was performed at optimal and stress temperatures in are , VF36, and the VF36 F3H overexpression line at multiple timepoints across pollen tube elongation. Differentially expressed gene numbers increased with duration of elevated temperature in all genotypes, with the largest number in are . These findings suggest potential agricultural interventions to combat the negative effects of heat-induced ROS in pollen that leads to reproductive failure. One sentence summary: Flavonol antioxidants reduce the negative impacts of elevated temperatures on pollen performance by reducing levels of heat induced reactive oxygen species and modulation of heat-induced changes in the pollen transcriptome.

Layered evolution of gene expression in "superfast" muscles for courtship.

Identifying the molecular process of complex trait evolution is a core goal of biology. However, pinpointing the specific context and timing of trait-associated changes within the molecular evolutionary history of an organism remains an elusive goal. We study this topic by exploring the molecular basis of elaborate courtship evolution, which represents an extraordinary example of trait innovation. Within the behaviorally diverse radiation of Central and South American manakin birds, species from two separate lineages beat their wings together using specialized "superfast" muscles to generate a "snap" that helps attract mates. Here, we develop an empirical approach to analyze phylogenetic lineage-specific shifts in gene expression in the key snap-performing muscle and then integrate these findings with comparative transcriptomic sequence analysis. We find that rapid wing displays are associated with changes to a wide range of molecular processes that underlie extreme muscle performance, including changes to calcium trafficking, myocyte homeostasis and metabolism, and hormone action. We furthermore show that these changes occur gradually in a layered manner across the species history, wherein which ancestral genetic changes to many of these molecular systems are built upon by later species-specific shifts that likely finalized the process of display performance adaptation. Our study demonstrates the potential for combining phylogenetic modeling of tissue-specific gene expression shifts with phylogenetic analysis of lineage-specific sequence changes to reveal holistic evolutionary histories of complex traits.

The Catalytic-Dependent and -Independent Roles of Lsd1 and Lsd2 Lysine Demethylases in Heterochromatin Formation in Schizosaccharomyces pombe.

In eukaryotes, heterochromatin plays a critical role in organismal development and cell fate acquisition, through regulating gene expression. The evolutionarily conserved lysine-specific demethylases, Lsd1 and Lsd2, remove mono- and dimethylation on histone H3, serving complex roles in gene expression. In the fission yeast Schizosaccharomyces pombe, null mutations of Lsd1 and Lsd2 result in either severe growth defects or inviability, while catalytic inactivation causes minimal defects, indicating that Lsd1 and Lsd2 have essential functions beyond their known demethylase activity. Here, we show that catalytic mutants of Lsd1 or Lsd2 partially assemble functional heterochromatin at centromeres in RNAi-deficient cells, while the C-terminal truncated alleles of Lsd1 or Lsd2 exacerbate heterochromatin formation at all major heterochromatic regions, suggesting that Lsd1 and Lsd2 repress heterochromatic transcripts through mechanisms both dependent on and independent of their catalytic activities. Lsd1 and Lsd2 are also involved in the establishment and maintenance of heterochromatin. At constitutive heterochromatic regions, Lsd1 and Lsd2 regulate one another and cooperate with other histone modifiers, including the class II HDAC Clr3 and the Sirtuin family protein Sir2 for gene silencing, but not with the class I HDAC Clr6. Our findings explore the roles of lysine-specific demethylases in epigenetic gene silencing at heterochromatic regions.

A Deep-sequencing-assisted, Spontaneous Suppressor Screen in the Fission Yeast Schizosaccharomyces pombe.

A genetic screen for mutant alleles that suppress phenotypic defects caused by a mutation is a powerful approach to identify genes that belong to closely related biochemical pathways. Previous methods such as the Synthetic Genetic Array (SGA) analysis, and random mutagenesis techniques using ultraviolet (UV) or chemicals like ethyl methanesulfonate (EMS) or N-ethyl-N- nitrosourea (ENU), have been widely used but are often costly and laborious. Also, these mutagen-based screening methods are frequently associated with severe side effects on the organism, inducing multiple mutations that add to the complexity of isolating the suppressors. Here, we present a simple and effective protocol to identify suppressor mutations in mutants which confer a growth defect in Schizosaccharomyces pombe. The fitness of cells with a growth deficiency in standard rich liquid media or synthetic liquid media can be monitored for recovery using an automated 96-well plate reader over an extended period. Once a cell acquires a suppressor mutation in the culture, its descendants outcompete those of the parental cells. The recovered cells that have a competitive growth advantage over the parental cells can then be isolated and backcrossed with the parental cells. The suppressor mutations are then identified using whole-genome sequencing. Using this approach, we have successfully isolated multiple suppressors that alleviate the severe growth defects caused by loss of Elf1, an AAA+ family ATPase that is important in nuclear mRNA transport and maintenance of genomic stability. There are currently over 400 genes in S. pombe with mutants conferring a growth defect. As many of these genes are uncharacterized, we propose that our method will hasten the identification of novel functional interactions with this user-friendly, high-throughput approach.

Group A Rotavirus VP1 Polymerase and VP2 Core Shell Proteins: Intergenotypic Sequence Variation and In Vitro Functional Compatibility.

Group A rotaviruses (RVAs) are classified according to a nucleotide sequence-based system that assigns a genotype to each of the 11 double-stranded RNA (dsRNA) genome segments. For the segment encoding the VP1 polymerase, 22 genotypes (R1 to R22) are defined with an 83% nucleotide identity cutoff value. For the segment encoding the VP2 core shell protein, which is a functional VP1-binding partner, 20 genotypes (C1 to C20) are defined with an 84% nucleotide identity cutoff value. However, the extent to which the VP1 and VP2 proteins encoded by these genotypes differ in their sequences or interactions has not been described. Here, we sought to (i) delineate the relationships and sites of variation for VP1 and VP2 proteins belonging to the known RVA genotypes and (ii) correlate intergenotypic sequence diversity with functional VP1-VP2 interaction(s) during dsRNA synthesis. Using bioinformatic approaches, we revealed which VP1 and VP2 genotypes encode divergent proteins and identified the positional locations of amino acid changes in the context of known structural domains/subdomains. We then employed an in vitro dsRNA synthesis assay to test whether genotype R1, R2, R4, and R7 VP1 polymerases could be enzymatically activated by genotype C1, C2, C4, C5, and C7 VP2 core shell proteins. Genotype combinations that were incompatible informed the rational design and in vitro testing of chimeric mutant VP1 and VP2 proteins. The results of this study connect VP1 and VP2 nucleotide-level diversity to protein-level diversity for the first time, and they provide new insights into regions/residues critical for VP1-VP2 interaction(s) during viral genome replication.IMPORTANCE Group A rotaviruses (RVAs) are widespread in nature, infecting numerous mammalian and avian hosts and causing severe gastroenteritis in human children. RVAs are classified using a system that assigns a genotype to each viral gene according to its nucleotide sequence. To date, 22 genotypes have been described for the gene encoding the viral polymerase (VP1), and 20 genotypes have been described for the gene encoding the core shell protein (VP2). Here, we analyzed if/how the VP1 and VP2 proteins encoded by the known RVA genotypes differ from each other in their sequences. We also used a biochemical approach to test whether the intergenotypic sequence differences influenced how VP1 and VP2 functionally engage each other to mediate RNA synthesis in a test tube. This work is important because it increases our understanding of RVA protein-level diversity and raises new ideas about the VP1-VP2 binding interface(s) that is important for viral replication.

Loss of Elongation-Like Factor 1 Spontaneously Induces Diverse, RNase H-Related Suppressor Mutations in Schizosaccharomyces pombe.

A healthy individual may carry a detrimental genetic trait that is masked by another genetic mutation. Such suppressive genetic interactions, in which a mutant allele either partially or completely restores the fitness defect of a particular mutant, tend to occur between genes that have a confined functional connection. Here we investigate a self-recovery phenotype in Schizosaccharomyces pombe, mediated by suppressive genetic interactions that can be amplified during cell culture. Cells without Elf1, an AAA+ family ATPase, have severe growth defects initially, but quickly recover growth rates near to those of wild-type strains by acquiring suppressor mutations. elf1Δ cells accumulate RNAs within the nucleus and display effects of genome instability such as sensitivity to DNA damage, increased incidence of lagging chromosomes, and mini-chromosome loss. Notably, the rate of phenotypic recovery was further enhanced in elf1Δ cells when RNase H activities were abolished and significantly reduced upon overexpression of RNase H1, suggesting that loss of Elf1-related genome instability can be resolved by RNase H activities, likely through eliminating the potentially mutagenic DNA-RNA hybrids caused by RNA nuclear accumulation. Using whole genome sequencing, we mapped a few consistent suppressors of elf1Δ including mutated Cue2, Rpl2702, and SPBPJ4664.02, suggesting previously unknown functional connections between Elf1 and these proteins. Our findings describe a mechanism by which cells bearing mutations that cause fitness defects and genome instability may accelerate the fitness recovery of their population through quickly acquiring suppressors. We propose that this mechanism may be universally applicable to all microorganisms in large-population cultures.

Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of life.

PREMISE OF THE STUDY: Phylogenetic support has been difficult to evaluate within the green plant tree of life partly due to a lack of specificity between conflicted versus poorly informed branches. As data sets continue to expand in both breadth and depth, new support measures are needed that are more efficient and informative. METHODS: We describe the Quartet Sampling (QS) method, a quartet-based evaluation system that synthesizes several phylogenetic and genomic analytical approaches. QS characterizes discordance in large-sparse and genome-wide data sets, overcoming issues of alignment sparsity and distinguishing strong conflict from weak support. We tested QS with simulations and recent plant phylogenies inferred from variously sized data sets. KEY RESULTS: QS scores demonstrated convergence with increasing replicates and were not strongly affected by branch depth. Patterns of QS support from different phylogenies led to a coherent understanding of ancestral branches defining key disagreements, including the relationships of Ginkgo to cycads, magnoliids to monocots and eudicots, and mosses to liverworts. The relationships of ANA-grade angiosperms (Amborella, Nymphaeales, Austrobaileyales), major monocot groups, bryophytes, and fern families are likely highly discordant in their evolutionary histories, rather than poorly informed. QS can also detect discordance due to introgression in phylogenomic data. CONCLUSIONS: Quartet Sampling is an efficient synthesis of phylogenetic tests that offers more comprehensive and specific information on branch support than conventional measures. The QS method corroborates growing evidence that phylogenomic investigations that incorporate discordance testing are warranted when reconstructing complex evolutionary histories, in particular those surrounding ANA-grade, monocots, and nonvascular plants.

Encoding Data Using Biological Principles: The Multisample Variant Format for Phylogenomics and Population Genomics.

Rapid progress in the fields of phylogenomics and population genomics has driven increases in both the size of multi-genomic datasets and the number and complexity of genome-wide analyses. We present the Multisample Variant Format, specifically designed to store multiple sequence alignments for phylogenomics and population genomic analysis. The signature feature of MVF is a distinctive encoding of aligned sites with specific biological information content (e.g., invariant, low-coverage). This biological pattern-based encoding of sequence data allows for rapid filtering and quality control of data and speeds up computation for many analyses. Similar to other modern formats, MVF has a simple data structure and flexible header structure to accommodate project metadata, allowing to also serve as an effective data publication and sharing format. We also propose several variants of the MVF format to accommodate protein and codon alignments, quality scores, and a mix of de novo and reference-aligned data. Using the MVFtools package, MVF files can be converted from other common sequence formats. MVFtools completes tasks ranging from simple transformation and filtering operations to complex genome-wide visualizations in only a few minutes, even on large datasets. In addition to presentation of MVF and MVFtools, we also discuss the application both in MVF and other existing data formats of the broader concept of using biological principles and patterns to inform sequence data encoding.

Heterogeneous molecular processes among the causes of how sequence similarity scores can fail to recapitulate phylogeny.

Sequence similarity tools like Basic Local Alignment Search Tool (BLAST) are essential components of many functional genetic, genomic, phylogenetic and bioinformatic studies. Many modern analysis pipelines use significant sequence similarity scores (p- or E-values) and the ranked order of BLAST matches to test a wide range of hypotheses concerning homology, orthology, the timing of de novo gene birth/death and gene family expansion/contraction. Despite significant contrary findings, many of these tests still implicitly assume that stronger or higher-ranked E-value scores imply closer phylogenetic relationships between sequences. Here, we demonstrate that even though a general relationship does exist between the phylogenetic distance of two sequences and their E-value, significant and misleading errors occur in both the completeness and the order of results under realistic evolutionary scenarios. These results provide additional details to past evidence showing that studies should avoid drawing direct inferences of evolutionary relatedness from measures of sequence similarity alone, and should instead, where possible, use more rigorous phylogeny-based methods.

Molecular mechanisms of postmating prezygotic reproductive isolation uncovered by transcriptome analysis.

Little is known about the physiological responses and genetic mutations associated with reproductive isolation between species, especially for postmating prezygotic isolating barriers. Here, we examine changes in gene expression that accompany the expression of 'unilateral incompatibility' (UI)-a postmating prezygotic barrier in which fertilization is prevented by gamete rejection in the reproductive tract [in this case of pollen tubes (male gametophytes)] in one direction of a species cross, but is successful in the reciprocal crossing direction. We use whole-transcriptome sequencing of multiple developmental stages of male and female tissues in two Solanum species that exhibit UI to: (i) identify transcript differences between UI-competent and UI noncompetent tissues; (ii) characterize transcriptional changes specifically associated with the phenotypic expression of UI; and (iii) using these comparisons, evaluate the behaviour of a priori candidate loci for UI and identify new candidates for future manipulative work. In addition to describing transcriptome-wide changes in gene expression that accompany this isolating barrier, we identify at least five strong candidates for involvement in postmating prezygotic incompatibility between species. These include three novel candidates and two candidates that are strongly supported by prior developmental, functional, and quantitative trait locus mapping studies. These latter genes are known molecular players in the intraspecific expression of mate choice via genetic self-incompatibility, and our study supports prior evidence that these inter- and intraspecific postmating prezygotic reproductive behaviours share specific genetic and molecular mechanisms.

Powerful methods for detecting introgressed regions from population genomic data.

Understanding the types and functions of genes that are able to cross species boundaries-and those that are not-is an important step in understanding the forces maintaining species as largely independent lineages across the remainder of the genome. With large next-generation sequencing data sets we are now able to ask whether introgression has occurred across the genome, and multiple methods have been proposed to detect the signature of such events. Here, we introduce a new summary statistic that can be used to test for introgression, RNDmin , that makes use of the minimum pairwise sequence distance between two population samples relative to divergence to an outgroup. We find that our method offers a modest increase in power over other, related tests, but that all such tests have high power to detect introgressed loci when migration is recent and strong. RNDmin is robust to variation in the mutation rate, and remains reliable even when estimates of the divergence time between sister species are inaccurate. We apply RNDmin to population genomic data from the African mosquitoes Anopheles quadriannulatus and A. arabiensis, identifying three novel candidate regions for introgression. Interestingly, one of the introgressed loci is on the X chromosome, but outside of an inversion separating these two species. Our results suggest that significant, but rare, sharing of alleles is occurring between species that diverged more than 1 million years ago, and that application of these methods to additional systems are likely to reveal similar results.

Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation.

Speciation events often occur in rapid bursts of diversification, but the ecological and genetic factors that promote these radiations are still much debated. Using whole transcriptomes from all 13 species in the ecologically and reproductively diverse wild tomato clade (Solanum sect. Lycopersicon), we infer the species phylogeny and patterns of genetic diversity in this group. Despite widespread phylogenetic discordance due to the sorting of ancestral variation, we date the origin of this radiation to approximately 2.5 million years ago and find evidence for at least three sources of adaptive genetic variation that fuel diversification. First, we detect introgression both historically between early-branching lineages and recently between individual populations, at specific loci whose functions indicate likely adaptive benefits. Second, we find evidence of lineage-specific de novo evolution for many genes, including loci involved in the production of red fruit color. Finally, using a "PhyloGWAS" approach, we detect environment-specific sorting of ancestral variation among populations that come from different species but share common environmental conditions. Estimated across the whole clade, small but substantial and approximately equal fractions of the euchromatic portion of the genome are inferred to contribute to each of these three sources of adaptive genetic variation. These results indicate that multiple genetic sources can promote rapid diversification and speciation in response to new ecological opportunity, in agreement with our emerging phylogenomic understanding of the complexity of both ancient and recent species radiations.

Detection and Polarization of Introgression in a Five-Taxon Phylogeny.

When multiple speciation events occur rapidly in succession, discordant genealogies due to incomplete lineage sorting (ILS) can complicate the detection of introgression. A variety of methods, including the [Formula: see text]-statistic (a.k.a. the "ABBA-BABA test"), have been proposed to infer introgression in the presence of ILS for a four-taxon clade. However, no integrated method exists to detect introgression using allelic patterns for more complex phylogenies. Here we explore the issues associated with previous systems of applying [Formula: see text]-statistics to a larger tree topology, and propose new [Formula: see text] tests as an integrated framework to infer both the taxa involved in and the direction of introgression for a symmetric five-taxon phylogeny. Using theory and simulations, we show that the [Formula: see text] statistics correctly identify the introgression donor and recipient lineages, even at low rates of introgression. [Formula: see text] is also shown to have extremely low false-positive rates. The [Formula: see text] tests are computationally inexpensive to calculate and can easily be applied to phylogenomic data sets, both genome-wide and in windows of the genome. In addition, we explore both the principles and problems of introgression detection in even more complex phylogenies.

Mosquito genomics. Extensive introgression in a malaria vector species complex revealed by phylogenomics.

Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here, we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between nonsister species.

Quantitative genetic analysis indicates natural selection on leaf phenotypes across wild tomato species (Solanum sect. Lycopersicon; Solanaceae).

Adaptive evolution requires both raw genetic material and an accessible path of high fitness from one fitness peak to another. In this study, we used an introgression line (IL) population to map quantitative trait loci (QTL) for leaf traits thought to be associated with adaptation to precipitation in wild tomatoes (Solanum sect. Lycopersicon; Solanaceae). A QTL sign test showed that several traits likely evolved under directional natural selection. Leaf traits correlated across species do not share a common genetic basis, consistent with a scenario in which selection maintains trait covariation unconstrained by pleiotropy or linkage disequilibrium. Two large effect QTL for stomatal distribution colocalized with key genes in the stomatal development pathway, suggesting promising candidates for the molecular bases of adaptation in these species. Furthermore, macroevolutionary transitions between vastly different stomatal distributions may not be constrained when such large-effect mutations are available. Finally, genetic correlations between stomatal traits measured in this study and data on carbon isotope discrimination from the same ILs support a functional hypothesis that the distribution of stomata affects the resistance to CO2 diffusion inside the leaf, a trait implicated in climatic adaptation in wild tomatoes. Along with evidence from previous comparative and experimental studies, this analysis indicates that leaf traits are an important component of climatic niche adaptation in wild tomatoes and demonstrates that some trait transitions between species could have involved few, large-effect genetic changes, allowing rapid responses to new environmental conditions.

More accurate phylogenies inferred from low-recombination regions in the presence of incomplete lineage sorting.

When speciation events occur in rapid succession, incomplete lineage sorting (ILS) can cause disagreement among individual gene trees. The probability that ILS affects a given locus is directly related to its effective population size (Ne ), which in turn is proportional to the recombination rate if there is strong selection across the genome. Based on these expectations, we hypothesized that low-recombination regions of the genome, as well as sex chromosomes and nonrecombining chromosomes, should exhibit lower levels of ILS. We tested this hypothesis in phylogenomic datasets from primates, the Drosophila melanogaster clade, and the Drosophila simulans clade. In all three cases, regions of the genome with low or no recombination showed significantly stronger support for the putative species tree, although results from the X chromosome differed among clades. Our results suggest that recurrent selection is acting in these low-recombination regions, such that current levels of diversity also reflect past decreases in the effective population size at these same loci. The results also demonstrate how considering the genomic context of a gene tree can assist in more accurate determination of the true species phylogeny, especially in cases where a whole-genome phylogeny appears to be an unresolvable polytomy.

Sex chromosomes evolved from independent ancestral linkage groups in winged insects.

The evolution of a pair of chromosomes that differ in appearance between males and females (heteromorphic sex chromosomes) has occurred repeatedly across plants and animals. Recent work has shown that the male heterogametic (XY) and female heterogametic (ZW) sex chromosomes evolved independently from different pairs of homomorphic autosomes in the common ancestor of birds and mammals but also that X and Z chromosomes share many convergent molecular features. However, little is known about how often heteromorphic sex chromosomes have either evolved convergently from different autosomes or in parallel from the same pair of autosomes and how universal patterns of molecular evolution on sex chromosomes really are. Among winged insects with sequenced genomes, there are male heterogametic species in both the Diptera (e.g., Drosophila melanogaster) and the Coleoptera (Tribolium castaneum), female heterogametic species in the Lepidoptera (Bombyx mori), and haplodiploid species in the Hymenoptera (e.g., Nasonia vitripennis). By determining orthologous relationships among genes on the X and Z chromosomes of insects with sequenced genomes, we are able to show that these chromosomes are not homologous to one another but are homologous to autosomes in each of the other species. These results strongly imply that heteromorphic sex chromosomes have evolved independently from different pairs of ancestral chromosomes in each of the insect orders studied. We also find that the convergently evolved X chromosomes of Diptera and Coleoptera share genomic features with each other and with vertebrate X chromosomes, including excess gene movement from the X to the autosomes. However, other patterns of molecular evolution--such as increased codon bias, decreased gene density, and the paucity of male-biased genes on the X--differ among the insect X and Z chromosomes. Our results provide evidence for both differences and nearly universal similarities in patterns of evolution among independently derived sex chromosomes.

No excess gene movement is detected off the avian or lepidopteran Z chromosome.

Most of our knowledge of sex-chromosome evolution comes from male heterogametic (XX/XY) taxa. With the genome sequencing of multiple female heterogametic (ZZ/ZW) taxa, we can now ask whether there are patterns of evolution common to both sex chromosome systems. In all XX/XY systems examined to date, there is an excess of testis-biased retrogenes moving from the X chromosome to the autosomes, which is hypothesized to result from either sexually antagonistic selection or escape from meiotic sex chromosome inactivation (MSCI). We examined RNA-mediated (retrotransposed) and DNA-mediated gene movement in two independently evolved ZZ/ZW systems, birds (chicken and zebra finch) and lepidopterans (silkworm). Even with sexually antagonistic selection likely operating in both taxa and MSCI having been identified in the chicken, we find no evidence for an excess of genes moving from the Z chromosome to the autosomes in either lineage. We detected no excess for either RNA- or DNA-mediated duplicates, across a range of approaches and methods. We offer some potential explanations for this difference between XX/XY and ZZ/ZW sex chromosome systems, but further work is needed to distinguish among these hypotheses. Regardless of the root causes, we have identified an additional, potentially inherent, difference between XX/XY and ZZ/ZW systems.

On the possible role of robustness in the evolution of infectious diseases.

Robustness describes the capacity for a biological system to remain canalized despite perturbation. Genetic robustness affords maintenance of phenotype despite mutational input, necessarily involving the role of epistasis. Environmental robustness is phenotypic constancy in the face of environmental variation, where epistasis may be uninvolved. Here we discuss genetic and environmental robustness, from the standpoint of infectious disease evolution, and suggest that robustness may be a unifying principle for understanding how different disease agents evolve. We focus especially on viruses with RNA genomes due to their importance in the evolution of emerging diseases and as model systems to test robustness theory. We present new data on adaptive constraints for a model RNA virus challenged to evolve in response to UV radiation. We also draw attention to other infectious disease systems where robustness theory may prove useful for bridging evolutionary biology and biomedicine, especially the evolution of antibiotic resistance in bacteria, immune evasion by influenza, and malaria parasite infections.