Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

Molecular evolution of male reproduction across species with highly divergent sperm morphology in diverse murine rodents

Sperm competition can drive rapid evolution of male reproductive traits, but it remains unclear how variation in sperm competition intensity shapes phenotypic and molecular diversity across clades. Old World mice and rats (subfamily Murinae) comprise a rapid radiation and exhibit incredible diversity in sperm morphology and production. We combined phenotype and sequence data to model the evolution of reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass, a trait reflective of reduced sperm competition. Several sperm traits were associated with relative testes mass, suggesting that mating system evolution likely selects for convergent traits related to sperm competitive ability. Molecular evolutionary rates of spermatogenesis proteins also correlated with relative testes mass, but in an unexpected direction. We predicted that sperm competition would result in rapid divergence among species with large relative testes mass, but instead found that many spermatogenesis genes evolve more rapidly in species with smaller relative testes mass due to relaxed purifying selection. While some reproductive genes evolved under positive selection, relaxed selection played a greater role underlying rapid evolution in small testes species. Our work demonstrates that sexual selection can impose strong purifying selection shaping the evolution of male reproduction.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage

Exposure to ultraviolet radiation (UVR) is harmful to living cells, leading organisms to evolve protective mechanisms against UVR-induced cellular damage and stress.1,2 UVR, particularly UVB (280-320 nm), can damage proteins and DNA, leading to errors during DNA repair and replication. Excessive UVR can induce cellular death. Aquatic organisms face risk of UV exposure as biologically harmful levels of UVB can penetrate >10 m in clear water.3 While melanin is the only known sunscreen in vertebrates, it often emerges late in embryonic development, rendering embryos of many species vulnerable during the earlier stages. Algae and microbes produce a class of sunscreening compounds known as mycosporine-like amino acids (MAAs).4 Fish eggs contain a similar compound called gadusol, whose role as a sunscreen has yet to be tested despite its discovery over 40 years ago.5 The recent finding that many vertebrate genomes contain a biosynthetic pathway for gadusol suggests that many fish may produce and use this molecule as a sunscreen.6 We generated a gadusol-deficient mutant zebrafish to investigate the role of gadusol in protecting fish embryos and larvae from UVR. Our results demonstrate that maternally provided gadusol is the primary sunscreen in embryonic and larval development, while melanin provides modest secondary protection. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Highly Dynamic Gene Family Evolution Suggests Changing Roles for PON Genes Within Metazoa

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage

Ultraviolet radiation (UVR) and its deleterious effects on living cells selects for UVR-protective mechanisms. Organisms across the tree of life evolved a variety of natural sunscreens to prevent UVR-induced cellular damage and stress. However, in vertebrates, only melanin is known to act as a sunscreen. Here we demonstrate that gadusol, a transparent compound discovered over 40 years ago in fish eggs, is a maternally provided sunscreen required for survival of embryonic and larval zebrafish exposed to UVR. Mutating an enzyme involved in gadusol biosynthesis increases the formation of cyclobutane pyrimidine dimers, a hallmark of UVB-induced DNA damage. Compared to the contributions of melanin and the chorion, gadusol is the primary sunscreening mechanism in embryonic and larval fish. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Multiple 9-1-1 complexes promote homolog synapsis, DSB repair, and ATR signaling during mammalian meiosis

DNA damage response mechanisms have meiotic roles that ensure successful gamete formation. While completion of meiotic double-strand break (DSB) repair requires the canonical RAD9A-RAD1-HUS1 (9A-1-1) complex, mammalian meiocytes also express RAD9A and HUS1 paralogs, RAD9B and HUS1B, predicted to form alternative 9-1-1 complexes. The RAD1 subunit is shared by all predicted 9-1-1 complexes and localizes to meiotic chromosomes even in the absence of HUS1 and RAD9A. Here, we report that testis-specific disruption of RAD1 in mice resulted in impaired DSB repair, germ cell depletion, and infertility. Unlike Hus1 or Rad9a disruption, Rad1 loss in meiocytes also caused severe defects in homolog synapsis, impaired phosphorylation of ATR targets such as H2AX, CHK1, and HORMAD2, and compromised meiotic sex chromosome inactivation. Together, these results establish critical roles for both canonical and alternative 9-1-1 complexes in meiotic ATR activation and successful prophase I completion.

Activation by cleavage of the epithelial Na+ channel α and γ subunits independently coevolved with the vertebrate terrestrial migration

Vertebrates evolved mechanisms for sodium conservation and gas exchange in conjunction with migration from aquatic to terrestrial habitats. Epithelial Na⁺ channel (ENaC) function is critical to systems responsible for extracellular fluid homeostasis and gas exchange. ENaC is activated by cleavage at multiple specific extracellular polybasic sites, releasing inhibitory tracts from the channel’s α and γ subunits. We found that proximal and distal polybasic tracts in ENaC subunits coevolved, consistent with the dual cleavage requirement for activation observed in mammals. Polybasic tract pairs evolved with the terrestrial migration and the appearance of lungs, coincident with the ENaC activator aldosterone, and appeared independently in the α and γ subunits. In summary, sites within ENaC for protease activation developed in vertebrates when renal Na⁺ conservation and alveolar gas exchange were required for terrestrial survival.

Experimental exchange of paralogous domains in the MLH family provides evidence of sub-functionalization after gene duplication

Baker's yeast contains a large number of duplicated genes; some function redundantly, whereas others have more specialized roles. We used the MLH family of DNA mismatch repair (MMR) proteins as a model to better understand the steps that lead to gene specialization following a gene duplication event. We focused on two highly conserved yeast MLH proteins, Pms1 and Mlh3, with Pms1 having a major role in the repair of misincorporation events during DNA replication and Mlh3 acting to resolve recombination intermediates in meiosis to form crossovers. The baker's yeast Mlh3 and Pms1 proteins are significantly diverged (19% overall identity), suggesting that an extensive number of evolutionary steps, some major, others involving subtle refinements, took place to diversify the MLH proteins. Using phylogenetic and molecular approaches, we provide evidence that all three domains (N-terminal ATP binding, linker, C-terminal endonuclease/MLH interaction) in the MLH protein family are critical for conferring pathway specificity. Importantly, mlh3 alleles in the ATP binding and endonuclease domains improved MMR functions in strains lacking the Pms1 protein and did not disrupt Mlh3 meiotic functions. This ability for mlh3 alleles to complement the loss of Pms1 suggests that an ancestral Pms1/Mlh3 protein was capable of performing both MMR and crossover functions. Our strategy for analyzing MLH pathway specificity provides an approach to understand how paralogs have evolved to support distinct cellular processes.

Evolutionary, proteomic, and experimental investigations suggest the extracellular matrix of cumulus cells mediates fertilization outcomes

Studies of fertilization biology often focus on sperm and egg interactions. However, before gametes interact, mammalian sperm must pass through the cumulus layer; in mice, this consists of several thousand cells tightly glued together with hyaluronic acid and other proteins. To better understand the role of cumulus cells and their surrounding matrix, we perform proteomic experiments on cumulus oophorus complexes (COCs) in house mice (Mus musculus), producing over 24,000 mass spectra to identify 711 proteins. Seven proteins known to stabilize hyaluronic acid and the extracellular matrix were especially abundant (using spectral counts as an indirect proxy for abundance). Through comparative evolutionary analyses, we show that three of these evolve rapidly, a classic signature of genes that influence fertilization rate. Some of the selected sites overlap regions of the protein known to impact function. In a follow-up experiment, we compared COCs from females raised in two different social environments. Female mice raised in the presence of multiple males produced COCs that were smaller and more resistant to sperm-derived hyaluronidase compared to females raised in the presence of a single male, consistent with a previous study that demonstrated such females produced COCs that were more resistant to fertilization. Although cumulus cells are often thought of as enhancers of fertilization, our evolutionary, proteomic, and experimental investigations implicate their extracellular matrix as a potential mediator of fertilization outcomes.

A Drosophila screen identifies NKCC1 as a modifier of NGLY1 deficiency

N-Glycanase 1 (NGLY1) is a cytoplasmic deglycosylating enzyme. Loss-of-function mutations in the NGLY1 gene cause NGLY1 deficiency, which is characterized by developmental delay, seizures, and a lack of sweat and tears. To model the phenotypic variability observed among patients, we crossed a Drosophila model of NGLY1 deficiency onto a panel of genetically diverse strains. The resulting progeny showed a phenotypic spectrum from 0 to 100% lethality. Association analysis on the lethality phenotype, as well as an evolutionary rate covariation analysis, generated lists of modifying genes, providing insight into NGLY1 function and disease. The top association hit was Ncc69 (human NKCC1/2), a conserved ion transporter. Analyses in NGLY1-/- mouse cells demonstrated that NKCC1 has an altered average molecular weight and reduced function. The misregulation of this ion transporter may explain the observed defects in secretory epithelium function in NGLY1 deficiency patients.

RERconverge: an R package for associating evolutionary rates with convergent traits

MOTIVATION

When different lineages of organisms independently adapt to similar environments, selection often acts repeatedly upon the same genes, leading to signatures of convergent evolutionary rate shifts at these genes. With the increasing availability of genome sequences for organisms displaying a variety of convergent traits, the ability to identify genes with such convergent rate signatures would enable new insights into the molecular basis of these traits.

RESULTS

Here we present the R package RERconverge, which tests for association between relative evolutionary rates of genes and the evolution of traits across a phylogeny. RERconverge can perform associations with binary and continuous traits, and it contains tools for visualization and enrichment analyses of association results.

AVAILABILITY AND IMPLEMENTATION

RERconverge source code, documentation and a detailed usage walk-through are freely available at https://github.com/nclark-lab/RERconverge. Datasets for mammals, Drosophila and yeast are available at https://bit.ly/2J2QBnj.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MCM8IP activates the MCM8-9 helicase to promote DNA synthesis and homologous recombination upon DNA damage

Homologous recombination (HR) mediates the error-free repair of DNA double-strand breaks to maintain genomic stability. Here we characterize C17orf53/MCM8IP, an OB-fold containing protein that binds ssDNA, as a DNA repair factor involved in HR. MCM8IP-deficient cells exhibit HR defects, especially in long-tract gene conversion, occurring downstream of RAD51 loading, consistent with a role for MCM8IP in HR-dependent DNA synthesis. Moreover, loss of MCM8IP confers cellular sensitivity to crosslinking agents and PARP inhibition. Importantly, we report that MCM8IP directly associates with MCM8-9, a helicase complex mutated in primary ovarian insufficiency, and RPA1. We additionally show that the interactions of MCM8IP with MCM8-9 and RPA facilitate HR and promote replication fork progression and cellular viability in response to treatment with crosslinking agents. Mechanistically, MCM8IP stimulates the helicase activity of MCM8-9. Collectively, our work identifies MCM8IP as a key regulator of MCM8-9-dependent DNA synthesis during DNA recombination and replication.

Pan-mammalian analysis of molecular constraints underlying extended lifespan

Although lifespan in mammals varies over 100-fold, the precise evolutionary mechanisms underlying variation in longevity remain unknown. Species-specific genetic changes have been observed in long-lived species including the naked mole-rat, bats, and the bowhead whale, but these adaptations do not generalize to other mammals. We present a novel method to identify associations between rates of protein evolution and continuous phenotypes across the entire mammalian phylogeny. Unlike previous analyses that focused on individual species, we treat absolute and relative longevity as quantitative traits and demonstrate that these lifespan traits affect the evolutionary constraint on hundreds of genes. Specifically, we find that genes related to cell cycle, DNA repair, cell death, the IGF1 pathway, and immunity are under increased evolutionary constraint in large and long-lived mammals. For mammals exceptionally long-lived for their body size, we find increased constraint in inflammation, DNA repair, and NFKB-related pathways. Strikingly, these pathways have considerable overlap with those that have been previously reported to have potentially adaptive changes in single-species studies, and thus would be expected to show decreased constraint in our analysis. This unexpected finding of increased constraint in many longevity-associated pathways underscores the power of our quantitative approach to detect patterns that generalize across the mammalian phylogeny.

Robust Method for Detecting Convergent Shifts in Evolutionary Rates

Identifying genomic elements underlying phenotypic adaptations is an important problem in evolutionary biology. Comparative analyses learning from convergent evolution of traits are gaining momentum in accurately detecting such elements. We previously developed a method for predicting phenotypic associations of genetic elements by contrasting patterns of sequence evolution in species showing a phenotype with those that do not. Using this method, we successfully demonstrated convergent evolutionary rate shifts in genetic elements associated with two phenotypic adaptations, namely the independent subterranean and marine transitions of terrestrial mammalian lineages. Our original method calculates gene-specific rates of evolution on branches of phylogenetic trees using linear regression. These rates represent the extent of sequence divergence on a branch after removing the expected divergence on the branch due to background factors. The rates calculated using this regression analysis exhibit an important statistical limitation, namely heteroscedasticity. We observe that the rates on branches that are longer on average show higher variance, and describe how this problem adversely affects the confidence with which we can make inferences about rate shifts. Using a combination of data transformation and weighted regression, we have developed an updated method that corrects this heteroscedasticity in the rates. We additionally illustrate the improved performance offered by the updated method at robust detection of convergent rate shifts in phylogenetic trees of protein-coding genes across mammals, as well as using simulated tree data sets. Overall, we present an important extension to our evolutionary-rates-based method that performs more robustly and consistently at detecting convergent shifts in evolutionary rates.

Evolutionary rate covariation analysis of E-cadherin identifies Raskol as a regulator of cell adhesion and actin dynamics in Drosophila

The adherens junction couples the actin cytoskeletons of neighboring cells to provide the foundation for multicellular organization. The core of the adherens junction is the cadherin-catenin complex that arose early in the evolution of multicellularity to link actin to intercellular adhesions. Over time, evolutionary pressures have shaped the signaling and mechanical functions of the adherens junction to meet specific developmental and physiological demands. Evolutionary rate covariation (ERC) identifies proteins with correlated fluctuations in evolutionary rate that can reflect shared selective pressures and functions. Here we use ERC to identify proteins with evolutionary histories similar to the Drosophila E-cadherin (DE-cad) ortholog. Core adherens junction components α-catenin and p120-catenin displayed positive ERC correlations with DE-cad, indicating that they evolved under similar selective pressures during evolution between Drosophila species. Further analysis of the DE-cad ERC profile revealed a collection of proteins not previously associated with DE-cad function or cadherin-mediated adhesion. We then analyzed the function of a subset of ERC-identified candidates by RNAi during border cell (BC) migration and identified novel genes that function to regulate DE-cad. Among these, we found that the gene CG42684, which encodes a putative GTPase activating protein (GAP), regulates BC migration and adhesion. We named CG42684 raskol (“to split” in Russian) and show that it regulates DE-cad levels and actin protrusions in BCs. We propose that Raskol functions with DE-cad to restrict Ras/Rho signaling and help guide BC migration. Our results demonstrate that a coordinated selective pressure has shaped the adherens junction and this can be leveraged to identify novel components of the complexes and signaling pathways that regulate cadherin-mediated adhesion.

The establishment of intercellular adhesions facilitated the genesis of multicellular organisms. The adherens junction, which links the actin cytoskeletons of neighboring cells, arose early in the evolution of multicellularity and selective pressures have shaped its function and molecular composition over time. In this study, we used evolutionary rate covariation (ERC) analysis to examine the evolutionary history of the adherens junction and to identify proteins that coevolved with the core adherens junction protein Drosophila E-cadherin (DE-cad). ERC analysis of DE-cad revealed a collection of proteins with similar evolutionary histories. We then tested the role of ERC-identified candidates in border cell migration in the fly egg chamber, a process that requires the coordinated regulation of cell-cell adhesion and cell motility. Among these, we found that a previously uncharacterized gene CG42684, which encodes a putative GTPase activating protein (GAP), regulates the collective cell migration of border cells, stabilizes cell-cell adhesions and regulates the actin dynamics. Our results demonstrate that components of the adherens junction share an evolutionary history and that ERC analysis is a powerful method to identify novel components of cell adhesion complexes in Drosophila.

Ancient convergent losses of Paraoxonase 1 yield potential risks for modern marine mammals

Mammals diversified by colonizing drastically different environments, with each transition yielding numerous molecular changes, including losses of protein function. Though not initially deleterious, these losses could subsequently carry deleterious pleiotropic consequences. We have used phylogenetic methods to identify convergent functional losses across independent marine mammal lineages. In one extreme case, Paraoxonase 1 (PON1) accrued lesions in all marine lineages, while remaining intact in all terrestrial mammals. These lesions coincide with PON1 enzymatic activity loss in marine species' blood plasma. This convergent loss is likely explained by parallel shifts in marine ancestors' lipid metabolism and/or bloodstream oxidative environment affecting PON1's role in fatty acid oxidation. PON1 loss also eliminates marine mammals' main defense against neurotoxicity from specific man-made organophosphorus compounds, implying potential risks in modern environments.

The molecular genetic basis of herbivory between butterflies and their host plants

Interactions between herbivorous insects and their host-plants are a central component of terrestrial food webs and a critical topic in agriculture, where a substantial fraction of potential crop yield is lost annually to pests. Important insights into plant-insect interactions have come from research on specific plant defenses and insect detoxification mechanisms. Yet, much remains unknown about the molecular mechanisms that mediate plant-insect interactions. Here we use multiple genome-wide approaches to map the molecular basis of herbivory from both plant and insect perspectives, focusing on butterflies and their larval host-plants. Parallel genome-wide association studies in the Cabbage White butterfly, Pieris rapae, and its host-plant, Arabidopsis thaliana, pinpointed a small number of butterfly and plant genes that influenced herbivory. These genes, along with much of the genome, were regulated in a dynamic way over the time course of the feeding interaction. Comparative analyses, including diverse butterfly/plant systems, showed a variety of genome-wide responses to herbivory, yet a core set of highly conserved genes in butterflies as well as their host-plants. These results greatly expand our understanding of the genomic causes and evolutionary consequences of ecological interactions across two of nature’s most diverse taxa, butterflies and flowering plants.

Subterranean mammals show convergent regression in ocular genes and enhancers, along with adaptation to tunneling

The underground environment imposes unique demands on life that have led subterranean species to evolve specialized traits, many of which evolved convergently. We studied convergence in evolutionary rate in subterranean mammals in order to associate phenotypic evolution with specific genetic regions. We identified a strong excess of vision- and skin-related genes that changed at accelerated rates in the subterranean environment due to relaxed constraint and adaptive evolution. We also demonstrate that ocular-specific transcriptional enhancers were convergently accelerated, whereas enhancers active outside the eye were not. Furthermore, several uncharacterized genes and regulatory sequences demonstrated convergence and thus constitute novel candidate sequences for congenital ocular disorders. The strong evidence of convergence in these species indicates that evolution in this environment is recurrent and predictable and can be used to gain insights into phenotype-genotype relationships.

Structural complexity and molecular heterogeneity of a butterfly ejaculate reflect a complex history of selection

Male ejaculates are often structurally complex, and this complexity is likely to influence key reproductive interactions between males and females. However, despite its potential evolutionary significance, the molecular underpinnings of ejaculate structural complexity have received little empirical attention. To address this knowledge gap, we sought to understand the biochemical and functional properties of the structurally complex ejaculates of Pieris rapae butterflies. Males in this species produce large ejaculates called spermatophores composed of an outer envelope, an inner matrix, and a bolus of sperm. Females are thought to benefit from the nutrition contained in the soluble inner matrix through increases in longevity and fecundity. However, the indigestible outer envelope of the spermatophore delays female remating, allowing males to monopolize paternity for longer. Here, we show that these two nonsperm-containing spermatophore regions, the inner matrix and the outer envelope, differ in their protein composition and functional properties. We also reveal how these divergent protein mixtures are separately stored in the male reproductive tract and sequentially transferred to the female reproductive tract during spermatophore assembly. Intriguingly, we discovered large quantities of female-derived proteases in both spermatophore regions shortly after mating, which may contribute to spermatophore digestion and hence, female control over remating rate. Finally, we report evidence of past selection on these spermatophore proteins and female proteases, indicating a complex evolutionary history. Our findings illustrate how structural complexity of ejaculates may allow functionally and/or spatially associated suites of proteins to respond rapidly to divergent selective pressures, such as sexual conflict or reproductive cooperation.

Hundreds of Genes Experienced Convergent Shifts in Selective Pressure in Marine Mammals

Mammal species have made the transition to the marine environment several times, and their lineages represent one of the classical examples of convergent evolution in morphological and physiological traits. Nevertheless, the genetic mechanisms of their phenotypic transition are poorly understood, and investigations into convergence at the molecular level have been inconclusive. While past studies have searched for convergent changes at specific amino acid sites, we propose an alternative strategy to identify those genes that experienced convergent changes in their selective pressures, visible as changes in evolutionary rate specifically in the marine lineages. We present evidence of widespread convergence at the gene level by identifying parallel shifts in evolutionary rate during three independent episodes of mammalian adaptation to the marine environment. Hundreds of genes accelerated their evolutionary rates in all three marine mammal lineages during their transition to aquatic life. These marine-accelerated genes are highly enriched for pathways that control recognized functional adaptations in marine mammals, including muscle physiology, lipid-metabolism, sensory systems, and skin and connective tissue. The accelerations resulted from both adaptive evolution as seen in skin and lung genes, and loss of function as in gustatory and olfactory genes. In regard to sensory systems, this finding provides further evidence that reduced senses of taste and smell are ubiquitous in marine mammals. Our analysis demonstrates the feasibility of identifying genes underlying convergent organism-level characteristics on a genome-wide scale and without prior knowledge of adaptations, and provides a powerful approach for investigating the physiological functions of mammalian genes.

Dynamic digestive physiology of a female reproductive organ in a polyandrous butterfly

Reproductive traits experience high levels of selection because of their direct ties to fitness, often resulting in rapid adaptive evolution. Much of the work in this area has focused on male reproductive traits. However, a more comprehensive understanding of female reproductive adaptations and their relationship to male characters is crucial to uncover the relative roles of sexual cooperation and conflict in driving co-evolutionary dynamics between the sexes. We focus on the physiology of a complex female reproductive adaptation in butterflies and moths: a stomach-like organ in the female reproductive tract called the bursa copulatrix that digests the male ejaculate (spermatophore). Little is known about how the bursa digests the spermatophore. We characterized bursa proteolytic capacity in relation to female state in the polyandrous butterfly Pieris rapae. We found that the virgin bursa exhibits extremely high levels of proteolytic activity. Furthermore, in virgin females, bursal proteolytic capacity increases with time since eclosion and ambient temperature, but is not sensitive to the pre-mating social environment. Post copulation, bursal proteolytic activity decreases rapidly before rebounding toward the end of a mating cycle, suggesting active female regulation of proteolysis and/or potential quenching of proteolysis by male ejaculate constituents. Using transcriptomic and proteomic approaches, we report identities for nine proteases actively transcribed by bursal tissue and/or expressed in the bursal lumen that may contribute to observed bursal proteolysis. We discuss how these dynamic physiological characteristics may function as female adaptations resulting from sexual conflict over female remating rate in this polyandrous butterfly.

Candida albicans Transcriptional Profiling Within Biliary Fluid From a Patient With Cholangitis, Before and After Antifungal Treatment and Surgical Drainage

We used ribonucleic acid sequencing to profile Candida albicans transcription within biliary fluid from a patient with cholangitis; samples were collected before and after treatment with fluconazole and drainage. Candida albicans transcriptomes at the infection site distinguished treated from untreated cholangitis. After treatment, 1131 C. albicans genes were differentially expressed in biliary fluid. Up-regulated genes were enriched in hyphal growth, cell wall organization, adhesion, oxidation reduction, biofilm, and fatty acid and ergosterol biosynthesis. This is the first study to define Candida global gene expression during deep-seated human infection. Successful treatment of cholangitis induced C. albicans genes involved in fluconazole responses and pathogenesis.

The Budding Yeast Ubiquitin Protease Ubp7 Is a Novel Component Involved in S Phase Progression

DNA damage must be repaired in an accurate and timely fashion to preserve genome stability. Cellular mechanisms preventing genome instability are crucial to human health because genome instability is considered a hallmark of cancer. Collectively referred to as the DNA damage response, conserved pathways ensure proper DNA damage recognition and repair. The function of numerous DNA damage response components is fine-tuned by posttranslational modifications, including ubiquitination. This not only involves the enzyme cascade responsible for conjugating ubiquitin to substrates but also requires enzymes that mediate directed removal of ubiquitin. Deubiquitinases remove ubiquitin from substrates to prevent degradation or to mediate signaling functions. The Saccharomyces cerevisiae deubiquitinase Ubp7 has been characterized previously as an endocytic factor. However, here we identify Ubp7 as a novel factor affecting S phase progression after hydroxyurea treatment and demonstrate an evolutionary and genetic interaction of Ubp7 with DNA damage repair pathways of homologous recombination and nucleotide excision repair. We find that deletion of UBP7 sensitizes cells to hydroxyurea and cisplatin and demonstrate that factors that stabilize replication forks are critical under these conditions. Furthermore, ubp7Δ cells exhibit an S phase progression defect upon checkpoint activation by hydroxyurea treatment. ubp7Δ mutants are epistatic to factors involved in histone maintenance and modification, and we find that a subset of Ubp7 is chromatin-associated. In summary, our results suggest that Ubp7 contributes to S phase progression by affecting the chromatin state at replication forks, and we propose histone H2B ubiquitination as a potential substrate of Ubp7.

Genetic and phenotypic influences on copulatory plug survival in mice

Across a diversity of animals, male seminal fluid coagulates upon ejaculation to form a hardened structure known as a copulatory plug. Previous studies suggest that copulatory plugs evolved as a mechanism for males to impede remating by females, but detailed investigations into the time course over which plugs survive in the female's reproductive tract are lacking. Here, we cross males from eight inbred strains to females from two inbred strains of house mice (Mus musculus domesticus). Plug survival was significantly affected by male genotype. Against intuition, plug survival time was negatively correlated with plug size: long-lasting plugs were small and relatively more susceptible to proteolysis. Plug size was associated with divergence in major protein composition of seminal vesicle fluid, suggesting that changes in gene expression may play an important role in plug dynamics. In contrast, we found no correlation to genetic variation in the protein-coding regions of five genes thought to be important in copulatory plug formation (Tgm4, Svs1, Svs2, Svs4 and Svs5). Our study demonstrates a complex relationship between copulatory plug characteristics and survival. We discuss several models to explain unexpected variation in plug phenotypes.

ERC analysis: web-based inference of gene function via evolutionary rate covariation

The recent explosion of comparative genomics data presents an unprecedented opportunity to construct gene networks via the evolutionary rate covariation (ERC) signature. ERC is used to identify genes that experienced similar evolutionary histories, and thereby draws functional associations between them. The ERC Analysis website allows researchers to exploit genome-wide datasets to infer novel genes in any biological function and to explore deep evolutionary connections between distinct pathways and complexes. The website provides five analytical methods, graphical output, statistical support and access to an increasing number of taxonomic groups.

AVAILABILITY AND IMPLEMENTATION

Analyses and data at http://csb.pitt.edu/erc_analysis/

CONTACT

nclark@pitt.edu.

Overlapping Patterns of Rapid Evolution in the Nucleic Acid Sensors cGAS and OAS1 Suggest a Common Mechanism of Pathogen Antagonism and Escape

A diverse subset of pattern recognition receptors (PRRs) detects pathogen-associated nucleic acids to initiate crucial innate immune responses in host organisms. Reflecting their importance for host defense, pathogens encode various countermeasures to evade or inhibit these immune effectors. PRRs directly engaged by pathogen inhibitors often evolve under recurrent bouts of positive selection that have been described as molecular ‘arms races.’ Cyclic GMP-AMP synthase (cGAS) was recently identified as a key PRR. Upon binding cytoplasmic double-stranded DNA (dsDNA) from various viruses, cGAS generates the small nucleotide secondary messenger cGAMP to signal activation of innate defenses. Here we report an evolutionary history of cGAS with recurrent positive selection in the primate lineage. Recent studies indicate a high degree of structural similarity between cGAS and 2’-5’-oligoadenylate synthase 1 (OAS1), a PRR that detects double-stranded RNA (dsRNA), despite low sequence identity between the respective genes. We present comprehensive comparative evolutionary analysis of cGAS and OAS1 primate sequences and observe positive selection at nucleic acid binding interfaces and distributed throughout both genes. Our data revealed homologous regions with strong signatures of positive selection, suggesting common mechanisms employed by unknown pathogen encoded inhibitors and similar modes of evasion from antagonism. Our analysis of cGAS diversification also identified alternately spliced forms missing multiple sites under positive selection. Further analysis of selection on the OAS family in primates, which comprises OAS1, OAS2, OAS3 and OASL, suggests a hypothesis where gene duplications and domain fusion events result in paralogs that provide another means of escaping pathogen inhibitors. Together our comparative evolutionary analysis of cGAS and OAS provides new insights into distinct mechanisms by which key molecular sentinels of the innate immune system have adapted to circumvent viral-encoded inhibitors.

A pathogen’s ability to infect new individuals within and across species is largely driven by its capacity to hijack cellular machinery and overcome the immune system. Pathogens have evolved multiple means to evade and shut down host immunity. Typically, mechanisms of inactivation involve direct interactions between host and pathogen factors. To escape inhibition over the course of generations, host factors frequently evolve in a manner that disrupts interactions at specific interfaces with pathogen factors. Likewise, pathogens adapt to restore such interactions, and these genetic tug-of-wars have been described as “molecular-arms races.” Here we focus on the adaptation of two critical host immune factors, cGAS and OAS that share identity in protein structures despite very limited genetic similarity. Our analysis identifies a variety of ways, including amino acid changes on protein surfaces, by which these host factors appear to escape pathogen-mediated inhibition. Surprisingly, some amino acid substitutions are located at equivalent sites suggesting that cGAS and OAS may have adapted to evade common pathogen encoded inhibitors. These data also identify protein surfaces that are targeted by viruses to inhibit host immunity. Taken together our results indicate the existence of critical, yet-to-be identified viral antagonists of cGAS and OAS.

Digestive organ in the female reproductive tract borrows genes from multiple organ systems to adopt critical functions

Persistent adaptive challenges are often met with the evolution of novel physiological traits. Although there are specific examples of single genes providing new physiological functions, studies on the origin of complex organ functions are lacking. One such derived set of complex functions is found in the Lepidopteran bursa copulatrix, an organ within the female reproductive tract that digests nutrients from the male ejaculate or spermatophore. Here, we characterized bursa physiology and the evolutionary mechanisms by which it was equipped with digestive and absorptive functionality. By studying the transcriptome of the bursa and eight other tissues, we revealed a suite of highly expressed and secreted gene products providing the bursa with a combination of stomach-like traits for mechanical and enzymatic digestion of the male spermatophore. By subsequently placing these bursa genes in an evolutionary framework, we found that the vast majority of their novel digestive functions were co-opted by borrowing genes that continue to be expressed in nonreproductive tissues. However, a number of bursa-specific genes have also arisen, some of which represent unique gene families restricted to Lepidoptera and may provide novel bursa-specific functions. This pattern of promiscuous gene borrowing and relatively infrequent evolution of tissue-specific duplicates stands in contrast to studies of the evolution of novelty via single gene co-option. Our results suggest that the evolution of complex organ-level phenotypes may often be enabled (and subsequently constrained) by changes in tissue specificity that allow expression of existing genes in novel contexts, such as reproduction. The extent to which the selective pressures encountered in these novel roles require resolution via duplication and sub/neofunctionalization is likely to be determined by the need for specialized reproductive functionality. Thus, complex physiological phenotypes such as that found in the bursa offer important opportunities for understanding the relative role of pleiotropy and specialization in adaptive evolution.

Evolutionary signatures amongst disease genes permit novel methods for gene prioritization and construction of informative gene-based networks

Genes involved in the same function tend to have similar evolutionary histories, in that their rates of evolution covary over time. This coevolutionary signature, termed Evolutionary Rate Covariation (ERC), is calculated using only gene sequences from a set of closely related species and has demonstrated potential as a computational tool for inferring functional relationships between genes. To further define applications of ERC, we first established that roughly 55% of genetic diseases posses an ERC signature between their contributing genes. At a false discovery rate of 5% we report 40 such diseases including cancers, developmental disorders and mitochondrial diseases. Given these coevolutionary signatures between disease genes, we then assessed ERC's ability to prioritize known disease genes out of a list of unrelated candidates. We found that in the presence of an ERC signature, the true disease gene is effectively prioritized to the top 6% of candidates on average. We then apply this strategy to a melanoma-associated region on chromosome 1 and identify MCL1 as a potential causative gene. Furthermore, to gain global insight into disease mechanisms, we used ERC to predict molecular connections between 310 nominally distinct diseases. The resulting “disease map” network associates several diseases with related pathogenic mechanisms and unveils many novel relationships between clinically distinct diseases, such as between Hirschsprung's disease and melanoma. Taken together, these results demonstrate the utility of molecular evolution as a gene discovery platform and show that evolutionary signatures can be used to build informative gene-based networks.

Molecular evolution has informed our understanding of gene function; however, classical methods have largely been static in their implementation, focusing on single genes. Here, we present and prove the utility of a dynamic, network-based understanding of molecular evolution to infer relationships between genes associated with human diseases. We have shown previously that groups of genes within functional niches tend to share similar evolutionary histories. Exploiting the availability of whole genomes from multiple species, these histories can be numerically scored and dynamically compared to one another using a sequence-based signature termed Evolutionary Rate Covariation (ERC). To explore potential applications, we characterized ERC amongst disease genes and found that many diseases contain significant ERC signatures between their contributing genes. We show that ERC can also prioritize “true” disease genes amongst unrelated gene candidates. Lastly, these signatures can serve as a foundation for creating instructive gene-based networks, unveiling novel relationships between diseases thought to be clinically distinct. Our hope is that this study will add to the increasing evidence that advancing our understanding of molecular evolution can be a crucial asset in large-scale gene discovery pursuits (Link to our webserver that provides intuitive ERC analysis tools: http://csb.pitt.edu/erc_analysis/).

Evolutionary rate covariation reveals shared functionality and coexpression of genes

Evolutionary rate covariation (ERC) is a phylogenetic signature that reflects the covariation of a pair of proteins over evolutionary time. ERC is typically elevated between interacting proteins and so is a promising signature to characterize molecular and functional interactions across the genome. ERC is often assumed to result from compensatory changes at interaction interfaces (i.e., intermolecular coevolution); however, its origin is still unclear and is likely to be complex. Here, we determine the biological factors responsible for ERC in a proteome-wide data set of 4459 proteins in 18 budding yeast species. We show that direct physical interaction is not required to produce ERC, because we observe strong correlations between noninteracting but cofunctional enzymes. We also demonstrate that ERC is uniformly distributed along the protein primary sequence, suggesting that intermolecular coevolution is not generally responsible for ERC between physically interacting proteins. Using multivariate analysis, we show that a pair of proteins is likely to exhibit ERC if they share a biological function or if their expression levels coevolve between species. Thus, ERC indicates shared function and coexpression of protein pairs and not necessarily coevolution between sites, as has been assumed in previous studies. This full interpretation of ERC now provides us with a powerful tool to assign uncharacterized proteins to functional groups and to determine the interconnectedness between entire genetic pathways.

An assessment of lead controls for torch cutting and rivet removal on steel structures

The use of engineering and work practice controls to protect workers from lead-containing dusts and fumes generated during rehabilitation of steel structures is mandated by the Occupational Safety and Health Administration (OSHA) Lead in Construction Standard (1993). Because the implementation and assessment of controls can be problematic in the rugged and dynamic construction environment, industrial hygienists should understand the effectiveness and limitations of controls adopted. The present investigation assesses the efficacy of two controls to reduce lead exposure: paint removal prior to oxy-acetylene torch cutting of steel, and encapsulation of rivets prior to their removal. A task-based exposure assessment approach was used to evaluate these tasks at three sites. Exposures at one site without controls were compared to exposures at sites with controls. Comparison of the results via an analysis of variance (0.05 significance level) indicates that, for torch cutting, exposures at the control site were not significantly different from those at an uncontrolled site (p = 0.14). The results for rivet busting show no significant differences in exposures at the control site compared to the uncontrolled site (p = 0.08). Results are also presented from two control sites where work was done in enclosed spaces. Two main difficulties in applying the controls are explored: technical and managerial. Technical problems during torch cutting included the penetration of paint into the steel profile and the configuration of the structures. For rivet busting, working within an enclosure was an important factor. Management problems arose both from a lack of coordination among different contractors, and from a failure to provide day-to-day guidance and assessment of the control. Important components of a program to implement controls are preplanning and coordination of control implementation, frequent testing of control efficacy, and a method for timely intervention to correct deficiencies.

Less common orbital fracture patterns: the role of computed tomography in the management of depression of the inferior oblique origin and lateral rectus involvement in blow-in fractures

During the past decade, advances in radiographic imaging have made it possible for the surgeon managing orbital fractures to adopt a rational therapeutic strategy based on a knowledge of alterations in surgical anatomy secondary to traumatic injury. To illustrate the value of computed tomography in the surgeon's armamentarium for management of orbital fractures, cases are presented in which imaging proved decisive in planning a course of therapy. Two patients presented with two types of isolated lateral blow-in fracture, an uncommon fracture pattern. The other cases underscore the value of defining involvement of the inferior oblique origin and lateral rectus muscles in imaging complex orbital fractures, issues not emphasized in earlier literature. Although diplopia alone does not always warrant surgical intervention, diplopia in the context of computed tomography-defined muscle entrapment or muscle origin displacement justifies operative therapy. These cases demonstrate the value of computed tomography in directing surgical therapy with resolution of diplopia and prevention and correction of enophthalmos.

Patterns of sealant delivery under variable third party requirements

Despite evidence indicating dental sealants to be effective in preventing caries on the pit and fissure surfaces of teeth, only 8 percent of school-aged US children had received sealants as of 1986-87. While many rationales have been suggested and scientifically answered for this low level of utilization, issues of cost and cost effectiveness remain considerable barriers to many insurers and public programs. This study reports dentist behavior when sealants were added with few restrictions as a new benefit to an existing third party system of care. Results are compared to another third party program with stricter reimbursement policies. Overall, sealants were provided in a reasonable fashion to both groups of patients, relative to teeth selected for sealing and costs. Only a small proportion of patients receiving sealants were under age six or above age 19. Even under the most liberal program, sealants were predominantly targeted to teeth at highest risk to decay. Still, for both groups, a significant proportion of the sealants were placed in teeth at points in time quite distant from expected eruption patterns and, thus, their periods of highest risk to decay. Overall, dentists appeared to use sealants only minimally in their practices. These data suggest that sealants can be added to third party dental programs with little overall risk of inappropriate use or abuse. They also suggest that specific efforts are merited to educate providers as to the most effective times at which to provide sealants for preventing the maximum amount of decay in a population.

Managing organizational stress in nursing

The strategic management of a nursing organization in a changing health care market creates pressure on nursing administrators and on their organizations. The effective management of organizational stress potentiates quality nursing care and reduces the costs related to stress and burnout. This article describes a two-tiered program designed to manage organizational stress. The first component focuses on educating nursing leadership about the nature and impact of organizational stress and how to minimize it. The second component is the instrument used to measure the perceived effectiveness of our interventions.

The human resource consultant in nursing

The authors have developed a role that supports nurse managers in developing their management skills and effectiveness. Nursing administrators will find the role useful because it provides a different avenue for the strategic planning of change. The role and examples of its use in the support of nursing are described. Considerations in operationalizing this role in a nursing organization also are discussed.