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/.

Chlamydial infections of fish: Diverse pathogens and emerging causes of disease in aquaculture species

Chlamydial infections of fish are emerging as an important cause of disease in new and established aquaculture industries. To date, epitheliocystis, a skin and gill disease associated with infection by these obligate intracellular pathogens, has been described in over 90 fish species, including hosts from marine and fresh water environments. Aided by advances in molecular detection and typing, recent years have seen an explosion in the description of these epitheliocystis-related chlamydial pathogens of fish, significantly broadening our knowledge of the genetic diversity of the order Chlamydiales. Remarkably, in most cases, it seems that each new piscine host studied has revealed the presence of a phylogenetically unique and novel chlamydial pathogen, providing researchers with a fascinating opportunity to understand the origin, evolution and adaptation of their traditional terrestrial chlamydial relatives. Despite the advances in this area, much still needs to be learnt about the epidemiology of chlamydial infections in fish if these pathogens are to be controlled in farmed environments. The lack of in vitro methods for culturing of chlamydial pathogens of fish is a major hindrance to this field. This review provides an update on our current knowledge of the taxonomy and diversity of chlamydial pathogens of fish, discusses the impact of these infections on the health, and highlights further areas of research required to understand the biology and epidemiology of this important emerging group of fish pathogens of aquaculture species.

Evaluation of the relationship between Chlamydia pecorum sequence types and disease using a species-specific multi-locus sequence typing scheme (MLST)

Chlamydia pecorum is globally associated with several ovine diseases including keratoconjunctivitis and polyarthritis. The exact relationship between the variety of C. pecorum strains reported and the diseases described in sheep remains unclear, challenging efforts to accurately diagnose and manage infected flocks. In the present study, we applied C. pecorum multi-locus sequence typing (MLST) to C. pecorum positive samples collected from sympatric flocks of Australian sheep presenting with conjunctivitis, conjunctivitis with polyarthritis, or polyarthritis only and with no clinical disease (NCD) in order to elucidate the exact relationships between the infecting strains and the range of diseases. Using Bayesian phylogenetic and cluster analyses on 62 C. pecorum positive ocular, vaginal and rectal swab samples from sheep presenting with a range of diseases and in a comparison to C. pecorum sequence types (STs) from other hosts, one ST (ST 23) was recognised as a globally distributed strain associated with ovine and bovine diseases such as polyarthritis and encephalomyelitis. A second ST (ST 69) presently only described in Australian animals, was detected in association with ovine as well as koala chlamydial infections. The majority of vaginal and rectal C. pecorum STs from animals with NCD and/or anatomical sites with no clinical signs of disease in diseased animals, clustered together in a separate group, by both analyses. Furthermore, 8/13 detected STs were novel. This study provides a platform for strain selection for further research into the pathogenic potential of C. pecorum in animals and highlights targets for potential strain-specific diagnostic test development.

Experimental and theoretical characterization of 5,10- diminoporphodimethenes: Dearomatized porphyrinoids from palladium-catalyzed hydrazinations of 5,10-diarylporphyrins

A new dearomatized porphyrinoid, 5,10-diiminoporphodimethene (5,10-DIPD), has been prepared by palladium-catalyzed hydrazination of 5,10-dibromo-15,20-bis(3,5-di-tert-butylphenyl)porphyrin and its nickel(II) complex, by using ethyl and 4-methoxybenzyl carbazates. The oxidative dearomatization of the porphyrin ring occurs in high yield. Further oxidation with 2,3-dichloro-5,6-dicyanobenzoquinone forms the corresponding 5,10-bis(azocarboxylates), thereby restoring the porphyrin aromaticity. The UV/visible spectra of the NiII DIPDs exhibit remarkable redshifts of the lowest-energy bands to 780 nm, and differential pulse voltammetry reveals a contracted electrochemical HOMO–LUMO gap of 1.44 V. Density functional theory (DFT) was used to calculate the optimized geometries and frontier molecular orbitals of model 5,10-DIPD Ni7c and 5,10-bis(azocarboxylate) Ni8c. The conformations of the carbamate groups and the configurations of the CNZ unit were considered in conjunction with the NOESY spectra, to generate the global minimum geometry and two other structures with slightly higher energies. In the absence of solution data regarding conformations, ten possible local minimum conformations were considered for Ni8c. Partition of the porphyrin macrocycle into tri- and monopyrrole fragments in Ni7c and the inclusion of terminal conjugating functional groups generate unique frontier molecular orbital distributions and a HOMO–LUMO transition with a strong element of charge transfer from the monopyrrole ring. Time-dependent DFT calculations were performed for the three lowest-energy structures of Ni7c and Ni8c, and weighting according to their energies allowed the prediction of the electronic spectra. The calculations reproduce the lower-energy regions of the spectra and the overall forms of the spectra with high accuracy, but agreement is not as good in the Soret region below 450 nm.

Combining the diaza-Diels-Alder reaction and palladium-catalyzed aminations to prepare amino-substituted porphyrins

A strategy to tackle the synthesis of azoporphyrins with unsubstituted terminal meso positions was investigated. It comprised the combination of diaza-Diels–Alder (DADA) reaction of 1,3-dienes with dialkyl azodicarboxylates, decarboxylative hydrolysis of the bis(carbamates), palladium-catalyzed amination of bromoporphyrin precursors, and retro-DADA reactions to release the ultimate targets. The somewhat confused historical results on the DADA reactions of 1,3-cyclohexadiene were clarified, but the hydrolyses yielded extremely air-sensitive amines which decomposed completely in minutes via autooxidation and retro-DADA reaction. With anthracene or 2,3-dimethyl-1,3-butadiene as the diene, the synthesis of azoporphyrin was not achieved but three amino-substituted porphyrins were obtained in moderate yields under mild conditions. The X-ray crystal structures of several of the intermediates and the final aminoanthracene-porphyrin nickel(II) complex were determined.

Assistive classification for improving the efficiency of avian species richness surveys

Avian species richness surveys, which measure the total number of unique avian species, can be conducted via remote acoustic sensors. An immense quantity of data can be collected, which, although rich in useful information, places a great workload on the scientists who manually inspect the audio. To deal with this big data problem, we calculated acoustic indices from audio data at a one-minute resolution and used them to classify one-minute recordings into five classes. By filtering out the non-avian minutes, we can reduce the amount of data by about 50% and improve the efficiency of determining avian species richness. The experimental results show that, given 60 one-minute samples, our approach enables to direct ecologists to find about 10% more avian species.

Computer-assisted sampling of acoustic data for more efficient determination of bird species richness

Bird species richness survey is one of the most intriguing ecological topics for evaluating environmental health. Here, bird species richness denotes the number of unique bird species in a particular area. Factors affecting the investigation of bird species richness include weather, observation bias, and most importantly, the prohibitive costs of conducting surveys at large spatiotemporal scales. Thanks to advances in recording techniques, these problems have been alleviated by deploying sensors for acoustic data collection. Although automated detection techniques have been introduced to identify various bird species, the innate complexity of bird vocalizations, the background noise present in the recording and the escalating volumes of acoustic data pose a challenging task on determination of bird species richness. In this paper we proposed a two-step computer-assisted sampling approach for determining bird species richness in one-day acoustic data. First, a classification model is built based on acoustic indices for filtering out minutes that contain few bird species. Then the classified bird minutes are ordered by an acoustic index and the redundant temporal minutes are removed from the ranked minute sequence. The experimental results show that our method is more efficient in directing experts for determination of bird species compared with the previous methods.

Amidinato- and triazenido-bridged binuclear complexes of palladium and platinum

The reaction of [M2Cl2(mu-Cl)(2)(PR3)(2)] (M=Pd or Pt; PR3=PEt3, PBu3, PMe2Ph, PMePh2) with lithium amidinate or sodium triazenide gave binuclear complexes containing amidinato- or triazenido-bridges, [M2Cl2(mu-ArNENAr)(2)(PR3)(2)] (E=CH, CMe or N). These complexes were characterized by elemental analysis and NMR (H-1, P-31 or Pt-195) data. The structures of two complexes, [(PdCl2)-Cl-2(mu-PhNC(Me)NPh)(2)(PMe2Ph)(2)] (10) and [Pt2Cl2(mu-PhNNNPh)(2)(PEt3)(2)] (11) were established by single crystal X-ray structural analyses. The Pt-195 NMR data Show coupling between two metal centers in the cis triazenido-bridged complex. The corresponding amidinate bridged complex does not show coupling. The role of the bridging ligand in mediating interaction between the metal centers is probed through Extended Huckel Theory (EHT) calculations. It is suggested that M-M interactions are primarily affected by the bridging ligands

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

trans-Dichloro-bis-(arylazoimidazole)palladium(II): Synthesis, Structure, Photoisomerization, and DFT Calculation

Reaction between PdCl2 and 1-alkyl-2-(arylazo)imidazole (RaaiR') or 1-alkyl-2-(naphthyl-alpha/beta-azo)imidazole (alpha/beta-NaiR') under reflux in ethanol has isolated complexes of compositions Pd(RaaiR')(2)Cl-2 (5, 6) and Pd(alpha/beta-NaiR')(2)Cl-2 (7, 8). The X-ray structure determination of one of the molecules, Pd(alpha-NaiBz)(2)Cl-2 (7c), has reported a trans-PdCl2 configuration, and alpha-NaiBz acts as monodentate N(imidazole) donor ligand. The spectral (IR, UV-vis, H-1 NMR) data support the structure. UV light irradiation (light source: Perkin-Elmer LS 55 spectrofluorimeter, Xenon discharge lamp, lambda = 360-396 nm) in a MeCN solution of the complexes shows E-to-Z isomerization of the coordinated azoimidazole unit. The reverse transformation, Z-to-E, is very slow with visible light irradiation. Quantum yields (phi(E-Z)) of E-to-Z isomerization are calculated, and phi is lower than that of the free ligand but comparable with those of Cd(II) and Hg(II) complexes of the same ligand. The Z-to-E isomerization is a thermally induced process. The activation energy (E-a) of Z-to-E isomerization is calculated by controlled-temperature experimentation. cis-Pd(azoimidazole)Cl-2 complexes (azomidazole acts as N(imidazole) and N(azo) Chelating ligand) do not respond upon light irradiation, which supports the idea that the presence of noncoordinated azo-N to make free azo (-N=N-) function is important to reveal photochromic activity. DFT calculation of Pd(alpha-NaiBz)(2)Cl-2 (7c) has suggested that the HOMO of the molecule is constituted of Pd (32%) and Cl (66%), and hence photo excitation may use the energy of Pd and Cl instead of that of the photofunctional -N=N-Ar motif; thus, the rate of photoisomerization and quantum yield decrease versus the free ligand values.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected

Subsampling multi-species trawl catches from tropical northern Australia: Does it matter which part of the catch is sampled?

Subsampling is a common method for estimating the abundance of species in trawl catches. However, the accuracy of subsampling in representing the total catch has not been assessed. To estimate one possible source of bias due to subsampling, we tested whether the position on trawler sorting trays from which subsamples were taken affected their ability to represent species in catches. This was done by sorting catches into 10 kg subsamples and comparing subsamples taken from different positions on the sorting tray. Comparisons were made after species were grouped into three categories of abundance, either 'rare', 'common' or 'abundant'. A generalised linear model analysis showed that taking subsamples from different positions on the sorting tray had no major effect on estimating the total numbers or weights of fish or invertebrates, or the total number of fish or invertebrate taxa, recorded in each position. Some individual taxa showed differences between positions on the sorting tray (11.5% of taxa ina three-position design; 25% in a five-position design). But consistent and meaningful patterns in the position of these taxa on the sorting tray could only be seen for the pony fish Leiognathus moretoniensis and the saucer scallop Amusium pleuronectes. Because most bycatch laxa are well mixed throughout the catch, subsamples can be taken from any position on trawler sorting trays without introducing bias.