The Global Invertebrate Genomics Alliance (GIGA): developing community resources to study diverse invertebrate genomes.

Over 95% of all metazoan (animal) species comprise the "invertebrates," but very few genomes from these organisms have been sequenced. We have, therefore, formed a "Global Invertebrate Genomics Alliance" (GIGA). Our intent is to build a collaborative network of diverse scientists to tackle major challenges (e.g., species selection, sample collection and storage, sequence assembly, annotation, analytical tools) associated with genome/transcriptome sequencing across a large taxonomic spectrum. We aim to promote standards that will facilitate comparative approaches to invertebrate genomics and collaborations across the international scientific community. Candidate study taxa include species from Porifera, Ctenophora, Cnidaria, Placozoa, Mollusca, Arthropoda, Echinodermata, Annelida, Bryozoa, and Platyhelminthes, among others. GIGA will target 7000 noninsect/nonnematode species, with an emphasis on marine taxa because of the unrivaled phyletic diversity in the oceans. Priorities for selecting invertebrates for sequencing will include, but are not restricted to, their phylogenetic placement; relevance to organismal, ecological, and conservation research; and their importance to fisheries and human health. We highlight benefits of sequencing both whole genomes (DNA) and transcriptomes and also suggest policies for genomic-level data access and sharing based on transparency and inclusiveness. The GIGA Web site (http://giga.nova.edu) has been launched to facilitate this collaborative venture.

Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing.

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the "rare biosphere" than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.

Ironing out the wrinkles in the rare biosphere through improved OTU clustering.

Deep sequencing of PCR amplicon libraries facilitates the detection of low-abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low-frequency, error-prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete-linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single-linkage preclustering methodology followed by an average-linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30-60% but does not reduce the fraction of OTUs in long-tailed rank abundance curves that defines the rare biosphere.

Genetic exchange within and between assemblages of Giardia duodenalis.

Meiotic sex evolved early in the history of eukaryotes. Giardia duodenalis (syn. Giardia lamblia, Giardia intestinalis), a parasitic protist belonging to an early diverging lineage of eukaryotes, shows no cytological or physiological evidence of meiotic or sexual processes. Recent molecular analyses challenge the idea that G. duodenalis is a strictly clonal organism by providing evidence of recombination between homologous chromosomes within one subgroup (Assemblage A) of this species as well as genetic transfer from one subgroup to another (Assemblage A-B). Because recombination is not well documented and because it is not known whether the observed inter-assemblage transfer represents true reciprocal genetic exchange or a non-sexual process, we analyzed genic sequences from all major subgroups (Assemblages A-G) of this species. For all assemblages, we detected molecular signatures consistent with meiotic sex or genetic exchange, including low levels of heterozygosity, as indicated by allelic sequence divergence within isolates, and intra- and inter-assemblage recombination. The identification of recombination between assemblages suggests a shared gene pool and calls into question whether it is appropriate to divide the genetically distinct assemblages of G. duodenalis into a species complex.

Molecular characterization of Giardia intestinalis haplotypes in marine animals: variation and zoonotic potential.

Giardia intestinalis is a microbial eukaryotic parasite that causes diarrheal disease in humans and other vertebrates worldwide. The negative effect on quality of life and economics caused by G. intestinalis may be increased by its potential status as a zoonosis, or a disease that can be transmitted from animals to humans. The zoonotic potential of G. intestinalis has been implied for over 2 decades, with human-infecting genotypes (belonging to the 2 major subgroups, Assemblages A and B) occurring in wildlife and domesticated animals. There are recent reports of G. intestinalis in shellfish, seals, sea lions and whales, suggesting that marine animals are also potential reservoirs of human disease. However, the prevalence, genetic diversity and effect of G. intestinalis in marine environments and the role that marine animals play in transmission of this parasite to humans are relatively unexplored. Here, we provide the first thorough molecular characterization of G. intestinalis in marine vertebrates. Using a multi-locus sequencing approach, we identify human-infecting G. intestinalis haplotypes of both Assemblages A and B in the fecal material of dolphins, porpoises, seals, herring gulls Larus argentatus, common eiders Somateria mollissima and a thresher shark Alopias vulpinus. Our results indicate that G. intestinalis is prevalent in marine ecosystems, and a wide range of marine hosts capable of harboring zoonotic forms of this parasite exist. The presence of G. intestinalis in marine ecosystems raises concerns about how this disease might be transmitted among different host species.

Long term seasonal dynamics of synechococcus population structure in the gulf of aqaba, northern red sea.

Spatial patterns of marine Synechococcus diversity across ocean domains have been reported on extensively. However, much less is known of seasonal and multiannual patterns of change in Synechococcus community composition. Here we report on the genotypic diversity of Synechococcus populations in the Gulf of Aqaba, Northern Red Sea, over seven annual cycles of deep mixing and stabile stratification, using ntcA as a phylogenetic marker. Synechococcus clone libraries were dominated by clade II and XII genotypes and a total of eight different clades were identified. Inclusion of ntcA sequences from the Global Ocean Sampling database in our analyses identified members of clade XII from beyond the Gulf of Aqaba, extending its known distribution. Most of the Synechococcus diversity was attributed to members of clade II during the spring bloom, while clade III contributed significantly to diversity during summer stratification. Clade XII diversity was most prevalent in fall and winter. Clade abundances were estimated from pyrosequencing of the V6 hypervariable region of 16S rRNA. Members of clade II dominated Synechococcus communities throughout the year, whereas the less frequent genotypes showed a pattern of seasonal succession. Based on the prevailing nutritional conditions we observed that clade I members thrive at higher nutrient concentrations during winter mixing. Clades V, VI and X became apparent during the transition periods between mixing and stratification. Clade III became prominent during sumeer stratification. We propose that members of clades V, VI, and X, and clade III are Synechococcus ecotypes that are adapted to intermediate and low nutrient levels respectively. This is the first time that molecular analyses have correlated population dynamics of Synechococcus genotypes with temporal fluctuations in nutrient regimes. Since these Synechococcus genotypes are routinely observed in the Gulf of Aqaba we suggest that seasonal fluctuations in nutrient levels create temporal niches that sustain their coexistence.

The identification of a new Giardia duodenalis assemblage in marine vertebrates and a preliminary analysis of G. duodenalis population biology in marine systems.

Giardia duodenalis is an intestinal parasite of many vertebrates. The presence of G. duodenalis in the marine environment due to anthropogenic and wildlife activity is well documented, including the contributions from untreated sewage and storm water, agricultural run-off and droppings from terrestrial animals. Recently, studies have detected this protistan parasite in the faeces of marine vertebrates such as whales, dolphins, seals and shore birds. To explore the population biology of G. duodenalis in marine life, we determined the prevalence of G. duodenalis in two species of seal (Halichoerus grypus, Phoca vitulina vitulina and Phoca vitulina richardsi) from the east and west coasts of the USA, sequenced two loci from G. duodenalis-positive samples to assess molecular diversity and examined G. duodenalis distribution amongst these seals and other marine vertebrates along the east coast. We found a significant difference in the presence of G. duodenalis between east and west coast seal species. Only the zoonotic lineages of G. duodenalis, Assemblages A and B and a novel lineage, which we designated as Assemblage H, were identified in marine vertebrates. Assemblages A and B are broadly distributed geographically and show a lack of host specificity. Only grey seal (Halichoerus grypus) samples and one gull sample (Larus argentatus) from a northern location of Cape Cod, Massachusetts, USA, showed the presence of Assemblage H haplotypes; only one other study of harbour seals from the Puget Sound region of Washington, USA previously recorded the presence of an Assemblage H haplotype. Assemblage H sequences form a monophyletic clade that appears as divergent from the other seven Assemblages of G. duodenalis as these assemblages are from each other. The discovery of a previously uncharacterised lineage of G. duodenalis suggests that this parasite has more genetic diversity and perhaps a larger host range than previously believed.

Analysis of expressed sequence tags of the cyclically parthenogenetic rotifer Brachionus plicatilis.

BACKGROUND: Rotifers are among the most common non-arthropod animals and are the most experimentally tractable members of the basal assemblage of metazoan phyla known as Gnathifera. The monogonont rotifer Brachionus plicatilis is a developing model system for ecotoxicology, aquatic ecology, cryptic speciation, and the evolution of sex, and is an important food source for finfish aquaculture. However, basic knowledge of the genome and transcriptome of any rotifer species has been lacking. METHODOLOGY/PRINCIPAL FINDINGS: We generated and partially sequenced a cDNA library from B. plicatilis and constructed a database of over 2300 expressed sequence tags corresponding to more than 450 transcripts. About 20% of the transcripts had no significant similarity to database sequences by BLAST; most of these contained open reading frames of significant length but few had recognized Pfam motifs. Sixteen transcripts accounted for 25% of the ESTs; four of these had no significant similarity to BLAST or Pfam databases. Putative up- and downstream untranslated regions are relatively short and AT rich. In contrast to bdelloid rotifers, there was no evidence of a conserved trans-spliced leader sequence among the transcripts and most genes were single-copy. CONCLUSIONS/SIGNIFICANCE: Despite the small size of this EST project it revealed several important features of the rotifer transcriptome and of individual monogonont genes. Because there is little genomic data for Gnathifera, the transcripts we found with no known function may represent genes that are species-, class-, phylum- or even superphylum-specific; the fact that some are among the most highly expressed indicates their importance. The absence of trans-spliced leader exons in this monogonont species contrasts with their abundance in bdelloid rotifers and indicates that the presence of this phenomenon can vary at the subphylum level. Our EST database provides a relatively large quantity of transcript-level data for B. plicatilis, and more generally of rotifers and other gnathiferan phyla, and can be browsed and searched at gmod.mbl.edu.

Accuracy and quality of massively parallel DNA pyrosequencing.

BACKGROUND: Massively parallel pyrosequencing systems have increased the efficiency of DNA sequencing, although the published per-base accuracy of a Roche GS20 is only 96%. In genome projects, highly redundant consensus assemblies can compensate for sequencing errors. In contrast, studies of microbial diversity that catalogue differences between PCR amplicons of ribosomal RNA genes (rDNA) or other conserved gene families cannot take advantage of consensus assemblies to detect and minimize incorrect base calls. RESULTS: We performed an empirical study of the per-base error rate for the Roche GS20 system using sequences of the V6 hypervariable region from cloned microbial ribosomal DNA (tag sequencing). We calculated a 99.5% accuracy rate in unassembled sequences, and identified several factors that can be used to remove a small percentage of low-quality reads, improving the accuracy to 99.75% or better. CONCLUSION: By using objective criteria to eliminate low quality data, the quality of individual GS20 sequence reads in molecular ecological applications can surpass the accuracy of traditional capillary methods.

Transmission dynamics of Cryptosporidium infection in a natural population of non-human primates at Polonnaruwa, Sri Lanka.

Infections from Cryptosporidium parvum are of interest not only to public health, but also to wildlife conservation, particularly when humans and livestock encroach on nature and thereby increase the risk of cross-species transmissions. To clarify this risk, we used polymerase chain reaction to examine the hypervariable region of the C. parvum 18S rRNA gene in feces from three monkey species. Samples were isolated from regions where disease transmission between monkeys, livestock, and humans was likely (soiled habitat) or unlikely (clean habitat). Monkey individuals, their social groups, and different species shared multiple genotypes/isolates of C. parvum. Ecological and molecular analyses suggested that Cryptosporidium infection among Toque macaques in soiled habitats was mainly the bovine genotype C. parvum. Monkeys inhabiting clean habitat, particularly gray and purple-faced langurs, lacked Cryptosporidium species/types associated with bovines. Livestock apparently was a main source of infection for wild primates.