Repair of 3-methylthymine and 1-methylguanine lesions by bacterial and human AlkB proteins

The Escherichia coli AlkB protein repairs 1-methyladenine (1-meA) and 3-methylcytosine (3-meC) lesions in DNA and RNA by oxidative demethylation, a reaction requiring ferrous iron and 2-oxoglutarate as cofactor and co-substrate, respectively. Here, we have studied the activity of AlkB proteins on 3-methylthymine (3-meT) and 1-methylguanine (1-meG), two minor lesions which are structurally analogous to 1-meA and 3-meC. AlkB as well as the human AlkB homologues, hABH2 and hABH3, were all able to demethylate 3-meT in a DNA oligonucleotide containing a single 3-meT residue. Also, 1-meG lesions introduced by chemical methylation of tRNA were efficiently removed by AlkB. Unlike 1-meA and 3-meC, nucleosides or bases corresponding to 1-meG or 3-meT did not stimulate the uncoupled, AlkB-mediated decarboxylation of 2-oxoglutarate. Our data show that 3-meT and 1-meG are repaired by AlkB, but indicate that the recognition of these substrates is different from that in the case of 1-meA and 3-meC.     

Bioinformatics and functional analysis define four distinct groups of AlkB DNA-dioxygenases in bacteria

The iron(II)- and 2-oxoglutarate (2OG)-dependent dioxygenase AlkB from Escherichia coli (EcAlkB) repairs alkylation damage in DNA by direct reversal. EcAlkB substrates include methylated bases, such as 1-methyladenine (m¹A) and 3-methylcytosine (m³C), as well as certain bulkier lesions, for example the exocyclic adduct 1,N⁶-ethenoadenine (εA). EcAlkB is the only bacterial AlkB protein characterized to date, and we here present an extensive bioinformatics and functional analysis of bacterial AlkB proteins. Based on sequence phylogeny, we show that these proteins can be subdivided into four groups: denoted 1A, 1B, 2A and 2B; each characterized by the presence of specific conserved amino acid residues in the putative nucleotide-recognizing domain. A scattered distribution of AlkB proteins from the four different groups across the bacterial kingdom indicates a substantial degree of horizontal transfer of AlkB genes. DNA repair activity was associated with all tested recombinant AlkB proteins. Notably, both a group 2B protein from Xanthomonas campestris and a group 2A protein from Rhizobium etli repaired etheno adducts, but had negligible activity on methylated bases. Our data indicate that the majority, if not all, of the bacterial AlkB proteins are DNA repair enzymes, and that some of these proteins do not primarily target methylated bases.     

Detection of Methylated CDO1 in Plasma of Colorectal Cancer; A PCR Study

Background

Cysteine biology is important for the chemosensitivity of cancer cells. Our research has focused on the epigenetic silencing of cysteine dioxygenase type 1 (CDO1) in colorectal cancer (CRC). In this study, we describe detection of CDO1 methylation in the plasma of CRC patients using methylation specific PCR (Q-MSP) and extensive analysis of the PCR reaction.

Methods

DNA was extracted from plasma, and analysed for methylation of the CDO1 gene using Q-MSP. The detection rate of CDO1 gene methylation was calculated and compared with that of diluted DNA extracted from “positive control” DLD1 cells. CDO1 gene methylation in the plasma of 40 CRC patients that were clinicopathologically analysed was then determined.

Results

(1) The cloned sequence analysis detected 93.3% methylation of the promoter CpG islands of the CDO1 gene of positive control DLD1 cells and 4.7% methylation of the negative control HepG2 CDO1 gene. (2) DLD1 CDO1 DNA could not be detected in this assay if the extracted DNA was diluted ∼1000 fold. The more DNA that was used for the PCR reaction, the more effectively it was amplified in Q-MSP. (3) By increasing the amount of DNA used, methylated CDO1 could be clearly detected in the plasma of 8 (20%) of the CRC patients. However, the percentage of CRC patients detected by methylated CDO1 in plasma was lower than that detected by CEA (35.9%) or CA19-9 (23.1%) in preoperative serum. Combination of CEA/CA19-9 plus plasma methylated CDO1 could increase the rate of detection of curable CRC patients (39.3%) as compared to CEA/CA19-9 (25%).

Conclusion

We have described detection of CDO1 methylation in the plasma of CRC patients. Although CDO1 methylation was not detected as frequently as conventional tumor markers, analysis of plasma CDO1 methylation in combination with CEA/CA19-9 levels increases the detection rate of curable CRC patients.     

Bimodal signatures of germline methylation are linked with gene expression plasticity in the coral Acropora millepora

Background

In invertebrates, genes belonging to dynamically regulated functional categories appear to be less methylated than “housekeeping” genes, suggesting that DNA methylation may modulate gene expression plasticity. To date, however, experimental evidence to support this hypothesis across different natural habitats has been lacking.

Results

Gene expression profiles were generated from 30 pairs of genetically identical fragments of coral Acropora millepora reciprocally transplanted between distinct natural habitats for 3 months. Gene expression was analyzed in the context of normalized CpG content, a well-established signature of historical germline DNA methylation. Genes with weak methylation signatures were more likely to demonstrate differential expression based on both transplant environment and population of origin than genes with strong methylation signatures. Moreover, the magnitude of expression differences due to environment and population were greater for genes with weak methylation signatures.

Conclusions

Our results support a connection between differential germline methylation and gene expression flexibility across environments and populations. Studies of phylogenetically basal invertebrates such as corals will further elucidate the fundamental functional aspects of gene body methylation in Metazoa.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1109) contains supplementary material, which is available to authorized users.     

A Two-Gene Blood Test for Methylated DNA Sensitive for Colorectal Cancer

Background

Specific genes are methylated with high frequency in colorectal neoplasia, and may leak into blood. Detection of multiple methylated DNA biomarkers in blood may improve assay sensitivity for colorectal cancer (CRC) relative to a single marker. We undertook a case-control study evaluating the presence of two methylation DNA markers, BCAT1 and IKZF1, in circulation to determine if they were complementary for detection of CRC.

Methods

Methylation-specific PCR assays were developed to measure the level of methylated BCAT1 and IKZF1 in DNA extracted from plasma obtained from colonoscopy-confirmed 144 healthy controls and 74 CRC cases.

Results

DNA yields ranged from 2 to 730 ng/mL plasma (mean 18.6ng/mL; 95% CI 11-26 ng/mL) and did not correlate with gender, age or CRC status. Methylated BCAT1 and IKZF1 DNA were detected in respectively 48 (65%) and 50 (68%) of the 74 cancers. In contrast, only 5 (4%) and 7 (5%) controls were positive for BCAT1 and IKZF1 DNA methylation, respectively. A two-gene classifier model (“either or” rule) improved segregation of CRC from controls, with 57 of 74 cancers (77%) compared to only 11 of 144 (7.6%) controls being positive for BCAT1 and/or IKZF1 DNA methylation. Increasing levels of methylated DNA were observed as CRC stage progressed.

Conclusions

Detection of methylated BCAT1 and/or IKZF1 DNA in plasma may have clinical application as a novel blood test for CRC. Combining the results from the two methylation-specific PCR assays improved CRC detection with minimal change in specificity. Further validation of this two-gene blood test with a view to application in screening is now indicated.     

Methylation-capture and Next-Generation Sequencing of free circulating DNA from human plasma

Background

Free circulating DNA (fcDNA) has many potential clinical applications, due to the non-invasive way in which it is collected. However, because of the low concentration of fcDNA in blood, genome-wide analysis carries many technical challenges that must be overcome before fcDNA studies can reach their full potential. There are currently no definitive standards for fcDNA collection, processing and whole-genome sequencing. We report novel detailed methodology for the capture of high-quality methylated fcDNA, library preparation and downstream genome-wide Next-Generation Sequencing. We also describe the effects of sample storage, processing and scaling on fcDNA recovery and quality.

Results

Use of serum versus plasma, and storage of blood prior to separation resulted in genomic DNA contamination, likely due to leukocyte lysis. Methylated fcDNA fragments were isolated from 5 donors using a methyl-binding protein-based protocol and appear as a discrete band of ~180 bases. This discrete band allows minimal sample loss at the size restriction step in library preparation for Next-Generation Sequencing, allowing for high-quality sequencing from minimal amounts of fcDNA. Following sequencing, we obtained 37×10⁶-86×10⁶ unique mappable reads, representing more than 50% of total mappable reads. The methylation status of 9 genomic regions as determined by DNA capture and sequencing was independently validated by clonal bisulphite sequencing.

Conclusions

Our optimized methods provide high-quality methylated fcDNA suitable for whole-genome sequencing, and allow good library complexity and accurate sequencing, despite using less than half of the recommended minimum input DNA.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-476) contains supplementary material, which is available to authorized users.     

Quilt Plots: A Simple Tool for the Visualisation of Large Epidemiological Data

Background

Graphical representation of data is one of the most easily comprehended forms of explanation. The current study describes a simple visualization tool which may allow greater understanding of medical and epidemiological data.

Method

We propose a simple tool for visualization of data, known as a “quilt plot”, that provides an alternative to presenting large volumes of data as frequency tables. Data from the Australian Needle and Syringe Program survey are used to illustrate “quilt plots”.

Conclusion

Visualization of large volumes of data using “quilt plots” enhances interpretation of medical and epidemiological data. Such intuitive presentations are particularly useful for the rapid assessment of problems in the data which cannot be readily identified by manual review. We recommend that, where possible, “quilt plots” be used along with traditional quantitative assessments of the data as an explanatory data analysis tool.     

Highly Mutagenic Exocyclic DNA Adducts Are Substrates for the Human Nucleotide Incision Repair Pathway

Background

Oxygen free radicals induce lipid peroxidation (LPO) that damages and breaks polyunsaturated fatty acids in cell membranes. LPO-derived aldehydes and hydroxyalkenals react with DNA leading to the formation of etheno(ε)-bases including 1,N ⁶-ethenoadenine (εA) and 3,N ⁴-ethenocytosine (εC). The εA and εC residues are highly mutagenic in mammalian cells and eliminated in the base excision repair (BER) pathway and/or by AlkB family proteins in the direct damage reversal process. BER initiated by DNA glycosylases is thought to be the major pathway for the removal of non-bulky endogenous base damage. Alternatively, in the nucleotide incision repair (NIR) pathway, the apurinic/apyrimidinic (AP) endonucleases can directly incise DNA duplex 5′ to a damaged base in a DNA glycosylase-independent manner.

Methodology/Principal Findings

Here we have characterized the substrate specificity of human major AP endonuclease 1, APE1, towards εA, εC, thymine glycol (Tg) and 7,8-dihydro-8-oxoguanine (8oxoG) residues when present in duplex DNA. APE1 cleaves oligonucleotide duplexes containing εA, εC and Tg, but not those containing 8oxoG. Activity depends strongly on sequence context. The apparent kinetic parameters of the reactions suggest that APE1 has a high affinity for DNA containing ε-bases but cleaves DNA duplexes at an extremely slow rate. Consistent with this observation, oligonucleotide duplexes containing an ε-base strongly inhibit AP site nicking activity of APE1 with IC₅₀ values in the range of 5–10 nM. MALDI-TOF MS analysis of the reaction products demonstrated that APE1-catalyzed cleavage of εA•T and εC•G duplexes generates, as expected, DNA fragments containing 5′-terminal ε-base residue.

Conclusions/Significance

The fact that ε-bases and Tg in duplex DNA are recognized and cleaved by APE1 in vitro, suggests that NIR may act as a backup pathway to BER to remove a large variety of genotoxic base lesions in human cells.     

Cost-Effectiveness of a Community Pharmacist-Led Sleep Apnea Screening Program – A Markov Model

Background

Despite the high prevalence and major public health ramifications, obstructive sleep apnea syndrome (OSAS) remains underdiagnosed. In many developed countries, because community pharmacists (CP) are easily accessible, they have been developing additional clinical services that integrate the services of and collaborate with other healthcare providers (general practitioners (GPs), nurses, etc.). Alternative strategies for primary care screening programs for OSAS involving the CP are discussed.

Objective

To estimate the quality of life, costs, and cost-effectiveness of three screening strategies among patients who are at risk of having moderate to severe OSAS in primary care.

Design

Markov decision model.

Data Sources

Published data.

Target Population

Hypothetical cohort of 50-year-old male patients with symptoms highly evocative of OSAS.

Time Horizon

The 5 years after initial evaluation for OSAS.

Perspective

Societal.

Interventions

Screening strategy with CP (CP-GP collaboration), screening strategy without CP (GP alone) and no screening.

Outcomes measures

Quality of life, survival and costs for each screening strategy.

Results of base-case analysis

Under almost all modeled conditions, the involvement of CPs in OSAS screening was cost effective. The maximal incremental cost for “screening strategy with CP” was about 455€ per QALY gained.

Results of sensitivity analysis

Our results were robust but primarily sensitive to the treatment costs by continuous positive airway pressure, and the costs of untreated OSAS. The probabilistic sensitivity analysis showed that the “screening strategy with CP” was dominant in 80% of cases. It was more effective and less costly in 47% of cases, and within the cost-effective range (maximum incremental cost effectiveness ratio at €6186.67/QALY) in 33% of cases.

Conclusions

CP involvement in OSAS screening is a cost-effective strategy. This proposal is consistent with the trend in Europe and the United States to extend the practices and responsibilities of the pharmacist in primary care.     

Extremely high Tp53 mutation load in esophageal squamous cell carcinoma in Golestan Province, Iran

Background

Golestan Province in northeastern Iran has one of the highest incidences of esophageal squamous cell carcinoma (ESCC) in the world with rates over 50 per 100,000 person-years in both sexes. We have analyzed TP53 mutation patterns in tumors from this high-risk geographic area in search of clues to the mutagenic processes involved in causing ESCC.

Methodology/Principal Findings

Biopsies of 119 confirmed ESCC tumor tissue from subjects enrolled in a case-control study conducted in Golestan Province were analyzed by direct sequencing of TP53 exons 2 through 11. Immunohistochemical staining for p53 was carried out using two monoclonal antibodies, DO7 and 1801. A total of 120 TP53 mutations were detected in 107/119 cases (89.9%), including 11 patients with double or triple mutations. The mutation pattern was heterogeneous with infrequent mutations at common TP53 “hotspots” but frequent transversions potentially attributable to environmental carcinogens forming bulky DNA adducts, including 40% at bases known as site of mutagenesis by polycyclic aromatic hydrocarbons (PAHs). Mutations showed different patterns according to the reported temperature of tea consumption, but no variation was observed in relation to ethnicity, tobacco or opium use, and alcoholic beverage consumption or urban versus rural residence.

Conclusion/Significance

ESCC tumors in people from Golestan Province show the highest rate of TP53 mutations ever reported in any cancer anywhere. The heterogeneous mutation pattern is highly suggestive of a causative role for multiple environmental carcinogens, including PAHs. The temperature and composition of tea may also influence mutagenesis.     

Rhomboids of Mycobacteria: Characterization Using an aarA Mutant of Providencia stuartii and Gene Deletion in Mycobacterium smegmatis

Background

Rhomboids are ubiquitous proteins with unknown roles in mycobacteria. However, bioinformatics suggested putative roles in DNA replication pathways and metabolite transport. Here, mycobacterial rhomboid-encoding genes were characterized; first, using the Providencia stuartii null-rhomboid mutant and then deleted from Mycobacterium smegmatis for additional insight in mycobacteria.

Methodology/Principal Findings

Using in silico analysis we identified in M. tuberculosis genome the genes encoding two putative rhomboid proteins; Rv0110 (referred to as “rhomboid protease 1”) and Rv1337 (“rhomboid protease 2”). Genes encoding orthologs of these proteins are widely represented in all mycobacterial species. When transformed into P. stuartii null-rhomboid mutant (ΔaarA), genes encoding mycobacterial orthologs of “rhomboid protease 2” fully restored AarA activity (AarA is the rhomboid protein of P. stuartii). However, most genes encoding mycobacterial “rhomboid protease 1” orthologs did not. Furthermore, upon gene deletion in M. smegmatis, the ΔMSMEG_4904 single mutant (which lost the gene encoding MSMEG_4904, orthologous to Rv1337, “rhomboid protease 2”) formed the least biofilms and was also more susceptible to ciprofloxacin and novobiocin, antimicrobials that inhibit DNA gyrase. However, the ΔMSMEG_5036 single mutant (which lost the gene encoding MSMEG_5036, orthologous to Rv0110, “rhomboid protease 1”) was not as susceptible. Surprisingly, the double rhomboid mutant ΔMSMEG_4904–ΔMSMEG_5036 (which lost genes encoding both homologs) was also not as susceptible suggesting compensatory effects following deletion of both rhomboid-encoding genes. Indeed, transforming the double mutant with a plasmid encoding MSMEG_5036 produced phenotypes of the ΔMSMEG_4904 single mutant (i.e. susceptibility to ciprofloxacin and novobiocin).

Conclusions/Significance

Mycobacterial rhomboid-encoding genes exhibit differences in complementing aarA whereby it''s only genes encoding “rhomboid protease 2” orthologs that fully restore AarA activity. Additionally, gene deletion data suggests inhibition of DNA gyrase by MSMEG_4904; however, the ameliorated effect in the double mutant suggests occurrence of compensatory mechanisms following deletion of genes encoding both rhomboids.     

Novel AlkB dioxygenases--alternative models for in silico and in vivo studies

Background

ALKBH proteins, the homologs of Escherichia coli AlkB dioxygenase, constitute a direct, single-protein repair system, protecting cellular DNA and RNA against the cytotoxic and mutagenic activity of alkylating agents, chemicals significantly contributing to tumor formation and used in cancer therapy. In silico analysis and in vivo studies have shown the existence of AlkB homologs in almost all organisms. Nine AlkB homologs (ALKBH1–8 and FTO) have been identified in humans. High ALKBH levels have been found to encourage tumor development, questioning the use of alkylating agents in chemotherapy. The aim of this work was to assign biological significance to multiple AlkB homologs by characterizing their activity in the repair of nucleic acids in prokaryotes and their subcellular localization in eukaryotes.

Methodology and Findings

Bioinformatic analysis of protein sequence databases identified 1943 AlkB sequences with eight new AlkB subfamilies. Since Cyanobacteria and Arabidopsis thaliana contain multiple AlkB homologs, they were selected as model organisms for in vivo research. Using E. coli alkB ⁻ mutant and plasmids expressing cyanobacterial AlkBs, we studied the repair of methyl methanesulfonate (MMS) and chloroacetaldehyde (CAA) induced lesions in ssDNA, ssRNA, and genomic DNA. On the basis of GFP fusions, we investigated the subcellular localization of ALKBHs in A. thaliana and established its mostly nucleo-cytoplasmic distribution. Some of the ALKBH proteins were found to change their localization upon MMS treatment.

Conclusions

Our in vivo studies showed highly specific activity of cyanobacterial AlkB proteins towards lesions and nucleic acid type. Subcellular localization and translocation of ALKBHs in A. thaliana indicates a possible role for these proteins in the repair of alkyl lesions. We hypothesize that the multiplicity of ALKBHs is due to their involvement in the metabolism of nucleo-protein complexes; we find their repair by ALKBH proteins to be economical and effective alternative to degradation and de novo synthesis.     

Cofactor Requirement of HpyAV Restriction Endonuclease

Background

Helicobacter pylori is the etiologic agent of common gastritis and a risk factor for gastric cancer. It is also one of the richest sources of Type II restriction-modification (R-M) systems in microorganisms.

Principal Findings

We have cloned, expressed and purified a new restriction endonuclease HpyAV from H. pylori strain 26695. We determined the HpyAV DNA recognition sequence and cleavage site as CCTTC 6/5. In addition, we found that HpyAV has a unique metal ion requirement: its cleavage activity is higher with transition metal ions than in Mg⁺⁺. The special metal ion requirement of HpyAV can be attributed to the presence of a HNH catalytic site similar to ColE9 nuclease instead of the canonical PD-X-D/EXK catalytic site found in many other REases. Site-directed mutagenesis was carried out to verify the catalytic residues of HpyAV. Mutation of the conserved metal-binding Asn311 and His320 to alanine eliminated cleavage activity. HpyAV variant H295A displayed approximately 1% of wt activity.

Conclusions/Significance

Some HNH-type endonucleases have unique metal ion cofactor requirement for optimal activities. Homology modeling and site-directed mutagenesis confirmed that HpyAV is a member of the HNH nuclease family. The identification of catalytic residues in HpyAV paved the way for further engineering of the metal binding site. A survey of sequenced microbial genomes uncovered 10 putative R-M systems that show high sequence similarity to the HpyAV system, suggesting lateral transfer of a prototypic HpyAV-like R-M system among these microorganisms.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-480) contains supplementary material, which is available to authorized users.     

Genetic and Biochemical Characterization of Human AP Endonuclease 1 Mutants Deficient in Nucleotide Incision Repair Activity

Background

Human apurinic/apyrimidinic endonuclease 1 (APE1) is a key DNA repair enzyme involved in both base excision repair (BER) and nucleotide incision repair (NIR) pathways. In the BER pathway, APE1 cleaves DNA at AP sites and 3′-blocking moieties generated by DNA glycosylases. In the NIR pathway, APE1 incises DNA 5′ to a number of oxidatively damaged bases. At present, physiological relevance of the NIR pathway is fairly well established in E. coli, but has yet to be elucidated in human cells.

Methodology/Principal Finding

We identified amino acid residues in the APE1 protein that affect its function in either the BER or NIR pathway. Biochemical characterization of APE1 carrying single K98A, R185A, D308A and double K98A/R185A amino acid substitutions revealed that all mutants exhibited greatly reduced NIR and 3′→5′ exonuclease activities, but were capable of performing BER functions to some extent. Expression of the APE1 mutants deficient in the NIR and exonuclease activities reduced the sensitivity of AP endonuclease-deficient E. coli xth nfo strain to an alkylating agent, methylmethanesulfonate, suggesting that our APE1 mutants are able to repair AP sites. Finally, the human NIR pathway was fully reconstituted in vitro using the purified APE1, human flap endonuclease 1, DNA polymerase β and DNA ligase I proteins, thus establishing the minimal set of proteins required for a functional NIR pathway in human cells.

Conclusion/Significance

Taken together, these data further substantiate the role of NIR as a distinct and separable function of APE1 that is essential for processing of potentially lethal oxidative DNA lesions.     

Identifying Fishes through DNA Barcodes and Microarrays

Background

International fish trade reached an import value of 62.8 billion Euro in 2006, of which 44.6% are covered by the European Union. Species identification is a key problem throughout the life cycle of fishes: from eggs and larvae to adults in fisheries research and control, as well as processed fish products in consumer protection.

Methodology/Principal Findings

This study aims to evaluate the applicability of the three mitochondrial genes 16S rRNA (16S), cytochrome b (cyt b), and cytochrome oxidase subunit I (COI) for the identification of 50 European marine fish species by combining techniques of “DNA barcoding” and microarrays. In a DNA barcoding approach, neighbour Joining (NJ) phylogenetic trees of 369 16S, 212 cyt b, and 447 COI sequences indicated that cyt b and COI are suitable for unambiguous identification, whereas 16S failed to discriminate closely related flatfish and gurnard species. In course of probe design for DNA microarray development, each of the markers yielded a high number of potentially species-specific probes in silico, although many of them were rejected based on microarray hybridisation experiments. None of the markers provided probes to discriminate the sibling flatfish and gurnard species. However, since 16S-probes were less negatively influenced by the “position of label” effect and showed the lowest rejection rate and the highest mean signal intensity, 16S is more suitable for DNA microarray probe design than cty b and COI. The large portion of rejected COI-probes after hybridisation experiments (>90%) renders the DNA barcoding marker as rather unsuitable for this high-throughput technology.

Conclusions/Significance

Based on these data, a DNA microarray containing 64 functional oligonucleotide probes for the identification of 30 out of the 50 fish species investigated was developed. It represents the next step towards an automated and easy-to-handle method to identify fish, ichthyoplankton, and fish products.