A high frequency of chromosomal duplications in unicellular algae is compensated by translational regulation

While duplications have long been recognized as a fundamental process driving major evolutionary innovations, direct estimates of spontaneous chromosome duplication rates, leading to aneuploid karyotypes, are scarce. Here, from mutation accumulation (MA) experiments, we provide the first estimates of spontaneous chromosome duplication rates in six unicellular eukaryotic species, which range from 1 × 10-4 to 1 × 10-3 per genome per generation. Although this is ∼4 to ∼50 times less frequent than spontaneous point mutations per genome, chromosome duplication events can affect 1 to 7% of the total genome size. In duplicated chromosomes, mRNA levels reflected gene copy numbers, but the level of translation estimated by polysome profiling revealed that dosage compensation must be occurring. In particular, one duplicated chromosome showed a 2.1-fold increase of mRNA but translation rates were decreased to 0.7-fold. Altogether, our results support previous observations of chromosome-dependent dosage compensation effects, providing evidence that compensation occurs during translation. We hypothesize that an unknown post-transcriptional mechanism modulates the translation of hundreds of transcripts from genes located on duplicated regions in eukaryotes.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Features of the Opportunistic Behaviour of the Marine Bacterium Marinobacter algicola in the Microalga Ostreococcus tauri Phycosphere

Although interactions between microalgae and bacteria are observed in both natural environment and the laboratory, the modalities of coexistence of bacteria inside microalgae phycospheres in laboratory cultures are mostly unknown. Here, we focused on well-controlled cultures of the model green picoalga Ostreococcus tauri and the most abundant member of its phycosphere, Marinobacter algicola. The prevalence of M. algicola in O. tauri cultures raises questions about how this bacterium maintains itself under laboratory conditions in the microalga culture. The results showed that M. algicola did not promote O. tauri growth in the absence of vitamin B₁₂ while M. algicola depended on O. tauri to grow in synthetic medium, most likely to obtain organic carbon sources provided by the microalgae. M. algicola grew on a range of lipids, including triacylglycerols that are known to be produced by O. tauri in culture during abiotic stress. Genomic screening revealed the absence of genes of two particular modes of quorum-sensing in Marinobacter genomes which refutes the idea that these bacterial communication systems operate in this genus. To date, the ‘opportunistic’ behaviour of M. algicola in the laboratory is limited to several phytoplanktonic species including Chlorophyta such as O. tauri. This would indicate a preferential occurrence of M. algicola in association with these specific microalgae under optimum laboratory conditions.

Metabolomic Insights into Marine Phytoplankton Diversity

The democratization of sequencing technologies fostered a leap in our knowledge of the diversity of marine phytoplanktonic microalgae, revealing many previously unknown species and lineages. The evolutionary history of the diversification of microalgae can be inferred from the analysis of their genome sequences. However, the link between the DNA sequence and the associated phenotype is notoriously difficult to assess, all the more so for marine phytoplanktonic microalgae for which the lab culture and, thus, biological experimentation is very tedious. Here, we explore the potential of a high-throughput untargeted metabolomic approach to explore the phenotypic-genotypic gap in 12 marine microalgae encompassing 1.2 billion years of evolution. We identified species- and lineage-specific metabolites. We also provide evidence of a very good correlation between the molecular divergence, inferred from the DNA sequences, and the metabolomic divergence, inferred from the complete metabolomic profiles. These results provide novel insights into the potential of chemotaxonomy in marine phytoplankton and support the hypothesis of a metabolomic clock, suggesting that DNA and metabolomic profiles co-evolve.

Fish larval recruitment to reefs is a thyroid hormone-mediated metamorphosis sensitive to the pesticide chlorpyrifos

Larval recruitment, the transition of pelagic larvae into reef-associated juveniles, is a critical step for the resilience of marine fish populations but its molecular control is unknown. Here, we investigate whether thyroid-hormones (TH) and their receptors (TR) coordinate the larval recruitment of the coral-reef-fish Acanthurus triostegus. We demonstrate an increase of TH-levels and TR-expressions in pelagic-larvae, followed by a decrease in recruiting juveniles. We generalize these observations in four other coral reef-fish species. Treatments with TH or TR-antagonist, as well as relocation to the open-ocean, disturb A. triostegus larvae transformation and grazing activity. Likewise, chlorpyrifos, a pesticide often encountered in coral-reefs, impairs A. triostegus TH-levels, transformation, and grazing activity, hence diminishing this herbivore’s ability to control the spread of reef-algae. Larval recruitment therefore corresponds to a TH-controlled metamorphosis, sensitive to endocrine disruption. This provides a framework to understand how larval recruitment, critical to reef-ecosystems maintenance, is altered by anthropogenic stressors.

Many animals go through a larval phase before developing into an adult. This transformation is called metamorphosis, and it is regulated by hormones of the thyroid gland in vertebrates. For example, most fish found on coral reefs actually spend the first part of their life as free-swimming larvae out in the ocean. The larvae usually look very different from the juveniles and adults. When these fish become juveniles, the larvae undergo a range of physical and behavioral changes to prepare for their life on the reef. Yet, until now it was not known what hormones control metamorphosis in these fish.

To address this question, Holzer, Besson et al. studied the convict surgeonfish Acanthurus triostegus. This herbivorous coral-reef fish lives in the Indo-Pacific Ocean, and the results showed that thyroid hormones do indeed regulate the metamorphosis of its larvae. This includes changing how the larvae behave and how their adult features develop. Further, Holzer, Besson et al. found that this was also true for four other coral-reef fish, including the lagoon triggerfish and the raccoon butterflyfish. In A. triostegus, thyroid hormones controlled the changes that enabled the juveniles to efficiently graze on algae growing on the reef such as an elongated gut.

When the fish larvae were then exposed to a pesticide called chlorpyrifos, a well-known reef pollutant, their hormone production was disturbed. This in turn affected their grazing behavior and also their metamorphosis. These fish had shortened, underdeveloped guts and could not graze on algae as effectively.

Herbivorous fish such as A. triostegus play a major role in supporting coral reef ecosystems by reducing algal cover and therefore promoting coral recruitment. These new findings show that pollutants from human activities could disturb the metamorphosis of coral-reef fish and, as a consequence, their ability to maintain the reefs. A next step will be to test what other factors can disrupt the hormones in coral-reef fish and thus pose a threat for fish populations and the coral-reef ecosystem.

A Viral Immunity Chromosome in the Marine Picoeukaryote, Ostreococcus tauri

Micro-algae of the genus Ostreococcus and related species of the order Mamiellales are globally distributed in the photic zone of world's oceans where they contribute to fixation of atmospheric carbon and production of oxygen, besides providing a primary source of nutrition in the food web. Their tiny size, simple cells, ease of culture, compact genomes and susceptibility to the most abundant large DNA viruses in the sea render them attractive as models for integrative marine biology. In culture, spontaneous resistance to viruses occurs frequently. Here, we show that virus-producing resistant cell lines arise in many independent cell lines during lytic infections, but over two years, more and more of these lines stop producing viruses. We observed sweeping over-expression of all genes in more than half of chromosome 19 in resistant lines, and karyotypic analyses showed physical rearrangements of this chromosome. Chromosome 19 has an unusual genetic structure whose equivalent is found in all of the sequenced genomes in this ecologically important group of green algae.

We propose that chromosome 19 of O. tauri is specialized in defence against viral attack, a constant threat for all planktonic life, and that the most likely cause of resistance is the over-expression of numerous predicted glycosyltransferase genes. O. tauri thus provides an amenable model for molecular analysis of genome evolution under environmental stress and for investigating glycan-mediated host-virus interactions, such as those seen in herpes, influenza, HIV, PBCV and mimivirus.

Marinobacter Dominates the Bacterial Community of the Ostreococcus tauri Phycosphere in Culture

Microalgal–bacterial interactions are commonly found in marine environments and are well known in diatom cultures maintained in laboratory. These interactions also exert strong effects on bacterial and algal diversity in the oceans. Small green eukaryote algae of the class Mamiellophyceae (Chlorophyta) are ubiquitous and some species, such as Ostreococcus spp., are particularly important in Mediterranean coastal lagoons, and are observed as dominant species during phytoplankton blooms in open sea. Despite this, little is known about the diversity of bacteria that might facilitate or hinder O. tauri growth. We show, using rDNA 16S sequences, that the bacterial community found in O. tauri RCC4221 laboratory cultures is dominated by γ-proteobacteria from the Marinobacter genus, regardless of the growth phase of O. tauri RCC4221, the photoperiod used, or the nutrient conditions (limited in nitrogen or phosphorous) tested. Several strains of Marinobacter algicola were detected, all closely related to strains found in association with taxonomically distinct organisms, particularly with dinoflagellates and coccolithophorids. These sequences were more distantly related to M. adhaerens, M. aquaeoli and bacteria usually associated to euglenoids. This is the first time, to our knowledge, that distinct Marinobacter strains have been found to be associated with a green alga in culture.

An improved genome of the model marine alga Ostreococcus tauri unfolds by assessing Illumina de novo assemblies

Background

Cost effective next generation sequencing technologies now enable the production of genomic datasets for many novel planktonic eukaryotes, representing an understudied reservoir of genetic diversity. O. tauri is the smallest free-living photosynthetic eukaryote known to date, a coccoid green alga that was first isolated in 1995 in a lagoon by the Mediterranean sea. Its simple features, ease of culture and the sequencing of its 13 Mb haploid nuclear genome have promoted this microalga as a new model organism for cell biology. Here, we investigated the quality of genome assemblies of Illumina GAIIx 75 bp paired-end reads from Ostreococcus tauri, thereby also improving the existing assembly and showing the genome to be stably maintained in culture.

Results

The 3 assemblers used, ABySS, CLCBio and Velvet, produced 95% complete genomes in 1402 to 2080 scaffolds with a very low rate of misassembly. Reciprocally, these assemblies improved the original genome assembly by filling in 930 gaps. Combined with additional analysis of raw reads and PCR sequencing effort, 1194 gaps have been solved in total adding up to 460 kb of sequence. Mapping of RNAseq Illumina data on this updated genome led to a twofold reduction in the proportion of multi-exon protein coding genes, representing 19% of the total 7699 protein coding genes. The comparison of the DNA extracted in 2001 and 2009 revealed the fixation of 8 single nucleotide substitutions and 2 deletions during the approximately 6000 generations in the lab. The deletions either knocked out or truncated two predicted transmembrane proteins, including a glutamate-receptor like gene.

Conclusion

High coverage (>80 fold) paired-end Illumina sequencing enables a high quality 95% complete genome assembly of a compact ~13 Mb haploid eukaryote. This genome sequence has remained stable for 6000 generations of lab culture.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1103) contains supplementary material, which is available to authorized users.

Evolution of codon usage in the smallest photosynthetic eukaryotes and their giant viruses

Prasinoviruses are among the largest viruses (>200 kb) and encode several hundreds of protein coding genes, including most genes of the DNA replication machinery and several genes involved in transcription and translation, as well as transfer RNAs (tRNAs). They can infect and lyse small eukaryotic planktonic marine green algae, thereby affecting global algal population dynamics. Here, we investigate the causes of codon usage bias (CUB) in one prasinovirus, OtV5, and its host Ostreococcus tauri, during a viral infection using microarray expression data. We show that 1) CUB in the host and in the viral genes increases with expression levels and 2) optimal codons use those tRNAs encoded by the most abundant host tRNA genes, supporting the notion of translational optimization by natural selection. We find evidence that viral tRNA genes complement the host tRNA pool for those viral amino acids whose host tRNAs are in short supply. We further discuss the coevolution of CUB in hosts and prasinoviruses by comparing optimal codons in three evolutionary diverged host–virus-specific pairs whose complete genome sequences are known.

Cryptic sex in the smallest eukaryotic marine green alga

Ostreococcus spp. are common worldwide oceanic picoeukaryotic pelagic algae. The complete genomes of three strains from different ecological niches revealed them to represent biologically distinct species despite their identical cellular morphologies (cryptic species). Their tiny genomes (13 Mb), with approximately 20 chromosomes, are colinear and densely packed with coding sequences, but no sexual life cycle has been described. Seventeen new strains of one of these species, Ostreococcus tauri, were isolated from 98 seawater samplings from the NW Mediterranean by filtering, culturing, cloning, and plating for single colonies and identification by sequencing their ribosomal 18S gene. In order to find the genetic markers for detection of polymorphisms and sexual recombination, we used an in silico approach to screen available genomic data. Intergenic regions of DNA likely to evolve neutrally were analyzed following polymerase chain reaction amplification of sequences using flanking primers from adjacent conserved coding sequences that were present as syntenic pairs in two different species of Ostreococcus. Analyses of such DNA regions from eight marker loci on two chromosomes from each strain revealed that the isolated O. tauri clones were haploid and that the overall level of polymorphism was approximately 0.01. Four different genetic tests for recombination showed that sexual exchanges must be inferred to account for the between-locus and between-chromosome marker combinations observed. However, our data suggest that sexual encounters are infrequent because we estimate the frequency of meioses/mitoses among the sampled strains to be 10(-6). Ostreococcus tauri and related species encode and express core genes for mitosis and meiosis, but their mechanisms of cell division and recombination, nevertheless, remain enigmatic because a classical eukaryotic spindle with 40 canonical microtubules would be much too large for the available approximately 0.9-microm(3) cellular volume.

Clues about the genetic basis of adaptation emerge from comparing the proteomes of two Ostreococcus ecotypes (Chlorophyta, Prasinophyceae)

We compared the proteomes of two picoplanktonic Ostreococcus unicellular green algal ecotypes to analyze the genetic basis of their adaptation with their ecological niches. We first investigated the function of the species-specific genes using Gene Ontology databases and similarity searches. Although most species-specific genes had no known function, we identified several species-specific functions involved in various cellular processes, which could be critical for environmental adaptations. Additionally, we investigated the rate of evolution of orthologous genes and its distribution across chromosomes. We show that faster evolving genes encode significantly more membrane or excreted proteins, consistent with the notion that selection acts on cell surface modifications that is driven by selection for resistance to viruses and grazers, keystone actors of phytoplankton evolution. The relationship between GC content and chromosome length also suggests that both strains have experienced recombination since their divergence and that lack of recombination on the two outlier chromosomes could explain part of their peculiar genomic features, including higher rates of evolution.

Life-cycle and genome of OtV5, a large DNA virus of the pelagic marine unicellular green alga Ostreococcus tauri

Large DNA viruses are ubiquitous, infecting diverse organisms ranging from algae to man, and have probably evolved from an ancient common ancestor. In aquatic environments, such algal viruses control blooms and shape the evolution of biodiversity in phytoplankton, but little is known about their biological functions. We show that Ostreococcus tauri, the smallest known marine photosynthetic eukaryote, whose genome is completely characterized, is a host for large DNA viruses, and present an analysis of the life-cycle and 186,234 bp long linear genome of OtV5. OtV5 is a lytic phycodnavirus which unexpectedly does not degrade its host chromosomes before the host cell bursts. Analysis of its complete genome sequence confirmed that it lacks expected site-specific endonucleases, and revealed the presence of 16 genes whose predicted functions are novel to this group of viruses. OtV5 carries at least one predicted gene whose protein closely resembles its host counterpart and several other host-like sequences, suggesting that horizontal gene transfers between host and viral genomes may occur frequently on an evolutionary scale. Fifty seven percent of the 268 predicted proteins present no similarities with any known protein in Genbank, underlining the wealth of undiscovered biological diversity present in oceanic viruses, which are estimated to harbour 200Mt of carbon.

Screening the Sargasso Sea metagenome for data to investigate genome evolution in Ostreococcus (Prasinophyceae, Chlorophyta)

The Sargasso Sea water shotgun sequencing unveiled an unprecedented glimpse of marine prokaryotic diversity and gene content. The sequence data was gathered from 0.8 microm filtered surface water extracts, and revealed picoeukaryotic (cell size<2 microm) sequences alongside the prokaryotic data. We used the available genome sequence of the picoeukaryote Ostreococcus tauri (Prasinophyceae, Chlorophyta) as a benchmark for the eukaryotic sequence content of the Sargasso Sea metagenome. Sequence data from at least two new Ostreococcus strains were identified and analyzed, and showed a bias towards higher coverage of the AT-rich organellar genomes. The Ostreococcus nuclear sequence data retrieved from the Sargasso metagenome is divided onto 731 scaffolds of average size 3917 bp, and covers 23% of the complete nuclear genome and 14% of the total number of protein coding genes in O. tauri. We used this environmental Ostreococcus sequence data to estimate the level of constraint on intronic and intergenic sequences in this compact genome.

A reanalysis of the indirect evidence for recombination in human mitochondrial DNA

In an attempt to resolve the controversy about whether recombination occurs in human mtDNA, we have analysed three recently published data sets of complete mtDNA sequences along with 10 RFLP data sets. We have analysed the relationship between linkage disequilibrium (LD) and distance between sites under a variety of conditions using two measures of LD, r2 and /D'/. We find that there is a negative correlation between r2 and distance in the majority of data sets, but no overall trend for /D'/. Five out of six mtDNA sequence data sets show an excess of homoplasy, but this could be due to either recombination or hypervariable sites. Two additional recombination detection methods used, Geneconv and Maximum Chi-Square, showed nonsignificant results. The overall significance of these findings is hard to quantify because of nonindependence, but our results suggest a lack of evidence for recombination in human mtDNA.

A broad survey of recombination in animal mitochondria

Recombination in mitochondrial DNA (mtDNA) remains a controversial topic. Here we present a survey of 279 animal mtDNA data sets, of which 12 were from asexual species. Using four separate tests, we show that there is widespread evidence of recombination; for one test as many as 14.2% of the data sets reject a model of clonal inheritance and in several data sets, including primates, the recombinants can be identified visually. We show that none of the tests give significant results for obligate clonal species (apomictic pathogens) and that the sexual species show significantly greater evidence of recombination than asexual species. For some data sets, such as Macaca nemestrina, additional data sets suggest that the recombinants are not artifacts. For others, it cannot be determined whether the recombinants are real or produced by laboratory error. Either way, the results have important implications for how mtDNA is sequenced and used.

Selection pressure-driven evolution of the Epstein-Barr virus-encoded oncogene LMP1 in virus isolates from Southeast Asia

The geographically constrained distribution of Epstein-Barr virus (EBV)-associated nasopharyngeal carcinoma (NPC) in southeast Asian populations suggests that both viral and host genetics may influence disease risk. Although susceptibility loci have been mapped within the human genome, the role of viral genetics in the focal distribution of NPC remains an enigma. Here we report a molecular phylogenetic analysis of an NPC-associated viral oncogene, LMP1, in a large panel of EBV isolates from southeast Asia and from Papua New Guinea, Africa, and Australia, regions of the world where NPC is and is not endemic, respectively. This analysis revealed that LMP1 sequences show a distinct geographic structure, indicating that the southeast Asian isolates have evolved as a lineage distinct from those of Papua New Guinea, African, and Australian isolates. Furthermore, a likelihood ratio test revealed that the C termini of the LMP1 sequences of the southeast Asian lineage are under significant positive selection pressure, particularly at some sites within the C-terminal activator regions. We also present evidence that although the N terminus and transmembrane region of LMP1 have undergone recombination, the C-terminal region of the gene has evolved without any history of recombination. Based on these observations, we speculate that selection pressure may be driving the LMP1 sequences in virus isolates from southeast Asia towards a more malignant phenotype, thereby influencing the endemic distribution of NPC in this region.

Estimating the distribution of fitness effects from DNA sequence data: implications for the molecular clock

We present a method for estimating the distribution of fitness effects of new amino acid mutations when those mutations can be assumed to be slightly advantageous, slightly deleterious, or strongly deleterious. We apply the method to mitochondrial data from several different species. In the majority of the data sets, the shape of the distribution is approximately exponential. Our results provide an estimate of the distribution of fitness effects of weakly selected mutations and provide a possible explanation for why the molecular clock is fairly constant across taxa and time.

Vanishing GC-rich isochores in mammalian genomes

To understand the origin and evolution of isochores-the peculiar spatial distribution of GC content within mammalian genomes-we analyzed the synonymous substitution pattern in coding sequences from closely related species in different mammalian orders. In primate and cetartiodactyls, GC-rich genes are undergoing a large excess of GC --> AT substitutions over AT --> GC substitutions: GC-rich isochores are slowly disappearing from the genome of these two mammalian orders. In rodents, our analyses suggest both a decrease in GC content of GC-rich isochores and an increase in GC-poor isochores, but more data will be necessary to assess the significance of this pattern. These observations question the conclusions of previous works that assumed that base composition was at equilibrium. Analysis of allele frequency in human polymorphism data, however, confirmed that in the GC-rich parts of the genome, GC alleles have a higher probability of fixation than AT alleles. This fixation bias appears not strong enough to overcome the large excess of GC --> AT mutations. Thus, whatever the evolutionary force (neutral or selective) at the origin of GC-rich isochores, this force is no longer effective in mammals. We propose a model based on the biased gene conversion hypothesis that accounts for the origin of GC-rich isochores in the ancestral amniote genome and for their decline in present-day mammals.

Hill-Robertson interference is a minor determinant of variations in codon bias across Drosophila melanogaster and Caenorhabditis elegans genomes

According to population genetics models, genomic regions with lower crossing-over rates are expected to experience less effective selection because of Hill-Robertson interference (HRi). The effect of genetic linkage is thought to be particularly important for a selection of weak intensity such as selection affecting codon usage. Consistent with this model, codon bias correlates positively with recombination rate in Drosophila melanogaster and Caenorhabditis elegans. However, in these species, the G+C content of both noncoding DNA and synonymous sites correlates positively with recombination, which suggests that mutation patterns and recombination are associated. To remove this effect of mutation patterns on codon bias, we used the synonymous sites of lowly expressed genes that are expected to be effectively neutral sites. We measured the differences between codon biases of highly expressed genes and their lowly expressed neighbors. In D. melanogaster we find that HRi weakly reduces selection on codon usage of genes located in regions of very low recombination; but these genes only comprise 4% of the total. In C. elegans we do not find any evidence for the effect of recombination on selection for codon bias. Computer simulations indicate that HRi poorly enhances codon bias if the local recombination rate is greater than the mutation rate. This prediction of the model is consistent with our data and with the current estimate of the mutation rate in D. melanogaster. The case of C. elegans, which is highly self-fertilizing, is discussed. Our results suggest that HRi is a minor determinant of variations in codon bias across the genome.

Expected relationship between the silent substitution rate and the GC content: implications for the evolution of isochores

The relationship between the silent substitution rate (Ks) and the GC content along the genome is a focal point of the debate about the origin of the isochore structure in vertebrates. Recent estimation of the silent substitution rate showed a positive correlation between Ks and GC content, in contradiction with the predictions of both the regional mutation bias model and the selection or biased gene conversion model. The aim of this paper is to help resolve this contradiction between theoretical studies and data. We analyzed the relationship between Ks and GC content under (1) uniform mutation bias, (2) a regional mutation bias, and (3) mutation bias and selection. We report that an increase in Ks with GC content is expected under mutation bias because of either nonequilibrium of the isochore structure or an increasing mutation rate from AT toward GC nucleotides in GC-richer isochores. We show by simulations that CpG deamination tends to increase the mutation rate with GC content in a regional mutation bias model. We also demonstrate that the relationship between Ks and GC under the selectionist or biased gene conversion model is positive under weak selection if the mutation selection equilibrium GC frequency is less than 0.5.

Multiplicative versus additive selection in relation to genome evolution: a simulation study

The evolution of molecular quantitative traits, such as codon usage bias or base frequencies, can be explained as the result of mutational biases alone, or as the result of mutation and selection. Whereas mutation models can be investigated easily, realistic modelling of selection-directed genome evolution is analytically intractable, and numerical calculations require substantial computer resources. We investigated the evolution of optimal codon frequency under additive and multiplicative effects of selected linked codons. We show that additive selective effects of many linked sites cannot be effective in genomes when the number of selected sites is greater than the effective population size, a realistic assumption according to current molecular data. We then discuss the implications of these results for isochore evolution in vertebrates.

High copy numbers of multiple transposable element families in an Australian population of Drosophila simulans

Sudden mobilization of transposable elements in Drosophila is a well-reported phenomenon but one that usually affects no more than a few elements (one to four). We report here the existence of a D. simulans natural population (Canberra) from Australia, which had high copy numbers for various transposable elements (transposons, LTR retrotransposons and non-LTR retrotransposons). The impact of transposable elements on the host genome and populations is discussed.