Pervasive contingency and entrenchment in a billion years of Hsp90 evolution

Interactions among mutations within a protein have the potential to make molecular evolution contingent and irreversible, but the extent to which epistasis actually shaped historical evolutionary trajectories is unclear. To address this question, we experimentally measured how the fitness effects of historical sequence substitutions changed during the billion-year evolutionary history of the heat shock protein 90 (Hsp90) ATPase domain beginning from a deep eukaryotic ancestor to modern Saccharomyces cerevisiae We found a pervasive influence of epistasis. Of 98 derived amino acid states that evolved along this lineage, about half compromise fitness when introduced into the reconstructed ancestral Hsp90. And the vast majority of ancestral states reduce fitness when introduced into the extant S. cerevisiae Hsp90. Overall, more than 75% of historical substitutions were contingent on permissive substitutions that rendered the derived state nondeleterious, became entrenched by subsequent restrictive substitutions that made the ancestral state deleterious, or both. This epistasis was primarily caused by specific interactions among sites rather than a general effect on the protein's tolerance to mutation. Our results show that epistasis continually opened and closed windows of mutational opportunity over evolutionary timescales, producing histories and biological states that reflect the transient internal constraints imposed by the protein's fleeting sequence states.

Energy demand and the context-dependent effects of genetic interactions underlying metabolism

Genetic effects are often context dependent, with the same genotype differentially affecting phenotypes across environments, life stages, and sexes. We used an environmental manipulation designed to increase energy demand during development to investigate energy demand as a general physiological explanation for context-dependent effects of mutations, particularly for those mutations that affect metabolism. We found that increasing the photoperiod during which Drosophila larvae are active during development phenocopies a temperature-dependent developmental delay in a mitochondrial-nuclear genotype with disrupted metabolism. This result indicates that the context-dependent fitness effects of this genotype are not specific to the effects of temperature and may generally result from variation in energy demand. The effects of this genotype also differ across life stages and between the sexes. The mitochondrial-nuclear genetic interaction disrupts metabolic rate in growing larvae, but not in adults, and compromises female, but not male, reproductive fitness. These patterns are consistent with a model where context-dependent genotype-phenotype relationships may generally arise from differences in energy demand experienced by individuals across environments, life stages, and sexes.

Association of NRG1 and AUTS2 genetic polymorphisms with Hirschsprung disease in a South Chinese population

Hirschsprung disease (HSCR) is a genetic disorder characterized by the absence of enteric ganglia. There are more than 15 genes identified as contributed to HSCR by family-based or population-based approaches. However, these findings were not fulfilled to explain the heritability of most sporadic cases. In this study, using 1470 HSCR and 1473 control subjects in South Chinese population, we replicated two variants in NRG1 (rs16879552, P = 1.05E-04 and rs7835688, P = 1.19E-07), and further clarified the two replicated SNPs were more essential for patients with short-segment aganglionosis (SHSCR) (P = 2.37E-05). We also tried to replicate the most prominent signal (rs7785360) in AUTS2, which was a potential susceptibility gene with HSCR. In our results, in terms of individual association, marginal effect was observed to affect the HSCR patients following recessive model (P = 0.089). Noteworthy, significant intergenic synergistic effect between rs16879552 (NRG1) and rs7785360 (AUTS2) was identified through cross-validation by logistic regression (P = 2.45E-03, OR = 1.53) and multifactor dimensionality reduction (MDR, P < 0.0001, OR = 1.77). Significant correlation was observed between expression of these two genes in the normal segments of the colons (P = 0.018), together with differential expression of these genes between aganglionic colonic segments and normal colonic segments of the HSCR patients (P value for AUTS2 <0.0001, P value for NRG1 = 0.0243). Although functional evaluation is required, we supply new evidence for the NRG1 to HSCR and raised up a new susceptibility gene AUTS2 to a specific symptom for the disease.

Mapping mutational effects along the evolutionary landscape of HIV envelope

The immediate evolutionary space accessible to HIV is largely determined by how single amino acid mutations affect fitness. These mutational effects can shift as the virus evolves. However, the prevalence of such shifts in mutational effects remains unclear. Here, we quantify the effects on viral growth of all amino acid mutations to two HIV envelope (Env) proteins that differ at >100 residues. Most mutations similarly affect both Envs, but the amino acid preferences of a minority of sites have clearly shifted. These shifted sites usually prefer a specific amino acid in one Env, but tolerate many amino acids in the other. Surprisingly, shifts are only slightly enriched at sites that have substituted between the Envs—and many occur at residues that do not even contact substitutions. Therefore, long-range epistasis can unpredictably shift Env’s mutational tolerance during HIV evolution, although the amino acid preferences of most sites are conserved between moderately diverged viral strains.

The virus that causes AIDS, or HIV, has a protein called Env on its surface, which is essential for the virus to infect cells. Env can also be recognized by the immune system, which then targets the virus for destruction or blocks it from infecting cells. Unfortunately, Env evolves very quickly, which means that HIV can evade our defenses. However, there are limits to how much this protein can change, since it still needs to perform its essential role in helping viruses enter cells.

In the century since HIV first appeared in human populations, the virus has evolved considerably. There are now many HIV strains that infect people, and they bear Env proteins with substantially different sequences. However, it is not clear if these changes in sequence have resulted in Envs from distinct strains being able to tolerate different mutations.

To examine this question, Haddox et al. compared how the Envs from two strains of HIV react to modifications in their sequences. They created all possible individual mutations in the proteins, and the resulting collections of mutated viruses were then tested for their ability to infect cells in the laboratory.

Most mutations had similar effects in both Env proteins. This allowed Haddox et al. to identify portions of the protein that easily accommodate changes, and portions that must remain unchanged for viruses to remain infectious—at least in the laboratory. Some of these mutations are under different types of pressures when the virus faces the immune system, and those were identified using computational approaches.

However, some mutations were tolerated differently by the two Env proteins. Therefore, viral strains differ in how their Env proteins can evolve. The parts of Env that showed differences in mutational tolerance between the strains were not necessarily the parts that differ in sequence. This shows that changes in sequence in one part of the protein can modify how other portions evolve.

It remains to be determined whether changes in tolerance to mutations translate into differences in how the virus can escape immunity. This is an important question given that the rapid evolution of Env is a major obstacle to creating a vaccine for HIV.

Interaction between rpsL and gyrA mutations affects the fitness and dual resistance of Mycobacterium tuberculosis clinical isolates against streptomycin and fluoroquinolones

Background

The interaction between different drug-resistant mutations is important to the development of drug resistance and its evolution. In this study, we aimed to reveal the potential relationships between mutations conferring resistance to two important antituberculosis drugs streptomycin (STR) and fluoroquinolones (FLQ).

Materials and methods

We used an in vitro competitive fitness assay to reveal the interactions between different mutations of rpsL and gyrA in drug-resistant Mycobacterium smegmatis, followed by the analysis of the frequency of rpsL and gyrA mutation combinations in 213 STR-FLQ dual-resistant clinical Mycobacterium tuberculosis isolates from Sichuan region, which was also investigated by the whole genome data from 3,056 global clinical M. tuberculosis isolates.

Results

The strains with K43R and K88R mutation in rpsL showed no difference in relative fitness compared with their susceptible ancestor, while K43N, K43M, K43T, and K88E exhibited a significantly lower relative fitness (P<0.05). For the FLQ-resistant mutants, all mutation types showed no difference in their relative fitness. Among STR-FLQ dual-resistant M. smegmatis strains, a lower fitness was detected in those with K43N/M/T and K88E instead of K43R and K88R mutations in rpsL. Among M. tuberculosis isolates harboring rpsL and gyrA dual mutations, the most two frequent combinatorial mutation types were K43R/D94G (n=37) and K43R/A90V (n=24), with the former being the most frequent one by both in vitro tests and clinical survey.

Conclusion

Our results suggest that the interaction between rpsL and gyrA mutations affects the fitness cost in STR-FLQ dual-resistant M. smegmatis and also the predilection of mutation combinations in clinical M. tuberculosis isolates.

Reconstructing a metazoan genetic pathway with transcriptome-wide epistasis measurements

RNA-sequencing (RNA-seq) is commonly used to identify genetic modules that respond to perturbations. In single cells, transcriptomes have been used as phenotypes, but this concept has not been applied to whole-organism RNA-seq. Also, quantifying and interpreting epistatic effects using expression profiles remains a challenge. We developed a single coefficient to quantify transcriptome-wide epistasis that reflects the underlying interactions and which can be interpreted intuitively. To demonstrate our approach, we sequenced four single and two double mutants of Caenorhabditis elegans From these mutants, we reconstructed the known hypoxia pathway. In addition, we uncovered a class of 56 genes with HIF-1-dependent expression that have opposite changes in expression in mutants of two genes that cooperate to negatively regulate HIF-1 abundance; however, the double mutant of these genes exhibits suppression epistasis. This class violates the classical model of HIF-1 regulation but can be explained by postulating a role of hydroxylated HIF-1 in transcriptional control.

Facilitating Anti-Cancer Combinatorial Drug Discovery by Targeting Epistatic Disease Genes

Due to synergistic effects, combinatorial drugs are widely used for treating complex diseases. However, combining drugs and making them synergetic remains a challenge. Genetic disease genes are considered a promising source of drug targets with important implications for navigating the drug space. Most diseases are not caused by a single pathogenic factor, but by multiple disease genes, in particular, interacting disease genes. Thus, it is reasonable to consider that targeting epistatic disease genes may enhance the therapeutic effects of combinatorial drugs. In this study, synthetic lethality gene pairs of tumors, similar to epistatic disease genes, were first targeted by combinatorial drugs, resulting in the enrichment of the combinatorial drugs with cancer treatment, which verified our hypothesis. Then, conventional epistasis detection software was used to identify epistatic disease genes from the genome wide association studies (GWAS) dataset. Furthermore, combinatorial drugs were predicted by targeting these epistatic disease genes, and five combinations were proven to have synergistic anti-cancer effects on MCF-7 cells through cell cytotoxicity assay. Combined with the three-dimensional (3D) genome-based method, the epistatic disease genes were filtered and were more closely related to disease. By targeting the filtered gene pairs, the efficiency of combinatorial drug discovery has been further improved.

Plasmid persistence: costs, benefits, and the plasmid paradox

Plasmids are extrachromosomal DNA elements that can be found throughout bacteria, as well as in other domains of life. Nonetheless, the evolutionary processes underlying the persistence of plasmids are incompletely understood. Bacterial plasmids may encode genes for traits that are sometimes beneficial to their hosts, such as antimicrobial resistance, virulence, heavy metal tolerance, and the catabolism of unique nutrient sources. In the absence of selection for these traits, however, plasmids generally impose a fitness cost on their hosts. As such, plasmid persistence presents a conundrum: models predict that costly plasmids will be lost over time or that beneficial plasmid genes will be integrated into the host genome. However, laboratory and comparative studies have shown that plasmids can persist for long periods, even in the absence of positive selection. Several hypotheses have been proposed to explain plasmid persistence, including host-plasmid co-adaptation, plasmid hitchhiking, cross-ecotype transfer, and high plasmid transfer rates, but there is no clear evidence that any one model adequately resolves the plasmid paradox.

Primer in Genetics and Genomics, Article 7-Multifactorial Concepts: Gene-Gene Interactions

Most common disorders affecting human health are not attributable to simple Mendelian (single-gene) inheritance patterns. Rather, the risk of developing a complex disease is often the result of interactions across genes, whereby one gene modifies the phenotype of another gene. These types of interactions can occur between two or more genes and are referred to as epistasis. There are five major types of epistatic interactions, but in human genetics, additive epistasis is most often discussed and includes both positive and negative subtypes. Detecting epistatic interactions can be quite difficult because seemingly unrelated genes can interact with and influence each other. As a result of this complexity, statistical geneticists are constantly developing new methods to enhance detection, but there are disadvantages to each proposed method. In this article, we explore the concept of epistasis, discuss different types of epistatic interactions, and provide a brief introduction to statistical methods researchers use to uncover sets of epistatic interactions. Then, we consider Alzheimer's disease as an exemplar for a disease with epistatic effects. Finally, we provide helpful resources, where nurses can learn more about epistasis in order to incorporate these methods into their own program of research.

Non-additive Effects in Genomic Selection

In the last decade, genomic selection has become a standard in the genetic evaluation of livestock populations. However, most procedures for the implementation of genomic selection only consider the additive effects associated with SNP (Single Nucleotide Polymorphism) markers used to calculate the prediction of the breeding values of candidates for selection. Nevertheless, the availability of estimates of non-additive effects is of interest because: (i) they contribute to an increase in the accuracy of the prediction of breeding values and the genetic response; (ii) they allow the definition of mate allocation procedures between candidates for selection; and (iii) they can be used to enhance non-additive genetic variation through the definition of appropriate crossbreeding or purebred breeding schemes. This study presents a review of methods for the incorporation of non-additive genetic effects into genomic selection procedures and their potential applications in the prediction of future performance, mate allocation, crossbreeding, and purebred selection. The work concludes with a brief outline of some ideas for future lines of that may help the standard inclusion of non-additive effects in genomic selection.

The genomic and ecological context of hybridization affects the probability that symmetrical incompatibilities drive hybrid speciation

Despite examples of homoploid hybrid species, theoretical work describing when, where, and how we expect homoploid hybrid speciation to occur remains relatively rare. Here, I explore the probability of homoploid hybrid speciation due to "symmetrical incompatibilities" under different selective and genetic scenarios. Through simulation, I test how genetic architecture and selection acting on traits that do not themselves generate incompatibilities interact to affect the probability that hybrids evolve symmetrical incompatibilities with their parent species. Unsurprisingly, selection against admixture at "adaptive" loci that are linked to loci that generate incompatibilities tends to reduce the probability of evolving symmetrical incompatibilities. By contrast, selection that favors admixed genotypes at adaptive loci can promote the evolution of symmetrical incompatibilities. The magnitude of these outcomes is affected by the strength of selection, aspects of genetic architecture such as linkage relationships and the linear arrangement of loci along a chromosome, and the amount of hybridization following the formation of a hybrid zone. These results highlight how understanding the nature of selection, aspects of the genetics of traits affecting fitness, and the strength of reproductive isolation between hybridizing taxa can all be used to inform when we expect to observe homoploid hybrid speciation due to symmetrical incompatibilities.

Weak Epistasis May Drive Adaptation in Recombining Bacteria

The impact of epistasis on the evolution of multi-locus traits depends on recombination. While sexually reproducing eukaryotes recombine so frequently that epistasis between polymorphisms is not considered to play a large role in short-term adaptation, many bacteria also recombine, some to the degree that their populations are described as "panmictic" or "freely recombining." However, whether this recombination is sufficient to limit the ability of selection to act on epistatic contributions to fitness is unknown. We quantify homologous recombination in five bacterial pathogens and use these parameter estimates in a multilocus model of bacterial evolution with additive and epistatic effects. We find that even for highly recombining species (e.g., Streptococcus pneumoniae or Helicobacter pylori), selection on weak interactions between distant mutations is nearly as efficient as for an asexual species, likely because homologous recombination typically transfers only short segments. However, for strong epistasis, bacterial recombination accelerates selection, with the dynamics dependent on the amount of recombination and the number of loci. Epistasis may thus play an important role in both the short- and long-term adaptive evolution of bacteria, and, unlike in eukaryotes, is not limited to strong effect sizes, closely linked loci, or other conditions that limit the impact of recombination.

Pharmacogenomics of Hypertension and Preeclampsia: Focus on Gene-Gene Interactions

Hypertension is a leading cause of cardiovascular mortality, but only about half of patients on antihypertensive therapy achieve blood pressure control. Preeclampsia is defined as pregnancy-induced hypertension and proteinuria, and is associated with increased maternal and perinatal mortality and morbidity. Similarly, a large number of patients with preeclampsia are non-responsive to antihypertensive therapy. Pharmacogenomics may help to guide the personalized treatment for non-responsive hypertensive patients. There is evidence for the association of genetic variants with variable response to the most commonly used antihypertensive drugs. However, further replication is needed to confirm these associations in different populations. The failure to replicate findings from single-locus association studies has prompted the search for novel statistical methods for data analysis, which are required to detect the complex effects from multiple genes to drug response phenotypes. Notably, gene–gene interaction analyses have been applied to pharmacogenetic studies, including antihypertensive drug response. In this perspective article, we present advances of considering the interactions among genetic polymorphisms of different candidate genes within pathways relevant to antihypertensive drug response, and we highlight recent findings related to gene–gene interactions on pharmacogenetics of hypertension and preeclampsia. Finally, we discuss the future directions that are needed to unravel additional genes and variants involved in the responsiveness to antihypertensive drugs.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

Toward an integrated map of genetic interactions in cancer cells

Cancer genomes often harbor hundreds of molecular aberrations. Such genetic variants can be drivers or passengers of tumorigenesis and create vulnerabilities for potential therapeutic exploitation. To identify genotype-dependent vulnerabilities, forward genetic screens in different genetic backgrounds have been conducted. We devised MINGLE, a computational framework to integrate CRISPR/Cas9 screens originating from different libraries building on approaches pioneered for genetic network discovery in model organisms. We applied this method to integrate and analyze data from 85 CRISPR/Cas9 screens in human cancer cells combining functional data with information on genetic variants to explore more than 2.1 million gene-background relationships. In addition to known dependencies, we identified new genotype-specific vulnerabilities of cancer cells. Experimental validation of predicted vulnerabilities identified GANAB and PRKCSH as new positive regulators of Wnt/β-catenin signaling. By clustering genes with similar genetic interaction profiles, we drew the largest genetic network in cancer cells to date. Our scalable approach highlights how diverse genetic screens can be integrated to systematically build informative maps of genetic interactions in cancer, which can grow dynamically as more data are included.

Ensemble survival tree models to reveal pairwise interactions of variables with time-to-events outcomes in low-dimensional setting

Unraveling interactions among variables such as genetic, clinical, demographic and environmental factors is essential to understand the development of common and complex diseases. To increase the power to detect such variables interactions associated with clinical time-to-events outcomes, we borrowed established concepts from random survival forest (RSF) models. We introduce a novel RSF-based pairwise interaction estimator and derive a randomization method with bootstrap confidence intervals for inferring interaction significance. Using various linear and nonlinear time-to-events survival models in simulation studies, we first show the efficiency of our approach: true pairwise interaction-effects between variables are uncovered, while they may not be accompanied with their corresponding main-effects, and may not be detected by standard semi-parametric regression modeling and test statistics used in survival analysis. Moreover, using a RSF-based cross-validation scheme for generating prediction estimators, we show that informative predictors may be inferred. We applied our approach to an HIV cohort study recording key host gene polymorphisms and their association with HIV change of tropism or AIDS progression. Altogether, this shows how linear or nonlinear pairwise statistical interactions of variables may be efficiently detected with a predictive value in observational studies with time-to-event outcomes.

Diversity of Functionally Permissive Sequences in the Receptor-Binding Site of Influenza Hemagglutinin

Influenza A virus hemagglutinin (HA) initiates viral entry by engaging host receptor sialylated glycans via its receptor-binding site (RBS). The amino acid sequence of the RBS naturally varies across avian and human influenza virus subtypes and is also evolvable. However, functional sequence diversity in the RBS has not been fully explored. Here, we performed a large-scale mutational analysis of the RBS of A/WSN/33 (H1N1) and A/Hong Kong/1/1968 (H3N2) HAs. Many replication-competent mutants not yet observed in nature were identified, including some that could escape from an RBS-targeted broadly neutralizing antibody. This functional sequence diversity is made possible by pervasive epistasis in the RBS 220-loop and can be buffered by avidity in viral receptor binding. Overall, our study reveals that the HA RBS can accommodate a much greater range of sequence diversity than previously thought, which has significant implications for the complex evolutionary interrelationships between receptor specificity and immune escape.

Conflict between heterozygote advantage and hybrid incompatibility in haplodiploids (and sex chromosomes)

In many diploid species, the sex chromosomes play a special role in mediating reproductive isolation. In haplodiploids, where females are diploid and males haploid, the whole genome behaves similarly to the X/Z chromosomes of diploids. Therefore, haplodiploid systems can serve as a model for the role of sex chromosomes in speciation and hybridization. A previously described population of Finnish Formica wood ants displays genome-wide signs of ploidally and sexually antagonistic selection resulting from hybridization. Here, hybrid females have increased survivorship but hybrid males are inviable. To understand how the unusual hybrid population may be maintained, we developed a mathematical model with hybrid incompatibility, female heterozygote advantage, recombination and assortative mating. The rugged fitness landscape resulting from the co-occurrence of heterozygote advantage and hybrid incompatibility results in a sexual conflict in haplodiploids, which is caused by the ploidy difference. Thus, whereas heterozygote advantage always promotes long-term polymorphism in diploids, we find various outcomes in haplodiploids in which the population stabilizes either in favour of males, females or via maximizing the number of introgressed individuals. We discuss these outcomes with respect to the potential long-term fate of the Finnish wood ant population and provide approximations for the extension of the model to multiple incompatibilities. Moreover, we highlight the general implications of our results for speciation and hybridization in haplodiploids versus diploids and how the described fitness relationships could contribute to the outstanding role of sex chromosomes as hotspots of sexual antagonism and genes involved in speciation.

Genetic characterization and fine mapping of S25, a hybrid male sterility gene, on rice chromosome 12

Hybrid male sterility genes are important factors in creating postzygotic reproductive isolation barriers in plants. One such gene, S25, is known to cause severe transmission ratio distortion in inter-subspecific progeny of cultivated rice Oryza sativa ssp. indica and japonica. To further characterize the S25 gene, we fine-mapped and genetically characterized the S25 gene using near-isogenic lines with reciprocal genetic backgrounds. We mapped the S25 locus within the 0.67-1.02 Mb region on rice chromosome 12. Further genetic analyses revealed that S25 substantially reduced male fertility in the japonica background, but not in the indica background. In first-generation hybrid progeny, S25 had a milder effect than it had in the japonica background. These results suggest that the expression of S25 is epistatically regulated by at least one partially dominant gene present in the indica genome. This finding supports our previous studies showing that hybrid male sterility due to pollen killer genes results from epistatic interaction with other genes that are hidden in the genetic background.

Gene-gene interactions in gastrointestinal cancer susceptibility

Cancer arises from complex, multi-layer interactions between diverse genetic and environmental factors. Genetic studies have identified multiple loci associated with tumor susceptibility. However, little is known about how germline polymorphisms interact with one another and with somatic mutations within a tumor to mediate acquisition of cancer traits. Here, we survey recent studies showing gene-gene interactions, also known as epistases, affecting genetic susceptibility in colorectal, gastric and esophageal cancers. We also catalog epistasis types and cancer hallmarks with respect to the interacting genes. A total of 22 gene variation pairs displayed all levels of statistical epistasis, including synergistic, redundant, suppressive and co-suppressive interactions. Five genes primarily involved in base excision repair formed a linear topology in the interaction network, MUTYH-OGG1-XRCC1-PARP1-MMP2, and three genes in mTOR cell-proliferation pathway formed another linear network, PRKAG2-RPS6KB1-PIK3CA. Discrete pairwise epistasis was also found in nucleotide excision repair, detoxification, proliferation, TP53, TGF-β and other pathways. We propose that three modes of biological interaction underlie the molecular mechanisms for statistical epistasis. The direct binding, linear pathway and convergence modes can exhibit any level of statistical epistasis in susceptibility to gastrointestinal cancers, and this is likely true for other complex diseases as well. This review highlights the link between cancer hallmarks and susceptibility genes.

Slow and temperature-mediated pathogen adaptation to a nonspecific fungicide in agricultural ecosystem

The spread of antimicrobial resistance and global change in air temperature represent two major phenomena that are exerting a disastrous impact on natural and social issues but investigation of the interaction between these phenomena in an evolutionary context is limited. In this study, a statistical genetic approach was used to investigate the evolution of antimicrobial resistance in agricultural ecosystem and its association with local air temperature, precipitation, and UV radiation. We found no resistance to mancozeb, a nonspecific fungicide widely used in agriculture for more than half a century, in 215 Alternaria alternata isolates sampled from geographic locations along a climatic gradient and cropping system representing diverse ecotypes in China, consistent with low resistance risk in many nonspecific fungicides. Genetic variance accounts for ~35% of phenotypic variation, while genotype-environment interaction is negligible, suggesting that heritability plays a more important role in the evolution of resistance to mancozeb in plant pathogens than phenotypic plasticity. We also found that tolerance to mancozeb in agricultural ecosystem is under constraining selection and significantly associated with local air temperature, possibly resulting from a pleiotropic effect of resistance with thermal and other ecological adaptations. The implication of these results for fungicide and other antimicrobial management in the context of global warming is discussed.

Fitness Tradeoffs of Antibiotic Resistance in Extraintestinal Pathogenic Escherichia coli

Evolutionary trade-offs occur when selection on one trait has detrimental effects on other traits. In pathogenic microbes, it has been hypothesized that antibiotic resistance trades off with fitness in the absence of antibiotic. Although studies of single resistance mutations support this hypothesis, it is unclear whether trade-offs are maintained over time, due to compensatory evolution and broader effects of genetic background. Here, we leverage natural variation in 39 extraintestinal clinical isolates of Escherichia coli to assess trade-offs between growth rates and resistance to fluoroquinolone and cephalosporin antibiotics. Whole-genome sequencing identifies a broad range of clinically relevant resistance determinants in these strains. We find evidence for a negative correlation between growth rate and antibiotic resistance, consistent with a persistent trade-off between resistance and growth. However, this relationship is sometimes weak and depends on the environment in which growth rates are measured. Using in vitro selection experiments, we find that compensatory evolution in one environment does not guarantee compensation in other environments. Thus, even in the face of compensatory evolution and other genetic background effects, resistance may be broadly costly, supporting the use of drug restriction protocols to limit the spread of resistance. Furthermore, our study demonstrates the power of using natural variation to study evolutionary trade-offs in microbes.

Coadaptation of mitochondrial and nuclear genes, and the cost of mother's curse

Strict maternal inheritance renders the mitochondrial genome susceptible to accumulating mutations that harm males, but are otherwise benign or beneficial for females. This 'mother's curse' effect can degrade male survival and fertility if unopposed by counteracting evolutionary processes. Coadaptation between nuclear and mitochondrial genomes-with nuclear genes evolving to compensate for male-harming mitochondrial substitutions-may ultimately resolve mother's curse. However, males are still expected to incur a transient fitness cost during mito-nuclear coevolution, and it remains unclear how severe such costs should be. We present a population genetic analysis of mito-nuclear coadaptation to resolve mother's curse effects, and show that the magnitude of the 'male mitochondrial load'-the negative impact of mitochondrial substitutions on male fitness components-may be large, even when genetic variation for compensatory evolution is abundant. We also find that the male load is surprisingly sensitive to population size: male fitness costs of mito-nuclear coevolution are particularly pronounced in both small and large populations, and minimized in populations of intermediate size. Our results reveal complex interactions between demography and genetic constraints during the resolution of mother's curse, suggesting potentially widespread species differences in susceptibility to mother's curse effects.

Epistatic interactions influence terrestrial-marine functional shifts in cetacean rhodopsin

Like many aquatic vertebrates, whales have blue-shifting spectral tuning substitutions in the dim-light visual pigment, rhodopsin, that are thought to increase photosensitivity in underwater environments. We have discovered that known spectral tuning substitutions also have surprising epistatic effects on another function of rhodopsin, the kinetic rates associated with light-activated intermediates. By using absorbance spectroscopy and fluorescence-based retinal release assays on heterologously expressed rhodopsin, we assessed both spectral and kinetic differences between cetaceans (killer whale) and terrestrial outgroups (hippo, bovine). Mutation experiments revealed that killer whale rhodopsin is unusually resilient to pleiotropic effects on retinal release from key blue-shifting substitutions (D83N and A292S), largely due to a surprisingly specific epistatic interaction between D83N and the background residue, S299. Ancestral sequence reconstruction indicated that S299 is an ancestral residue that predates the evolution of blue-shifting substitutions at the origins of Cetacea. Based on these results, we hypothesize that intramolecular epistasis helped to conserve rhodopsin's kinetic properties while enabling blue-shifting spectral tuning substitutions as cetaceans adapted to aquatic environments. Trade-offs between different aspects of molecular function are rarely considered in protein evolution, but in cetacean and other vertebrate rhodopsins, may underlie multiple evolutionary scenarios for the selection of specific amino acid substitutions.

Clinical Implications of Digenic Inheritance and Epistasis in Primary Immunodeficiency Disorders

The existence of epistasis in humans was first predicted by Bateson in 1909. Epistasis describes the non-linear, synergistic interaction of two or more genetic loci, which can substantially modify disease severity or result in entirely new phenotypes. The concept has remained controversial in human genetics because of the lack of well-characterized examples. In humans, it is only possible to demonstrate epistasis if two or more genes are mutated. In most cases of epistasis, the mutated gene products are likely to be constituents of the same physiological pathway leading to severe disruption of a cellular function such as antibody production. We have recently described a digenic family, who carry mutations of TNFRSF13B/TACI as well as TCF3 genes. Both genes lie in tandem along the immunoglobulin isotype switching and secretion pathway. We have shown they interact in an epistatic way causing severe immunodeficiency and autoimmunity in the digenic proband. With the advent of next generation sequencing, it is likely other families with digenic inheritance will be identified. Since digenic inheritance does not always cause epistasis, we propose an epistasis index which may help quantify the effects of the two mutations. We also discuss the clinical implications of digenic inheritance and epistasis in humans with primary immunodeficiency disorders.