Ecological niche modelling and DNA sequencing support a new, morphologically cryptic beetle species unveiled by DNA barcoding

<b>Background:</b> DNA sequencing techniques used to estimate biodiversity, such as DNA barcoding, may reveal cryptic species. However, disagreements between barcoding and morphological data have already led to controversy. Species delimitation should therefore not be based on mtDNA alone. Here, we explore the use of nDNA and bioclimatic modelling in a new species of aquatic beetle revealed by mtDNA sequence data. <br /><br /><b>Methodology/Principal Findings:</b> The aquatic beetle fauna of Australia is characterised by high degrees of endemism, including local radiations such as the genus Antiporus. Antiporus femoralis was previously considered to exist in two disjunct, but morphologically indistinguishable populations in south-western and south-eastern Australia. We constructed a phylogeny of Antiporus and detected a deep split between these populations. Diagnostic characters from the highly variable nuclear protein encoding arginine kinase gene confirmed the presence of two isolated populations. We then used ecological niche modelling to examine the climatic niche characteristics of the two populations. All results support the status of the two populations as distinct species. We describe the south-western species as Antiporus occidentalis sp.n. <br /><br /><b>Conclusion/Significance:</b> In addition to nDNA sequence data and extended use of mitochondrial sequences, ecological niche modelling has great potential for delineating morphologically cryptic species. &Acirc;&copy; 2011 Hawlitschek et al.

The development of 20 microsatellite loci for the Australian marine mollusk, Donax deltoides, through next generation DNA sequencing

&nbsp;Next Generation DNA sequencing was used to develop a suite of microsatellite markers for the marine mollusk, Donax deltoides. A total of 20 polymorphic loci were identified and 12 characterized using 30 individuals from a single population (Venus Bay) in south eastern Australia. We observed moderate to high genetic variation across most loci (mean number of alleles per locus = 7.3; mean heterozygosity = 0.633) with only a single locus (Ddel32) displaying significant deviation from Hardy&ndash;Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, however two loci were found to be influenced by null alleles. The 10 viable markers characterized in the present study provide a valuable resource for future population genetic assessments and fisheries management of D. deltoides in Australia.

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as 'endangered' under the Environment Protection and Biodiversity Conservation Act 1999, and 'vulnerable' under the International Union for Conservation of Nature's Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F ST = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.

Mitochondrial DNA analysis of field populations of Helicoverpa armigera (Lepidoptera : Noctuidae) and of its relationship to H. zea

<b>Background</b><br />Helicoverpa armigera and H. zea are amongst the most significant polyphagous pest lepidopteran species in the Old and New Worlds respectively. Separation of H. armigera and H. zea is difficult and is usually only achieved through morphological differences in the genitalia. They are capable of interbreeding to produce fertile offspring. The single species status of H. armigera has been doubted, due to its wide distribution and plant host range across the Old World. This study explores the global genetic diversity of H. armigera and its evolutionary relationship to H zea.<br /><br /><b>Results</b><br />We obtained partial (511 bp) mitochondrial DNA (mtDNA) Cytochrome Oxidase-I (COI) sequences for 249 individuals of H. armigera sampled from Australia, Burkina Faso, Uganda, China, India and Pakistan which were associated with various host plants. Single nucleotide polymorphisms (SNPs) within the partial COI gene differentiated H. armigera populations into 33 mtDNA haplotypes. Shared haplotypes between continents, low F-statistic values and low nucleotide diversity between countries (0.0017 &ndash; 0.0038) suggests high mobility in this pest. Phylogenetic analysis of four major Helicoverpa pest species indicates that H. punctigera is basal to H. assulta, which is in turn basal to H. armigera and H. zea. Samples from North and South America suggest that H. zea is also a single species across its distribution. Our data reveal short genetic distances between H. armigera and H. zea which seem to have been established via a founder event from H. armigera stock at around 1.5 million years ago.<br /><br /><b>Conclusion</b><br />Our mitochondrial DNA sequence data supports the single species status of H. armigera across Africa, Asia and Australia. The evidence for inter-continental gene flow observed in this study is consistent with published evidence of the capacity of this species to migrate over long distances. The finding of high genetic similarity between Old World H. armigera and New World H. zea emphasises the need to consider work on both pests when building pest management strategies for either. <br />

Geographical variation and population structure in the White-rumped Sandpiper Calidris fuscicollis as shown by morphology, mitochondrial DNA and carbon isotope ratios

We studied the population structure of a high arctic breeding wader bird species, the White-rumped Sandpiper Calidris fuscicollis. Breeding adults, chicks and juveniles were sampled at seven localities throughout the species' breeding range in arctic Canada in 1999. The mitochondrial control region was analysed by DNA sequencing, feathers were analysed for carbon isotope ratios (C13/C12) by isotope ratio mass spectrometry, and morphological measurements were analysed using principal component analyses, taking the effect of sex into account (identified by molecular genetic methods). In general, our results support the notion that the White-rumped Sandpiper is a monotypic species with no subspecies, and most of the morphological and genetic variation occurs within sites. Nevertheless, some differences between sites were found. Birds from the two northernmost sites (Ellesmere and Devon Islands) had relatively longer bill and wing and shorter tarsus than birds sampled further south, possibly reflecting genetic differences between populations. The carbon isotope ratios were higher at the easternmost site (Baffin Island), revealing differences in the isotope content of the food. The mtDNA sequences showed no significant differentiation between sites and no pattern of isolation-by-distance was found. Based on the mtDNA variation, the species was estimated to have a long-term effective population size of approximately 9,000 females. The species shows no clear evidence of any population expansion or decline. Our results indicate that carbon isotope ratios, and possibly also certain mtDNA haplotypes, may be useful as tools for identifying the breeding origin of White-rumped Sandpipers on migration and wintering sites.

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

Development and application of polymerase chain reaction-based assays for Rathayibacter Toxicus and a Bacteriophage associated with annual ryegrass (Lolium Rigidum) toxicity

Annual ryegrass toxicity (ARGT) is responsible for significant stock losses in South Australia and Western Australia. The toxicity is caused by corynetoxins produced by the bacterium <i>Rathayibacter toxicus</i> (with the possible involvement of a bacteriophage), which infects annual ryegrass (<i>Lolium</i> <i>rigidum</i>). Polymerase chain reaction (PCR)-based assays, compatible with an existing enzyme-linked immunosorbent assay for the corynetoxins, have been developed and used to screen <i>L. rigidum</i> for both the presence of <i>R</i>. <i>toxicus</i> and for the bacteriophage isolate NCPPB 3778. The results from analysing bacterially infected galls from toxic grain screenings showed a positive correlation between the presence of the bacterium and corynetoxins but not with the bacteriophage. Analysis of pasture-derived samples of annual ryegrass showed about a 50% correlation of corynetoxins with bacterial presence and about a 5% correlation of phage with the presence of the bacterium. These observations support the potential application of the PCR-based assays in providing a useful, complementary tool in the assessment of the likelihood of pasture and feed to cause ARGT and to enable a better understanding of the complex aetiology of ARGT.<br />

A parallel DNA fragment assembly algorithm based on eulerian superpath approach

Fragments assembly is among the core problems in the research of Genome. Although many assembly tools based on the &quot;overlap-layout-consensus&quot; paradigm are widely used such as in the Human Genome Project currently, they still can not resolve the &quot;repeats problem&quot; in the DNA sequencing. For the purpose of resolving such problem, Pevzner et al. put forward a new Euler Superpath assembly algorithm. But it needs a big and complex de Bruijin graph which consumes large amounts of memories i.e. becomes the bottleneck of the performance. We present a parallel DNA fragment assembly algorithm based on the Eularian Superpath theory and solve the bottleneck in the current assembly program. The experimental results demonstrate that our approach has a good scalability, and can be used in DNA assembly of middle and large size of eukaryote genome.

Extracting DNA from museum bird eggs, and whole genome amplification of archive DNA

We present a comprehensive protocol for extracting DNA from egg membranes and other internal debris recovered from the interior of blown museum bird eggs. A variety of commercially available DNA extraction methods were found to be applicable. DNA sequencing of polymerase chain reaction (PCR) products for a 176-bp fragment of mitochondrial DNA was successful for most egg samples (&gt; 78%) even though the amount of DNA extracted (mean = 14.71 &plusmn; 4.55 ng/&micro;L) was significantly less than that obtained for bird skin samples (mean = 67.88 &plusmn; 4.77 ng/&micro;L). For PCR and sequencing of snipe (Gallinago) DNA, we provide eight new primers for the &lsquo;DNA barcode&rsquo; region of COI mtDNA. In various combinations, the primers target a range of PCR products sized from 72 bp to the full &lsquo;barcode&rsquo; of 751 bp. Not all possible combinations were tested with archive snipe DNA, but we found a significantly better success rate of PCR amplification for a shorter 176-bp target compared with a larger 288-bp fragment (67% vs. 39%). Finally, we explored the feasibility of whole genome amplification (WGA) for extending the use of archive DNA in PCR and sequencing applications. Of two WGA approaches, a PCR-based method was found to be able to amplify whole genomic DNA from archive skins and eggs from museum bird collections. After WGA, significantly more archive egg samples produced visible PCR products on agarose (56.9% before WGA vs. 79.0% after WGA). However, overall sequencing success did not improve significantly (78.8% compared with 83.0%).

RNASeqBrowser: a genome browser for simultaneous visualization of raw strand specific RNAseq reads and UCSC genome browser custom tracks

BACKGROUND: Strand specific RNAseq data is now more common in RNAseq projects. Visualizing RNAseq data has become an important matter in Analysis of sequencing data. The most widely used visualization tool is the UCSC genome browser that introduced the custom track concept that enabled researchers to simultaneously visualize gene expression at a particular locus from multiple experiments. Our objective of the software tool is to provide friendly interface for visualization of RNAseq datasets. <br /><br />RESULTS: This paper introduces a visualization tool (RNASeqBrowser) that incorporates and extends the functionality of the UCSC genome browser. For example, RNASeqBrowser simultaneously displays read coverage, SNPs, InDels and raw read tracks with other BED and wiggle tracks -- all being dynamically built from the BAM file. Paired reads are also connected in the browser to enable easier identification of novel exon/intron borders and chimaeric transcripts. Strand specific RNAseq data is also supported by RNASeqBrowser that displays reads above (positive strand transcript) or below (negative strand transcripts) a central line. Finally, RNASeqBrowser was designed for ease of use for users with few bioinformatic skills, and incorporates the features of many genome browsers into one platform. <br /><br />CONCLUSIONS: The features of RNASeqBrowser: (1) RNASeqBrowser integrates UCSC genome browser and NGS visualization tools such as IGV. It extends the functionality of the UCSC genome browser by adding several new types of tracks to show NGS data such as individual raw reads, SNPs and InDels. (2) RNASeqBrowser can dynamically generate RNA secondary structure. It is useful for identifying non-coding RNA such as miRNA. (3) Overlaying NGS wiggle data is helpful in displaying differential expression and is simple to implement in RNASeqBrowser. (4) NGS data accumulates a lot of raw reads. Thus, RNASeqBrowser collapses exact duplicate reads to reduce visualization space. Normal PC's can show many windows of NGS individual raw reads without much delay. (5) Multiple popup windows of individual raw reads provide users with more viewing space. This avoids existing approaches (such as IGV) which squeeze all raw reads into one window. This will be helpful for visualizing multiple datasets simultaneously. RNASeqBrowser and its manual are freely available at http://www.australianprostatecentre.org/research/software/rnaseqbrowser or http://sourceforge.net/projects/rnaseqbrowser/.

Creatine and the creatine transporter: a review

The cellular role of creatine (Cr) and Cr phosphate (CrP) has been studied extensively in neural, cardiac and skeletal muscle. Several studies have demonstrated that alterations in the cellular total Cr (Cr + CrP) concentration in these tissues can produce marked functional and/or structural change. The primary aim of this review was to critically evaluate the literature that has examined the regulation of cellular total Cr content. In particular, the review focuses on the regulation of the activity and gene expression of the Cr transporter (CreaT), which is primarily responsible for cellular Cr uptake. Two CreaT genes (CreaT1 and CreaT2) have been identified and their chromosomal location and DNA sequencing have been completed. From these data, putative structures of the CreaT proteins have been formulated. Transcription products of the CreaT2 gene are expressed exclusively in the testes, whereas CreaT1 transcripts are found in a variety of tissues. Recent research has measured the expression of the CreaT1 protein in several tissues including neural, cardiac and skeletal muscle. There is very little information available about the factors regulating CreaT gene expression. There is some evidence that suggests the intracellular Cr concentration may be involved in the regulatory process but there is much more to learn before this process is understood. The activity of the CreaT protein is controlled by many factors. These include substrate concentration, transmembrane Na+ gradients, cellular location, and various hormones. It is also likely that transporter activity is influenced by its phosphorylation state and by its interaction with other plasma membrane proteins. The extent of CreaT protein glycosylation may vary within cells, the functional significance of which remains unclear.<br />

Functional and metabolic remodelling in GLUT4-deficient hearts confers hyper-responsiveness to substrate intervention

Impaired glucose uptake is associated with both cardiac hypertrophy and contractile dysfunction, but whether there are common underlying&nbsp; mechanisms linking these conditions is yet to be determined. Using a &lsquo;gene dose&rsquo; Cre-Lox GLUT4-deficient murine model, we examined the effect of suppressed glucose availability on global myocardial gene expression and glycolysis substrate bypass on the function of isolated perfused hearts. Performance of hearts from 22- to 60-week-old male GLUT4 knockout (KO, &gt; 95% reduction in GLUT4), GLUT4 knockdown (KD, 85% reduction in cardiac GLUT4) and C57Bl/6 wild-type (WT) controls was measured ex vivo in Langendorff mode perfusion. DNA microarray was used to profile mRNA expression differences between GLUT4-KO and GLUT4-KD hearts. At 22 weeks, GLUT4-KO hearts exhibited cardiac hypertrophy and impaired contractile function ex vivo, characterized by a 40% decrease in developed pressure. At 60 weeks, dysfunction was accentuated in GLUT4-KO hearts and evident in GLUT4-KD hearts. Exogenous pyruvate (5 mM) restored systolic pressure to a level equivalent to WT (GLUT4-KO, 176.8 &plusmn; 13.2 mmHg vs. WT, 146.4 &plusmn; 9.56 mmHg) in 22-week-old GLUT4-KO hearts but not in 60-week-old GLUT4-KO hearts. In GLUT4-KO, DNA microarray analysis detected downregulation of a number of genes centrally involved in mitochondrial oxidation and upregulation of other genes indicative of a shift to cytosolic &beta;-oxidation of long chain fatty acids. A direct link between cardiomyocyte GLUT4 deficiency, hypertrophy and contractile dysfunction is demonstrated. These data provide mechanistic insight into the myocardial metabolic adaptations associated with short and long-term insulin resistance and indicate a window of opportunity for substrate intervention and functional &lsquo;rescue&rsquo;.<br />

Molecular systematic studies of freshwater crayfish of the the genus Cherax

This study investigated the molecular systematics and zoogeography of freshwater crayfish species (yabbies) of the genus Cherax throughout Australia. Using DNA sequencing of up to five gene regions and phylogenetic methods of analyses, this study provided new insights into the taxonomy and evolution of this important group of Australian freshwater invertebrates.