Ana Pavasovic

One and the same: integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis

The invasive fruit fly Bactrocera invadens Drew, Tsuruta & White, and the Oriental fruit fly Bactrocera dorsalis (Hendel) are highly destructive horticultural pests of global significance. Bactrocera invadens originates from the Indian subcontinent and has recently invaded all of sub‐Saharan Africa, while B. dorsalis principally occurs from the Indian subcontinent towards southern China and South‐east Asia. High morphological and genetic similarity has cast doubt over whether B. invadens is a distinct species from B. dorsalis. Addressing this issue within an integrative taxonomic framework, we sampled from across the geographic distribution of both taxa and: (i) analysed morphological variation, including those characters considered diagnostic (scutum colour, length of aedeagus, width of postsutural lateral vittae, wing size, and wing shape); (ii) sequenced four loci (ITS1, ITS2, cox1 and nad4) for phylogenetic inference; and (iii) generated a cox1 haplotype network to examine population structure. Molecular analyses included the closely related species, Bactrocera kandiensis Drew & Hancock. Scutum colour varies from red‐brown to fully black for individuals from Africa and the Indian subcontinent. All individuals east of the Indian subcontinent are black except for a few red‐brown individuals from China. The postsutural lateral vittae width of B. invadens is narrower than B. dorsalis from eastern Asia, but the variation is clinal, with subcontinent B. dorsalis populations intermediate in size. Aedeagus length, wing shape and wing size cannot discriminate between the two taxa. Phylogenetic analyses failed to resolve B. invadens from B. dorsalis, but did resolve B. kandiensis. Bactrocera dorsalis and B. invadens shared cox1 haplotypes, yet the haplotype network pattern does not reflect current taxonomy or patterns in thoracic colour. Some individuals of B. dorsalis/B. invadens possessed haplotypes more closely related to B. kandiensis than to conspecifics, suggestive of mitochondrial introgression between these species. The combined evidence fails to support the delimitation of B. dorsalis and B. invadens as separate biological species. Consequently, existing biological data for B. dorsalis may be applied to the invasive population in Africa. Our recommendation, in line with other recent publications, is that B. invadens be synonymized with B. dorsalis.

A transcriptome resource for the koala ( Phascolarctos cinereus ): insights into koala retrovirus transcription and sequence diversity

BackgroundThe koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus.ResultsRNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene.Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population.ConclusionsThis transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org.

Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of De novo assemblers

BackgroundThe sequencing, de novo assembly and annotation of transcriptome datasets generated with next generation sequencing (NGS) has enabled biologists to answer genomic questions in non-model species with unprecedented ease. Reliable and accurate de novo assembly and annotation of transcriptomes, however, is a critically important step for transcriptome assemblies generated from short read sequences. Typical benchmarks for assembly and annotation reliability have been performed with model species. To address the reliability and accuracy of de novo transcriptome assembly in non-model species, we generated an RNAseq dataset for an intertidal gastropod mollusc species, Nerita melanotragus, and compared the assembly produced by four different de novo transcriptome assemblers; Velvet, Oases, Geneious and Trinity, for a number of quality metrics and redundancy.ResultsTranscriptome sequencing on the Ion Torrent PGM™ produced 1,883,624 raw reads with a mean length of 133 base pairs (bp). Both the Trinity and Oases de novo assemblers produced the best assemblies based on all quality metrics including fewer contigs, increased N50 and average contig length and contigs of greater length. Overall the BLAST and annotation success of our assemblies was not high with only 15-19% of contigs assigned a putative function.ConclusionsWe believe that any improvement in annotation success of gastropod species will require more gastropod genome sequences, but in particular an increase in mollusc protein sequences in public databases. Overall, this paper demonstrates that reliable and accurate de novo transcriptome assemblies can be generated from short read sequencers with the right assembly algorithms.

The koala immunological toolkit: sequence identification and comparison of key markers of the koala (Phascolarctos cinereus) immune response

Abstract. The koala (Phascolarctos cinereus) is an Australian marsupial that continues to experience significant population declines. Infectious diseases caused by pathogens such as Chlamydia are proposed to have a major role. Very few species-specific immunological reagents are available, severely hindering our ability to respond to the threat of infectious diseases in the koala. In this study, we utilise data from the sequencing of the koala transcriptome to identify key immunological markers of the koala adaptive immune response and cytokines known to be important in the host response to chlamydial infection in other species. This report describes the identification and preliminary sequence analysis of (1) T lymphocyte glycoprotein markers (CD4, CD8); (2) IL-4, a marker for the Th2 response; (3) cytokines such as IL-6, IL-12 and IL-1β, that have been shown to have a role in chlamydial clearance and pathology in other hosts; and (4) the sequences for the koala immunoglobulins, IgA, IgG, IgE and IgM. These sequences will enable the development of a range of immunological reagents for understanding the koala’s innate and adaptive immune responses, while also providing a resource that will enable continued investigations into the origin and evolution of the marsupial immune system.

Comparative Analysis and Distribution of Omega-3 lcPUFA Biosynthesis Genes in Marine Molluscs.

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.

Comparative analysis of gill transcriptomes of two freshwater crayfish, Cherax cainii and C. destructor.

We undertook deep sequencing of gill transcriptomes from two freshwater crayfish, Cherax cainii and Cherax destructor, in order to generate genomic resources for future genomics research. Over 83 and 100 million high quality (quality score (Q)≥30) paired-end Illumina reads (150 bp) were assembled into 147,101 and 136,622 contigs in C. cainii and C. destructor, respectively. A total of 24,630 and 23,623 contigs received significant BLASTx hits and allowed the identification of multiple gill expressed candidate genes associated with pH and salinity balance. These functionally annotated transcripts will provide a resource to facilitate comparative genomic research in the genus Cherax, and in particular allow insights into respiratory and osmoregulatory physiology of this group of animals.

Sea Anemones: Quiet Achievers in the Field of Peptide Toxins

Sea anemones have been understudied as a source of peptide and protein toxins, with relatively few examined as a source of new pharmacological tools or therapeutic leads. This is surprising given the success of some anemone peptides that have been tested, such as the potassium channel blocker from Stichodactyla helianthus known as ShK. An analogue of this peptide, ShK-186, which is now known as dalazatide, has successfully completed Phase 1 clinical trials and is about to enter Phase 2 trials for the treatment of autoimmune diseases. One of the impediments to the exploitation of sea anemone toxins in the pharmaceutical industry has been the difficulty associated with their high-throughput discovery and isolation. Recent developments in multiple ‘omic’ technologies, including genomics, transcriptomics and proteomics, coupled with advanced bioinformatics, have opened the way for large-scale discovery of novel sea anemone toxins from a range of species. Many of these toxins will be useful pharmacological tools and some will hopefully prove to be valuable therapeutic leads.

Transcriptomic investigation of wound healing and regeneration in the cnidarian Calliactis polypus.

Wound healing and regeneration in cnidarian species, a group that forms the sister phylum to Bilateria, remains poorly characterised despite the ability of many cnidarians to rapidly repair injuries, regenerate lost structures, or re-form whole organisms from small populations of somatic cells. Here we present results from a fully replicated RNA-Seq experiment to identify genes that are differentially expressed in the sea anemone Calliactis polypus following catastrophic injury. We find that a large-scale transcriptomic response is established in C. polypus, comprising an abundance of genes involved in tissue patterning, energy dynamics, immunity, cellular communication, and extracellular matrix remodelling. We also identified a substantial proportion of uncharacterised genes that were differentially expressed during regeneration, that appear to be restricted to cnidarians. Overall, our study serves to both identify the role that conserved genes play in eumetazoan wound healing and regeneration, as well as to highlight the lack of information regarding many genes involved in this process. We suggest that functional analysis of the large group of uncharacterised genes found in our study may contribute to better understanding of the regenerative capacity of cnidarians, as well as provide insight into how wound healing and regeneration has evolved in different lineages.

Draft Genome Sequence of Enterococcus faecalis Strain PF3, Isolated from Adélie Penguin Feces from Antarctica

ABSTRACT Enterococcus faecalis is one of the leading causes of nosocomial infections and is a common commensal organism in humans and other animals. In this study, we report a draft genome sequence for the E. faecalis strain PF3, isolated from Adélie penguin feces collected from Warriner Island, Antarctica.

Insights into the phylogenetic and molecular evolutionary histories of Fad and Elovl gene families in Actiniaria

Abstract The biosynthesis of long‐chain polyunsaturated fatty acids (LC‐PUFAs, ≥ C20) is reliant on the action of desaturase and elongase enzymes, which are encoded by the fatty acyl desaturase (Fad) and elongation of very long‐chain fatty acid (Elovl) gene families, respectively. In Metazoa, research investigating the distribution and evolution of these gene families has been restricted largely to Bilateria. Here, we provide insights into the phylogenetic and molecular evolutionary histories of the Fad and Elovl gene families in Cnidaria, the sister phylum to Bilateria. Four model cnidarian genomes and six actiniarian transcriptomes were interrogated. Analysis of the fatty acid composition of a candidate cnidarian species, Actinia tenebrosa, was performed to determine the baseline profile of this species. Phylogenetic analysis revealed lineage‐specific gene duplication in actiniarians for both the Fad and Elovl gene families. Two distinct cnidarian Fad clades clustered with functionally characterized Δ5 and Δ6 proteins from fungal and plant species, respectively. Alternatively, only a single cnidarian Elovl clade clustered with functionally characterized Elovl proteins (Elovl4), while two additional clades were identified, one actiniarian‐specific (Novel ElovlA) and the another cnidarian‐specific (Novel ElovlB). In actiniarians, selection analyses revealed pervasive purifying selection acting on both gene families. However, codons in the Elovl gene family show patterns of nucleotide variation consistent with the action of episodic diversifying selection following gene duplication events. Significantly, these codons may encode amino acid residues that are functionally important for Elovl proteins to target and elongate different precursor fatty acids. In A. tenebrosa, the fatty acid analysis revealed an absence of LC‐PUFAs > C20 molecules and implies that the Elovl enzymes are not actively contributing to the elongation of these LC‐PUFAs. Overall, this study has revealed that actiniarians possess Fad and Elovl genes required for the biosynthesis of some LC‐PUFAs, and that these genes appear to be distinct from bilaterians.