Positive selection of yeast nonhomologous end-joining genes and a retrotransposon conflict hypothesis

Genetic conflicts in budding yeast: The 2μ plasmid as a model selfish element

Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2μ) plasmid as a model selfish element. The 2μ is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2μ plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

Recurrent Duplication and Diversification of a Vital DNA Repair Gene Family Across Drosophila

Maintaining genome integrity is vital for organismal survival and reproduction. Essential, broadly conserved DNA repair pathways actively preserve genome integrity. However, many DNA repair proteins evolve adaptively. Ecological forces like UV exposure are classically cited drivers of DNA repair evolution. Intrinsic forces like repetitive DNA, which also imperil genome integrity, have received less attention. We recently reported that a Drosophila melanogaster-specific DNA satellite array triggered species-specific, adaptive evolution of a DNA repair protein called Spartan/MH. The Spartan family of proteases cleave hazardous, covalent crosslinks that form between DNA and proteins ("DNA-protein crosslink repair"). Appreciating that DNA satellites are both ubiquitous and universally fast-evolving, we hypothesized that satellite DNA turnover spurs adaptive evolution of DNA-protein crosslink repair beyond a single gene and beyond the D. melanogaster lineage. This hypothesis predicts pervasive Spartan gene family diversification across Drosophila species. To study the evolutionary history of the Drosophila Spartan gene family, we conducted population genetic, molecular evolution, phylogenomic, and tissue-specific expression analyses. We uncovered widespread signals of positive selection across multiple Spartan family genes and across multiple evolutionary timescales. We also detected recurrent Spartan family gene duplication, divergence, and gene loss. Finally, we found that ovary-enriched parent genes consistently birthed functionally diverged, testis-enriched daughter genes. To account for Spartan family diversification, we introduce a novel mechanistic model of antagonistic coevolution that links DNA satellite evolution and adaptive regulation of Spartan protease activity. This framework promises to accelerate our understanding of how DNA repeats drive recurrent evolutionary innovation to preserve genome integrity.

Diverse and mobile: eccDNA-based identification of carrot low-copy-number LTR retrotransposons active in callus cultures

Long terminal repeat retrotransposons (LTR-RTs) are mobilized via an RNA intermediate using a 'copy and paste' mechanism, and account for the majority of repetitive DNA in plant genomes. As a side effect of mobilization, the formation of LTR-RT-derived extrachromosomal circular DNAs (eccDNAs) occurs. Thus, high-throughput sequencing of eccDNA can be used to identify active LTR-RTs in plant genomes. Despite the release of a reference genome assembly, carrot LTR-RTs have not yet been thoroughly characterized. LTR-RTs are abundant and diverse in the carrot genome. We identified 5976 carrot LTR-RTs, 2053 and 1660 of which were attributed to Copia and Gypsy superfamilies, respectively. They were further classified into lineages, families and subfamilies. More diverse LTR-RT lineages, i.e. lineages comprising many low-copy-number subfamilies, were more frequently associated with genic regions. Certain LTR-RT lineages have been recently active in Daucus carota. In particular, low-copy-number LTR-RT subfamilies, e.g. those belonging to the DcAle lineage, have significantly contributed to carrot genome diversity as a result of continuing activity. We utilized eccDNA sequencing to identify and characterize two DcAle subfamilies, Alex1 and Alex3, active in carrot callus. We documented 14 and 32 de novo insertions of Alex1 and Alex3, respectively, which were positioned in non-repetitive regions.

DNA damage repair system in C57BL/6 J mice is evolutionarily stable

Background

DNA damage repair (DDR) system is vital in maintaining genome stability and survival. DDR consists of over 160 genes in 7 different pathways to repair specific type of DNA damage caused by external and internal damaging factors. The functional importance of DDR system implies that evolution could play important roles in maintaining its functional intactness to perform its function. Indeed, it has been observed that positive selection is present in BRCA1 and BRCA2 (BRCA), which are key genes in homologous recombination pathway of DDR system, in the humans and its close relatives of chimpanzee and bonobos. Efforts have been made to investigate whether the same selection could exist for BRCA in other mammals but found no evidence so far. However, as most of the studies in non-human mammals analyzed only a single or few individuals in the studied species, the observation may not reflect the true status in the given species. Furthermore, few studies have studied evolution selection in other DDR genes except BRCA. In current study, we used laboratory mouse C57BL/6 J as a model to address evolution selection on DDR genes in non-primate mammals by dynamically monitoring genetic variation across 30 generations in C57BL/6 J.

Results

Using exome sequencing, we collected coding sequences of 169 DDR genes from 44 C57BL/6 J individual genomes in 2018. We compared the coding sequences with the mouse reference genome sequences derived from 1998 C57BL/6 J DNA, and with the mouse Eve6B reference genome sequences derived from 2003 C57BL/6 J DNA, covering 30 generations of C57BL/6 J from 1998 to 2018. We didn’t identify meaningful coding variation in either Brca1 or Brca2, or in 167 other DDR genes across the 30 generations. In the meantime, we did identify 812 coding variants in 116 non-DNA damage repair genes during the same period, which served as a quality control to validate the reliability of our analytic pipeline and the negative results in DDR genes.

Conclusions

DDR genes in laboratory mouse strain C57BL/6 J were not under positive selection across its 30-generation period, highlighting the possibility that DDR system in rodents could be evolutionarily stable.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07983-7.

Nuku, a family of primate retrocopies derived from KU70

The gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70 have the hallmark of long interspersed element-1 mediated retrotransposition with evidence of target-site duplications, the poly-A tails, and the absence of introns. Evolutionary analysis of this expanded family of KU70-derived "NUKU" retrocopies reveals that these genes are both ancient and also actively being created in extant primate species. NUKU retrocopies show evidence of functional divergence away from KU70, as evinced by their altered pattern of tissue expression and possible tissue-specific translation. Molecular modeling predicted that amino acid changes in Nuku2p at the interaction interface with Ku80p would prevent the assembly of the Ku heterodimer. The lack of Nuku2p-Ku80p interaction was confirmed by yeast two-hybrid assay, which contrasts the robust interaction of Ku70p-Ku80p. While several NUKU retrocopies appear to have been degraded by mutation, NUKU2 shows evidence of positive natural selection, suggesting that this retrocopy is undergoing neofunctionalization. Although Nuku proteins do not appear to antagonize retrovirus transduction in cell culture, the observed expansion and rapid evolution of NUKUs could be being driven by alternative selective pressures related to infectious disease or an undefined role in primate physiology.

Adaptive evolution of an essential telomere protein restricts telomeric retrotransposons

Essential, conserved cellular processes depend not only on essential, strictly conserved proteins but also on essential proteins that evolve rapidly. To probe this poorly understood paradox, we exploited the rapidly evolving Drosophila telomere-binding protein, cav/HOAP, which protects chromosomes from lethal end-to-end fusions. We replaced the D. melanogaster HOAP with a highly diverged version from its close relative, D. yakuba. The D. yakuba HOAP ('HOAP[yak]') localizes to D. melanogaster telomeres and protects D. melanogaster chromosomes from fusions. However, HOAP[yak] fails to rescue a previously uncharacterized HOAP function: silencing of the specialized telomeric retrotransposons that, instead of telomerase, maintain chromosome length in Drosophila. Whole genome sequencing and cytogenetics of experimentally evolved populations revealed that HOAP[yak] triggers telomeric retrotransposon proliferation, resulting in aberrantly long telomeres. This evolution-generated, separation-of-function allele resolves the paradoxical observation that a fast-evolving essential gene directs an essential, strictly conserved function: telomeric retrotransposon containment, not end-protection, requires evolutionary innovation at HOAP.

Overexpressed coiled-coil domain containing protein 8 (CCDC8) mediates newly synthesized HIV-1 Gag lysosomal degradation

Normally, HIV-1 enters into CD4+ cells through membrane fusion, and newly synthesized HIV-1 viral proteins assemble on the plasma membrane to form viral particles and bud out. In the previous study, we found host factor coiled-coil domain containing protein 8 (CCDC8) can strongly inhibit HIV-1 production, but the underline mechanism is not clear. Here we show that overexpression of CCDC8 reverses the normal HIV-1 production process, and causes newly assembled HIV-1 Gag particles to be endocytosed on the plasma membrane, rather than budding out. Live-cell imaging system captured the moment of CCDC8-mediated Gag internalization on the plasma membrane, and the speed of Gag turnover is up to 1.53 μm/s, much faster than Gag assembly on the plasma membrane. After Gag internalization, it accumulates in the cellular organelle-lysosome for degradation, but not proteasome, autophagosome, endoplasmic reticulum, clathrin or recycling endosome. In addition, CCDC8 is a membrane-associated protein, and N-terminal of CCDC8 is very important for membrane binding, and also important for inhibition of Gag assembly. C-terminal deletion of CCDC8 has a little effect on anti-HIV-1 effect. Moreover, CCDC8 is phosphorylated at amino acid threonine T87 and serine S261, and mono-methylated at lysine K491. Alanine mutations of T87A, S261A and K491A singly or in combination do not affect CCDC8 anti-HIV activity. In conclusion, overexpression of CCDC8 can cause newly assembled HIV-1 Gag particles on the plasma membrane to be endocytosed and degraded in lysosome.

Yeast XRS2 and human NBN gene: Experimental evidence for homology using codon optimized cDNA

The genes, XRS2 in Saccharomyces cerevisiae and NBN in mammals, have little sequence identity at the amino acid level. Nevertheless, they are both found together with MRE11 and RAD50 in a highly conserved protein complex which functions in the repair of DNA double-strand breaks. Here, we have examined the evolutionary and functional relationship of these two genes by cross-complementation experiments. These experiments necessitated sequence correction for specific codon usage before they could be successfully conducted. We present evidence that despite extreme sequence divergence nibrin can, at least partially, replace Xrs2 in the cellular DNA damage response, and Xrs2 is able to promote nuclear localization of MRE11 in NBS cells. We discuss that the extreme sequence divergence reflects a unique adaptive pressure during evolution related to the specific eukaryotic role for both Xrs2 and nibrin in the subcellular localisation of the DNA repair complex. This, we suggest, is of particular relevance when cells are infected by viruses. The conflict hypothesis of co-evolution of DNA repair genes and DNA viruses may thus explain the very low sequence identity of these two homologous genes.

Control of yeast retrotransposons mediated through nucleoporin evolution

Yeasts serve as hosts to several types of genetic parasites. Few studies have addressed the evolutionary trajectory of yeast genes that control the stable co-existence of these parasites with their host cell. In Saccharomyces yeasts, the retrovirus-like Ty retrotransposons must access the nucleus. We show that several genes encoding components of the yeast nuclear pore complex have experienced natural selection for substitutions that change the encoded protein sequence. By replacing these S. cerevisiae genes with orthologs from other Saccharomyces species, we discovered that natural sequence changes have affected the mobility of Ty retrotransposons. Specifically, changing the genetic sequence of NUP84 or NUP82 to match that of other Saccharomyces species alters the mobility of S. cerevisiae Ty1 and Ty3. Importantly, all tested housekeeping functions of NUP84 and NUP82 remained equivalent across species. Signatures of natural selection, resulting in altered interactions with viruses and parasitic genetic elements, are common in host defense proteins. Yet, few instances have been documented in essential housekeeping proteins. The nuclear pore complex is the gatekeeper of the nucleus. This study shows how the evolution of this large, ubiquitous eukaryotic complex can alter the replication of a molecular parasite, but concurrently maintain essential host functionalities regarding nucleocytoplasmic trafficking.

The frenemies within: viruses, retrotransposons and plasmids that naturally infect Saccharomyces yeasts

Viruses are a major focus of current research efforts because of their detrimental impact on humanity and their ubiquity within the environment. Bacteriophages have long been used to study host-virus interactions within microbes, but it is often forgotten that the single-celled eukaryote Saccharomyces cerevisiae and related species are infected with double-stranded RNA viruses, single-stranded RNA viruses, LTR-retrotransposons and double-stranded DNA plasmids. These intracellular nucleic acid elements have some similarities to higher eukaryotic viruses, i.e. yeast retrotransposons have an analogous lifecycle to retroviruses, the particle structure of yeast totiviruses resembles the capsid of reoviruses and segregation of yeast plasmids is analogous to segregation strategies used by viral episomes. The powerful experimental tools available to study the genetics, cell biology and evolution of S. cerevisiae are well suited to further our understanding of how cellular processes are hijacked by eukaryotic viruses, retrotransposons and plasmids. This article has been written to briefly introduce viruses, retrotransposons and plasmids that infect Saccharomyces yeasts, emphasize some important cellular proteins and machineries with which they interact, and suggest the evolutionary consequences of these interactions. Copyright © 2017 John Wiley & Sons, Ltd.

Recurrent Innovation at Genes Required for Telomere Integrity in Drosophila

Telomeres are nucleoprotein complexes at the ends of linear chromosomes. These specialized structures ensure genome integrity and faithful chromosome inheritance. Recurrent addition of repetitive, telomere-specific DNA elements to chromosome ends combats end-attrition, while specialized telomere-associated proteins protect naked, double-stranded chromosome ends from promiscuous repair into end-to-end fusions. Although telomere length homeostasis and end-protection are ubiquitous across eukaryotes, there is sporadic but building evidence that the molecular machinery supporting these essential processes evolves rapidly. Nevertheless, no global analysis of the evolutionary forces that shape these fast-evolving proteins has been performed on any eukaryote. The abundant population and comparative genomic resources of Drosophila melanogaster and its close relatives offer us a unique opportunity to fill this gap. Here we leverage population genetics, molecular evolution, and phylogenomics to define the scope and evolutionary mechanisms driving fast evolution of genes required for telomere integrity. We uncover evidence of pervasive positive selection across multiple evolutionary timescales. We also document prolific expansion, turnover, and expression evolution in gene families founded by telomeric proteins. Motivated by the mutant phenotypes and molecular roles of these fast-evolving genes, we put forward four alternative, but not mutually exclusive, models of intra-genomic conflict that may play out at very termini of eukaryotic chromosomes. Our findings set the stage for investigating both the genetic causes and functional consequences of telomere protein evolution in Drosophila and beyond.

Nomadic genetic elements contribute to oncogenic translocations: Implications in carcinogenesis

Chromosomal translocations as molecular signatures have been reported in various malignancies but, the mechanism behind which is largely unknown. Swapping of chromosomal fragments occurs by induction of double strand breaks (DSBs), most of which were initially assumed de novo. However, decoding of human genome proved that transposable elements (TE) might have profound influence on genome integrity. TEs are highly conserved mobile genetic elements that generate DSBs, subsequently resulting in large chromosomal rearrangements. Previously TE insertions were thought to be harmless, but recently gains attention due to the origin of spectrum of post-insertional genomic alterations and subsequent transcriptional alterations leading to development of deleterious effects mainly carcinogenesis. Though the existing knowledge on the cancer-associated TE dynamics is very primitive, exploration of underlying mechanism promises better therapeutic strategies for cancer. Thus, this review focuses on the prevalence of TE in the genome, associated genomic instability upon transposition activation and impact on tumorigenesis.

The Ty1 Retrotransposon Restriction Factor p22 Targets Gag

A novel form of copy number control (CNC) helps maintain a low number of Ty1 retrovirus-like transposons in the Saccharomyces genome. Ty1 produces an alternative transcript that encodes p22, a trans-dominant negative inhibitor of Ty1 retrotransposition whose sequence is identical to the C-terminal half of Gag. The level of p22 increases with copy number and inhibits normal Ty1 virus-like particle (VLP) assembly and maturation through interactions with full length Gag. A forward genetic screen for CNC-resistant (CNC^R) mutations in Ty1 identified missense mutations in GAG that restore retrotransposition in the presence of p22. Some of these mutations map within a predicted UBN2 domain found throughout the Ty1/copia family of long terminal repeat retrotransposons, and others cluster within a central region of Gag that is referred to as the CNC^R domain. We generated multiple alignments of yeast Ty1-like Gag proteins and found that some Gag proteins, including those of the related Ty2 elements, contain non-Ty1 residues at multiple CNC^R sites. Interestingly, the Ty2-917 element is resistant to p22 and does not undergo a Ty1-like form of CNC. Substitutions conferring CNC^R map within predicted helices in Ty1 Gag that overlap with conserved sequence in Ty1/copia, suggesting that p22 disturbs a central function of the capsid during VLP assembly. When hydrophobic residues within predicted helices in Gag are mutated, Gag level remains unaffected in most cases yet VLP assembly and maturation is abnormal. Gag CNC^R mutations do not alter binding to p22 as determined by co-immunoprecipitation analyses, but instead, exclude p22 from Ty1 VLPs. These findings suggest that the CNC^R alleles enhance retrotransposition in the presence of p22 by allowing productive Gag-Gag interactions during VLP assembly. Our work also expands the strategies used by retroviruses for developing resistance to Gag-like restriction factors to now include retrotransposons.

The presence of transposable elements in the eukaryotic genome threatens genomic stability and normal gene function, thus various defense mechanisms exist to silence element expression and target integration to benign locations in the genome. Even though the budding yeast Saccharomyces lacks many of the defense systems present in other eukaryotes, including RNAi, DNA methylation, and APOBEC3 proteins, they maintain low numbers of mobile elements in their genome. In the case of the Saccharomyces retrotransposon Ty1, a system called copy number control (CNC) helps determine the number of elements in the genome. Recently, we demonstrated that the mechanism of CNC relies on a trans-acting protein inhibitor of Ty1 expressed from the element itself. This protein inhibitor, called p22, impacts the replication of Ty1 as its copy number increases. To identify a molecular target of p22, mutagenized Ty1 was subjected to a forward genetic screen for CNC-resistance. Mutations in specific domains of Gag, including the UBN2 Gag motif and a novel region we have named the CNC^R domain, confer CNC^R by preventing the incorporation of p22 into assembling virus-like particles (VLPs), which restores maturation and completion of the Ty1 life cycle. The mechanism of Ty1 inhibition by p22 is conceptually similar to Gag-like restriction factors in mammals since they inhibit normal particle function. In particular, resistance to p22 and the enJS56A1 restriction factor of sheep involves exclusion of the restriction factor during particle assembly, although Ty1 CNC^R achieves this in a way that is distinct from the Jaagsiekte retrovirus escape mutants. Our work introduces an intriguing variation on resistance mechanisms to retroviral restriction factors.

The Zinc-Finger Antiviral Protein ZAP Inhibits LINE and Alu Retrotransposition

Long INterspersed Element-1 (LINE-1 or L1) is the only active autonomous retrotransposon in the human genome. To investigate the interplay between the L1 retrotransposition machinery and the host cell, we used co-immunoprecipitation in conjunction with liquid chromatography and tandem mass spectrometry to identify cellular proteins that interact with the L1 first open reading frame-encoded protein, ORF1p. We identified 39 ORF1p-interacting candidate proteins including the zinc-finger antiviral protein (ZAP or ZC3HAV1). Here we show that the interaction between ZAP and ORF1p requires RNA and that ZAP overexpression in HeLa cells inhibits the retrotransposition of engineered human L1 and Alu elements, an engineered mouse L1, and an engineered zebrafish LINE-2 element. Consistently, siRNA-mediated depletion of endogenous ZAP in HeLa cells led to a ~2-fold increase in human L1 retrotransposition. Fluorescence microscopy in cultured human cells demonstrated that ZAP co-localizes with L1 RNA, ORF1p, and stress granule associated proteins in cytoplasmic foci. Finally, molecular genetic and biochemical analyses indicate that ZAP reduces the accumulation of full-length L1 RNA and the L1-encoded proteins, yielding mechanistic insight about how ZAP may inhibit L1 retrotransposition. Together, these data suggest that ZAP inhibits the retrotransposition of LINE and Alu elements.

Long INterspersed Element-1 (LINE-1 or L1) is the only active autonomous retrotransposon in the human genome. L1s comprise ~17% of human DNA and it is estimated that an average human genome has ~80–100 active L1s. L1 moves throughout the genome via a “copy-and-paste” mechanism known as retrotransposition. L1 retrotransposition is known to cause mutations; thus, it stands to reason that the host cell has evolved mechanisms to protect the cell from unabated retrotransposition. Here, we demonstrate that the zinc-finger antiviral protein (ZAP) inhibits the retrotransposition of human L1 and Alu retrotransposons, as well as related retrotransposons from mice and zebrafish. Biochemical and genetic data suggest that ZAP interacts with L1 RNA. Fluorescent microscopy demonstrates that ZAP associates with L1 in cytoplasmic foci that co-localize with stress granule proteins. Mechanistic analyses suggest that ZAP reduces the expression of full-length L1 RNA and the L1-encoded proteins, thereby providing mechanistic insight for how ZAP may restricts retrotransposition. Importantly, these data suggest that ZAP initially may have evolved to combat endogenous retrotransposons and subsequently was co-opted as a viral restriction factor.

The influence of selection for protein stability on dN/dS estimations

Understanding the relative contributions of various evolutionary processes-purifying selection, neutral drift, and adaptation-is fundamental to evolutionary biology. A common metric to distinguish these processes is the ratio of nonsynonymous to synonymous substitutions (i.e., dN/dS) interpreted from the neutral theory as a null model. However, from biophysical considerations, mutations have non-negligible effects on the biophysical properties of proteins such as folding stability. In this work, we investigated how stability affects the rate of protein evolution in phylogenetic trees by using simulations that combine explicit protein sequences with associated stability changes. We first simulated myoglobin evolution in phylogenetic trees with a biophysically realistic approach that accounts for 3D structural information and estimates of changes in stability upon mutation. We then compared evolutionary rates inferred directly from simulation to those estimated using maximum-likelihood (ML) methods. We found that the dN/dS estimated by ML methods (ωML) is highly predictive of the per gene dN/dS inferred from the simulated phylogenetic trees. This agreement is strong in the regime of high stability where protein evolution is neutral. At low folding stabilities and under mutation-selection balance, we observe deviations from neutrality (per gene dN/dS > 1 and dN/dS < 1). We showed that although per gene dN/dS is robust to these deviations, ML tests for positive selection detect statistically significant per site dN/dS > 1. Altogether, we show how protein biophysics affects the dN/dS estimations and its subsequent interpretation. These results are important for improving the current approaches for detecting positive selection.

Rapid evolution of BRCA1 and BRCA2 in humans and other primates

Background

The maintenance of chromosomal integrity is an essential task of every living organism and cellular repair mechanisms exist to guard against insults to DNA. Given the importance of this process, it is expected that DNA repair proteins would be evolutionarily conserved, exhibiting very minimal sequence change over time. However, BRCA1, an essential gene involved in DNA repair, has been reported to be evolving rapidly despite the fact that many protein-altering mutations within this gene convey a significantly elevated risk for breast and ovarian cancers.

Results

To obtain a deeper understanding of the evolutionary trajectory of BRCA1, we analyzed complete BRCA1 gene sequences from 23 primate species. We show that specific amino acid sites have experienced repeated selection for amino acid replacement over primate evolution. This selection has been focused specifically on humans and our closest living relatives, chimpanzees (Pan troglodytes) and bonobos (Pan paniscus). After examining BRCA1 polymorphisms in 7 bonobo, 44 chimpanzee, and 44 rhesus macaque (Macaca mulatta) individuals, we find considerable variation within each of these species and evidence for recent selection in chimpanzee populations. Finally, we also sequenced and analyzed BRCA2 from 24 primate species and find that this gene has also evolved under positive selection.

Conclusions

While mutations leading to truncated forms of BRCA1 are clearly linked to cancer phenotypes in humans, there is also an underlying selective pressure in favor of amino acid-altering substitutions in this gene. A hypothesis where viruses are the drivers of this natural selection is discussed.

Genome rearrangements and pervasive meiotic drive cause hybrid infertility in fission yeast

Hybrid sterility is one of the earliest postzygotic isolating mechanisms to evolve between two recently diverged species. Here we identify causes underlying hybrid infertility of two recently diverged fission yeast species Schizosaccharomyces pombe and S. kambucha, which mate to form viable hybrid diploids that efficiently complete meiosis, but generate few viable gametes. We find that chromosomal rearrangements and related recombination defects are major but not sole causes of hybrid infertility. At least three distinct meiotic drive alleles, one on each S. kambucha chromosome, independently contribute to hybrid infertility by causing nonrandom spore death. Two of these driving loci are linked by a chromosomal translocation and thus constitute a novel type of paired meiotic drive complex. Our study reveals how quickly multiple barriers to fertility can arise. In addition, it provides further support for models in which genetic conflicts, such as those caused by meiotic drive alleles, can drive speciation.

DOI:http://dx.doi.org/10.7554/eLife.02630.001

It is widely thought that all of the billions of species on Earth are descended from a common ancestor. New species are created via a process called speciation, and nature employs various ‘barriers’ to keep closely related species distinct from one another. One of these barriers is called hybrid sterility. Horses and donkeys, for example, can mate to produce hybrids called mules, but mules cannot produce offspring of their own because they are infertile.

Hybrid sterility can occur for a number of reasons. Mules are infertile because they inherit 32 chromosomes from their horse parent, but only 31 chromosomes from their donkey parent—and so have an odd chromosome that they cannot pair-off when they make sperm or egg cells. However, even if a hybrid inherits the same number of chromosomes from each parent, if the chromosomes from the two parents have different structures, the hybrid may still be infertile.

Zanders et al. have now looked at two species of fission yeast—S. pombe and S. kambucha—that share 99.5% of their DNA sequence. Although hybrids of these two species inherit three chromosomes from each parent, the majority of spores (the yeast equivalent of sperm) that these hybrids produce fail to develop into new yeast cells. Zanders et al. identified two causes of this infertility: one of these was chromosomal rearrangement; the other was due to three different sites in the DNA of S. kambucha that interfere with the development of the spores that inherit S. pombe chromosomes.

Since these two yeast species are so closely related, the findings of Zanders et al. reveal how quickly multiple barriers to fertility can arise. In addition, these findings provide further support for models in which conflicts between different genes in genomes can drive the process of speciation.

DOI:http://dx.doi.org/10.7554/eLife.02630.002

Genome-wide screens for sensitivity to ionizing radiation identify the fission yeast nonhomologous end joining factor Xrc4

Nonhomologous end joining (NHEJ) is the main means for repairing DNA double-strand breaks (DSBs) in human cells. Molecular understanding of NHEJ has benefited from analyses in the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. In human cells, the DNA ligation reaction of the classical NHEJ pathway is carried out by a protein complex composed of DNA ligase IV (LigIV) and XRCC4. In S. cerevisiae, this reaction is catalyzed by a homologous complex composed of Dnl4 and Lif1. Intriguingly, no homolog of XRCC4 has been found in S. pombe, raising the possibility that such a factor may not always be required for classical NHEJ. Here, through screening the ionizing radiation (IR) sensitivity phenotype of a genome-wide fission yeast deletion collection in both the vegetative growth state and the spore state, we identify Xrc4, a highly divergent homolog of human XRCC4. Like other fission yeast NHEJ factors, Xrc4 is critically important for IR resistance of spores, in which no homologous recombination templates are available. Using both extrachromosomal and chromosomal DSB repair assays, we show that Xrc4 is essential for classical NHEJ. Exogenously expressed Xrc4 colocalizes with the LigIV homolog Lig4 at the chromatin region of the nucleus in a mutually dependent manner. Furthermore, like their human counterparts, Xrc4 and Lig4 interact with each other and this interaction requires the inter-BRCT linker and the second BRCT domain of Lig4. Our discovery of Xrc4 suggests that an XRCC4 family protein is universally required for classical NHEJ in eukaryotes.

Conceptual and empirical challenges of ascribing functions to transposable elements

Media attention and the subsequent scientific backlash engendered by the claim by spokespeople for the Encyclopedia of DNA Elements (ENCODE) project that 80% of the human genome has a biochemical function highlight the need for a clearer understanding of function concepts in biology. This article provides an overview of two major function concepts that have been developed in the philosophy of science--the causal role concept and the selected effects concept--and their relevance to ENCODE. Unlike in some previous critiques, the ENCODE project is not considered problematic here because it employed a causal role definition of function (which is relatively common in genetics) but because of how this concept was misused. In addition, several unique challenges that arise when dealing with transposable elements (TEs) but that were ignored by ENCODE are highlighted. These include issues surrounding TE-level versus organism-level selection, the origins versus the persistence of elements, and accidental versus functional organism-level benefits. Finally, some key questions are presented that should be addressed in any study aiming to ascribe functions to major portions of large eukaryotic genomes, the majorities of which are made up of transposable elements.

Catalytic and noncatalytic roles of the CtIP endonuclease in double-strand break end resection

The carboxy-terminal binding protein (CtBP)-interacting protein (CtIP) is known to function in 5' strand resection during homologous recombination, similar to the budding yeast Sae2 protein, but its role in this process is unclear. Here, we characterize recombinant human CtIP and find that it exhibits 5' flap endonuclease activity on branched DNA structures, independent of the MRN complex. Phosphorylation of CtIP at known damage-dependent sites and other sites is essential for its catalytic activity, although the S327 and T847 phosphorylation sites are dispensable. A catalytic mutant of CtIP that is deficient in endonuclease activity exhibits wild-type levels of homologous recombination at restriction enzyme-generated breaks but is deficient in processing topoisomerase adducts and radiation-induced breaks in human cells, suggesting that the nuclease activity of CtIP is specifically required for the removal of DNA adducts at sites of DNA breaks.

Positive selection of primate genes that promote HIV-1 replication

Evolutionary analyses have revealed that most host-encoded restriction factors against HIV have experienced virus-driven selection during primate evolution. However, HIV also depends on the function of many human proteins, called host factors, for its replication. It is not clear whether virus-driven selection shapes the evolution of host factor genes to the extent that it is known to shape restriction factor genes. We show that five out of 40 HIV host factor genes (13%) analyzed do bear strong signatures of positive selection. Some of these genes (CD4, NUP153, RANBP2/NUP358) have been characterized with respect to the HIV lifecycle, while others (ANKRD30A/NY-BR-1 and MAP4) remain relatively uncharacterized. One of these, ANKRD30A, shows the most rapid evolution within this set of genes and is induced by interferon stimulation. We discuss how evolutionary analysis can aid the study of host factors for viral replication, just as it has the study of host immunity systems.

Dual host-virus arms races shape an essential housekeeping protein

Relentless selective pressures exerted by viruses trigger arms race dynamics that shape the evolution of even critical host genes like those involved in iron homeostasis.

Transferrin Receptor (TfR1) is the cell-surface receptor that regulates iron uptake into cells, a process that is fundamental to life. However, TfR1 also facilitates the cellular entry of multiple mammalian viruses. We use evolutionary and functional analyses of TfR1 in the rodent clade, where two families of viruses bind this receptor, to mechanistically dissect how essential housekeeping genes like TFR1 successfully balance the opposing selective pressures exerted by host and virus. We find that while the sequence of rodent TfR1 is generally conserved, a small set of TfR1 residue positions has evolved rapidly over the speciation of rodents. Remarkably, all of these residues correspond to the two virus binding surfaces of TfR1. We show that naturally occurring mutations at these positions block virus entry while simultaneously preserving iron-uptake functionalities, both in rodent and human TfR1. Thus, by constantly replacing the amino acids encoded at just a few residue positions, TFR1 divorces adaptation to ever-changing viruses from preservation of key cellular functions. These dynamics have driven genetic divergence at the TFR1 locus that now enforces species-specific barriers to virus transmission, limiting both the cross-species and zoonotic transmission of these viruses.

Genetic differences between mammalian species dictate the patterns of viral infection observed in nature. They also define how viruses must evolve in order to infect new mammalian hosts, giving rise to new and sometimes pandemic diseases. Because viruses must enter cells before they can replicate, new diseases often emerge when existing viruses evolve the ability to bind to the cell-surface receptor of a new species. At the same time, host cell receptors also evolve to counteract virus attacks. This back-and-forth evolution between virus and host can lead to an arms race that shapes the sequences of the proteins involved. In wild rodent populations, the retrovirus MMTV and New World arenaviruses both exploit Transferrin Receptor 1 (TfR1) to enter the cells of their hosts. Here we show that the physical interactions between these viruses and TfR1 have triggered evolutionary arms race dynamics that have directly modified the sequence of TfR1 and at least one of the viruses involved. Computational evolutionary analysis allowed us to identify specific residues in TfR1 that define patterns of viral infection in nature. The approach presented here can theoretically be applied to the study of any virus, through analysis of host genes known to be key to controlling viral infection. As such, this approach can expand our understanding of how viruses emerge from wildlife reservoirs, and how they drive the evolution of host genes.

Population-Sequencing as a Biomarker for Sample Characterization

Sequencing is accepted as the “gold” standard for genetic analysis and continues to be used as a validation and reference tool. The idea of using sequence analysis directly for sample characterization has been met with skepticism. However, herein, utility of direct use of sequencing to identify multiple genomes present in samples is presented and reviewed. All samples and “pure” isolates are populations of genomes. Population-Sequencing is the use of probabilistic matching tools in combination with large volumes of sequence information to identify genomes present, based on DNA analysis across entire genomes to determine genome assignments, to calculate confidence scores of major and minor genome content. Accurate genome identification from mixtures without culture purification steps can achieve phylogenetic classification by direct analysis of millions of DNA fragments. Genome sequencing data of mixtures can function as biomarkers for use to interrogate genetic content of samples and to establish a sample profile, inclusive of major and minor genome components, drill down to identify rare SNP and mutation events, compare relatedness of genetic content between samples, profile-to-profile, and provide a probabilistic or statistical scoring confidence for sample characterization and attribution. The application of Population-Sequencing will facilitate sample characterization and genome identification strategies.

Human immunodeficiency virus type 1 Vpu and cellular TASK proteins suppress transcription of unintegrated HIV-1 DNA

BACKGROUND

Unintegrated HIV-1 DNA serves as transcriptionally active templates in HIV-infected cells. Several host factors including NF-κβ enhance HIV-1 transcription. HIV-1 induced NF-κβ activation can be suppressed by viral protein U (Vpu). Interestingly HIV-1 Vpu shares amino acid homology with cellular Twik-related Acid Sensitive K+ (TASK) channel 1 and the proteins physically interact in cultured cells and AIDS lymphoid tissue. Furthermore, the first transmembrane domain of TASK-1 is functionally interchangeable with Vpu and like Vpu enhances HIV-1 release.

RESULTS

Here we further characterize the role of TASK channels and Vpu in HIV-1 replication. We demonstrate that both TASK channels and Vpu can preferentially inhibit transcription of unintegrated HIV-1 DNA. Interestingly, TASK-1 ion channel function is not required and suppression of HIV-1 transcription by TASK-1 and Vpu was reversed by overexpression of RelA (NF-κβ p65).

CONCLUSION

TASK proteins and Vpu suppress transcription of unintegrated HIV-1 DNA through an NF-κβ-dependent mechanism. Taken together these findings support a possible physiological role for HIV-1 Vpu and TASK proteins as modulators of transcription of unintegrated HIV-1 DNA genomes.

Divergence of the yeast transcription factor FZF1 affects sulfite resistance

Changes in gene expression are commonly observed during evolution. However, the phenotypic consequences of expression divergence are frequently unknown and difficult to measure. Transcriptional regulators provide a mechanism by which phenotypic divergence can occur through multiple, coordinated changes in gene expression during development or in response to environmental changes. Yet, some changes in transcriptional regulators may be constrained by their pleiotropic effects on gene expression. Here, we use a genome-wide screen for promoters that are likely to have diverged in function and identify a yeast transcription factor, FZF1, that has evolved substantial differences in its ability to confer resistance to sulfites. Chimeric alleles from four Saccharomyces species show that divergence in FZF1 activity is due to changes in both its coding and upstream noncoding sequence. Between the two closest species, noncoding changes affect the expression of FZF1, whereas coding changes affect the expression of SSU1, a sulfite efflux pump activated by FZF1. Both coding and noncoding changes also affect the expression of many other genes. Our results show how divergence in the coding and promoter region of a transcription factor alters the response to an environmental stress.

Changes in gene regulation are thought to play an important role in evolution. While variation in gene expression between species is common, it is hard to identify the phenotypic consequences of this variation since many changes in gene expression may have subtle or no phenotypic effects. In this study, we investigate changes in sulfite resistance and gene expression caused by the transcription factor, FZF1, that has evolved rapidly during the divergence of related yeast species. We find that divergence in the ability of FZF1 to confer sulfite resistance is mediated by changes in its expression as well as changes in its protein structure, both of which cause changes in the expression of other genes. Our results show how the combination of multiple changes within a transcription factor can produce substantial changes in phenotype and the expression of many genes.

The retrohoming of linear group II intron RNAs in Drosophila melanogaster occurs by both DNA ligase 4-dependent and -independent mechanisms

Mobile group II introns are bacterial retrotransposons that are thought to have invaded early eukaryotes and evolved into introns and retroelements in higher organisms. In bacteria, group II introns typically retrohome via full reverse splicing of an excised intron lariat RNA into a DNA site, where it is reverse transcribed by the intron-encoded protein. Recently, we showed that linear group II intron RNAs, which can result from hydrolytic splicing or debranching of lariat RNAs, can retrohome in eukaryotes by performing only the first step of reverse splicing, ligating their 3' end to the downstream DNA exon. Reverse transcription then yields an intron cDNA, whose free end is linked to the upstream DNA exon by an error-prone process that yields junctions similar to those formed by non-homologous end joining (NHEJ). Here, by using Drosophila melanogaster NHEJ mutants, we show that linear intron RNA retrohoming occurs by major Lig4-dependent and minor Lig4-independent mechanisms, which appear to be related to classical and alternate NHEJ, respectively. The DNA repair polymerase θ plays a crucial role in both pathways. Surprisingly, however, mutations in Ku70, which functions in capping chromosome ends during NHEJ, have only moderate, possibly indirect effects, suggesting that both Lig4 and the alternate end-joining ligase act in some retrohoming events independently of Ku. Another potential Lig4-independent mechanism, reverse transcriptase template switching from the intron RNA to the upstream exon DNA, occurs in vitro, but gives junctions differing from the majority in vivo. Our results show that group II introns can utilize cellular NHEJ enzymes for retromobility in higher organisms, possibly exploiting mechanisms that contribute to retrotransposition and mitigate DNA damage by resident retrotransposons. Additionally, our results reveal novel activities of group II intron reverse transcriptases, with implications for retrohoming mechanisms and potential biotechnological applications.

Birth, death, and replacement of karyopherins in Drosophila

Nucleocytoplasmic transport is a broadly conserved process across eukaryotes. Despite its essential function and conserved mechanism, components of the nuclear transport apparatus have been implicated in genetic conflicts in Drosophila, especially in the male germ line. The best understood case is represented by a truncated RanGAP gene duplication that is part of the segregation distorter system in Drosophila melanogaster. Consistent with the hypothesis that the nuclear transport pathway is at the heart of mediating genetic conflicts, both nucleoporins and directionality imposing components of nuclear transport have previously been shown to evolve under positive selection. Here, we present a comprehensive phylogenomic analysis of importins (karyopherins) in Drosophila evolution. Importins are adaptor molecules that physically mediate the transport of cargo molecules and comprise the third component of the nuclear transport apparatus. We find that importins have been repeatedly gained and lost throughout various stages of Drosophila evolution, including two intriguing examples of an apparently coincident loss and gain of nonorthologous and noncanonical importin-α. Although there are a few signatures of episodic positive selection, genetic innovation in importin evolution is more evident in patterns of recurrent gene birth and loss specifically for function in Drosophila testes, which is consistent with their role in supporting host genomes defense against segregation distortion.

The role of unintegrated DNA in HIV infection

Integration of the reverse transcribed viral genome into host chromatin is the hallmark of retroviral replication. Yet, during natural HIV infection, various unintegrated viral DNA forms exist in abundance. Though linear viral cDNA is the precursor to an integrated provirus, increasing evidence suggests that transcription and translation of unintegrated DNAs prior to integration may aid productive infection through the expression of early viral genes. Additionally, unintegrated DNA has the capacity to result in preintegration latency, or to be rescued and yield productive infection and so unintegrated DNA, in some circumstances, may be considered to be a viral reservoir. Recently, there has been interest in further defining the role and function of unintegrated viral DNAs, in part because the use of anti-HIV integrase inhibitors leads to an abundance of unintegrated DNA, but also because of the potential use of non-integrating lentiviral vectors in gene therapy and vaccines. There is now increased understanding that unintegrated viral DNA can either arise from, or be degraded through, interactions with host DNA repair enzymes that may represent a form of host antiviral defence. This review focuses on the role of unintegrated DNA in HIV infection and additionally considers the potential implications for antiviral therapy.

Evidence for divergent evolution of growth temperature preference in sympatric Saccharomyces species

The genus Saccharomyces currently includes eight species in addition to the model yeast Saccharomyces cerevisiae, most of which can be consistently isolated from tree bark and soil. We recently found sympatric pairs of Saccharomyces species, composed of one cryotolerant and one thermotolerant species in oak bark samples of various geographic origins. In order to contribute to explain the occurrence in sympatry of Saccharomyces species, we screened Saccharomyces genomic data for protein divergence that might be correlated to distinct growth temperature preferences of the species, using the dN/dS ratio as a measure of protein evolution rates and pair-wise species comparisons. In addition to proteins previously implicated in growth at suboptimal temperatures, we found that glycolytic enzymes were among the proteins exhibiting higher than expected divergence when one cryotolerant and one thermotolerant species are compared. By measuring glycolytic fluxes and glycolytic enzymatic activities in different species and at different temperatures, we subsequently show that the unusual divergence of glycolytic genes may be related to divergent evolution of the glycolytic pathway aligning its performance to the growth temperature profiles of the different species. In general, our results support the view that growth temperature preference is a trait that may have undergone divergent selection in the course of ecological speciation in Saccharomyces.

Genome-wide approaches to the study of adaptive gene expression evolution: systematic studies of evolutionary adaptations involving gene expression will allow many fundamental questions in evolutionary biology to be addressed

The role of gene expression in evolutionary adaptation has been a subject of debate for over 40 years. cis-regulation of transcription has been proposed to be the primary source of morphological novelty in evolution, though this is based on only a handful of examples. Recently the first genome-wide studies of gene expression adaptation have been published, giving us an initial global view of this process. Systematic studies such as these will allow a number of key questions currently facing the field of gene expression evolution to be addressed.

Genome-wide survey of natural selection on functional, structural, and network properties of polymorphic sites in Saccharomyces paradoxus

BACKGROUND

To characterize the genetic basis of phenotypic evolution, numerous studies have identified individual genes that have likely evolved under natural selection. However, phenotypic changes may represent the cumulative effect of similar evolutionary forces acting on functionally related groups of genes. Phylogenetic analyses of divergent yeast species have identified functional groups of genes that have evolved at significantly different rates, suggestive of differential selection on the functional properties. However, due to environmental heterogeneity over long evolutionary timescales, selection operating within a single lineage may be dramatically different, and it is not detectable via interspecific comparisons alone. Moreover, interspecific studies typically quantify selection on protein-coding regions using the D(n)/D(s) ratio, which cannot be extended easily to study selection on noncoding regions or synonymous sites. The population genetic-based analysis of selection operating within a single lineage ameliorates these limitations.

FINDINGS

We investigated selection on several properties associated with genes, promoters, or polymorphic sites, by analyzing the derived allele frequency spectrum of single nucleotide polymorphisms (SNPs) in 28 strains of Saccharomyces paradoxus. We found evidence for significant differential selection between many functionally relevant categories of SNPs, underscoring the utility of function-centric approaches for discovering signatures of natural selection. When comparable, our findings are largely consistent with previous studies based on interspecific comparisons, with one notable exception: our study finds that mutations from an ancient amino acid to a relatively new amino acid are selectively disfavored, whereas interspecific comparisons have found selection against ancient amino acids. Several of our findings have not been addressed through prior interspecific studies: we find that synonymous mutations from preferred to unpreferred codons are selected against and that synonymous SNPs in the linker regions of proteins are relatively less constrained than those within protein domains.

CONCLUSIONS

We present the first global survey of selection acting on various functional properties in S. paradoxus. We found that selection pressures previously detected over long evolutionary timescales have also shaped the evolution of S. paradoxus. Importantly, we also make novel discoveries untenable via conventional interspecific analyses.

Genomes in conflict: maintaining genome integrity during virus infection

The cellular surveillance network for sensing and repairing damaged DNA prevents an array of human diseases, and when compromised it can lead to genomic instability and cancer. The carefully maintained cellular response to DNA damage is challenged during viral infection, when foreign DNA is introduced into the cell. The battle between virus and host generates a genomic conflict. The host attempts to limit viral infection and protect its genome, while the virus deploys tactics to eliminate, evade, or exploit aspects of the cellular defense. Studying this conflict has revealed that the cellular DNA damage response machinery comprises part of the intrinsic cellular defense against viral infection. In this review we examine recent advances in this emerging field. We identify common themes used by viruses in their attempts to commandeer or circumvent the host cell's DNA repair machinery, and highlight potential outcomes of the conflict for both virus and host.

Ancient and recent adaptive evolution of primate non-homologous end joining genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.

Because all cells experience DNA damage, they must also have mechanisms for repairing DNA. When the proteins that repair DNA malfunction, mutation and disease often result. Based on their fundamental importance, DNA repair proteins would be expected to be well preserved over evolutionary time in order to ensure optimal DNA repair function. However, a previous genome-wide study of molecular evolution in Saccharomyces yeast identified the non-homologous end joining (NHEJ) DNA repair pathway as one of the two most rapidly evolving pathways in the yeast genome. In order to analyze the evolution of this pathway in humans, we have generated large evolutionary sequence sets of NHEJ genes from our primate relatives. Similar to the scenario in yeast, several genes in this pathway are evolving rapidly in primate genomes and in modern human populations. Thus, complex and seemingly opposite selective forces are shaping the evolution of these important DNA repair genes. The finding that NHEJ genes are rapidly evolving in species groups as diverse as yeasts and primates indicates a systematic perturbation of the NHEJ pathway, one that is potentially important to human health.

Protein evolution in yeast transcription factor subnetworks

When averaged over the full yeast protein-protein interaction and transcriptional regulatory networks, protein hubs with many interaction partners or regulators tend to evolve significantly more slowly due to increased negative selection. However, genome-wide analysis of protein evolution in the subnetworks of associations involving yeast transcription factors (TFs) reveals that TF hubs do not tend to evolve significantly more slowly than TF non-hubs. This result holds for all four major types of TF hubs: interaction hubs, regulatory in-degree and out-degree hubs, as well as co-regulatory hubs that jointly regulate target genes with many TFs. Furthermore, TF regulatory in-degree hubs tend to evolve significantly more quickly than TF non-hubs. Most importantly, the correlations between evolutionary rate (K(A)/K(S)) and degrees for TFs are significantly more positive than those for generic proteins within the same global protein-protein interaction and transcriptional regulatory networks. Compared to generic protein hubs, TF hubs operate at a higher level in the hierarchical structure of cellular networks, and hence experience additional evolutionary forces (relaxed negative selection or positive selection through network rewiring). The striking difference between the evolution of TF hubs and generic protein hubs demonstrates that components within the same global network can be governed by distinct organizational and evolutionary principles.

A novel method to detect proteins evolving at correlated rates: identifying new functional relationships between coevolving proteins

Interacting proteins evolve at correlated rates, possibly as the result of evolutionary pressures shared by functional groups and/or coevolution between interacting proteins. This evolutionary signature can be exploited to learn more about protein networks and to infer functional relationships between proteins on a genome-wide scale. Multiple methods have been introduced that detect correlated evolution using amino acid distances. One assumption made by these methods is that the neutral rate of nucleotide substitution is uniform over time; however, this is unlikely and such rate heterogeneity would adversely affect amino acid distance methods. We explored alternative methods that detect correlated rates using protein-coding nucleotide sequences in order to better estimate the rate of nonsynonymous substitution at each branch (d(N)) normalized by the underlying synonymous substitution rate (d(S)). Our novel likelihood method, which was robust to realistic simulation parameters, was tested on Drosophila nuclear pore proteins, which form a complex with well-documented physical interactions. The method revealed significantly correlated evolution between nuclear pore proteins, where members of a stable subcomplex showed stronger correlations compared with those proteins that interact transiently. Furthermore, our likelihood approach was better able to detect correlated evolution among closely related species than previous methods. Hence, these sequence-based methods are a complementary approach for detecting correlated evolution and could be applied genome-wide to provide candidate protein-protein interactions and functional group assignments using just coding sequences.

Segregating YKU80 and TLC1 alleles underlying natural variation in telomere properties in wild yeast

In yeast, as in humans, telomere length varies among individuals and is controlled by multiple loci. In a quest to define the extent of variation in telomere length, we screened 112 wild-type Saccharomyces sensu stricto isolates. We found extensive telomere length variation in S. paradoxus isolates. This phenotype correlated with their geographic origin: European strains were observed to have extremely short telomeres (<150 bp), whereas American isolates had telomeres approximately three times as long (>400 bp). Insertions of a URA3 gene near telomeres allowed accurate analysis of individual telomere lengths and telomere position effect (TPE). Crossing the American and European strains resulted in F1 spores with a continuum of telomere lengths consistent with what would be predicted if many quantitative trait loci (QTLs) were involved in length maintenance. Variation in TPE is similarly quantitative but only weakly correlated with telomere length. Genotyping F1 segregants indicated several QTLs associated with telomere length and silencing variation. These QTLs include likely candidate genes but also map to regions where there are no known genes involved in telomeric properties. We detected transgressive segregation for both phenotypes. We validated by reciprocal hemizygosity that YKU80 and TLC1 are telomere-length QTLs in the two S. paradoxus subpopulations. Furthermore, we propose that sequence divergence within the Ku heterodimer generates negative epistasis within one of the allelic combinations (American-YKU70 and European-YKU80) resulting in very short telomeres.

Telomere length is a complex trait that varies among individuals. Its regulation is critical to the process of aging, and altered length control can result in either senescence or immortalization. We detected extreme variation between different subpopulations of the wild yeast S. paradoxus, the closest relative to S. cerevisiae. By tagging individual telomeric ends in these two groups, we show that regardless of the total number of telomeric repeats, the critical length at which any telomere is replenished remains conserved. To detect the quantitative trait loci (QTLs) behind the length variation, we used the two sub-populations with the most polar distribution to generate progeny and perform linkage analysis. Further, we validated that naturally occurring sequence variations in YKU80 and TLC1, two genes previously shown to be important for telomere length maintenance, can explain part of the variation. We also identified other loci that influence both telomere length and gene silencing. Further investigation will provide more insights into the underlying genetic mechanism behind normal telomere regulation, potentially relevant in aging and aging-related disease such as cancer.

Repair-mediated duplication by capture of proximal chromosomal DNA has shaped vertebrate genome evolution

DNA double-strand breaks (DSBs) are a common form of cellular damage that can lead to cell death if not repaired promptly. Experimental systems have shown that DSB repair in eukaryotic cells is often imperfect and may result in the insertion of extra chromosomal DNA or the duplication of existing DNA at the breakpoint. These events are thought to be a source of genomic instability and human diseases, but it is unclear whether they have contributed significantly to genome evolution. Here we developed an innovative computational pipeline that takes advantage of the repetitive structure of genomes to detect repair-mediated duplication events (RDs) that occurred in the germline and created insertions of at least 50 bp of genomic DNA. Using this pipeline we identified over 1,000 probable RDs in the human genome. Of these, 824 were intra-chromosomal, closely linked duplications of up to 619 bp bearing the hallmarks of the synthesis-dependent strand-annealing repair pathway. This mechanism has duplicated hundreds of sequences predicted to be functional in the human genome, including exons, UTRs, intron splice sites and transcription factor binding sites. Dating of the duplication events using comparative genomics and experimental validation revealed that the mechanism has operated continuously but with decreasing intensity throughout primate evolution. The mechanism has produced species-specific duplications in all primate species surveyed and is contributing to genomic variation among humans. Finally, we show that RDs have also occurred, albeit at a lower frequency, in non-primate mammals and other vertebrates, indicating that this mechanism has been an important force shaping vertebrate genome evolution.

The repair of DNA double-strand breaks (DSBs) is essential for the maintenance of genome integrity. The mechanisms by which DSBs are repaired have been the subject of intense experimental investigations. It has emerged that several imperfect repair pathways exist in eukaryotes that have the potential to result in chromosomal alterations, including genomic duplications. However, it remains unclear to what extent these imperfect repair events have contributed to shaping genomes throughout evolution. Here we introduce an innovative computational approach that takes advantage of the repetitive nature of eukaryotic genomes to identify repair-mediated duplications (RD) that occurred during evolution. We discovered over one thousand RDs in the human genome, with two-thirds resulting from the capture of a chromosomal DNA segment located in close proximity to the presumed site of the DSB, giving rise to local genomic duplications. Comparative genomic analyses reveal that the mechanism has operated continuously, but with decreasing intensity during primate evolution, generating species-specific duplications in all primates surveyed and generating genomic variation among humans. Finally, we show that RDs have also occurred in non-primate mammals and other vertebrates, indicating that this is a previously under-appreciated force shaping vertebrate genomes.

Genetic evidence that the non-homologous end-joining repair pathway is involved in LINE retrotransposition

Long interspersed elements (LINEs) are transposable elements that proliferate within eukaryotic genomes, having a large impact on eukaryotic genome evolution. LINEs mobilize via a process called retrotransposition. Although the role of the LINE-encoded protein(s) in retrotransposition has been extensively investigated, the participation of host-encoded factors in retrotransposition remains unclear. To address this issue, we examined retrotransposition frequencies of two structurally different LINEs—zebrafish ZfL2-2 and human L1—in knockout chicken DT40 cell lines deficient in genes involved in the non-homologous end-joining (NHEJ) repair of DNA and in human HeLa cells treated with a drug that inhibits NHEJ. Deficiencies of NHEJ proteins decreased retrotransposition frequencies of both LINEs in these cells, suggesting that NHEJ is involved in LINE retrotransposition. More precise characterization of ZfL2-2 insertions in DT40 cells permitted us to consider the possibility of dual roles for NHEJ in LINE retrotransposition, namely to ensure efficient integration of LINEs and to restrict their full-length formation.

Long interspersed elements (LINEs) are transposable elements that mobilize and amplify their own copies within eukaryotic genomes. Although LINEs had been considered as “junk” DNA, recent studies have suggested that the LINE-induced alterations of host chromosomes are a major driving force for eukaryotic genome evolution. LINEs mobilize via a mechanism called retrotransposition, in which transcribed LINE RNA is reverse transcribed into DNA that is then integrated into the host chromosome. Although the role of LINE-encoded proteins in retrotransposition has been revealed, the participation of host-encoded proteins has not been well investigated. Here, using knockout chicken DT40 cell lines, we present genetic evidence that the host-encoded proteins involved in repair of DNA double-strand breaks participate in LINE retrotransposition. More precise characterization of LINE insertions in DT40 cells suggested dual roles for these host DNA repair proteins in LINE retrotransposition; one function is required for efficient integration of LINEs and the other restricts their full-length formation.

Protein kinase R reveals an evolutionary model for defeating viral mimicry

Distinguishing self from non-self is a fundamental biological challenge. Many pathogens exploit the challenge of self discrimination by employing mimicry to subvert key cellular processes including the cell cycle, apoptosis and cytoskeletal dynamics. Other mimics interfere with immunity. Poxviruses encode K3L, a mimic of eIF2alpha, which is the substrate of protein kinase R (PKR), an important component of innate immunity in vertebrates. The PKR-K3L interaction exemplifies the conundrum imposed by viral mimicry. To be effective, PKR must recognize a conserved substrate (eIF2alpha) while avoiding rapidly evolving substrate mimics such as K3L. Using the PKR-K3L system and a combination of phylogenetic and functional analyses, we uncover evolutionary strategies by which host proteins can overcome mimicry. We find that PKR has evolved under intense episodes of positive selection in primates. The ability of PKR to evade viral mimics is partly due to positive selection at sites most intimately involved in eIF2alpha recognition. We also find that adaptive changes on multiple surfaces of PKR produce combinations of substitutions that increase the odds of defeating mimicry. Thus, although it can seem that pathogens gain insurmountable advantages by mimicking cellular components, host factors such as PKR can compete in molecular 'arms races' with mimics because of evolutionary flexibility at protein interaction interfaces challenged by mimicry.

Network hubs buffer environmental variation in Saccharomyces cerevisiae

Regulatory and developmental systems produce phenotypes that are robust to environmental and genetic variation. A gene product that normally contributes to this robustness is termed a phenotypic capacitor. When a phenotypic capacitor fails, for example when challenged by a harsh environment or mutation, the system becomes less robust and thus produces greater phenotypic variation. A functional phenotypic capacitor provides a mechanism by which hidden polymorphism can accumulate, whereas its failure provides a mechanism by which evolutionary change might be promoted. The primary example to date of a phenotypic capacitor is Hsp90, a molecular chaperone that targets a large set of signal transduction proteins. In both Drosophila and Arabidopsis, compromised Hsp90 function results in pleiotropic phenotypic effects dependent on the underlying genotype. For some traits, Hsp90 also appears to buffer stochastic variation, yet the relationship between environmental and genetic buffering remains an important unresolved question. We previously used simulations of knockout mutations in transcriptional networks to predict that many gene products would act as phenotypic capacitors. To test this prediction, we use high-throughput morphological phenotyping of individual yeast cells from single-gene deletion strains to identify gene products that buffer environmental variation in Saccharomyces cerevisiae. We find more than 300 gene products that, when absent, increase morphological variation. Overrepresented among these capacitors are gene products that control chromosome organization and DNA integrity, RNA elongation, protein modification, cell cycle, and response to stimuli such as stress. Capacitors have a high number of synthetic-lethal interactions but knockouts of these genes do not tend to cause severe decreases in growth rate. Each capacitor can be classified based on whether or not it is encoded by a gene with a paralog in the genome. Capacitors with a duplicate are highly connected in the protein–protein interaction network and show considerable divergence in expression from their paralogs. In contrast, capacitors encoded by singleton genes are part of highly interconnected protein clusters whose other members also tend to affect phenotypic variability or fitness. These results suggest that buffering and release of variation is a widespread phenomenon that is caused by incomplete functional redundancy at multiple levels in the genetic architecture.

Most species maintain abundant genetic variation and experience a wide range of environmental conditions, yet phenotypic differences between individuals are usually small. This phenomenon, known as phenotypic robustness, presents an apparent contradiction: if biological systems are so resistant to variation, how do they diverge and adapt through evolutionary time? Here, we address this question by investigating the molecular mechanisms that underlie phenotypic robustness and how these mechanisms can be broken to produce phenotypic heterogeneity. We identify genes that contribute to phenotypic robustness in yeast by analyzing the variance of morphological phenotypes in a comprehensive collection of single-gene knockout strains. We find that ∼5% of yeast genes break phenotypic robustness when knocked out. The products of these genes tend to be involved in critical cellular processes, including maintaining DNA stability, processing RNA, modifying proteins, and responding to stressful environments. These genes tend to interact genetically with a large number of other genes, and their products tend to interact physically with a large number of other gene products. Our results suggest that loss of phenotypic robustness might be a common phenomenon during evolution that occurs when cellular networks are disrupted.

A genome-wide screen inSaccharomyces cerevisiae identifies over 300 gene products that buffer environmental variation--dubbed phenotypic capacitors--and function as hubs in protein-protein and synthetic-lethal interaction networks.

Different strategies to persist: the pogo-like Lemi1 transposon produces miniature inverted-repeat transposable elements or typical defective elements in different plant genomes

Miniature inverted-repeat transposable elements (MITEs) are a particular type of defective class II elements present in genomes as high-copy-number populations of small and highly homogeneous elements. While virtually all class II transposon families contain non-autonomous defective transposon copies, only a subset of them have a related MITE family. At present it is not known in which circumstances MITEs are generated instead of typical class II defective transposons. The ability to produce MITEs could be an exclusive characteristic of particular transposases, could be related to a particular structure of certain defective class II elements, or could be the consequence of particular constraints imposed by certain host genomes on transposon populations. We describe here a new family of pogo-like transposons from Medicago truncatula closely related to the Arabidopsis Lemi1 element that we have named MtLemi1. In contrast to the Arabidopsis Lemi1, present as a single-copy element and associated with hundreds of related Emigrant MITEs, MtLemi1 has attained >30 copies and has not generated MITEs. This shows that a particular transposon can adopt completely different strategies to colonize genomes. The comparison of AtLemi1 and MtLemi1 reveals transposase-specific domains and possible regulatory sequences that could be linked to the ability to produce MITEs.

HIV-1 accessory proteins--ensuring viral survival in a hostile environment

One of the features of primate immunodeficiency viruses (HIVs and SIVs) that distinguishes them from other retroviruses is the array of "accessory" proteins they encode. Here, we discuss recent advances in understanding the interactions of the HIV-1 Nef, Vif, Vpu, and Vpr proteins with factors and pathways expressed in cells of the immune system. In at least three instances, the principal activity of the accessory proteins appears to be evasion from various forms of cell-mediated (or intrinsic), antiviral resistance. Broadly speaking, the HIV-1 accessory proteins modify the local environment within infected cells to ensure viral persistence, replication, dissemination, and transmission.

Sae2 is an endonuclease that processes hairpin DNA cooperatively with the Mre11/Rad50/Xrs2 complex

Mre11/Rad50 complexes in all organisms function in the repair of DNA double-strand breaks. In budding yeast, genetic evidence suggests that the Sae2 protein is essential for the processing of hairpin DNA intermediates and meiotic double-strand breaks by Mre11/Rad50 complexes, but the biochemical basis of this functional relationship is not known. Here we demonstrate that recombinant Sae2 binds DNA and exhibits endonuclease activity on single-stranded DNA independently of Mre11/Rad50 complexes, but hairpin DNA structures are cleaved cooperatively in the presence of Mre11/Rad50 or Mre11/Rad50/Xrs2. Hairpin structures are not processed at the tip by Sae2 but rather at single-stranded DNA regions adjacent to the hairpin. Truncation and missense mutants of Sae2 inactivate this endonuclease activity in vitro and fail to complement Deltasae2 strains in vivo for meiosis and recombination involving hairpin intermediates, suggesting that the catalytic activities of Sae2 are important for its biological functions.

Restricting the ligation step of non-homologous end-joining

Non-homologous end-joining is an important pathway for repairing DNA double-strand breaks. The budding yeast Saccharomyces cerevisiae possesses two proteins, Nej1/Lif2 and Ntr1/Spp382, which play a role in restricting the activity of Dnl4-Lif1, the complex that executes the final ligation step of this process.