High mutational rates of large-scale duplication and deletion in Daphnia pulex.

Knowledge of the genome-wide rate and spectrum of mutations is necessary to understand the origin of disease and the genetic variation driving all evolutionary processes. Here, we provide a genome-wide analysis of the rate and spectrum of mutations obtained in two Daphnia pulex genotypes via separate mutation-accumulation (MA) experiments. Unlike most MA studies that utilize haploid, homozygous, or self-fertilizing lines, D. pulex can be propagated ameiotically while maintaining a naturally heterozygous, diploid genome, allowing the capture of the full spectrum of genomic changes that arise in a heterozygous state. While base-substitution mutation rates are similar to those in other multicellular eukaryotes (about 4 × 10(-9) per site per generation), we find that the rates of large-scale (>100 kb) de novo copy-number variants (CNVs) are significantly elevated relative to those seen in previous MA studies. The heterozygosity maintained in this experiment allowed for estimates of gene-conversion processes. While most of the conversion tract lengths we report are similar to those generated by meiotic processes, we also find larger tract lengths that are indicative of mitotic processes. Comparison of MA lines to natural isolates reveals that a majority of large-scale CNVs in natural populations are removed by purifying selection. The mutations observed here share similarities with disease-causing, complex, large-scale CNVs, thereby demonstrating that MA studies in D. pulex serve as a system for studying the processes leading to such alterations.

Common transcriptional mechanisms for visual photoreceptor cell differentiation among Pancrustaceans.

A hallmark of visual rhabdomeric photoreceptors is the expression of a rhabdomeric opsin and uniquely associated phototransduction molecules, which are incorporated into a specialized expanded apical membrane, the rhabdomere. Given the extensive utilization of rhabdomeric photoreceptors in the eyes of protostomes, here we address whether a common transcriptional mechanism exists for the differentiation of rhabdomeric photoreceptors. In Drosophila, the transcription factors Pph13 and Orthodenticle (Otd) direct both aspects of differentiation: rhabdomeric opsin transcription and rhabdomere morphogenesis. We demonstrate that the orthologs of both proteins are expressed in the visual systems of the distantly related arthropod species Tribolium castaneum and Daphnia magna and that their functional roles are similar in these species. In particular, we establish that the Pph13 homologs have the ability to bind a subset of Rhodopsin core sequence I sites and that these sites are present in key phototransduction genes of both Tribolium and Daphnia. Furthermore, Pph13 and Otd orthologs are capable of executing deeply conserved functions of photoreceptor differentiation as evidenced by the ability to rescue their respective Drosophila mutant phenotypes. Pph13 homologs are equivalent in their ability to direct both rhabdomere morphogenesis and opsin expression within Drosophila, whereas Otd paralogs demonstrate differential abilities to regulate photoreceptor differentiation. Finally, loss-of-function analyses in Tribolium confirm the conserved requirement of Pph13 and Otd in regulating both rhabdomeric opsin transcription and rhabdomere morphogenesis. Taken together, our data identify components of a regulatory framework for rhabdomeric photoreceptor differentiation in Pancrustaceans, providing a foundation for defining ancestral regulatory modules of rhabdomeric photoreceptor differentiation.

Population-genomic insights into the evolutionary origin and fate of obligately asexual Daphnia pulex.

Despite much theoretical work, the molecular-genetic causes and evolutionary consequences of asexuality remain largely undetermined. Asexual animal species are rare, evolutionarily short-lived, and thought to suffer mutational meltdown as a result of lack of recombination. Whole-genome analysis of 11 sexual and 11 asexual genotypes of Daphnia pulex indicates that current asexual lineages are in fact very young, exhibit no signs of purifying selection against accumulating mutations, and have extremely high rates of gene conversion and deletion. The reconstruction of chromosomal haplotypes in regions containing SNP markers associated with asexuality (chromosomes VIII and IX) indicates that introgression from a sister species, Daphnia pulicaria, underlies the origin of the asexual phenotype. Silent-site divergence of the shared chromosomal haplotypes of asexuals indicates that the spread of asexuality is as recent as 1,250 y, although the origin of the meiosis-suppressing element or elements could be substantially older. In addition, using previous estimates of the gene conversion rate from Daphnia mutation accumulation lines, we are able to age each asexual lineage. Although asexual lineages originate from wide crosses that introduce elevated individual heterozygosities on clone foundation, they also appear to be constrained by the inbreeding-like effect of loss of heterozygosity that accrues as gene conversion and hemizygous deletion expose preexisting recessive deleterious alleles of asexuals, limiting their evolutionary longevity. Our study implies that the buildup of newly introduced deleterious mutations (i.e., Muller's ratchet) may not be the dominant force imperiling nonrecombining populations of D. pulex, as previously proposed.

Extensive error in the number of genes inferred from draft genome assemblies.

Current sequencing methods produce large amounts of data, but genome assemblies based on these data are often woefully incomplete. These incomplete and error-filled assemblies result in many annotation errors, especially in the number of genes present in a genome. In this paper we investigate the magnitude of the problem, both in terms of total gene number and the number of copies of genes in specific families. To do this, we compare multiple draft assemblies against higher-quality versions of the same genomes, using several new assemblies of the chicken genome based on both traditional and next-generation sequencing technologies, as well as published draft assemblies of chimpanzee. We find that upwards of 40% of all gene families are inferred to have the wrong number of genes in draft assemblies, and that these incorrect assemblies both add and subtract genes. Using simulated genome assemblies of Drosophila melanogaster, we find that the major cause of increased gene numbers in draft genomes is the fragmentation of genes onto multiple individual contigs. Finally, we demonstrate the usefulness of RNA-Seq in improving the gene annotation of draft assemblies, largely by connecting genes that have been fragmented in the assembly process.

Extraordinary genome stability in the ciliate Paramecium tetraurelia.

Mutation plays a central role in all evolutionary processes and is also the basis of genetic disorders. Established base-substitution mutation rates in eukaryotes range between ∼5 × 10(-10) and 5 × 10(-8) per site per generation, but here we report a genome-wide estimate for Paramecium tetraurelia that is more than an order of magnitude lower than any previous eukaryotic estimate. Nevertheless, when the mutation rate per cell division is extrapolated to the length of the sexual cycle for this protist, the measure obtained is comparable to that for multicellular species with similar genome sizes. Because Paramecium has a transcriptionally silent germ-line nucleus, these results are consistent with the hypothesis that natural selection operates on the cumulative germ-line replication fidelity per episode of somatic gene expression, with the germ-line mutation rate per cell division evolving downward to the lower barrier imposed by random genetic drift. We observe ciliate-specific modifications of widely conserved amino acid sites in DNA polymerases as one potential explanation for unusually high levels of replication fidelity.

p53 Superfamily proteins in marine bivalve cancer and stress biology.

The human p53 tumour suppressor protein is inactivated in many cancers and is also a major player in apoptotic responses to cellular stress. The p53 protein and the two other members of this protein family (p63, p73) are encoded by distinct genes and their functions have been extensively documented for humans and some other vertebrates. The structure and relative expression levels for members of the p53 superfamily have also been reported for most major invertebrate taxa. The functions of homologous proteins have been investigated for only a few invertebrates (specifically, p53 in flies, nematodes and recently a sea anemone). These studies of classical model organisms all suggest that the gene family originally evolved to mediate apoptosis of damaged germ cells or to protect germ cells from genotoxic stress. Here, we have correlated data from a number of molluscan and other invertebrate sequencing projects to provide a framework for understanding p53 signalling pathways in marine bivalve cancer and stress biology. These data suggest that (a) the two identified p53 and p63/73-like proteins in soft shell clam (Mya arenaria), blue mussel (Mytilus edulis) and Northern European squid (Loligo forbesi) have identical core sequences and may be splice variants of a single gene, while some molluscs and most other invertebrates have two or more distinct genes expressing different p53 family members; (b) transcriptional activation domains (TADs) in bivalve p53 and p63/73-like protein sequences are 67-69% conserved with human p53, while those in ecdysozoan, cnidarian, placozoan and choanozoan eukaryotes are ≤33% conserved; (c) the Mdm2 binding site in the transcriptional activation domain is 100% conserved in all sequenced bivalve p53 proteins (e.g. Mya, Mytilus, Crassostrea and Spisula) but is not present in other non-deuterostome invertebrates; (d) an Mdm2 homologue has been cloned for Mytilus trossulus; (e) homologues for both human p53 upstream regulatory and transcriptional target genes exist in molluscan genomes (missing are ARF, CIP1 and BH3 only proteins) and (f) p53 is demonstrably involved in bivalve haemocyte and germinoma cancers. We usually do not know enough about the molecular biology of marine invertebrates to address molecular mechanisms that characterize particular diseases. Understanding the molecular basis of naturally occurring diseases in marine bivalves is a virtually unexplored aspect of toxicoproteomics and genomics and related drug discovery. Additionally, increases in coastal development and concomitant increases in aquatic pollutants have driven interest in developing models appropriate for evaluating potential hazardous compounds or conditions found in the aquatic environment. Data reviewed in this study are coupled with recent developments in our understanding the molecular biology of the marine bivalve p53 superfamily. Taken together, they suggest that both structurally and functionally, bivalve p53 family proteins are the most highly conserved members of this gene superfamily so far identified outside of higher vertebrates and invertebrate chordates. Marine bivalves provide some of the most relevant and best understood models currently available for experimental studies by biomedical and marine environmental researchers.

Extensive, recent intron gains in Daphnia populations.

Rates and mechanisms of intron gain and loss have traditionally been inferred from alignments of highly conserved genes sampled from phylogenetically distant taxa. We report a population-genomic approach that detected 24 discordant intron/exon boundaries between the whole-genome sequences of two Daphnia pulex isolates. Sequencing of presence/absence loci across a collection of D. pulex isolates and outgroup Daphnia species shows that most polymorphisms are a consequence of recent gains, with parallel gains often occurring at the same locations in independent allelic lineages. More than half of the recent gains are associated with short sequence repeats, suggesting an origin via repair of staggered double-strand breaks. By comparing the allele-frequency spectrum of intron-gain alleles with that for derived single-base substitutions, we also provide evidence that newly arisen introns are intrinsically deleterious and tend to accumulate in population-genetic settings where random genetic drift is a relatively strong force.

Evaluating high-throughput sequencing as a method for metagenomic analysis of nematode diversity.

Nematodes play an important role in ecosystem processes, yet the relevance of nematode species diversity to ecology is unknown. Because nematode identification of all individuals at the species level using standard techniques is difficult and time-consuming, nematode communities are not resolved down to the species level, leaving ecological analysis ambiguous. We assessed the suitability of massively parallel sequencing for analysis of nematode diversity from metagenomic samples. We set up four artificial metagenomic samples involving 41 diverse reference nematodes in known abundances. Two samples came from pooling polymerase chain reaction products amplified from single nematode species. Two additional metagenomic samples consisted of amplified products of DNA extracted from pooled nematode species. Amplified products involved two rapidly evolving ~400-bp sections coding for the small and large subunit of rRNA. The total number of reads ranged from 4159 to 14771 per metagenomic sample. Of these, 82% were > 199 bp in length. Among the reads > 199 bp, 86% matched the referenced species with less than three nucleotide differences from a reference sequence. Although neither rDNA section recovered all nematode species, the use of both loci improved the detection level of nematode species from 90 to 97%. Overall, results support the suitability of massively parallel sequencing for identification of nematodes. In contrast, the frequency of reads representing individual species did not correlate with the number of individuals in the metagenomic samples, suggesting that further methodological work is necessary before it will be justified for inferring the relative abundances of species within a nematode community.