First report of molecular characterization and phylogeny of Trichuris fossor Hall, 1916 (Nematoda: Trichuridae).

Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.

Genome Announcement: The Draft Genomes of Two Radopholus similis populations from Costa Rica.

Radopholus similis is an economically important pest of both banana and citrus in tropical regions. Here we present draft genomes from two populations of R. similis from Costa Rica that were created and assembled using short read libraries from Illumina HiSeq technology.

Transcription through the eye of a needle: daily and annual cyclic gene expression variation in Douglas-fir needles.

BACKGROUND: Perennial growth in plants is the product of interdependent cycles of daily and annual stimuli that induce cycles of growth and dormancy. In conifers, needles are the key perennial organ that integrates daily and seasonal signals from light, temperature, and water availability. To understand the relationship between seasonal cycles and seasonal gene expression responses in conifers, we examined diurnal and circannual needle mRNA accumulation in Douglas-fir (Pseudotsuga menziesii) needles at diurnal and circannual scales. Using mRNA sequencing, we sampled 6.1 × 109 reads from 19 trees and constructed a de novo pan-transcriptome reference that includes 173,882 tree-derived transcripts. Using this reference, we mapped RNA-Seq reads from 179 samples that capture daily and annual variation. RESULTS: We identified 12,042 diurnally-cyclic transcripts, 9299 of which showed homology to annotated genes from other plant genomes, including angiosperm core clock genes. Annual analysis revealed 21,225 circannual transcripts, 17,335 of which showed homology to annotated genes from other plant genomes. The timing of maximum gene expression is associated with light intensity at diurnal scales and photoperiod at annual scales, with approximately half of transcripts reaching maximum expression +/- 2 h from sunrise and sunset, and +/- 20 days from winter and summer solstices. Comparisons with published studies from other conifers shows congruent behavior in clock genes with Japanese cedar (Cryptomeria), and a significant preservation of gene expression patterns for 2278 putative orthologs from Douglas-fir during the summer growing season, and 760 putative orthologs from spruce (Picea) during the transition from fall to winter. CONCLUSIONS: Our study highlight the extensive diurnal and circannual transcriptome variability demonstrated in conifer needles. At these temporal scales, 29% of expressed transcripts show a significant diurnal cycle, and 58.7% show a significant circannual cycle. Remarkably, thousands of genes reach their annual peak activity during winter dormancy. Our study establishes the fine-scale timing of daily and annual maximum gene expression for diverse needle genes in Douglas-fir, and it highlights the potential for using this information for evaluating hypotheses concerning the daily or seasonal timing of gene activity in temperate-zone conifers, and for identifying cyclic transcriptome components in other conifer species.

Paternal Mitochondrial Transmission in Intra-Species Caenorhabditis briggsae Hybrids.

To study mitochondrial-nuclear genetic interactions in the nematode Caenorhabditis briggsae, our three laboratories independently created 38 intra-species cytoplasmic-nuclear hybrid (cybrid) lines. Although the cross design combines maternal mitotypes with paternal nuclear genotypes, eight lines (21%) unexpectedly contained paternal mitotypes. All eight share in common ancestry of one of two genetically related strains. This unexpected parallel observation of paternal mitochondrial transmission, undesirable given our intent of creating cybrids, provides a serendipitous experimental model and framework to study the molecular and evolutionary basis of uniparental mitochondrial inheritance.

The Draft Genome of Globodera ellingtonae.

Globodera ellingtonae is a newly described potato cyst nematode (PCN) found in Idaho, Oregon, and Argentina. Here, we present a genome assembly for G. ellingtonae, a relative of the quarantine nematodes G. pallida and G. rostochiensis, produced using data from Illumina and Pacific Biosciences DNA sequencing technologies.

Genome-wide polymorphism and signatures of selection in the symbiotic sea anemone Aiptasia.

BACKGROUND: Coral reef ecosystems are declining in response to global climate change and anthropogenic impacts. Yet patterns of standing genetic variation within cnidarian species, a major determinant of adaptive potential, are virtually unknown at genome-scale resolution. We explore patterns of genome-wide polymorphism and identify candidate loci under selection in the sea anemone Aiptasia, an important laboratory model system for studying the symbiosis between corals and dinoflagellate algae of the genus Symbiodinium. RESULTS: Low coverage genome sequencing revealed large genetic distances among globally widespread lineages, novel candidate targets of selection, and considerably higher heterozygosity than previously reported for Aiptasia. More than 670,000 single nucleotide polymorphisms were identified among 10 Aiptasia individuals including two pairs of genetic clones. Evolutionary relationships based on genome-wide polymorphism supported the current paradigm of a genetically distinct population from the US South Atlantic that harbors diverse Symbiodinium clades. However, anemones from the US South Atlantic demonstrated a striking lack of shared derived polymorphism. Heterozygosity was an important feature shaping nucleotide diversity patterns: at any given SNP site, more than a third of individuals genotyped were heterozygotes, and heterozygosity within individual genomes ranged from 0.37-0.58%. Analysis of nonsynonymous and synonymous sites suggested that highly heterozygous regions are evolving under relaxed purifying selection compared to the rest of the Aiptasia genome. Genes previously identified as having elevated evolutionary rates in Aiptasia compared to other cnidarians were found in our study to be under strong purifying selection within Aiptasia. Candidate targets of selection, including lectins and genes involved in Rho GTPase signalling, were identified based on unusual signatures of nucleotide diversity, Tajima's D, and heterozygosity compared to genome-wide averages. CONCLUSIONS: This study represents the first genome-wide analysis of Tajima's D in a cnidarian. Our results shed light on patterns of intraspecific genome-wide polymorphism in a model for studies of coral-algae symbiosis and present genetic targets for future research on evolutionary and cellular processes in early-diverging metazoans.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers.

BACKGROUND: The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae. RESULTS: The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles. CONCLUSIONS: The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.

Endogenous ROS levels in C. elegans under exogenous stress support revision of oxidative stress theory of life-history tradeoffs.

BACKGROUND: The oxidative stress theory of life-history tradeoffs states that oxidative stress caused by damaging free radicals directly underpins tradeoffs between reproduction and longevity by altering the allocation of energetic resources between these tasks. We test this theory by characterizing the effects of exogenous oxidative insult and its interaction with thermal stress and diet quality on a suite of life-history traits and correlations in Caenorhabditis elegans nematodes. We also quantify demographic aging rates and endogenous reactive oxygen species (ROS) levels in live animals. RESULTS: Our findings indicate a tradeoff between investment in reproduction and antioxidant defense (somatic maintenance) consistent with theoretical predictions, but correlations between standard life-history traits yield little evidence that oxidative stress generates strict tradeoffs. Increasing oxidative insult, however, shows a strong tendency to uncouple positive phenotypic correlations and, in particular, to reduce the correlation between reproduction and lifespan. We also found that mild oxidative insult results in lower levels of endogenous ROS accompanied by hormetic changes in lifespan, demographic aging, and reproduction that disappear in combined-stress treatments--consistent with the oxidative stress theory of aging. CONCLUSIONS: Our findings demonstrate that oxidative stress is a direct contributor to life-history trait variation and that traditional tradeoffs are not necessary to invoke oxidative stress as a mediator of relationships between life-history traits, supporting previous calls for revisions to theory.

Sea anemones possess dynamic mitogenome structures.

A notable feature of hexacoral mitogenomes is the presence of complex self-catalytic group I introns. We investigated mitogenome structural variations and evolutionary mechanisms in actiniarian sea anemones based on the complete mitogenome sequence of the cold-water sea anemone species Urticina eques, Bolocera tuediae, Hormathia digitata and Metridium senile, and two isolates of the sub-tropical Aiptasia pulchella. Whole genome sequencing at 50 times coverage of B. tuediae and H. digitata indicated low mtDNA copy number of per haploid nuclear genome and presence of rare haplotypes. A group I intron inserted in ND5 was found to host essential mitochondrial protein genes in all species, and an additional truncated copy of ND5 in B. tuediae. A second group I intron (inserted in COI) that contained a homing endonuclease gene (HEG) was present in all mtDNA examined. Different variants of HEGs were observed, and included expressed elements fused in-frame with upstream exons and free-standing HEGs embedded within the intron. A notable hallmark of HEGs was a high extent of overlap with ribozyme structural elements; the U. eques HEG overlapped with the entire intron. We reconstructed the evolutionary history of the COI intron from insertion at unoccupied cognate sites, through HEG degradation, to intron loss. We also identified a novel insertion element in U. eques that contained two expressed protein-coding genes. An evolutionary analysis of the sea anemone mtDNA genes revealed higher substitution rates in the HEG and the insertion sequence as compared to the other loci, indicating relaxed selective pressures in these elements. We conclude that sea anemone mitogenomes are surprisingly dynamic in structure despite the economical organization and low sequence mutation rate.

The transcriptional consequences of mutation and natural selection in Caenorhabditis elegans.

The evolutionary importance of gene-expression divergence is unclear: some studies suggest that it is an important mechanism for evolution by natural selection, whereas others claim that most between-species regulatory changes are neutral or nearly neutral. We examined global transcriptional divergence patterns in a set of Caenorhabditis elegans mutation-accumulation lines and natural isolate lines to provide insights into the evolutionary importance of transcriptional variation and to discriminate between the forces of mutation and natural selection in shaping the evolution of gene expression. We detected the effects of selection on transcriptional divergence patterns and characterized them with respect to coexpressed gene sets, chromosomal clustering of expression changes and functional gene categories. We directly compared observed transcriptional variation patterns in the mutation-accumulation and natural isolate lines to a neutral model of transcriptome evolution to show that strong stabilizing selection dominates the evolution of transcriptional change for thousands of C. elegans expressed sequences.

Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans.

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.

The worm in the world and the world in the worm.

Caenorhabditis elegans is a preeminent model organism, but the natural ecology of this nematode has been elusive. A four-year survey of French orchards published in BMC Biology reveals thriving populations of C. elegans (and Caenorhabditis briggsae) in rotting fruit and plant stems. Rather than being simply a 'soil nematode', C. elegans appears to be a 'plant-rot nematode'. These studies signal a growing interest in the integrated genomics and ecology of these tractable animals.

High mitogenomic evolutionary rates and time dependency.

Using entire modern and ancient mitochondrial genomes of Adélie penguins (Pygoscelis adeliae) that are up to 44000 years old, we show that the rates of evolution of the mitochondrial genome are two to six times greater than those estimated from phylogenetic comparisons. Although the rate of evolution at constrained sites, including nonsynonymous positions and RNAs, varies more than twofold with time (between shallow and deep nodes), the rate of evolution at synonymous sites remains the same. The time-independent neutral evolutionary rates reported here would be useful for the study of recent evolutionary events.

A genome-wide view of Caenorhabditis elegans base-substitution mutation processes.

Knowledge of mutation processes is central to understanding virtually all evolutionary phenomena and the underlying nature of genetic disorders and cancers. However, the limitations of standard molecular mutation detection methods have historically precluded a genome-wide understanding of mutation rates and spectra in the nuclear genomes of multicellular organisms. We applied two high-throughput DNA sequencing technologies to identify and characterize hundreds of spontaneously arising base-substitution mutations in 10 Caenorhabditis elegans mutation-accumulation (MA)-line nuclear genomes. C. elegans mutation rate estimates were similar to previous calculations based on smaller numbers of mutations. Mutations were distributed uniformly within and among chromosomes and were not associated with recombination rate variation in the MA lines, suggesting that intragenomic variation in genetic hitchhiking and/or background selection are primarily responsible for the chromosomal distribution patterns of polymorphic nucleotides in C. elegans natural populations. A strong mutational bias from G/C to A/T nucleotides was detected in the MA lines, implicating oxidative DNA damage as a major endogenous mutagenic force in C. elegans. The observed mutational bias also suggests that the C. elegans nuclear genome cannot be at equilibrium because of mutation alone. Transversions dominate the spectrum of spontaneous mutations observed here, whereas transitions dominate patterns of allegedly neutral polymorphism in natural populations of C. elegans and many other animal species; this observation challenges the assumption that natural patterns of molecular variation in noncoding regions of the nuclear genome accurately reflect underlying mutation processes.

Natural variation in life history and aging phenotypes is associated with mitochondrial DNA deletion frequency in Caenorhabditis briggsae.

BACKGROUND: Mutations that impair mitochondrial functioning are associated with a variety of metabolic and age-related disorders. A barrier to rigorous tests of the role of mitochondrial dysfunction in aging processes has been the lack of model systems with relevant, naturally occurring mitochondrial genetic variation. Toward the goal of developing such a model system, we studied natural variation in life history, metabolic, and aging phenotypes as it relates to levels of a naturally-occurring heteroplasmic mitochondrial ND5 deletion recently discovered to segregate among wild populations of the soil nematode, Caenorhabditis briggsae. The normal product of ND5 is a central component of the mitochondrial electron transport chain and integral to cellular energy metabolism. RESULTS: We quantified significant variation among C. briggsae isolates for all phenotypes measured, only some of which was statistically associated with isolate-specific ND5 deletion frequency. We found that fecundity-related traits and pharyngeal pumping rate were strongly inversely related to ND5 deletion level and that C. briggsae isolates with high ND5 deletion levels experienced a tradeoff between early fecundity and lifespan. Conversely, oxidative stress resistance was only weakly associated with ND5 deletion level while ATP content was unrelated to deletion level. Finally, mean levels of reactive oxygen species measured in vivo showed a significant non-linear relationship with ND5 deletion level, a pattern that may be driven by among-isolate variation in antioxidant or other compensatory mechanisms. CONCLUSIONS: Our findings suggest that the ND5 deletion may adversely affect fitness and mitochondrial functioning while promoting aging in natural populations, and help to further establish this species as a useful model for explicit tests of hypotheses in aging biology and mitochondrial genetics.

Evolution of Caenorhabditis mitochondrial genome pseudogenes and Caenorhabditis briggsae natural isolates.

Although most metazoan mitochondrial genomes are highly streamlined and encode little noncoding DNA outside of the "AT" region, the accumulation of mitochondrial pseudogenes and other types of noncoding DNA has been observed in a growing number of animal groups. The nematode species Caenorhabditis briggsae harbors two mitochondrial DNA (mtDNA) pseudogenes, named Psinad5-1 and Psinad5-2, presumably derived from the nad5 protein-coding gene. Here, we provide an in-depth analysis of mtDNA pseudogene evolution in C. briggsae natural isolates and related Caenorhabditis species. Mapping the observed presence and absence of the pseudogenes onto phylogenies suggests that Psinad5-1 originated in the ancestor to C. briggsae and its recently discovered outcrossing relative species Caenorhabditis sp. 5 and Caenorhabditis sp. 9. However, Psinad5-1 was not detected in Caenorhabditis sp. 9 natural isolates, suggesting a lineage-specific loss of this pseudogene in this species. Our results corroborated the previous finding that Psinad5-2 originated within C. briggsae. The observed pattern of mitochondrial pseudogene gain and loss in Caenorhabditis was inconsistent with predictions of the tandem duplication-random loss model of mitochondrial genome evolution and suggests that intralineage recombination-like mechanisms might play a major role in Caenorhabditis mtDNA evolution. Natural variation was analyzed at the pseudogenes and flanking mtDNA sequences in 141 geographically diverse C. briggsae natural isolates. Although phylogenetic analysis placed the majority of isolates into the three previously established major intraspecific clades of C. briggsae, two new and unexpected haplotypes fell outside of these conventional groupings. Psinad5-2 copy number variation was observed among C. briggsae isolates collected from the same geographic site. Patterns of nucleotide diversity were analyzed in Psinad5-1 and Psinad5-2, and confidence intervals were found to overlap values from synonymous sites in protein-coding genes, consistent with neutral expectations. Our findings provide new insights into the mode and tempo of mitochondrial genome and pseudogene evolution both within and between Caenorhabditis nematode species.

Genomic Analyses of Globodera pallida, A Quarantine Agricultural Pathogen in Idaho.

Globodera pallida is among the most significant plant-parasitic nematodes worldwide, causing major damage to potato production. Since it was discovered in Idaho in 2006, eradication efforts have aimed to contain and eradicate G. pallida through phytosanitary action and soil fumigation. In this study, we investigated genome-wide patterns of G. pallida genetic variation across Idaho fields to evaluate whether the infestation resulted from a single or multiple introduction(s) and to investigate potential evolutionary responses since the time of infestation. A total of 53 G. pallida samples (~1,042,000 individuals) were collected and analyzed, representing five different fields in Idaho, a greenhouse population, and a field in Scotland that was used for external comparison. According to genome-wide allele frequency and fixation index (Fst) analyses, most of the genetic variation was shared among the G. pallida populations in Idaho fields pre-fumigation, indicating that the infestation likely resulted from a single introduction. Temporal patterns of genome-wide polymorphisms involving (1) pre-fumigation field samples collected in 2007 and 2014 and (2) pre- and post-fumigation samples revealed nucleotide variants (SNPs, single-nucleotide polymorphisms) with significantly differentiated allele frequencies indicating genetic differentiation. This study provides insights into the genetic origins and adaptive potential of G. pallida invading new environments.