VIP Barcoding: Composition vector-based software for rapid species identification based on DNA barcoding

Species identification based on short sequences of DNA markers, that is, DNA barcoding, has emerged as an integral part of modern taxonomy. However, software for the analysis of large and multilocus barcoding data sets is scarce. The Basic Local Alignment Search Tool (BLAST) is currently the fastest tool capable of handling large databases (e.g. >5000 sequences), but its accuracy is a concern and has been criticized for its local optimization. However, current more accurate software requires sequence alignment or complex calculations, which are time-consuming when dealing with large data sets during data preprocessing or during the search stage. Therefore, it is imperative to develop a practical program for both accurate and scalable species identification for DNA barcoding. In this context, we present VIP Barcoding: a user-friendly software in graphical user interface for rapid DNA barcoding. It adopts a hybrid, two-stage algorithm. First, an alignment-free composition vector (CV) method is utilized to reduce searching space by screening a reference database. The alignment-based K2P distance nearest-neighbour method is then employed to analyse the smaller data set generated in the first stage. In comparison with other software, we demonstrate that VIP Barcoding has (i) higher accuracy than Blastn and several alignment-free methods and (ii) higher scalability than alignment-based distance methods and character-based methods. These results suggest that this platform is able to deal with both large-scale and multilocus barcoding data with accuracy and can contribute to DNA barcoding for modern taxonomy. VIP Barcoding is free and available at http://msl.sls.cuhk.edu.hk/vipbarcoding/.

Copper(II) Complexes of l-Arginine as Netropsin Mimics Showing DNA Cleavage Activity in Red Light

Copper(II) complexes [Cu(L-arg)(2)](NO3)(2) (1) and [Cu(L-arg)(B)Cl]Cl (2-5), where B is a heterocyclic base, namely, 2,2'-bipyridine (bpy, 2), 1,10-phenanthroline (phen, 3), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq, 4), and dipyrido[3,2-a:2',3'-c)phenazine (dppz, 5), are prepared and their DNA binding and photoinduced DNA cleavage activity studied. Ternary complex 3, structurally characterized using X-ray crystallography, shows a square-pyramidal (4 + 1) coordination geometry in which the N,O-donor L-arginine and N,N-donor 1,10-phenanthroline form the basal plane with one chloride at the elongated axial site. The complex has a pendant cationic guanidinium moiety. The one-electron paramagnetic complexes display a metal-centered d-d band in the range of 590-690 nm in aqueous DMF They show quasireversible cyclic voltammetric response due to the Cu(II)/Cu(I) couple in the range of -0.1 to -0.3 V versus a saturated calomel electrode in a DMF-Tris HCl buffer (pH 7.2). The DNA binding propensity of the complexes is studied using various techniques. Copper(II) bis-arginate 1 mimics the minor groove binder netropsin by showing preferential binding to the AT-rich sequence of double-strand (ds) DNA. DNA binding study using calf thymus DNA gives an order: 5 (L-arg-dppz) >= 1 (biS-L-arg) > 4 (L-arg-dpq) > 3 (L-arg-phen) >> 2 (L-arg-bpy). Molecular docking calculations reveal that the complexes bind through extensive hydrogen bonding and electrostatic interactions with ds-DNA. The complexes cleave supercoiled pUC19 DNA in the presence of 3-mercaptopropionic acid as a reducing agent forming hydroxyl ((OH)-O-center dot) radicals. The complexes show oxidative photoinduced DNA cleavage activity in UV-A light of 365 nm and red light of 647.1 nm (Ar-Kr mixed-gas-ion laser) in a metal-assisted photoexcitation process forming singlet oxygen (O-1(2)) species in a type-II pathway. All of the complexes, barring complex 2, show efficient DNA photocleavage activity. Complexes 4 and 5 exhibit significant double-strand breaks of DNA in red light of 647.1 nm due to the presence of two photosensitizers, namely, L-arginine and dpq or dppz in the molecules.

Arsenic trioxide promotes mitochondrial DNA mutation and cell apoptosis in primary APL cells and NB4 cell line

This study aimed to investigate the effects of arsenic trioxide (As2O3) on the mitochondrial DNA (mtDNA) of acute promyelocytic leukemia (APL) cells. The NB4 cell line was treated with 2.0 μmol/L As2O3in vitro, and the primary APL cells were treated with 2.0 μmol/L As2O3in vitro and 0.16 mg kg-1 d-1 As2O3in vivo. The mitochondrial DNA of all the cells above was amplified by PCR, directly sequenced and analyzed by Sequence Navigatore and Factura software. The apoptosis rates were assayed by flow cytometry. Mitochondrial DNA mutation in the D-loop region was found in NB4 and APL cells before As2O3 use, but the mutation spots were remarkably increased after As2O3 treatment, which was positively correlated to the rates of cellular apoptosis, the correlation coefficient: rNB4-As2O3=0.973818, and rAPL-As2O3=0.934703. The mutation types include transition, transversion, codon insertion or deletion, and the mutation spots in all samples were not constant and regular. It is revealed that As2O3 aggravates mtDNA mutation in the D-loop region of acute promyelocytic leukemia cells both in vitro and in vivo. Mitochondrial DNA might be one of the targets of As2O3 in APL treatment.

Crystalline A-DNA: The X-Ray Analysis of the Fragment d(G-G-T-A-T-A-C-C)

An A-DNA type double helical conformation was observed in the single crystal X-ray structure of the octamer d(G-G-T-A-T-A-C-C), 1, and its 5-bromouracil-containing analogue, 2. The structure of the isomorphous crystals (space group P61) was solved by a search technique based on packing criteria and R-factor calculations, with use of only low order data. At the present stage of refinement the R factors are 31 % for 1 and 28 % for 2 at a resolution of 2.25 A (0.225 nm). The molecules interact through their minor grooves by hydrogen bonding and base to sugar van der Waals contacts. The stable A conformation observed in the crystal may have some structural relevance to promoter regions where the T-A-T-A sequence is frequently found.

Identification of small juvenile scombrids from northwest tropical Australia using mitochondrial DNA cytochrome b sequences

Small juveniles of the nine species of scombrids in Australian waters are morphologically similar to one another and, consequently, difficult to identify to species level. We show that the sequence of the mitochondrial DNA cytochrome b gene region is a powerful tool for identification of these young fish. Using this method, we identified 50 juvenile scombrids collected from Exmouth Bay, Western Australia. Six species of scombrids were apparent in this sample of fish: narrow-barred Spanish mackerel (Scomberomorus commerson), Indian mackerel (Rastrelliger kanagurta), frigate tuna (Auxis thazard), bullet tuna (Auxis rochei), leaping bonito (Cybiosarda elegans), and kawakawa (Euthynnus affinis). The presence of Indian mackerel, frigate tuna, leaping bonito, and kawakawa is the first indication that coastal waters may be an important spawning habitat for these species, although offshore spawning may also occur. The occurrence of small juvenile S. commerson was predicted from the known spawning patterns of that species, but other mackerel species (Scomberomorus munroi, Scomberomorus queenslandicus, Scomberomorus semifasiciatus) likely to be spawning during the sampling period were not detected among the 50 small juveniles analyzed here.

DNA structural variability as a factor in gene expression and evolution

Redundant DNA can buffer sequence dependent structural deviations from an ideal double helix. Buffering serves a mechanistic function by reducing extraneous conformational effects which could interfere with readout or which would impose energetic constraints on evolution. It also serves an evolutionary function by allowing for gradual variations in conformation-dependent regulation of gene expression. Such gradualism is critical for the rate of evolution. The buffer structure concept provides a new interpretation for repetitive DNA and for exons and introns.

The assessment of TGGE for the detection of interspecific and intergeneric DNA-marker polymorphism within Solanaceae

RFLP markers are currently the most appropriate marker system for the identification of uncharacterised polymorphism at the interspecific and intergeneric level. Given the benefits of a PCR-based marker system and the availability of sequence information for many Solanaceous cDNA clones, it is now possible to target conserved fragments, for primer development, that flank sequences possessing interspecific polymorphism. The potential outcome is the development of a suite of markers that amplify widely in Solanaceae. Temperature gradient gel electrophoresis (TGGE) is a relatively inexpensive gel-based system that is suitable for the detection of most single-base changes. TGGE can be used to screen for both known and unknown polymorphisms, and has been assessed here, for the development of PCR-based markers that are useful for the detection of interspecific variation within Solanaceae. Fifteen markers are presented where differences between Lycopersicon esculentum and L. pennellii have been detected by TGGE. The markers were assessed on a wider selection of plant species and found to be potentially useful for the identification of interspecific and intergeneric polymorphism in Solanaceous plants.

DNA bar-coding reveals source and patterns of Thaumastocoris peregrinus invasions in South Africa and South America.

Thaumastocoris peregrinus is a recently introduced invertebrate pest of non-native Eucalyptus plantations in the Southern Hemisphere. It was first reported from South Africa in 2003 and in Argentina in 2005. Since then, populations have grown explosively and it has attained an almost ubiquitous distribution over several regions in South Africa on 26 Eucalyptus species. Here we address three key questions regarding this invasion, namely whether only one species has been introduced, whether there were single or multiple introductions into South Africa and South America and what the source of the introduction might have been. To answer these questions, bar-coding using mitochondrial DNA (COI) sequence diversity was used to characterise the populations of this insect from Australia, Argentina, Brazil, South Africa and Uruguay. Analyses revealed three cryptic species in Australia, of which only T. peregrinus is represented in South Africa and South America. Thaumastocoris peregrinus populations contained eight haplotypes, with a pairwise nucleotide distance of 0.2-0.9% from seventeen locations in Australia. Three of these haplotypes are shared with populations in South America and South Africa, but the latter regions do not share haplotypes. These data, together with the current distribution of the haplotypes and the known direction of original spread in these regions, suggest that at least three distinct introductions of the insect occurred in South Africa and South America before 2005. The two most common haplotypes in Sydney, one of which was also found in Brisbane, are shared with the non-native regions. Sydney populations of T. peregrinus, which have regularly reached outbreak levels in recent years, might thus have served as source of these three distinct introductions into other regions of the Southern Hemisphere.

DNA-based identification of Quambalaria pitereka causing severe leaf blight of Corymbia citriodora in China

Quambalaria spp. include serious plant pathogens, causing leaf and shoot blight of Corymbia and Eucalyptus spp. In this study, a disease resembling Quambalaria leaf blight was observed on young Corymbia citriodora trees in a plantation in the Guangdong Province of China. Comparisons of rDNA sequence data showed that the causal agent of the disease is Q. pitereka. This study provides the first report of Quambalaria leaf blight from China, and it is also the first time that this pathogen has been found on trees outside the native range of Eucalypts.

Sequences that adopt non-B-DNA conformation in form V DNA as probed by enzymic methylation

pBR322 form V DNA is a highly torsionally strained molecule with a linking number of zero. We have used sequence- specific DNA methylases as probes for B-DNA in this molecule, exploiting the inability of methylases to methylate single-stranded DNA and Z-DNA, both of which are known to occur in form V DNA. Some sequences in form V DNA were shown to be totally in the B-form, others were totally in an altered, unmethylatable conformation, while still other sites appeared to exist partly in altered and partly in normal B-conformation. Some potential Z-forming sequences (alternating pyrimidine/purine) of less than seven base-pairs were not in the Z conformation in form V DNA, whereas others did adopt an altered structure, indicating a modulating influence of flanking sequences. Furthermore, regions of imperfect alternating pyrimidine/purine structure were sometimes capable of adopting an altered structure. In addition, some regions of altered structure had no apparent Z-forming sequences, nor were they in polypurine stretches, which have also been proposed to form left-handed DNA. These non-B-DNA conformations may represent novel left-handed helical structures or sequences that become single stranded under torsional strain. Long regions of either altered (unmethylatable) DNA or B-DNA were not always observed. In fact, one region showed three transitions between B-like DNA and altered structure within 26 base-pairs.

DNA markers linked to the R(2) rust resistance gene in sunflower (Helianthus annuus L.) facilitate anticipatory breeding for this disease variant.

Pre-emptive breeding for host disease resistance is an effective strategy for combating and managing devastating incursions of plant pathogens. Comprehensive, long-term studies have revealed that virulence to the R (2) sunflower (Helianthus annuus L.) rust resistance gene in the line MC29 does not exist in the Australian rust (Puccinia helianthi) population. We report in this study the identification of molecular markers linked to this gene. The three simple sequence repeat (SSR) markers ORS795, ORS882, and ORS938 were linked in coupling to the gene, while the SSR marker ORS333 was linked in repulsion. Reliable selection for homozygous-resistant individuals was efficient when the three markers, ORS795, ORS882, and ORS333, were used in combination. Phenotyping for this resistance gene is not possible in Australia without introducing a quarantinable race of the pathogen. Therefore, the availability of reliable and heritable DNA-based markers will enable the efficient deployment of this gene, permitting a more effective strategy for generating sustainable commercial cultivars containing this rust resistance gene.

Comparative analysis of Tritrichomonas foetus (Riedmuller, 1928) cat genotype, T. foetus (Riedmuller, 1928) cattle genotype and Tritrichomonas suis (Davaine, 1875) at 10 DNA loci

The parasitic protists in the genus Tritrichomonas cause significant disease in domestic cattle and cats. To assess the genetic diversity of feline and bovine isolates of Tritrichomonas foetus (Riedmuller, 1928) Wenrich and Emmerson, 1933, we used 10 different genetic regions, namely the protein coding genes of cysteine proteases 1,2 and 4-9 (CP1, 2, 4-9) involved in the pathogenesis of the disease caused by the parasite. The cytosolic malate dehydrogenase 1 (MDH1) and internal transcribed spacer region 2 of the rDNA unit (ITS2) were included as additional markers. The gene sequences were compared with those of Tritrichomonas suis (Davaine. 1875) Morgan and Hawkins, 1948 and Tritrichomonas mobilensis Culberson et al., 1986. The study revealed 100% identity for all 10 genes among all feline isolates (=T. foetus cat genotype), 100% identity among all bovine isolates (=T. foetus cattle genotype) and a genetic distinctness of 1% between the cat and cattle genotypes of T. foetus. The cattle genotype of T. foetus was 100% identical to T. suis at nine loci (CP1, 2,4-8, ITS2, MDH1). At CP9, three out of four T. suis isolates were identical to the T. foetus cattle genotype, while the T. suis isolate SUI-H3B sequence contained a single unique nucleotide substitution. Tritrichomonas mobilensis was 0.4% and 0.7% distinct from the cat and cattle genotypes of T. foetus, respectively. The genetic differences resulted in amino acid changes in the CP genes, most pronouncedly in CP2, potentially providing a platform for elucidation of genotype-specific host-pathogen interactions of T. foetus. On the basis of this data we judge T. suis and T. foetus to be subjective synonyms. For the first time, on objective nomenclatural grounds, the authority of T. suis is given to Davaine, 1875, rather than the commonly cited Gruby and Delafond, 1843. To maintain prevailing usage of T. foetus, we are suppressing the senior synomym T. suis Davaine, 1875 according to Article 23.9, because it has never been used as a valid name after 1899 and T. foetus is widely discussed as the cause of bovine trichomonosis. Thus bovine, feline and porcine isolates should all be given the name T. foetus. This promotes the stability of T. foetus for the veterinary and economically significant venereal parasite causing bovine trichomonosis. (C) 2012 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Mitochondrial DNA diversity of Cleruchoides noackae (Hymenoptera: Mymaridae): a potential biological control agent for Thaumastocoris peregrinus (Hemiptera: Thaumastocoridae)

Carpintero and Dellap, (Hemiptera: Thaumastocoridae) is a native Australian sap-feeding insect that has become invasive and seriously damaging to commercially grown in the Southern Hemisphere. Lin and Huber (Hymenoptera: Mymaridae) was recently discovered as an egg parasitoid of the Thaumastocoridae in Australia. Mitochondrial DNA (mtDNA; cytochrome oxidase subunit I, COI) sequence diversity amongst 104 individuals from these native populations revealed 24 sequence haplotypes. The COI haplotypes of individuals collected from the Sydney and Southeast Queensland clustered in distinct groups, indicating limited spread of the insect between the regions. Individuals collected from Perth in Western Australia were represented by four COI haplotypes. Although this population is geographically more isolated from other populations, two COI haplotypes were identical to haplotypes found in the Sydney region. The results suggest that has recently been introduced into Perth, possibly from the Sydney area. The high mtDNA diversity and limited spread that is suggested for is in contrast to the lack of geographic associated mtDNA diversity and extensive spread of . If implemented as a biological control agent, this factor will need to be considered in collecting and releasing .

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: - 1) priming site polymorphism in the template leads to inferior or erratic amplification; - 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions, and; - 3) at least occasionally, a PCR primer straddles an exon–intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3′ end of coding sequence and 3′ UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

Characterization of Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease, reveals a unique mode of DNA-binding, helical distortion and cleavage, compared to a canonical LAGLIDADG homing endonuclease

Mycobacterium leprae, which has undergone reductive evolution leaving behind a minimal set of essential genes, has retained intervening sequences in four of its genes implicating a vital role for them in the survival of the leprosy bacillus. A single in-frame intervening sequence has been found embedded within its recA gene. Comparison of M. leprae recA intervening sequence with the known intervening sequences indicated that it has the consensus amino acid sequence necessary for being a LAGLIDADG-type homing endonuclease. In light of massive gene decay and function loss in the leprosy bacillus, we sought to investigate whether its recA intervening sequence encodes a catalytically active homing endonuclease. Here we show that the purified M. leprae RecA intein (PI-MleI) binds to cognate DNA and displays endonuclease activity in the presence of alternative divalent cations, Mg2+ or Mn2+. A combination of approaches including four complementary footprinting assays such as DNase I, Cu/phenanthroline, methylation protection and KMnO4, enhancement of 2-aminopurine fluorescence and mapping of the cleavage site revealed that PI-MleI binds to cognate DNA flanking its insertion site, induces helical distortion at the cleavage site and generates two staggered double-strand breaks. Taken together, these results implicate that PI-MleI possess a modular structure with separate domains for DNA target recognition and cleavage, each with distinct sequence preferences. From a biological standpoint, it is tempting to speculate that our findings have implications for understanding the evolution of LAGLIDADG family of homing endonucleases