A blackberry (Rubus L.) expressed sequence tag library for the development of simple sequence repeat markers

Background

The recent development of novel repeat-fruiting types of blackberry (Rubus L.) cultivars, combined with a long history of morphological marker-assisted selection for thornlessness by blackberry breeders, has given rise to increased interest in using molecular markers to facilitate blackberry breeding. Yet no genetic maps, molecular markers, or even sequences exist specifically for cultivated blackberry. The purpose of this study is to begin development of these tools by generating and annotating the first blackberry expressed sequence tag (EST) library, designing primers from the ESTs to amplify regions containing simple sequence repeats (SSR), and testing the usefulness of a subset of the EST-SSRs with two blackberry cultivars.

Results

A cDNA library of 18,432 clones was generated from expanding leaf tissue of the cultivar Merton Thornless, a progenitor of many thornless commercial cultivars. Among the most abundantly expressed of the 3,000 genes annotated were those involved with energy, cell structure, and defense. From individual sequences containing SSRs, 673 primer pairs were designed. Of a randomly chosen set of 33 primer pairs tested with two blackberry cultivars, 10 detected an average of 1.9 polymorphic PCR products.

Conclusion

This rate predicts that this library may yield as many as 940 SSR primer pairs detecting 1,786 polymorphisms. This may be sufficient to generate a genetic map that can be used to associate molecular markers with phenotypic traits, making possible molecular marker-assisted breeding to compliment existing morphological marker-assisted breeding in blackberry.     

GONOME: measuring correlations between GO terms and genomic positions

Background

The design of oligonucleotides and PCR primers for studying large genomes is complicated by the redundancy of sequences. The eukaryotic genomes are particularly difficult to study due to abundant repeats. The speed of most existing primer evaluation programs is not sufficient for large-scale experiments.

Results

In order to improve the efficiency and success rate of automatic primer/oligo design, we created a novel method which allows rapid masking of repeats in large sequence files, for example in eukaryotic genomes. It also allows the detection of all alternative binding sites of PCR primers and the prediction of PCR products. The new method was implemented in a collection of efficient programs, the GENOMEMASKER package. The performance of the programs was compared to other similar programs. We also modified the PRIMER3 program, to be able to design primers from lowercase-masked sequences.

Conclusion

The GENOMEMASKER package is able to mask the entire human genome for non-unique primers within 6 hours and find locations of all binding sites for 10 000 designed primer pairs within 10 minutes. Additionally, it predicts all alternative PCR products from large genomes for given primer pairs.     

OpenADAM: an open source genome-wide association data management system for Affymetrix SNP arrays

Background

The goal of DNA barcoding is to develop a species-specific sequence library for all eukaryotes. A 650 bp fragment of the cytochrome c oxidase 1 (CO1) gene has been used successfully for species-level identification in several animal groups. It may be difficult in practice, however, to retrieve a 650 bp fragment from archival specimens, (because of DNA degradation) or from environmental samples (where universal primers are needed).

Results

We used a bioinformatics analysis using all CO1 barcode sequences from GenBank and calculated the probability of having species-specific barcodes for varied size fragments. This analysis established the potential of much smaller fragments, mini-barcodes, for identifying unknown specimens. We then developed a universal primer set for the amplification of mini-barcodes. We further successfully tested the utility of this primer set on a comprehensive set of taxa from all major eukaryotic groups as well as archival specimens.

Conclusion

In this study we address the important issue of minimum amount of sequence information required for identifying species in DNA barcoding. We establish a novel approach based on a much shorter barcode sequence and demonstrate its effectiveness in archival specimens. This approach will significantly broaden the application of DNA barcoding in biodiversity studies.     

Rapid single step subcloning procedure by combined action of type II and type IIs endonucleases with ligase

Background

The subcloning of a DNA fragment from an entry vector into a destination vector is a routinely performed task in molecular biology labs.

Results

We here present a novel benchtop procedure to achieve rapid recombination into any destination vector of choice with the sole requirement of an endonuclease recognition site. The method relies on a specifically designed entry vector and the combined action of type II and type IIs endonucleases with ligase. The formulation leads to accumulation of a single stable cloning product representing the desired insert carrying destination vector.

Conclusion

The described method provides a fast single step procedure for routine subcloning from an entry vector into a series of destination vectors with the same restriction enzyme recognition site.     

Boolean network simulations for life scientists

Background

We present a C++ class library for Monte Carlo simulation of molecular systems, including proteins in solution. The design is generic and highly modular, enabling multiple developers to easily implement additional features. The statistical mechanical methods are documented by extensive use of code comments that – subsequently – are collected to automatically build a web-based manual.

Results

We show how an object oriented design can be used to create an intuitively appealing coding framework for molecular simulation. This is exemplified in a minimalistic C++ program that can calculate protein protonation states. We further discuss performance issues related to high level coding abstraction.

Conclusion

C++ and the Standard Template Library (STL) provide a high-performance platform for generic molecular modeling. Automatic generation of code documentation from inline comments has proven particularly useful in that no separate manual needs to be maintained.     

Recodon: Coalescent simulation of coding DNA sequences with recombination, migration and demography

Background

Coalescent simulations have proven very useful in many population genetics studies. In order to arrive to meaningful conclusions, it is important that these simulations resemble the process of molecular evolution as much as possible. To date, no single coalescent program is able to simulate codon sequences sampled from populations with recombination, migration and growth.

Results

We introduce a new coalescent program, called Recodon, which is able to simulate samples of coding DNA sequences under complex scenarios in which several evolutionary forces can interact simultaneously (namely, recombination, migration and demography). The basic codon model implemented is an extension to the general time-reversible model of nucleotide substitution with a proportion of invariable sites and among-site rate variation. In addition, the program implements non-reversible processes and mixtures of different codon models.

Conclusion

Recodon is a flexible tool for the simulation of coding DNA sequences under realistic evolutionary models. These simulations can be used to build parameter distributions for testing evolutionary hypotheses using experimental data. Recodon is written in C, can run in parallel, and is freely available from http://darwin.uvigo.es/.     

Effect of X-ray irradiation on olive shoot culture evaluated by morphological measurements, nuclear DNA content and SSR and AFLP markers

Key message

This article presents the first investigation of X-ray irradiation treatments of olive shoot culture followed by detailed genetic analysis of induced mutation events. This was achieved by comparing morphological measurements, studies of nuclear DNA content and molecular marker (SSR and AFLP) analysis.

Abstract

Traditional olive cultivars often need minor phenotypic corrections, which can be assessed using mutation breeding. We investigated the effects of X-ray irradiation at 0, 10, 30 and 60 Gy delivered to in vitro grown shoots of the olive variety ‘Canino’. Morphological measurements showed that growth parameters did not decrease at a 10-Gy dose, that rooting ability was diminished at a 30-Gy dose and that shoot growth was inhibited at a 60-Gy dose. Measurements of acclimatized plants showed that nuclear DNA content decreased with irradiation dose from 2.972 to 2.963 and 2.935 pg/2C nucleus of control, for 10 and 30 Gy treated plants, respectively. SSR profiling using 6 microsatellite loci revealed no polymorphism. AFLP analysis generated in total 630 scorable fragments within 13 primer combinations, of which 179 (28.4 %) were polymorphic. Although no major morphological changes were observed at a 10-Gy dose, a high number of absent/novel fragments were detected by AFLP analysis, indicating that mutation events can also be detected for doses at which growth parameters are not affected.     

Blocking primers to enhance PCR amplification of rare sequences in mixed samples – a case study on prey DNA in Antarctic krill stomachs

Background

Identification of DNA sequence diversity is a powerful means for assessing the species present in environmental samples. The most common molecular strategies for estimating taxonomic composition depend upon PCR with universal primers that amplify an orthologous DNA region from a range of species. The diversity of sequences within a sample that can be detected by universal primers is often compromised by high concentrations of some DNA templates. If the DNA within the sample contains a small number of sequences in relatively high concentrations, then less concentrated sequences are often not amplified because the PCR favours the dominant DNA types. This is a particular problem in molecular diet studies, where predator DNA is often present in great excess of food-derived DNA.

Results

We have developed a strategy where a universal PCR simultaneously amplifies DNA from food items present in DNA purified from stomach samples, while the predator's own DNA is blocked from amplification by the addition of a modified predator-specific blocking primer. Three different types of modified primers were tested out; one annealing inhibiting primer overlapping with the 3' end of one of the universal primers, another annealing inhibiting primer also having an internal modification of five dI molecules making it a dual priming oligo, and a third elongation arrest primer located between the two universal primers. All blocking primers were modified with a C3 spacer. In artificial PCR mixtures, annealing inhibiting primers proved to be the most efficient ones and this method reduced predator amplicons to undetectable levels even when predator template was present in 1000 fold excess of the prey template. The prey template then showed strong PCR amplification where none was detectable without the addition of blocking primer. Our method was applied to identifying the winter food of one of the most abundant animals in the world, the Antarctic krill, Euphausia superba. Dietary item DNA was PCR amplified from a range of species in krill stomachs for which we had no prior sequence knowledge.

Conclusion

We present a simple, robust and cheap method that is easily adaptable to many situations where a rare DNA template is to be PCR amplified in the presence of a higher concentration template with identical PCR primer binding sites.     

An algorithm and program for finding sequence specific oligo-nucleotide probes for species identification

Background

The identification of species or species groups with specific oligo-nucleotides as molecular signatures is becoming increasingly popular for bacterial samples. However, it shows also great promise for other small organisms that are taxonomically difficult to tract.

Results

We have devised here an algorithm that aims to find the optimal probes for any given set of sequences. The program requires only a crude alignment of these sequences as input and is optimized for performance to deal also with very large datasets. The algorithm is designed such that the position of mismatches in the probes influences the selection and makes provision of single nucleotide outloops. Program implementations are available for Linux and Windows.     

Redasoft Plasmid 1.1: software for easy, efficient cloning and map drawing

A computer program is described which creates circular and linear genetic maps with or without sequences, performs restriction analysis, and simulates basic molecular cloning operations. Redasoft Plasmid 1.1 is a user-friendly program for Windows 95/98/NT/2000, designed to generate high quality genetic maps for presentation and publication and to aid the molecular biologist in selecting restriction endonucleases for common molecular cloning experiments. The program incorporates an integrated web browser and can automatically generate complete, labeled maps from sequences on the Internet.     

Design of a novel quantitative PCR (QPCR)-based protocol for genotyping mice carrying the neuroprotective Wallerian degeneration slow (Wld s ) gene

Background

The strain of MeCP2^tm1.1Bird mice is a broadly used model for Rett syndrome. Because males carrying the invalidated MeCP2 locus are sterile, this strain has to be maintained in a heterozygous state. All animals therefore have to be genotyped at every generation to discriminate those carrying the invalidated allele (+/- females and y/- males) from those that do not. This is conveniently carried out by PCR on tail genomic DNA but because the primer pairs described initially for this purpose yield very similar size DNA bands on the WT and the KO alleles, this requires to carry out two independent PCR reactions on tail DNA preparations from all animals.

Results

After cloning and sequencing the PCR fragment amplified on the KO allele, we tested several sets of primers that were designed to yield PCR fragments of different sizes on the KO and WT alleles.

Conclusion

We have thus identified a set of three primers that allows for efficient genotyping of the animals by a single PCR reaction. Furthermore, using of this set of primers also resolves a recurrent problem related to the tendency of one of the initial primers to give rise to a non specific band because of its capacity to anneal at both ends of a repeated genomic element which we have identified as MurvyLTR.     

Genome-wide characterization of the biggest grass, bamboo, based on 10,608 putative full-length cDNA sequences

Background

With the availability of rice and sorghum genome sequences and ongoing efforts to sequence genomes of other cereal and energy crops, the grass family (Poaceae) has become a model system for comparative genomics and for better understanding gene and genome evolution that underlies phenotypic and ecological divergence of plants. While the genomic resources have accumulated rapidly for almost all major lineages of grasses, bamboo remains the only large subfamily of Poaceae with little genomic information available in databases, which seriously hampers our ability to take a full advantage of the wealth of grass genomic data for effective comparative studies.

Results

Here we report the cloning and sequencing of 10,608 putative full length cDNAs (FL-cDNAs) primarily from Moso bamboo, Phyllostachys heterocycla cv. pubescens, a large woody bamboo with the highest ecological and economic values of all bamboos. This represents the third largest FL-cDNA collection to date of all plant species, and provides the first insight into the gene and genome structures of bamboos. We developed a Moso bamboo genomic resource database that so far contained the sequences of 10,608 putative FL-cDNAs and nearly 38,000 expressed sequence tags (ESTs) generated in this study.

Conclusion

Analysis of FL-cDNA sequences show that bamboo diverged from its close relatives such as rice, wheat, and barley through an adaptive radiation. A comparative analysis of the lignin biosynthesis pathway between bamboo and rice suggested that genes encoding caffeoyl-CoA O-methyltransferase may serve as targets for genetic manipulation of lignin content to reduce pollutants generated from bamboo pulping.     

A BAC-based integrated linkage map of the silkworm Bombyx mori

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.     

A diploid wheat TILLING resource for wheat functional genomics

Background

Cultivated peanut (Arachis hypogaea L.) is an important crop worldwide, valued for its edible oil and digestible protein. It has a very narrow genetic base that may well derive from a relatively recent single polyploidization event. Accordingly molecular markers have low levels of polymorphism and the number of polymorphic molecular markers available for cultivated peanut is still limiting.

Results

Here, we report a large set of BAC-end sequences (BES), use them for developing SSR (BES-SSR) markers, and apply them in genetic linkage mapping. The majority of BESs had no detectable homology to known genes (49.5%) followed by sequences with similarity to known genes (44.3%), and miscellaneous sequences (6.2%) such as transposable element, retroelement, and organelle sequences. A total of 1,424 SSRs were identified from 36,435 BESs. Among these identified SSRs, dinucleotide (47.4%) and trinucleotide (37.1%) SSRs were predominant. The new set of 1,152 SSRs as well as about 4,000 published or unpublished SSRs were screened against two parents of a mapping population, generating 385 polymorphic loci. A genetic linkage map was constructed, consisting of 318 loci onto 21 linkage groups and covering a total of 1,674.4 cM, with an average distance of 5.3 cM between adjacent loci. Two markers related to resistance gene homologs (RGH) were mapped to two different groups, thus anchoring 1 RGH-BAC contig and 1 singleton.

Conclusions

The SSRs mined from BESs will be of use in further molecular analysis of the peanut genome, providing a novel set of markers, genetically anchoring BAC clones, and incorporating gene sequences into a linkage map. This will aid in the identification of markers linked to genes of interest and map-based cloning.     

FastGroup: A program to dereplicate libraries of 16S rDNA sequences

Background

Ribosomal 16S DNA sequences are an essential tool for identifying and classifying microbes. High-throughput DNA sequencing now makes it economically possible to produce very large datasets of 16S rDNA sequences in short time periods, necessitating new computer tools for analyses. Here we describe FastGroup, a Java program designed to dereplicate libraries of 16S rDNA sequences. By dereplication we mean to: 1) compare all the sequences in a data set to each other, 2) group similar sequences together, and 3) output a representative sequence from each group. In this way, duplicate sequences are removed from a library.

Results

FastGroup was tested using a library of single-pass, bacterial 16S rDNA sequences cloned from coral-associated bacteria. We found that the optimal strategy for dereplicating these sequences was to: 1) trim ambiguous bases from the 5' end of the sequences and all sequence 3' of the conserved Bact517 site, 2) match the sequences from the 3' end, and 3) group sequences >=97% identical to each other.

Conclusions

The FastGroup program simplifies the dereplication of 16S rDNA sequence libraries and prepares the raw sequences for subsequent analyses.     

A novel method for cloning of coding sequences of highly toxic proteins

Background

During standard gene cloning, the recombinant protein appearing in bacteria as the result of expression leakage very often inhibits cell proliferation leading to blocking of the cloning procedure. Although different approaches can reduce transgene basal expression, the recombinant proteins, which even in trace amounts inhibit bacterial growth, can completely prevent the cloning process.