Mutational and transcriptional landscape of spontaneous gene duplications and deletions in Caenorhabditis elegans.

Gene duplication and deletion are pivotal processes shaping the structural and functional repertoire of genomes, with implications for disease, adaptation, and evolution. We employed a mutation accumulation (MA) framework partnered with high-throughput genomics to assess the molecular and transcriptional characteristics of newly arisen gene copy-number variants (CNVs) in Caenorhabditis elegans populations subjected to varying intensity of selection. Here, we report a direct spontaneous genome-wide rate of gene duplication of 2.9 × 10-5/gene per generation in C. elegans, the highest for any species to date. The rate of gene deletion is sixfold lower (5 × 10-6/gene per generation). Deletions of highly expressed genes are particularly deleterious, given their paucity in even the N = 1 lines with minimal efficacy of selection. The increase in average transcript abundance of new duplicates arising under minimal selection is significantly greater than twofold compared with single copies of the same gene, suggesting that genes in segmental duplications are frequently overactive at inception. The average increase in transcriptional activity of gene duplicates is greater in the N = 1 MA lines than in MA lines with larger population bottlenecks. There is an inverse relationship between the ancestral transcription levels of new gene duplicates and population size, with duplicate copies of highly expressed genes less likely to accumulate in larger populations. Our results demonstrate a fitness cost of increased transcription following duplication, which results in purifying selection against new gene duplicates. However, on average, duplications also provide a significant increase in gene expression that can facilitate adaptation to novel environmental challenges.

Long-term experimental evolution reveals purifying selection on piRNA-mediated control of transposable element expression.

BACKGROUND: Transposable elements (TEs) are an almost universal constituent of eukaryotic genomes. In animals, Piwi-interacting small RNAs (piRNAs) and repressive chromatin often play crucial roles in preventing TE transcription and thus restricting TE activity. Nevertheless, TE content varies widely across eukaryotes and the dynamics of TE activity and TE silencing across evolutionary time is poorly understood. RESULTS: Here, we used experimentally evolved populations of C. elegans to study the dynamics of TE expression over 409 generations. The experimental populations were evolved at population sizes of 1, 10 and 100 individuals to manipulate the efficiency of natural selection versus genetic drift. We demonstrate increased TE expression relative to the ancestral population, with the largest increases occurring in the smallest populations. We show that the transcriptional activation of TEs within active regions of the genome is associated with failure of piRNA-mediated silencing, whilst desilenced TEs in repressed chromatin domains retain small RNAs. Additionally, we find that the sequence context of the surrounding region influences the propensity of TEs to lose silencing through failure of small RNA-mediated silencing. CONCLUSIONS: Our results show that natural selection in C. elegans is responsible for maintaining low levels of TE expression, and provide new insights into the epigenomic features responsible.

Mitochondrial Mutation Rate, Spectrum and Heteroplasmy in Caenorhabditis elegans Spontaneous Mutation Accumulation Lines of Differing Population Size.

Mitochondrial genomes of metazoans, given their elevated rates of evolution, have served as pivotal markers for phylogeographic studies and recent phylogenetic events. In order to determine the dynamics of spontaneous mitochondrial mutations in small populations in the absence and presence of selection, we evolved mutation accumulation (MA) lines of Caenorhabditis elegans in parallel over 409 consecutive generations at three varying population sizes of N = 1, 10, and 100 hermaphrodites. The N =1 populations should have a minimal influence of natural selection to provide the spontaneous mutation rate and the expected rate of neutral evolution, whereas larger population sizes should experience increasing intensity of selection. New mutations were identified by Illumina paired-end sequencing of 86 mtDNA genomes across 35 experimental lines and compared with published genomes of natural isolates. The spontaneous mitochondrial mutation rate was estimated at 1.05 × 10-7/site/generation. A strong G/C→A/T mutational bias was observed in both the MA lines and the natural isolates. This suggests that the low G + C content at synonymous sites is the product of mutation bias rather than selection as previously proposed. The mitochondrial effective population size per worm generation was estimated to be 62. Although it was previously concluded that heteroplasmy was rare in C. elegans, the vast majority of mutations in this study were heteroplasmic despite an experimental regime exceeding 400 generations. The frequencies of frameshift and nonsynonymous mutations were negatively correlated with population size, which suggests their deleterious effects on fitness and a potent role for selection in their eradication.

The phylogenetic distribution and evolution of enzymes within the thymidine kinase 2-like gene family in metazoa.

Deoxyribonucleoside kinases (dNKs) carry out the rate-determining step in the nucleoside salvage pathway within all domains of life where the pathway is present, and, hence, are an indication on whether or not a species/genus retains the ability to salvage deoxyribonucleosides. Here, a phylogenetic tree is constructed for the thymidine kinase 2-like dNK gene family in metazoa. Each enzyme class (deoxycytidine, deoxyguanosine, and deoxythymidine kinases, as well as the multisubstrate dNKs) falls into a monophyletic clade. However, in vertebrates, dCK contains an apparent duplication with one paralog lost in mammals, and a number of crustacean genomes (like Caligus rogercresseyi and Lepeophtheirus salmonis) unexpectedly contain not only the multisubstrate dNKs, related to Drosophila multisubstrate dNK, but also a TK2-like kinase. Additionally, crustaceans (Daphnia, Caligus, and Lepeophtheirus) and some insects (Tribolium, Danaus, Pediculus, and Acyrthosiphon) contain several multisubstrate dNK-like enzymes which group paraphyletically within the arthropod clade. This might suggest that the multisubstrate dNKs underwent multiple rounds of duplications with differential retention of duplicate copies between insect families and more complete retention within some crustaceans and insects. Genomes of several basal animalia contain more than one dNK-like sequence, some of which group outside the remaining eukaryotes (both plants and animals) and/or with bacterial dNKs. Within the vertebrates, the mammalian genomes do not contain the second dCK, while birds, fish, and amphibians do retain it. Phasianidae (chicken and turkey) have lost dGK, while it has been retained in other bird lineages, like zebra finch. Reconstruction of the ancestral sequence between the multisubstrate arthropod dNKs and the TK2 clade of vertebrates followed by homology modeling and discrete molecular dynamics calculations on this sequence were performed to examine the evolutionary path which led to the two different enzyme classes. The structural models showed that the carboxyl terminus of the ancestral sequence is more helical than dNK, in common with TK2, although any implications of this for enzyme specificity will require biochemical validation. Finally, rate-shift and conservation-shift analysis between clades with different specificities uncovered candidate residues outside the active site pocket which may have contributed to differentiation in substrate specificity between enzyme clades.

Mutation rate and spectrum in obligately outcrossing Caenorhabditis elegans mutation accumulation lines subjected to RNAi-induced knockdown of the mismatch repair gene msh-2.

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutSα homolog, msh-2, in Caenorhabditis elegans male-female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male-female background substantially increased the frequency of nuclear base substitutions (∼23×) and small indels (∼328×) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

Mutation accumulation and horizontal gene transfer in Escherichia coli colonizing the gut of old mice.

The ecology and environment of the microbes that inhabit the mammalian intestine undergoes several changes as the host ages. Here, we ask if the selection pressure experienced by a new strain colonizing the aging gut differs from that in the gut of young adults. Using experimental evolution in mice after a short antibiotic treatment, as a model for a common clinical situation, we show that a new colonizing E. coli strain rapidly adapts to the aging gut via both mutation accumulation and bacteriophage-mediated horizontal gene transfer (HGT). The pattern of evolution of E. coli in aging mice is characterized by a larger number of transposable element insertions and intergenic mutations compared to that in young mice, which is consistent with the gut of aging hosts harboring a stressful and iron limiting environment.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size.

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

The evolution of catalytic residues and enzyme mechanism within the bacterial nucleoside phosphorylase superfamily 1.

Nucleoside phosphorylases are essential for the salvage and catabolism of nucleotides in bacteria and other organisms, and members of this enzyme superfamily have been of interest for the development of antimicrobial and cancer therapies. The nucleotide phosphorylase superfamily 1 encompasses a number of different enzymes which share a general superfold and catalytic mechanism, while they differ in the nature of the nucleophiles used and in the nature of characteristic active site residues. Recently, one subfamily, the uridine phosphorylases, has been subdivided into two types which differ with respect to the mechanism of transition state stabilization, as dictated by differences in critical amino acid residues. Little is known about the phylogenetic distribution and relationship of the two different types, as well as the relationship to other NP-1 superfamily members. Here comparative genomic analysis illustrates that UP-1s and UP-2s fall into monophyletic groups and are biased with respect to species representation. UP-1 evolved in Gram negative bacteria, while Gram positive species tend to predominantly contain UP-2. PNP (a sister clade to all UPs) contains both Gram positive and Gram negative species. The findings imply that the nucleoside phosphorylase superfamily 1 evolved through a series of three important duplications, leading to the separate, monophyletic enzyme families, coupled to individual lateral transfer events. Extensive horizontal transfer explains the occurrence of unexpected uridine phosphorylases in some genomes. This study provides a basis for understanding the evolution of uridine and purine nucleoside phosphorylases with respect to DNA/RNA metabolism and with potential utility in the design of antimicrobial and anti-tumor drugs.

The global distribution and evolution of deoxyribonucleoside kinases in bacteria.

Deoxyribonucleoside kinases (dNKs) are important to DNA metabolism, especially in environments where nucleosides are freely available to be absorbed and used for the salvage pathway. Little has previously been known about the complement of dNKs in different bacterial genomes. However, it was believed that Gram-negative bacteria had a single dNK, while Gram-positive bacteria possessed several. An analysis of 992 fully sequenced bacterial genomes, including both Gram-positive and Gram-negative organisms, was conducted to investigate the phylogenetic relationship of all TK1-like and non-TK1-like dNKs. It was illustrated that both gene families evolved through a number of duplications and horizontal gene transfers, leading to the presence of multiple dNKs in different types of bacteria. The findings of this study provide a backbone for further studies into the evolution of the interplay between the de novo and salvage pathways in DNA synthesis with respect to environmental availability of deoxyribonucleosides and metabolic processes generating the provisions of different dNTPs.

Deoxyribonucleoside kinases in two aquatic bacteria with high specificity for thymidine and deoxyadenosine.

Deoxyribonucleoside kinases (dNKs) are essential in the mammalian cell but their 'importance' in bacteria, especially aquatic ones, is less clear. We studied two aquatic bacteria, Gram-negative Flavobacterium psychrophilum JIP02/86 and Polaribacter sp. MED152, for their ability to salvage deoxyribonucleosides (dNs). Both had a Gram-positive-type thymidine kinase (TK1), which could phosphorylate thymidine, and one non-TK1 dNK, which could efficiently phosphorylate deoxyadenosine and slightly also deoxycytosine. Surprisingly, the four tested dNKs could not phosphorylate deoxyguanosine, and apparently, these two bacteria are missing this activity. When tens of available aquatic bacteria genomes were examined for the presence of dNKs, a majority had at least a TK1-like gene, but several lacked any dNKs. Apparently, among aquatic bacteria, the role of the dN salvage varies.

Toward a general model for the evolutionary dynamics of gene duplicates.

Gene duplication is an important process in the functional divergence of genes and genomes. Several processes have been described that lead to duplicate gene retention over different timescales after both smaller-scale events and whole-genome duplication, including neofunctionalization, subfunctionalization, and dosage balance. Two common modes of duplicate gene loss include nonfunctionalization and loss due to population dynamics (failed fixation). Previous work has characterized expectations of duplicate gene retention under the neofunctionalization and subfunctionalization models. Here, that work is extended to dosage balance using simulations. A general model for duplicate gene loss/retention is then presented that is capable of fitting expectations under the different models, is defined at t = 0, and decays to an orthologous asymptotic rate rather than zero, based upon a modified Weibull hazard function. The model in a maximum likelihood framework shows the property of identifiability, recovering the evolutionary mechanism and parameters of simulation. This model is also capable of recovering the evolutionary mechanism of simulation from data generated using an unrelated network population genetic model. Lastly, the general model is applied as part of a mixture model to recent gene duplicates from the Oikopleura dioica genome, suggesting that neofunctionalization may be an important process leading to duplicate gene retention in that organism.

Plasticity of animal genome architecture unmasked by rapid evolution of a pelagic tunicate.

Genomes of animals as different as sponges and humans show conservation of global architecture. Here we show that multiple genomic features including transposon diversity, developmental gene repertoire, physical gene order, and intron-exon organization are shattered in the tunicate Oikopleura, belonging to the sister group of vertebrates and retaining chordate morphology. Ancestral architecture of animal genomes can be deeply modified and may therefore be largely nonadaptive. This rapidly evolving animal lineage thus offers unique perspectives on the level of genome plasticity. It also illuminates issues as fundamental as the mechanisms of intron gain.

Growth of keratinocytes on porous films of poly(3-hydroxybutyrate) and poly(4-hydroxybutyrate) blended with hyaluronic acid and chitosan.

The objective of this study was to develop novel absorbable films suitable for use as a tissue-engineering scaffold for keratinocytes as a therapy for replacement of damaged skin. Poly(4-hydroxybutyrate) (P(4HB)) and poly (3-hydroxybutyrate) (P(3HB)) were blended with small amounts of the polysaccharides hyaluronic acid (HA), chitosan (CH), pectin and alginic acid, and were solution cast to produce porous films. The resulting composites had favorable mechanical properties, and these films were compared with two commercially available implantable films made of poly(L-lactide-co-D,L-lactide) (PLA copolymer) and HA benzyl ester. Tensile testing demonstrated that a high level of flexibility of P(4HB) was retained in the P(4HB)-polysaccharide composite films, whereas the P(3HB) film and its polysaccharide composites were stiffer and more brittle. The proliferation kinetics of adherent HaCaT keratinocytes on the films was examined in vitro. The porous surface of the P(4HB) and P(3HB) films blended with HA or CH promoted the growth of keratinocytes significantly. The order of maximum cell numbers on these films was P(4HB)/HA > P(4HB)/CH > P(3HB)/HA > P(3HB)/CH > P(3HB)/pectin > P(3HB)/alginic acid. Scanning electron microscopy and confocal laser scanning microscopy revealed differences in cell growth. Cells formed clusters on P(3HB) and its composites, while the cells grew as a confluent layer on P(4HB) and its composites. HaCaT cells formed large numbers of filaments only on P(4HB) films, indicating the excellent biocompatibility of this material. For the nonporous PHB films, the proliferation rate of cells was found to increase with decreasing hydrophobicity in the order: P(4HB) > P(3HB)/P(4HB) blend > P(3HB).