Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

PCR and real-time PCR primers developed for detection and identification of Bifidobacterium thermophilum in faeces

BACKGROUND: Culture-independent methods based on the 16S ribosomal RNA molecule are nowadays widely used for assessment of the composition of the intestinal microbiota, in relation to host health or probiotic efficacy. Because Bifidobacterium thermophilum was only recently isolated from human faeces until now, no specific real-time PCR (qPCR) assay has been developed for detection of this species as component of the bifidobacterial community of the human intestinal flora. RESULTS: Design of specific primers and probe was achieved based on comparison of 108 published bifidobacterial 16S rDNA sequences with the recently published sequence of the human faecal isolate B. thermophilum RBL67. Specificity of the primer was tested in silico by similarity search against the sequence database and confirmed experimentally by PCR amplification on 17 Bifidobacterium strains, representing 12 different species, and two Lactobacillus strains. The qPCR assay developed was linear for B. thermophilum RBL67 DNA quantities ranging from 0.02 ng/microl to 200 ng/microl and showed a detection limit of 10(5) cells per gram faeces. The application of this new qPCR assay allowed to detect the presence of B. thermophilum in one sample from a 6-month old breast-fed baby among 17 human faecal samples tested. Additionally, the specific qPCR primers in combination with selective plating experiments led to the isolation of F9K9, a faecal isolate from a 4-month old breast-fed baby. The 16S rDNA sequence of this isolate is 99.93% similar to that of B. thermophilum RBL67 and confirmed the applicability of the new qPCR assay in faecal samples. CONCLUSION: A new B. thermophilum-specific qPCR assay was developed based on species-specific target nucleotides in the 16S rDNA. It can be used to further characterize the composition of the bifidobacterial community in the human gastrointestinal tract. Until recently, B. thermophilum was considered as a species of animal origin, but here we confirm with the application of this new PCR assay the presence of B. thermophilum strains in the human gut.

Intestinal Tritrichomonas foetus infection in cats in Switzerland detected by in vitro cultivation and PCR

Tritrichomonas foetus, a parasite well known for its significance as venereally transmitted pathogen in cattle, has recently been identified as a cause of chronic large-bowel diarrhea in domestic cats in the US, UK, and, more recently, also in Norway. In a period of 3 months (October to December 2007), 45 cats of Switzerland suffering from chronic diarrhea were investigated for intestinal infections, including a search for trichomonads. A commercially available in vitro culture system was used to screen for infection, complemented with a PCR and subsequent amplicon sequencing to support speciation. The PCR is based upon amplification of a sequence derived from the internal transcribed spacer region 1 (ITS1) on the ribosomal RNA gene (rRNA) using primers designed to detect a broad range of genera and species belonging to the family of Trichomonadidae. The method was furthermore adapted to the uracil DNA glycosylase (UDG) system in order to prevent carry-over contamination and it included a recombinant internal control to track for inhibitory reactions. Eleven out of the 45 cats were culture-positive, as revealed by microscopic identification of trichomonadid organisms. One of the isolates was subjected to scanning electron microscopy and findings revealed the presence of three flagella, thus placing the isolate into the gender Tritrichomonas sp. PCR and subsequent amplicon sequencing were carried out with ten of the 11 isolates. A total homology with published T. foetus sequences was confirmed in all of the cases. T. foetus therefore appears to range among those organisms that can cause chronic diarrhea in cats in Switzerland.

The small ribosomal subunit from Thermus thermophilus at 4.5 Å resolution: Pattern fittings and the identification of a functional site

The electron density map of the small ribosomal subunit from Thermus thermophilus, constructed at 4.5 Å resolution, shows the recognizable morphology of this particle, as well as structural features that were interpreted as ribosomal RNA and proteins. Unbiased assignments, carried out by quantitative covalent binding of heavy atom compounds at predetermined sites, led to the localization of the surface of the ribosomal protein S13 at a position compatible with previous assignments, whereas the surface of S11 was localized at a distance of about twice its diameter from the site suggested for its center by neutron scattering. Proteins S5 and S7, whose structures have been determined crystallographically, were visually placed in the map with no alterations in their conformations. Regions suitable to host the fold of protein S15 were detected in several positions, all at a significant distance from the location of this protein in the neutron scattering map. Targeting the 16S RNA region, where mRNA docks to allow the formation of the initiation complex by a mercurated mRNA analog, led to the characterization of its vicinity.

Cell cycle-dependent regulation of RNA polymerase I transcription: The nucleolar transcription factor UBF is inactive in mitosis and early G1

Transcription of ribosomal RNA genes by RNA polymerase (pol) I oscillates during the cell cycle, being maximal in S and G2 phase, repressed during mitosis, and gradually recovering during G1 progression. We have shown that transcription initiation factor (TIF)-IB/SL1 is inactivated during mitosis by cdc2/cyclin B-directed phosphorylation of TAFI110. In this study, we have monitored reactivation of transcription after exit from mitosis. We demonstrate that the pol I factor UBF is also inactivated by phosphorylation but recovers with different kinetics than TIF-IB/SL1. Whereas TIF-IB/SL1 activity is rapidly regained on entry into G1, UBF is reactivated later in G1, concomitant with the onset of pol I transcription. Repression of pol I transcription in mitosis and early G1 can be reproduced with either extracts from cells synchronized in M or G1 phase or with purified TIF-IB/SL1 and UBF isolated in the presence of phosphatase inhibitors. The results suggest that two basal transcription factors, e.g., TIF-IB/SL1 and UBF, are inactivated at mitosis and reactivated by dephosphorylation at the exit from mitosis and during G1 progression, respectively.

Association of RNase mitochondrial RNA processing enzyme with ribonuclease P in higher ordered structures in the nucleolus: a possible coordinate role in ribosome biogenesis.

RNase mitochondrial RNA processing enzyme (MRP) is a nucleolar ribonucleoprotein particle that participates in 5.8S ribosomal RNA maturation in eukaryotes. This enzyme shares a polypeptide and an RNA structural motif with ribonuclease P (RNase P), a nuclear endoribonuclease originally described in the nucleus that processes RNA transcripts to generate their mature 5' termini. Both enzymes are also located in mitochondria. This report further characterizes the relationship between RNase MRP and RNase P. Antisense affinity selection with biotinylated 2'-O-methyl oligoribonucleotides and glycerol gradient fractionation experiments demonstrated that small subpopulations of RNase MRP and RNase P associate with each other in vivo in macromolecular complex, possibly 60-80S preribosomes. This latter notion was supported by fluorescence in situ hybridization experiments with antisense oligonucleotides that localized that RNA components of RNase MRP and RNase P to the nucleolus and to discrete cytoplasmic structures. These findings suggest that small subpopulations of RNase MRP and RNase P are physically associated, and that both may function in ribosomal RNA maturation or ribosome assembly.

Application of nox-restriction fragment length polymorphism for the differentiation of Brachyspira intestinal spirochetes isolated from pigs and poultry in Australia

Sixty-nine intestinal spirochetes isolated from pigs and poultry in eastern Australia were selected to evaluate the effectiveness of a species-specific PCR-based restriction fragment length polymorphism (RFLP) analysis of the Brachyspira nox gene. For comparative purposes, all isolates were subjected to species-specific PCRs for the pathogenic species Brachyspira hyodysenteriae and Brachyspira pilosicoli, and selected isolates were examined further by sequence analysis of the nox and 16S ribosomal RNA genes. Modifications to the original nox-RFLP method included direct inoculation of bacterial cells into the amplification mixture and purification of the PCR product, which further optimized the nox-RFLP for use in a veterinary diagnostic laboratory, producing sufficient product for both species identification and future comparisons. Although some novel profiles that prevented definitive identification were observed, the nox-RFLP method successfully classified 45 of 51 (88%) porcine and 15 of 18 (83%) avian isolates into 5 of the 6 recognized species of Brachyspira. This protocol represents a significant improvement over conventional methods currently used in veterinary diagnostic laboratories for rapid specific identification of Brachyspira spp. isolated from both pigs and poultry.

Molecular characterization of Thelastomatoidea (Nematoda : Oxyurida) from cockroaches in Australia

A molecular approach was used to genetically characterize 5 species (Aoruroides queenslandensis. Blattophila sphaerolaima, Cordonicola gibsoni, Desmicola ornato and Leidynemella fusiformis) belonging to the superfamily. Thelastomatoidea fi (Nematoda: Oxyurida), a group of pinworms that parasitizes terrestrial arthropods. The D3 domain of the large subunit Of nuclear ribosomal RNA (LSU) was sequenced for individual specimens, and the analysis of the sequence data allowed the genetic relationships of the 5 species to be studied dagger. The sequence variation in the D3 domain within individual species (0-1-8%) was significantly less than the differences among species (4(.)3-12(.)4%). Phylogenetic analyses, Using maximum parsimony, maximum likelihood, and neighbour-joining, tree-building methods, established relationships among the 5 species of Thelastomatoidea and Oxyuris equi (a species of the order Oxyurida). The molecular approach employed provides the prospect for developing DNA tools for the specific identification of the Thelastomatoidea, irrespective of developmental stage and sex, as a basis for systematic, ecological and/or population genetic investigations of members within this superfamily.

Nucleolar localization of aprataxin is dependent on interaction with nucleolin and on active ribosomal DNA transcription

The APTX gene, mutated in patients with the neurological disorder ataxia with oculomotor apraxia type 1 (AOA1), encodes a novel protein aprataxin. We describe here, the interaction and interdependence between aprataxin and several nucleolar proteins, including nucleolin, nucleophosmin and upstream binding factor-1 (UBF-1), involved in ribosomal RNA (rRNA) synthesis and cellular stress signalling. Interaction between aprataxin and nucleolin occurred through their respective N-terminal regions. In AOA1 cells lacking aprataxin, the stability of nucleolin was significantly reduced. On the other hand, down-regulation of nucleolin by RNA interference did not affect aprataxin protein levels but abolished its nucleolar localization suggesting that the interaction with nucleolin is involved in its nucleolar targeting. GFP-aprataxin fusion protein co-localized with nucleolin, nucleophosmin and UBF-1 in nucleoli and inhibition of ribosomal DNA transcription altered the distribution of aprataxin in the nucleolus, suggesting that the nature of the nucleolar localization of aprataxin is also dependent on ongoing rRNA synthesis. In vivo rRNA synthesis analysis showed only a minor decrease in AOA1 cells when compared with controls cells. These results demonstrate a cross-dependence between aprataxin and nucleolin in the nucleolus and while aprataxin does not appear to be directly involved in rRNA synthesis its nucleolar localization is dependent on this synthesis.

Ribosomal frameshifting and transcriptional slippage: from genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

Influence of historical landscapes, drainage evolution and ecological traits on patterns of genetic diversity in Southeast Asian freshwater snakehead fishes

Snakehead fishes in the family Channidae are obligate freshwater fishes represented by two extant genera, the African Parachannna and the Asian Channa. These species prefer still or slow flowing water bodies, where they are top predators that exercise high levels of parental care, have the ability to breathe air, can tolerate poor water quality, and interestingly, can aestivate or traverse terrestrial habitat in response to seasonal changes in freshwater habitat availability. These attributes suggest that snakehead fishes may possess high dispersal potential, irrespective of the terrestrial barriers that would otherwise constrain the distribution of most freshwater fishes. A number of biogeographical hypotheses have been developed to account for the modern distributions of snakehead fishes across two continents, including ancient vicariance during Gondwanan break-up, or recent colonisation tracking the formation of suitable climatic conditions. Taxonomic uncertainty also surrounds some members of the Channa genus, as geographical distributions for some taxa across southern and Southeast (SE) Asia are very large, and in one case is highly disjunct. The current study adopted a molecular genetics approach to gain an understanding of the evolution of this group of fishes, and in particular how the phylogeography of two Asian species may have been influenced by contemporary versus historical levels of dispersal and vicariance. First, a molecular phylogeny was constructed based on multiple DNA loci and calibrated with fossil evidence to provide a dated chronology of divergence events among extant species, and also within species with widespread geographical distributions. The data provide strong evidence that trans-continental distribution of the Channidae arose as a result of dispersal out of Asia and into Africa in the mid–Eocene. Among Asian Channa, deep divergence among lineages indicates that the Oligocene-Miocene boundary was a time of significant species radiation, potentially associated with historical changes in climate and drainage geomorphology. Mid-Miocene divergence among lineages suggests that a taxonomic revision is warranted for two taxa. Deep intra-specific divergence (~8Mya) was also detected between C. striata lineages that occur sympatrically in the Mekong River Basin. The study then examined the phylogeography and population structure of two major taxa, Channa striata (the chevron snakehead) and the C. micropeltes (the giant snakehead), across SE Asia. Species specific microsatellite loci were developed and used in addition to a mitochondrial DNA marker (Cyt b) to screen neutral genetic variation within and among wild populations. C. striata individuals were sampled across SE Asia (n=988), with the major focus being the Mekong Basin, which is the largest drainage basin in the region. The distributions of two divergent lineages were identified and admixture analysis showed that where they co-occur they are interbreeding, indicating that after long periods of evolution in isolation, divergence has not resulted in reproductive isolation. One lineage is predominantly confined to upland areas of northern Lao PDR to the north of the Khorat Plateau, while the other, which is more closely related to individuals from southern India, has a widespread distribution across mainland SE Asian and Sumatra. The phylogeographical pattern recovered is associated with past river networks, and high diversity and divergence among all populations sampled reveal that contemporary dispersal is very low for this taxon, even where populations occur in contiguous freshwater habitats. C. micropeltes (n=280) were also sampled from across the Mekong River Basin, focusing on the lower basin where it constitutes an important wild fishery resource. In comparison with C. striata, allelic diversity and genetic divergence among populations were extremely low, suggesting very recent colonisation of the greater Mekong region. Populations were significantly structured into at least three discrete populations in the lower Mekong. Results of this study have implications for establishing effective conservation plans for managing both species, that represent economically important wild fishery resources for the region. For C. micropeltes, it is likely that a single fisheries stock in the Tonle Sap Great Lake is being exploited by multiple fisheries operations, and future management initiatives for this species in this region will need to account for this. For C. striata, conservation of natural levels of genetic variation will require management initiatives designed to promote population persistence at very localised spatial scales, as the high level of population structuring uncovered for this species indicates that significant unique diversity is present at this fine spatial scale.

The complete mitochondrial genome of the gall-forming fly, Fergusonina taylori(Nelson and Yeates) (Diptera: Fergusoninidae)

The monogeneric family Fergusoninidae consists of gall-forming flies that, together with Fergusobia (Tylenchida: Neotylenchidae) nematodes, form the only known mutualistic association between insects and nematodes. In this study, the entire 16,000 bp mitochondrial genome of Fergusonina taylori Nelson and Yeates was sequenced. The circular genome contains one encoding region including 27 genes and one non-coding A þT-rich region. The arrangement of the proteincoding, ribosomal RNA (rRNA) and transfer RNA (tRNA) genes was the same as that found in the ancestral insect. Nucleotide composition is highly A þ T biased. All of the protein initiation codons are ATN, except for nad1 which begins with TTT. All 22 tRNA anticodons of F. taylori match those observed in Drosophila yakuba, and all form the typical cloverleaf structure except for tRNA-Ser (AGN) which lacks a dihydrouridine (DHU) arm. Secondary structural features of the rRNA genes of Fergusonina are similar to those proposed for other insects, with minor modifications. The mitochondrial genome of Fergusonina presented here may prove valuable for resolving the sister group to the Fergusoninidae, and expands the available mtDNA data sources for acalyptrates overall.

Ureaplasma parvum : understanding the complexities of intra-amniotic infection in an ovine model

The human Ureaplasma species are the most frequently isolated bacteria from the upper genital tract of pregnant women and can cause clinically asymptomatic, intra-uterine infections, which are difficult to treat with antimicrobials. Ureaplasma infection of the upper genital tract during pregnancy has been associated with numerous adverse outcomes including preterm birth, chorioamnionitis and neonatal respiratory diseases. The mechanisms by which ureaplasmas are able to chronically colonise the amniotic fluid and avoid eradication by (i) the host immune response and (ii) maternally-administered antimicrobials, remain virtually unexplored. To address this gap within the literature, this study investigated potential mechanisms by which ureaplasmas are able to cause chronic, intra-amniotic infections in an established ovine model. In this PhD program of research the effectiveness of standard, maternal erythromycin for the treatment of chronic, intra-amniotic ureaplasma infections was evaluated. At 55 days of gestation pregnant ewes received an intra-amniotic injection of either: a clinical Ureaplasma parvum serovar 3 isolate that was sensitive to macrolide antibiotics (n = 16); or 10B medium (n = 16). At 100 days of gestation, ewes were then randomised to receive either maternal erythromycin treatment (30 mg/kg/day for four days) or no treatment. Ureaplasmas were isolated from amniotic fluid, chorioamnion, umbilical cord and fetal lung specimens, which were collected at the time of preterm delivery of the fetus (125 days of gestation). Surprisingly, the numbers of ureaplasmas colonising the amniotic fluid and fetal tissues were not different between experimentally-infected animals that received erythromycin treatment or infected animals that did not receive treatment (p > 0.05), nor were there any differences in fetal inflammation and histological chorioamnionitis between these groups (p > 0.05). These data demonstrate the inability of maternal erythromycin to eradicate intra-uterine ureaplasma infections. Erythromycin was detected in the amniotic fluid of animals that received antimicrobial treatment (but not in those that did not receive treatment) by liquid chromatography-mass spectrometry; however, the concentrations were below therapeutic levels (<10 – 76 ng/mL). These findings indicate that the ineffectiveness of standard, maternal erythromycin treatment of intra-amniotic ureaplasma infections may be due to the poor placental transfer of this drug. Subsequently, the phenotypic and genotypic characteristics of ureaplasmas isolated from the amniotic fluid and chorioamnion of pregnant sheep after chronic, intra-amniotic infection and low-level exposure to erythromycin were investigated. At 55 days of gestation twelve pregnant ewes received an intra-amniotic injection of a clinical U. parvum serovar 3 isolate, which was sensitive to macrolide antibiotics. At 100 days of gestation, ewes received standard maternal erythromycin treatment (30 mg/kg/day for four days, n = 6) or saline (n = 6). Preterm fetuses were surgically delivered at 125 days of gestation and ureaplasmas were cultured from the amniotic fluid and the chorioamnion. The minimum inhibitory concentrations (MICs) of erythromycin, azithromycin and roxithromycin were determined for cultured ureaplasma isolates, and antimicrobial susceptibilities were different between ureaplasmas isolated from the amniotic fluid (MIC range = 0.08 – 1.0 mg/L) and chorioamnion (MIC range = 0.06 – 5.33 mg/L). However, the increased resistance to macrolide antibiotics observed in chorioamnion ureaplasma isolates occurred independently of exposure to erythromycin in vivo. Remarkably, domain V of the 23S ribosomal RNA gene (which is the target site of macrolide antimicrobials) of chorioamnion ureaplasmas demonstrated significant variability (125 polymorphisms out of 422 sequenced nucleotides, 29.6%) when compared to the amniotic fluid ureaplasma isolates and the inoculum strain. This sequence variability did not occur as a consequence of exposure to erythromycin, as the nucleotide substitutions were identical between chorioamnion ureaplasmas isolated from different animals, including those that did not receive erythromycin treatment. We propose that these mosaic-like 23S ribosomal RNA gene sequences may represent gene fragments transferred via horizontal gene transfer. The significant differences observed in (i) susceptibility to macrolide antimicrobials and (ii) 23S ribosomal RNA sequences of ureaplasmas isolated from the amniotic fluid and chorioamnion suggests that the anatomical site from which they were isolated may exert selective pressures that alter the socio-microbiological structure of the bacterial population, by selecting for genetic changes and altered antimicrobial susceptibility profiles. The final experiment for this PhD examined antigenic size variation of the multiple banded antigen (MBA, a surface-exposed lipoprotein and predicted ureaplasmal virulence factor) in chronic, intra-amniotic ureaplasma infections. Previously defined ‘virulent-derived’ and ‘avirulent-derived’ clonal U. parvum serovar 6 isolates (each expressing a single MBA protein) were injected into the amniotic fluid of pregnant ewes (n = 20) at 55 days of gestation, and amniotic fluid was collected by amniocentesis every two weeks until the time of near-term delivery of the fetus (at 140 days of gestation). Both the avirulent and virulent clonal ureaplasma strains generated MBA size variants (ranging in size from 32 – 170 kDa) within the amniotic fluid of pregnant ewes. The mean number of MBA size variants produced within the amniotic fluid was not different between the virulent (mean = 4.2 MBA variants) and avirulent (mean = 4.6 MBA variants) ureaplasma strains (p = 0.87). Intra-amniotic infection with the virulent strain was significantly associated with the presence of meconium-stained amniotic fluid (p = 0.01), which is an indicator of fetal distress in utero. However, the severity of histological chorioamnionitis was not different between the avirulent and virulent groups. We demonstrated that ureaplasmas were able to persist within the amniotic fluid of pregnant sheep for 85 days, despite the host mounting an innate and adaptive immune response. Pro-inflammatory cytokines (interleukin (IL)-1â, IL-6 and IL-8) were elevated within the chorioamnion tissue of pregnant sheep from both the avirulent and virulent treatment groups, and this was significantly associated with the production of anti-ureaplasma IgG antibodies within maternal sera (p < 0.05). These findings suggested that the inability of the host immune response to eradicate ureaplasmas from the amniotic cavity may be due to continual size variation of MBA surface-exposed epitopes. Taken together, these data confirm that ureaplasmas are able to cause long-term in utero infections in a sheep model, despite standard antimicrobial treatment and the development of a host immune response. The overall findings of this PhD project suggest that ureaplasmas are able to cause chronic, intra-amniotic infections due to (i) the limited placental transfer of erythromycin, which prevents the accumulation of therapeutic concentrations within the amniotic fluid; (ii) the ability of ureaplasmas to undergo rapid selection and genetic variation in vivo, resulting in ureaplasma isolates with variable MICs to macrolide antimicrobials colonising the amniotic fluid and chorioamnion; and (iii) antigenic size variation of the MBA, which may prevent eradication of ureaplasmas by the host immune response and account for differences in neonatal outcomes. The outcomes of this program of study have improved our understanding of the biology and pathogenesis of this highly adapted microorganism.

Beyond barcoding : a mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.

The complete mitochondrial genome of the flesh fly, Sarcophaga impatiens Walker (Diptera: Sarcophagidae)

Approximately 2500 fly species comprise the Sarcophagidae family worldwide. The complete mitochondrial genome of the carrion-breeding, forensically important Sarcophaga impatiens Walker (Diptera: Sarcophagidae) from Australia was sequenced. The 15,169 bp circular genome contains the 37 genes found in a typical Metazoan genome: 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. It also contains one non-coding A+T-rich region. The arrangement of the genes was the same as that found in the ancestral insect. All the protein initiation codons are ATN, except for cox1 that begins with TCG (encoding S). The 22 tRNA anticodons of S. impatiens are consistent with those observed in Drosophila yakuba, and all form the typical cloverleaf structure, except for tRNA-Ser(AGN) that lacks the DHU arm. The mitochondrial genome of Sarcophaga presented will be valuable for resolving phylogenetic relationships within the family Sarcophagidae and the order Diptera, and could be used to identify favourable genetic markers for species identifications for forensic purposes.