De Novo Transcriptome Analysis of the Common New Zealand Stick Insect Clitarchus hookeri (Phasmatodea) Reveals Genes Involved in Olfaction, Digestion and Sexual Reproduction.

Phasmatodea, more commonly known as stick insects, have been poorly studied at the molecular level for several key traits, such as components of the sensory system and regulators of reproduction and development, impeding a deeper understanding of their functional biology. Here, we employ de novo transcriptome analysis to identify genes with primary functions related to female odour reception, digestion, and male sexual traits in the New Zealand common stick insect Clitarchus hookeri (White). The female olfactory gene repertoire revealed ten odorant binding proteins with three recently duplicated, 12 chemosensory proteins, 16 odorant receptors, and 17 ionotropic receptors. The majority of these olfactory genes were over-expressed in female antennae and have the inferred function of odorant reception. Others that were predominantly expressed in male terminalia (n = 3) and female midgut (n = 1) suggest they have a role in sexual reproduction and digestion, respectively. Over-represented transcripts in the midgut were enriched with digestive enzyme gene families. Clitarchus hookeri is likely to harbour nine members of an endogenous cellulase family (glycoside hydrolase family 9), two of which appear to be specific to the C. hookeri lineage. All of these cellulase sequences fall into four main phasmid clades and show gene duplication events occurred early in the diversification of Phasmatodea. In addition, C. hookeri genome is likely to express γ-proteobacteria pectinase transcripts that have recently been shown to be the result of horizontal transfer. We also predicted 711 male terminalia-enriched transcripts that are candidate accessory gland proteins, 28 of which were annotated to have molecular functions of peptidase activity and peptidase inhibitor activity, two groups being widely reported to regulate female reproduction through proteolytic cascades. Our study has yielded new insights into the genetic basis of odour detection, nutrient digestion, and male sexual traits in stick insects. The C. hookeri reference transcriptome, together with identified gene families, provides a comprehensive resource for studying the evolution of sensory perception, digestive systems, and reproductive success in phasmids.

A new virus discovered by immunocapture of double-stranded RNA, a rapid method for virus enrichment in metagenomic studies.

Next-generation sequencing technologies enable the rapid identification of viral infection of diseased organisms. However, despite a consistent decrease in sequencing costs, it is difficult to justify their use in large-scale surveys without a virus sequence enrichment technique. As the majority of plant viruses have an RNA genome, a common approach is to extract the double-stranded RNA (dsRNA) replicative form, to enrich the replicating virus genetic material over the host background. The traditional dsRNA extraction is time-consuming and labour-intensive. We present an alternative method to enrich dsRNA from plant extracts using anti-dsRNA monoclonal antibodies in a pull-down assay. The extracted dsRNA can be amplified by reverse transcriptase-polymerase chain reaction and sequenced by next-generation sequencing. In our study, we have selected three distinct plant hosts: Maori potato (Solanum tuberosum), rengarenga (Arthropodium cirratum) and broadleaved dock (Rumex obtusifolius) representing a cultivated crop, a New Zealand-native ornamental plant and a weed, respectively. Of the sequence data obtained, 31-74% of the reads were of viral origin, and we identified five viruses including Potato virus Y and Potato virus S in potato; Turnip mosaic virus in rengarenga (a new host record); and in the dock sample Cherry leaf roll virus and a novel virus belonging to the genus Macluravirus. We believe that this new assay represents a significant opportunity to upscale virus ecology studies from environmental, primary industry and/or medical samples.

Evaluating a multigene environmental DNA approach for biodiversity assessment.

BACKGROUND: There is an increasing demand for rapid biodiversity assessment tools that have a broad taxonomic coverage. Here we evaluate a suite of environmental DNA (eDNA) markers coupled with next generation sequencing (NGS) that span the tree of life, comparing them with traditional biodiversity monitoring tools within ten 20×20 meter plots along a 700 meter elevational gradient. RESULTS: From six eDNA datasets (one from each of 16S, 18S, ITS, trnL and two from COI) we identified sequences from 109 NCBI taxonomy-defined phyla or equivalent, ranging from 31 to 60 for a given eDNA marker. Estimates of alpha and gamma diversity were sensitive to the number of sequence reads, whereas beta diversity estimates were less sensitive. The average within-plot beta diversity was lower than between plots for all markers. The soil beta diversity of COI and 18S markers showed the strongest response to the elevational variation of the eDNA markers (COI: r=0.49, p<0.001; 18S: r=0.48, p<0.001). Furthermore pairwise beta diversities for these two markers were strongly correlated with those calculated from traditional vegetation and invertebrate biodiversity measures. CONCLUSIONS: Using a soil-based eDNA approach, we demonstrate that standard phylogenetic markers are capable of recovering sequences from a broad diversity of eukaryotes, in addition to prokaryotes by 16S. The COI and 18S eDNA markers are the best proxies for aboveground biodiversity based on the high correlation between the pairwise beta diversities of these markers and those obtained using traditional methods.

Integration over song classification replicates: song variant analysis in the hihi.

Human expert analyses are commonly used in bioacoustic studies and can potentially limit the reproducibility of these results. In this paper, a machine learning method is presented to statistically classify avian vocalizations. Automated approaches were applied to isolate bird songs from long field recordings, assess song similarities, and classify songs into distinct variants. Because no positive controls were available to assess the true classification of variants, multiple replicates of automatic classification of song variants were analyzed to investigate clustering uncertainty. The automatic classifications were more similar to the expert classifications than expected by chance. Application of these methods demonstrated the presence of discrete song variants in an island population of the New Zealand hihi (Notiomystis cincta). The geographic patterns of song variation were then revealed by integrating over classification replicates. Because this automated approach considers variation in song variant classification, it reduces potential human bias and facilitates the reproducibility of the results.

The incidence of species-level paraphyly in animals: a re-assessment.

Species-level paraphyly was found by Funk and Omland (2003) to occur in 23% of animal species on the basis of a meta-analysis of published mitochondrial gene trees. Given the potential for bias in the selection of study organisms and the subsequent publication of their gene trees, I re-estimated the incidence of paraphyly in an independent dataset of publicly accessible COI sequences from the Barcode of Life Data System. Among 7368 animal species represented by two or more sequences, 19% were paraphyletic, slightly less than in the previous study. Rates within major taxonomic groups mirrored, but were slightly lower than, that observed earlier. Tests were made for operational factors that could inflate, and sampling effects that could underestimate, the rate of paraphyly. Overall the previous findings are confirmed. The observed incidence suggests that on average animal species diverged 2-3 Ne generations in the past, far short of the predicted 5 Ne generations required for complete monophyly.

Pyrosequencing reveals regional differences in fruit-associated fungal communities.

We know relatively little of the distribution of microbial communities generally. Significant work has examined a range of bacterial communities, but the distribution of microbial eukaryotes is less well characterized. Humans have an ancient association with grape vines (Vitis vinifera) and have been making wine since the dawn of civilization, and fungi drive this natural process. While the molecular biology of certain fungi naturally associated with vines and wines is well characterized, complementary investigations into the ecology of fungi associated with fruiting plants is largely lacking. DNA sequencing technologies allow the direct estimation of microbial diversity from a given sample, avoiding culture-based biases. Here, we use deep community pyrosequencing approaches, targeted at the 26S rRNA gene, to examine the richness and composition of fungal communities associated with grapevines and test for geographical community structure among four major regions in New Zealand (NZ). We find over 200 taxa using this approach, which is 10-fold more than previously recovered using culture-based methods. Our analyses allow us to reject the null hypothesis of homogeneity in fungal species richness and community composition across NZ and reveal significant differences between major areas.

Selection pressure on the hepatitis B virus pre-S/S and P open reading frames in Tongan subjects with a chronic hepatitis B virus infection.

Identification of the full repertoire of hepatitis B virus (HBV) peptides that are presented to CD8+ T cells by common HLA class I alleles will be useful for designing immunotherapies for chronic hepatitis B. One hundred and seventy five cloned sequences containing the pre-S/S and P open reading frames (ORF) of the HBV were obtained from serum HBV-DNA of HBeAg-positive (n=4) and HBeAg-negative (inactive healthy carriers (IHC), n=16) Tongan subjects with an inactive chronic HBV infection. In addition, 34 and 32 sequences were obtained 5.2±1.4 (mean±SD) years apart from eight subjects. PAML was used to identify codons in the pre-S/S and P ORFs that were under positive selection pressure (ω>1). The number of non-synonymous substitutions in these codons was compared in IHC who were homozygous for either HLA-B∗4001 (n=9) or HLA-B*5602 (n=7), and who were either positive (n=6) or negative (n=10) for HLA-A*02. 34 codons in the pre-S/S and 11 codons in the P ORFs were under positive selection pressure. There was a higher number of non-synonymous substitutions in these codons in HBeAg-negative versus HBeAg-positive subjects in the P (p=0.02) but not the pre-S/S (p=0.64) ORF. There was no association between any HLA class I allele and non-synonymous substitutions in these codons. There was no increase in positive selection pressure on the pre-S/S and P ORFs with time. In conclusion, we could not find HLA class I-restricted selection pressure on any pre-S/S or P ORF amino acid; raising the possibility that peptide-based immunotherapies for chronic hepatitis B may not require peptides from these ORFs.

Species Delimitation--a Geneious plugin for the exploration of species boundaries.

Species Delimitation is a plugin to the Geneious software to support the exploration of species boundaries in a gene tree. The user assigns taxa to putative species and the plugin computes statistics relating to the probability of the observed monophyly or exclusivity having occurred by chance in a coalescent process. It also assesses the within and between species genetic distances to infer the probability with which members of a putative species might be identified successfully with tree-based methods.

Selective sweeps at the organophosphorus insecticide resistance locus, Rop-1, have affected variation across and beyond the α-esterase gene cluster in the Australian sheep blowfly, Lucilia cuprina.

A major theoretical consequence of selection at a locus is the genetic hitchhiking of linked sites (selective sweep). The extent of hitchhiking around a gene is related to the strength of selection and the rate of recombination, with its impact diminishing with distance from the selected site. At the Rop-1 locus of the sheep blowfly, Lucilia cuprina, polymorphisms at two different sites within the LcαE7 gene encode forms of the protein that confer organophosphorus insecticide resistance. To assess the impact of selection at these two sites on variation around LcαE7, we sequenced regions within six other genes along chromosome IV across isogenic (IV) strains of L. cuprina. High levels of linkage disequilibrium, characterized by low haplotype number (K) and diversity (H), and significant R(2) values were observed for two genes, LcαE1 and LcαE10, both members of the same α-esterase gene cluster as LcαE7. A significant R(2) value was also observed for a gene predicted to be the next closest to LcαE7, AL03, but not for any of the other genes, LcRpL13a, Lcdsx, or LcAce. Skews in the site frequency spectra toward high-frequency variants were significant for LcαE1 (Fay and Wu's H = -2.91), LcαE10 (H = -1.85), and Lcdsx (H = -2.00). Since the selective sweeps, two forms of likely returning variation were observed, including variation in microsatellites in an intron of LcαE10 and a recombination event between LcαE7 and LcαE10. These data suggest that two incomplete soft sweeps have occurred at LcαE7 that have significantly affected variation across, and beyond, the α-esterase gene cluster of L. cuprina. The speed and impact of these selective sweeps on surrounding genomic variation and the ability of L. cuprina to respond to future environmental challenges are discussed.

Evolutionary relationships and divergence times among the native rats of Australia.

BACKGROUND: The genus Rattus is highly speciose and has a complex taxonomy that is not fully resolved. As shown previously there are two major groups within the genus, an Asian and an Australo-Papuan group. This study focuses on the Australo-Papuan group and particularly on the Australian rats. There are uncertainties regarding the number of species within the group and the relationships among them. We analysed 16 mitochondrial genomes, including seven novel genomes from six species, to help elucidate the evolutionary history of the Australian rats. We also demonstrate, from a larger dataset, the usefulness of short regions of the mitochondrial genome in identifying these rats at the species level. RESULTS: Analyses of 16 mitochondrial genomes representing species sampled from Australo-Papuan and Asian clades of Rattus indicate divergence of these two groups ~2.7 million years ago (Mya). Subsequent diversification of at least 4 lineages within the Australo-Papuan clade was rapid and occurred over the period from ~ 0.9-1.7 Mya, a finding that explains the difficulty in resolving some relationships within this clade. Phylogenetic analyses of our 126 taxon, but shorter sequence (1952 nucleotides long), Rattus database generally give well supported species clades. CONCLUSIONS: Our whole mitochondrial genome analyses are concordant with a taxonomic division that places the native Australian rats into the Rattus fuscipes species group. We suggest the following order of divergence of the Australian species. R. fuscipes is the oldest lineage among the Australian rats and is not part of a New Guinean radiation. R. lutreolus is also within this Australian clade and shallower than R. tunneyi while the R. sordidus group is the shallowest lineage in the clade. The divergences within the R. sordidus and R. leucopus lineages occurring about half a million years ago support the hypotheses of more recent interchanges of rats between Australia and New Guinea. While problematic for inference of deeper divergences, we report that the analysis of shorter mitochondrial sequences is very useful for species identification in rats.

Latitude, elevation and the tempo of molecular evolution in mammals.

Faster rates of microevolution have been recorded for plants and marine foraminifera occupying warmer low latitude environments relative to those occurring at higher latitudes. By contrast, because this rate heterogeneity has been attributed to a relationship between thermal habit and mutagenesis via a body temperature linkage, it has been assumed that microevolution in mammals should not also vary systematically with environmental temperature. However, this assumption has not previously been empirically examined. In this study, we tested for a thermally mediated influence on the tempo of microevolution among mammals using a comprehensive global dataset that included 260 mammal species, from 10 orders and 29 families. In contrast to theoretical predictions, we found that substitution rates in the cytochrome b gene have been substantially faster for species living in warmer latitudes and elevations relative to sister species living in cooler habitats. These results could not be attributed to factors otherwise thought to influence rates of microevolution, such as body mass differentials or genetic drift. Instead, the results indicate that the tempo of microevolution among mammals is either responding directly to the thermal environment or indirectly via an ecological mechanism such as the 'Red Queen' effect.

Slower tempo of microevolution in island birds: implications for conservation biology.

Whether microevolution on small islands differs from that on larger landmasses is a key question in biology with substantial implications for species conservation. However, due to the difficulties faced in producing adequately replicated samples and in controlling for confounding variables, prior attempts to examine evolutionary questions relating to habitat area and population size have produced equivocal results. Here we show, using experimental design criteria that reduce the potential for such confounding, that bird species on larger landmasses have higher rates of molecular evolution. The study involves a global dataset of 48 independent contrasts for the cytochrome b gene encompassing all possible paired sister species comparisons (from seven orders and 17 families) that were available at the time of dataset assembly. A more rapid evolutionary tempo in larger areas has important ramifications for biodiversity conservation because it indicates a new imperative, beyond that of simply maintaining preexisting genetic diversity, for securing large areas for threatened species. This result suggests that the trend of confining species to limited refugia is likely to be slowing the tempo of microevolution. That effect might constrain the potential for adaptive shifts in response to changing environments such as those associated with global warming.

Dating of divergences within the Rattus genus phylogeny using whole mitochondrial genomes.

The timing and order of divergences within the genus Rattus have, to date, been quite speculative. In order to address these important issues we sequenced six new whole mitochondrial genomes from wild-caught specimens from four species, Rattus exulans, Rattus praetor, Rattus rattus and Rattus tanezumi. The only rat whole mitochondrial genomes available previously were all from Rattus norvegicus specimens. Our phylogenetic and dating analyses place the deepest divergence within Rattus at approximately 3.5 million years ago (Mya). This divergence separates the New Guinean endemic R. praetor lineage from the Asian lineages. Within the Asian/Island Southeast Asian clade R. norvegicus diverged earliest at approximately 2.9Mya. R. exulans and the ancestor of the sister species R. rattus and R. tanezumi subsequently diverged at approximately 2.2Mya, with R. rattus and R. tanezumi separating as recently as approximately 0.4Mya. Our results give both a better resolved species divergence order and diversification dates within Rattus than previous studies.

Unsupervised bird song syllable classification using evolving neural networks.

Evolution of bird vocalizations is subjected to selection pressure related to their functions. Passerine bird songs are also under a neutral model of evolution because of the learning process supporting their transmission; thus they contain signals of individual, population, and species relationships. In order to retrieve this information, large amounts of data need to be processed. From vocalization recordings, songs are extracted and encoded as sequences of syllables before being compared. Encoding songs in such a way can be done either by ear and spectrogram visual analysis or by specific algorithms permitting reproducible studies. Here, a specific automatic method is presented to compute a syllable distance measure allowing an unsupervised classification of song syllables. Results obtained from the encoding of White-crowned Sparrow (Zonotrichia leucophrys pugetensis) songs are compared to human-based analysis.

Testing the reliability of genetic methods of species identification via simulation.

Although genetic methods of species identification, especially DNA barcoding, are strongly debated, tests of these methods have been restricted to a few empirical cases for pragmatic reasons. Here we use simulation to test the performance of methods based on sequence comparison (BLAST and genetic distance) and tree topology over a wide range of evolutionary scenarios. Sequences were simulated on a range of gene trees spanning almost three orders of magnitude in tree depth and in coalescent depth; that is, deep or shallow trees with deep or shallow coalescences. When the query's conspecific sequences were included in the reference alignment, the rate of positive identification was related to the degree to which different species were genetically differentiated. The BLAST, distance, and liberal tree-based methods returned higher rates of correct identification than did the strict tree-based requirement that the query was within, but not sister to, a single-species clade. Under this more conservative approach, ambiguous outcomes occurred in inverse proportion to the number of reference sequences per species. When the query's conspecific sequences were not in the reference alignment, only the strict tree-based approach was relatively immune to making false-positive identifications. Thresholds affected the rates at which false-positive identifications were made when the query's species was unrepresented in the reference alignment but did not otherwise influence outcomes. A conservative approach using the strict tree-based method should be used initially in large-scale identification systems, with effort made to maximize sequence sampling within species. Once the genetic variation within a taxonomic group is well characterized and the taxonomy resolved, then the choice of method used should be dictated by considerations of computational efficiency. The requirement for extensive genetic sampling may render these techniques inappropriate in some circumstances.

Recombination in feline lentiviral genomes during experimental cross-species infection.

Domestic cats develop an asymptomatic, productive infection with a feline immunodeficiency virus (PLV) derived from a naturally infected cougar (P. concolor). We previously demonstrated that there are extensive G to A substitutions, characteristic of host cytidine deaminase editing, and positive selection on reverse transcriptase in the PLV genome during this cross-species infection. In this study, we evaluated full-length viral genomes from each of four cats infected with PLV to determine if viral recombination occurred during this single source infection. Recombination rates were measurable in three of the four infected cats. In two of these animals, a single site in reverse transcriptase was under positive selection and there was significant topological incongruence among individual genes in the 3' half of the genomes. The break point was proximate to a splice site used for accessory gene expression. Our data indicate that recombination can facilitate lentivirus persistence in unfavorable environments such as a new host species.

Sources of variation in ancestral sequence reconstruction for HIV-1 envelope genes.

We characterized the variation in the reconstructed ancestor of 118 HIV-1 envelope gene sequences arising from the methods used for (a) estimating and (b) rooting the phylogenetic tree, and (c) reconstructing the ancestor on that tree, from (d) the sequence format, and from (e) the number of input sequences. The method of rooting the tree was responsible for most of the sequence variation both among the reconstructed ancestral sequences and between the ancestral and observed sequences. Variation in predicted 3-D structural properties of the ancestors mirrored their sequence variation. The observed sequence consensus and ancestral sequences from center-rooted trees were most similar in all predicted attributes. Only for the predicted number of N-glycosylation sites was there a difference between MP and ML methods of reconstruction. Taxon sampling effects were observed only for outgroup-rooted trees, not center-rooted, reflecting the occurrence of several divergent basal sequences. Thus, for sequences exhibiting a radial phylogenetic tree, as does HIV-1, most of the variation in the estimated ancestor arises from the method of rooting the phylogenetic tree. Those investigating the ancestors of genes exhibiting such a radial tree should pay particular attention to alternate rooting methods in order to obtain a representative sample of ancestors.

Mosaic structure of human coronavirus NL63, one thousand years of evolution.

Before the SARS outbreak only two human coronaviruses (HCoV) were known: HCoV-OC43 and HCoV-229E. With the discovery of SARS-CoV in 2003, a third family member was identified. Soon thereafter, we described the fourth human coronavirus (HCoV-NL63), a virus that has spread worldwide and is associated with croup in children. We report here the complete genome sequence of two HCoV-NL63 clinical isolates, designated Amsterdam 57 and Amsterdam 496. The genomes are 27,538 and 27,550 nucleotides long, respectively, and share the same genome organization. We identified two variable regions, one within the 1a and one within the S gene, whereas the 1b and N genes were most conserved. Phylogenetic analysis revealed that HCoV-NL63 genomes have a mosaic structure with multiple recombination sites. Additionally, employing three different algorithms, we assessed the evolutionary rate for the S gene of group Ib coronaviruses to be approximately 3 x 10(-4) substitutions per site per year. Using this evolutionary rate we determined that HCoV-NL63 diverged in the 11th century from its closest relative HCoV-229E.

Taxonomy: Sus bucculentus revisited.

In 1997, the rediscovery of Sus bucculentus in Laos was announced by Groves et al.--this wild pig species had gone unrecorded since first being described in 1892. Although the identification of the new specimen was based initially on morphology, the authors also used a 7% sequence divergence from the common Eurasian pig S. scrofa (based on their analysis of 327 base pairs of the gene encoding mitochondrial 12S ribosomal RNA) as support for the species status of S. bucculentus. Concerned about the large divergence reported for a relatively conserved gene, and the absence of the sequence in any public database, we analysed an additional tissue sample from the specimen and found only 0.6% divergence from S. scrofa. Our more extensive analysis places the sample within the S. scrofa clade, calling into question the species status of S. bucculentus and demonstrating the need for both phylogenetic and morphological evidence in defining species.