BatchPrimer3: a high throughput web application for PCR and sequencing primer design

PrimerCE: designing primers for cloning and gene expression

A number of primer design programs have been developed for diverse applications. However, none of these programs can be used to design primers for gene cloning aimed at expressing protein. Here we report the design of PrimerCE, which can be used to cover the whole process of gene cloning and expression. The main features of PrimerCE include inspection of restriction enzyme recognition sequence, open reading frame verification, stop codon inspection, base adjustment, primer optimization, sequence assembly and protein analysis. In addition to this, the program can be modified based on the different needs of users, e.g. new vector sequence and restriction enzyme recognition sequences can be integrated. With the use of PrimerCE, a pair of primers can be designed within minutes. The program has been proven to be efficient in designing primers in our high throughput cloning and gene expression experiments. The software is freely available at http://tch.hebau.edu.cn/shengm/download/down.html.

Forensic and genetic characterization of mtDNA from Pathans of Pakistan

Complete mitochondrial control region data were generated for 230 unrelated Pathans from North West Frontier Province and Federally Administered Tribal Areas of Pakistan. To confirm data quality and to explore the genetic structure of Pathans, mitochondrial DNA haplogroup affiliation was determined by shared haplogroup-specific polymorphisms in the control region and by the analysis of diagnostic coding region single-nucleotide polymorphisms using a multiplex system for the assignment of eight haplogroups: M, N1'5, W, R, R0, T, J, and U. Sequence comparison revealed that 193 haplotypes were defined by 215 variable sites when major insertions were ignored at nucleotide positions 16193, 309, and 573. From a phylogenetic perspective, Pathans have a heterogeneous origin, displaying a high percentage of West Eurasian haplogroups followed by haplogroups native to South Asia and a small fraction from East Asian lineages. In population comparisons, this ethnic group differed significantly from several other ethnic groups from Pakistan and surrounding countries. These results suggest that frequency estimates for mtDNA haplotypes should be determined for endogamous ethnic groups individually instead of pooling data for these subpopulations into a single dataset for the Pakistani population. Data presented here may contribute to the accuracy of forensic mtDNA comparisons in the Pathans of Pakistan.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

A comparison of restriction fragment length polymorphism, tetra primer amplification refractory mutation system PCR and unlabeled probe melting analysis for LTA+252 C>T SNP genotyping

BACKGROUND

From the wide range of methods currently available for genotyping, we wished to identify a quick, reliable and affordable approach for routine use in our laboratory for LTA+252 C>T SNP screening.

METHODS

We set up and compared three genotyping methods for SNP detection: restriction fragment length polymorphism (RFLP), tetra primer amplification refractory mutation system PCR (TPAP) and unlabeled probe melting analysis (UPMA). The SNP model used was LTA+252 C>T, a cytokine gene polymorphism that has been associated with response to treatment in rheumatoid arthritis. The study was performed using 46 samples from healthy Caucasian volunteers.

RESULTS

Allele and genotype distribution was similar to that previously described in the same population. All three genotyping methods showed good reproducibility and are suitable for a medium scale throughput molecular platform. UPMA was the most cost effective, reliable and safe method since it required the shortest technician time, could be performed in a single closed tube and involved automatic data analysis.

CONCLUSION

This work is the first to compare these three genotyping techniques and provides evidence for UPMA being the method of choice for LTA+252 C>T SNP genotyping.

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.

Male sex interspecies divergence and down regulation of expression of spermatogenesis genes in Drosophila sterile hybrids

Male sex genes have shown a pattern of rapid interspecies divergence at both the coding and gene expression level. A common outcome from crosses between closely-related species is hybrid male sterility. Phenotypic and genetic studies in Drosophila sterile hybrid males have shown that spermatogenesis arrest is postmeiotic with few exceptions, and that most misregulated genes are involved in late stages of spermatogenesis. Comparative studies of gene regulation in sterile hybrids and parental species have mainly used microarrays providing a whole genome representation of regulatory problems in sterile hybrids. Real-time PCR studies can reject or reveal differences not observed in microarray assays. Moreover, differences in gene expression between samples can be dependant on the source of RNA (e.g., whole body vs. tissue). Here we survey expression in D. simulans, D. mauritiana and both intra and interspecies hybrids using a real-time PCR approach for eight genes expressed at the four main stages of sperm development. We find that all genes show a trend toward under expression in the testes of sterile hybrids relative to parental species with only the two proliferation genes (bam and bgcn) and the two meiotic class genes (can and sa) showing significant down regulation. The observed pattern of down regulation for the genes tested can not fully explain hybrid male sterility. We discuss the down regulation of spermatogenesis genes in hybrids between closely-related species within the contest of rapid divergence experienced by the male genome, hybrid sterility and possible allometric changes due to subtle testes-specific developmental abnormalities.

Rapid direct PCR for ABO blood typing

Many different molecular typing methods have been reported to complement routine serological ABO blood typing in forensics. However, these ABO genotyping methods are often time-consuming and call for an initial DNA isolation step that requires the use of expensive kits or reagents. We report here a rapid direct ABO genotyping method that eliminates the need for DNA extraction from fresh blood, hair, and body fluid stains before PCR. Using a fast PCR instrument and an optimized polymerase, the genotyping method-which employs a multiplex allele-specific primer set for the simultaneous detection of three single-nucleotide polymorphism (SNP) sites (nucleotides 261, 526, and 803)-identifies A, B, O01/O02, O03, and cis-AB01 alleles in around 70 min from sample collection to electropherogram. Not only will this ABO genotyping method be efficiently used in forensic practice for rapid screening of samples before full-blown multilocus short tandem repeat profiling, but it will also demonstrate an example of rapid direct genotyping of SNPs that offers the advantages of time- and cost-efficiency, convenience, and reduced contamination during DNA analysis.

A single primer pair gives a specific ortholog amplicon in a wide range of Cyanobacteria and plastid-bearing organisms: applicability in inventory of reference material from collections and phylogenetic analysis

The scarcity of universally applied molecular markers for algae has resulted in the development of multiple, independent and not easily comparable systems. The goal of this work is to increase the number of available molecular markers and to generate easily comparable systems. Thereby, we have designed a primer pair capable of amplifying a broad range of organisms: Cyanobacteria, Chlorophyta, Chlorarachniophyta, Cryptophyta, Euglenida, Glaucophyta, Rhodophyta, Stramenopiles and Streptophyta including plants. This primer pair can amplify a portion of the 23S rRNA gene with sufficient variability to identify reference material form collections across a broad range of taxa and perform phylogenetic studies alongside other available markers.

Temporal expression and steroidal regulation of piRNA pathway genes (mael, piwi, vasa) during Silurana (Xenopus) tropicalis embryogenesis and early larval development

It has been extensively documented that exposure of amphibians and teleost fish to exogenous steroid hormones like estrogen, androgen, xenoestrogen or steroid biosynthesis inhibitors can impair their gonadal development or induce sex reversal against genotypic sex. However, the molecular pathways underlying sexual development and the effects of sex steroids or other exogenous hormones in these aquatic vertebrates remain elusive. Recently, a germ plasm-associated piRNA (piwi-interacting RNA) pathway has been shown to be a determinant in the development of animal gonadal germline cells. In the current study, we examined whether this piRNA pathway is involved in the regulation of sex steroid hormones in gonadal development. We firstly established developmental expression patterns of three key piRNA pathway genes (mael, piwi and vasa), during Silurana (Xenopus) tropicalis embryogenesis and early larval development. All three genes exhibit high expression at early developmental stages and have significantly decreased expression thereafter, indicating a very active involvement of piRNA pathway at the beginning of embryogenesis. We further examined gene expression changes of those genes in frog larvae exposed to two sex steroid biosynthesis inhibitors, fadrozole and finasteride, both of which are known to result in male-biased or female-biased phenotypes, respectively. We found that fadrozole and finasteride exposures increased the expression of piRNA pathway genes such as mael and vasa at the larval stage when the expression of piRNA pathway genes is programmed to be very low. Therefore, our results indicate that the piRNA pathway is likely a common pathway by which different sex steroid hormones regulate gonadal sex differentiation.

Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers

BACKGROUND

Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption and the goal of this study was to integrate a recently constructed physical map for the species with a microsatellite based genetic map using a BAC library from the genotype G19833 and the recombinant inbred line population DOR364 x G19833.

RESULTS

We searched for simple sequence repeats (SSRs) in the 89,017 BAC-end sequences (BES) from the physical map and genetically mapped any polymorphic BES-SSRs onto the genetic map. Among the BES it was possible to identify 623 contig-linked SSRs, most of which were highly AT-rich. A subgroup of 230 di-nucleotide and tri-nucleotide based SSR primer pairs from these BACs was tested on the mapping parents with 176 single copy loci and 114 found to be polymorphic markers. Of these, 99 were successfully integrated into the genetic map. The 99 linkages between the genetic and physical maps corresponded to an equal number of contigs containing a total of 5,055 BAC clones.

CONCLUSIONS

Class II microsatellites were more common in the BES than longer class I microsatellites. Both types of markers proved to be valuable for linking BAC clones to the genetic map and were successfully placed across all 11 linkage groups. The integration of common bean physical and genetic maps is an important part of comparative genome analysis and a prelude to positional cloning of agronomically important genes for this crop.

Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers

Background

Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption and the goal of this study was to integrate a recently constructed physical map for the species with a microsatellite based genetic map using a BAC library from the genotype G19833 and the recombinant inbred line population DOR364 × G19833.

Results

We searched for simple sequence repeats (SSRs) in the 89,017 BAC-end sequences (BES) from the physical map and genetically mapped any polymorphic BES-SSRs onto the genetic map. Among the BES it was possible to identify 623 contig-linked SSRs, most of which were highly AT-rich. A subgroup of 230 di-nucleotide and tri-nucleotide based SSR primer pairs from these BACs was tested on the mapping parents with 176 single copy loci and 114 found to be polymorphic markers. Of these, 99 were successfully integrated into the genetic map. The 99 linkages between the genetic and physical maps corresponded to an equal number of contigs containing a total of 5,055 BAC clones.

Conclusions

Class II microsatellites were more common in the BES than longer class I microsatellites. Both types of markers proved to be valuable for linking BAC clones to the genetic map and were successfully placed across all 11 linkage groups. The integration of common bean physical and genetic maps is an important part of comparative genome analysis and a prelude to positional cloning of agronomically important genes for this crop.

Genetic mapping of microsatellite markers around the arcelin bruchid resistance locus in common bean

The deployment in common beans (Phaseolus vulgaris L.) of arcelin-based bruchid resistance could help reduce post-harvest storage losses to the Mexican bean weevil [(Zabrotes subfasciatus (Boheman)]. Arcelin is a member of the arcelin-phytohemagglutinin-alpha-amylase inhibitor (APA) family of seed proteins, which has been extensively studied but not widely used in bean breeding programs. The purpose of this study was to evaluate microsatellite markers for genetic analysis of arcelin-based bruchid resistance and to determine the orientation of markers and the rate of recombination around the APA locus. A total of 10 previously developed microsatellites and 22 newly developed markers based on a sequenced BAC from the APA locus were screened for polymorphism and of these 15 were mapped with an F(2) population of 157 individuals resulting from a susceptible x resistant cross of SEQ1006 x RAZ106 that segregated for both the arcelin 1 allele and resistance to the bruchid, Z. subfasciatus. Microsatellites derived from APA gene sequences were linked within 0.8 cM of each other and were placed relative to the rest of the b04 linkage group. In a comparison of genetic to physical distance on the BAC sequence, recombination was found to be moderate with a ratio of 125 kb/cM, but repressed within the APA locus itself. Several markers were predicted to be very effective for genetic studies or marker-assisted selection, based on their significant associations with bruchid resistance and on low adult insect emergence and positions flanking the arcelin and phytohemagglutinin genes.

PCRTiler: automated design of tiled and specific PCR primer pairs

Efficiency and specificity of PCR amplification is dependent on several parameters, such as amplicon length, as well as hybridization specificity and melting temperature of primer oligonucleotides. Primer design is thus of critical importance for the success of PCR experiments, but can be a time-consuming and repetitive task, for example when large genomic regions are to be scanned for the presence of a protein of interest by chromatin immunoprecipitation experiments. We present here a webserver that allows the automated design of tiled primer pairs for any number of genomic loci. PCRTiler splits the target DNA sequences into smaller regions, and identifies candidate primers for each sub-region by running the well-known program Primer3 followed by the elimination of primers with a high cross-hybridization potential via BLAST. Tiling density and primer characteristics are specified by the user via a simple and user-friendly interface. The webserver can be accessed at http://pcrtiler.alaingervais.org:8080/PCRTiler. Additionally, users may download a standalone Java-based implementation of this software. Experimental validation of PCRTiler has demonstrated that it produces correct results. We have tiled a region of the human genome, in which 96 of 123 primer pairs worked in the first attempt, and 105 of 123 (85%) could be made to work by optimizing the conditions of the PCR assay.

RJPrimers: unique transposable element insertion junction discovery and PCR primer design for marker development

Transposable elements (TE) exist in the genomes of nearly all eukaryotes. TE mobilization through ‘cut-and-paste’ or ‘copy-and-paste’ mechanisms causes their insertions into other repetitive sequences, gene loci and other DNA. An insertion of a TE commonly creates a unique TE junction in the genome. TE junctions are also randomly distributed along chromosomes and therefore useful for genome-wide marker development. Several TE-based marker systems have been developed and applied to genetic diversity assays, and to genetic and physical mapping. A software tool ‘RJPrimers’ reported here allows for accurate identification of unique repeat junctions using BLASTN against annotated repeat databases and a repeat junction finding algorithm, and then for fully automated high-throughput repeat junction-based primer design using Primer3 and BatchPrimer3. The software was tested using the rice genome and genomic sequences of Aegilops tauschii. Over 90% of repeat junction primers designed by RJPrimers were unique. At least one RJM marker per 10 Kb sequence of A. tauschii was expected with an estimate of over 0.45 million such markers in a genome of 4.02 Gb, providing an almost unlimited source of molecular markers for mapping large and complex genomes. A web-based server and a command line-based pipeline for RJPrimers are both available at http://wheat.pw.usda.gov/demos/RJPrimers/.

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

Background

Massively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.

Results

We obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.

Conclusions

This sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

BACKGROUND

Massively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.

RESULTS

We obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.

CONCLUSIONS

This sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Evidence of genetic instability in tumors and normal nearby tissues

BACKGROUND

Comprehensive analyses have recently been performed on many human cancer tissues, leading to the identification of a number of mutated genes but providing no information on the variety of mutations present in each of them. This information is of interest to understand the possible origin of gene mutations that cause tumors.

METHODOLOGY/PRINCIPAL FINDINGS

We have analyzed the sequence heterogeneity of the transcripts of the human HPRT and G6PD single copy genes that are not considered tumor markers. Analyses have been performed on different colon cancers and on the nearby histologically normal tissues of two male patients. Several copies of each cDNA, which were produced by cloning the RT-PCR-amplified fragments of the specific mRNA, have been sequenced. Similar analyses have been performed on blood samples of two ostensibly healthy males as reference controls. The sequence heterogeneity of the HPRT and G6PD genes was also determined on DNA from tumor tissues. The employed analytical approach revealed the presence of low-frequency mutations not detectable by other procedures. The results show that genetic heterogeneity is detectable in HPRT and G6PD transcripts in both tumors and nearby healthy tissues of the two studied colon tumors. Similar frequencies of mutations are observed in patient genomic DNA, indicating that mutations have a somatic origin. HPRT transcripts show genetic heterogeneity also in healthy individuals, in agreement with previous results on human T-cells, while G6PD transcript heterogeneity is a characteristic of the patient tissues. Interestingly, data on TP53 show little, if any, heterogeneity in the same tissues.

CONCLUSIONS/SIGNIFICANCE

These findings show that genetic heterogeneity is a peculiarity not only of cancer cells but also of the normal tissue where a tumor arises.

Oxytocin stimulates expression of a noncoding RNA tumor marker in a human neuroblastoma cell line

AIMS

A noncoding RNA, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), is upregulated in several malignant tumors. Its expression in neuroblastoma, however, is not known, and the regulatory mechanisms of MALAT1 gene expression have not been elucidated. The aim of this study is to clarify how MALAT1 gene expression is altered by extracellular signals in the SK-N-SH neuroblastoma cell line and to define its proximal promoter in order to study the mechanism of MALAT1 gene expression.

METHODS

Transcript amounts were analyzed by real-time semiquantitative polymerase chain reaction (qPCR). Genes coregulated with MALAT1 were identified by DNA microarray analysis. The structure of the MALAT1 transcript was delineated using a tiling microarray, and the 5'-end was determined using the rapid amplification of cDNA ends (RACE) method. We investigated binding of the cyclic AMP-responsive element binding (CREB) transcription factor to the MALAT1 promoter by using chromatin immunoprecipitation (ChIP) followed by tiling array analysis, and the results were confirmed using ChIP-qPCR.

KEY FINDINGS

The posterior pituitary hormone oxytocin increased the levels of MALAT1 and immediate early gene transcripts as early as 15 min after stimulation. Although the expression of immediate early genes returned to basal levels after 3h, MALAT1 transcript levels peaked 6-24h after stimulation. We identified a shorter transcriptional initiation site and found that CREB binds to the defined proximal promoter of the MALAT1 gene.

SIGNIFICANCE

The expression of the tumor marker MALAT1 ncRNA is sensitive to cell surface receptor activation by oxytocin in a neuroblastoma cell line.

Inhibition of heat shock protein 90 attenuates adenylate cyclase sensitization after chronic morphine treatment

Cellular adaptations to chronic opioid treatment result in enhanced responsiveness of adenylate cyclase and an increase in forskolin- or agonist-stimulated cAMP production. It is, however, not known whether chaperone molecules such as heat shock proteins contribute to this adenylate cyclase sensitization. Here, we report that treatment of cells with geldanamycin, an inhibitor of heat shock protein 90 (Hsp90), led to effective attenuation of morphine-induced adenylate cyclase sensitization. In SK-N-SH human neuroblastoma cells, morphine significantly increased RNA transcript and protein levels of type I adenylate cyclase, leading to sensitization. Whole-genome tiling array analysis revealed that cAMP response element-binding protein, an important mediator for cellular adaptation to morphine, associated with the proximal promoter of Hsp90AB1 not only in SK-N-SH cells but also in rat PC12 and human embryonic kidney cells. Hsp90AB1 transcript and protein levels increased significantly during morphine treatment, and co-application of geldanamycin (0.1-10 nM) effectively suppressed the increase in forskolin-activated adenylate cyclase activation by 56%. Type I adenylate cyclase, but not Hsp90AB1, underwent significant degradation during geldanamycin treatment. These results indicate that Hsp90 is a new pharmacological target for the suppression of adenylate cyclase sensitization induced by chronic morphine treatment.

Microsatellites in Brassica unigenes: relative abundance, marker design, and use in comparative physical mapping and genome analysis

Microsatellites present in the transcribed regions of the genome have the potential to reveal functional diversity. Unigene sequence databases are the sources of such genic microsatellites with unique flanking sequences and genomic locations even in complex polyploids. The present study was designed to assay the unigenes of Brassica napus and B. rapa for various microsatellite repeats, and to design markers and use them in comparative genome analysis and study of evolution. The average frequency of microsatellites in Brassica unigenes was one in every 7.25 kb of sequence, as compared with one in every 8.57 kb of sequence in Arabidopsis thaliana . Trinucleotide motifs coding for serine and the dinucleotide motif GA were most abundant. We designed 2374 and 347 unigene-based microsatellite (UGMS) markers including 541 and 58 class I types in B. napus and B. rapa, respectively, and evaluated their use across diverse species and genera. Most of these markers (93.3%) gave successful amplification of target microsatellite motifs, which was confirmed by sequencing. Interspecific polymorphism between B. napus and B. rapa detected in silico for the UGMS markers was 4.16 times higher in 5′ untranslated regions than in coding sequences. Physical anchoring of Brassica UGMS markers on the A. thaliana genome indicated their significance in studying the evolutionary history of A. thaliana genomic duplications in relation to speciation. Comparative physical mapping identified 85% of Brassica unigenes as single copy and gave clues for the presence of conserved primordial gene order. Complex chromosomal rearrangements such as inversions, tandem and segmental duplications, and insertions/deletions were evident between A. thaliana and B. rapa genomes. The results obtained have encouraging implications for the use of UGMS markers in comparative genome analysis and for understanding evolutionary complexities in the family Brassicaceae.

PrimerBank: a resource of human and mouse PCR primer pairs for gene expression detection and quantification

PrimerBank (http://pga.mgh.harvard.edu/primerbank/) is a public resource for the retrieval of human and mouse primer pairs for gene expression analysis by PCR and Quantitative PCR (QPCR). A total of 306,800 primers covering most known human and mouse genes can be accessed from the PrimerBank database, together with information on these primers such as T(m), location on the transcript and amplicon size. For each gene, at least one primer pair has been designed and in many cases alternative primer pairs exist. Primers have been designed to work under the same PCR conditions, thus facilitating high-throughput QPCR. There are several ways to search for primers for the gene(s) of interest, such as by: GenBank accession number, NCBI protein accession number, NCBI gene ID, PrimerBank ID, NCBI gene symbol or gene description (keyword). In all, 26,855 primer pairs covering most known mouse genes have been experimentally validated by QPCR, agarose gel analysis, sequencing and BLAST, and all validation data can be freely accessed from the PrimerBank web site.

ConservedPrimers 2.0: a high-throughput pipeline for comparative genome referenced intron-flanking PCR primer design and its application in wheat SNP discovery

Background

In some genomic applications it is necessary to design large numbers of PCR primers in exons flanking one or several introns on the basis of orthologous gene sequences in related species. The primer pairs designed by this target gene approach are called "intron-flanking primers" or because they are located in exonic sequences which are usually conserved between related species, "conserved primers". They are useful for large-scale single nucleotide polymorphism (SNP) discovery and marker development, especially in species, such as wheat, for which a large number of ESTs are available but for which genome sequences and intron/exon boundaries are not available. To date, no suitable high-throughput tool is available for this purpose.

Results

We have developed, the ConservedPrimers 2.0 pipeline, for designing intron-flanking primers for large-scale SNP discovery and marker development, and demonstrated its utility in wheat. This tool uses non-redundant wheat EST sequences, such as wheat contigs and singleton ESTs, and related genomic sequences, such as those of rice, as inputs. It aligns the ESTs to the genomic sequences to identify unique colinear exon blocks and predicts intron lengths. Intron-flanking primers are then designed based on the intron/exon information using the Primer3 core program or BatchPrimer3. Finally, a tab-delimited file containing intron-flanking primer pair sequences and their primer properties is generated for primer ordering and their PCR applications. Using this tool, 1,922 bin-mapped wheat ESTs (31.8% of the 6,045 in total) were found to have unique colinear exon blocks suitable for primer design and 1,821 primer pairs were designed from these single- or low-copy genes for PCR amplification and SNP discovery. With these primers and subsequently designed genome-specific primers, a total of 1,527 loci were found to contain one or more genome-specific SNPs.

Conclusion

The ConservedPrimers 2.0 pipeline for designing intron-flanking primers was developed and its utility demonstrated. The tool can be used for SNP discovery, genetic variation assays and marker development for any target genome that has abundant ESTs and a related reference genome that has been fully sequenced. The ConservedPrimers 2.0 pipeline has been implemented as a command-line tool as well as a web application. Both versions are freely available at .

Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding

BACKGROUND

Cultivated tomato (Solanum lycopersicum L.) has narrow genetic diversity that makes it difficult to identify polymorphisms between elite germplasm. We explored array-based single feature polymorphism (SFP) discovery as a high-throughput approach for marker development in cultivated tomato.

RESULTS

Three varieties, FL7600 (fresh-market), OH9242 (processing), and PI114490 (cherry) were used as a source of genomic DNA for hybridization to oligonucleotide arrays. Identification of SFPs was based on outlier detection using regression analysis of normalized hybridization data within a probe set for each gene. A subset of 189 putative SFPs was sequenced for validation. The rate of validation depended on the desired level of significance (alpha) used to define the confidence interval (CI), and ranged from 76% for polymorphisms identified at alpha <or= 10-6 to 60% for those identified at alpha <or= 10-2. Validation percentage reached a plateau between alpha <or= 10-4 and alpha <or= 10-7, but failure to identify known SFPs (Type II error) increased dramatically at alpha <or= 10-6. Trough sequence validation, we identified 279 SNPs and 27 InDels in 111 loci. Sixty loci contained >or= 2 SNPs per locus. We used a subset of validated SNPs for genetic diversity analysis of 92 tomato varieties and accessions. Pairwise estimation of theta (Fst) suggested significant differentiation between collections of fresh-market, processing, vintage, Latin American (landrace), and S. pimpinellifolium accessions. The fresh-market and processing groups displayed high genetic diversity relative to vintage and landrace groups. Furthermore, the patterns of SNP variation indicated that domestication and early breeding practices have led to progressive genetic bottlenecks while modern breeding practices have reintroduced genetic variation into the crop from wild species. Finally, we examined the ratio of non-synonymous (Ka) to synonymous substitutions (Ks) for 20 loci with multiple SNPs (>or= 4 per locus). Six of 20 loci showed ratios of Ka/Ks >or= 0.9.

CONCLUSION

Array-based SFP discovery was an efficient method to identify a large number of molecular markers for genetics and breeding in elite tomato germplasm. Patterns of sequence variation across five major tomato groups provided insight into to the effect of human selection on genetic variation.

Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding

BACKGROUND

Cultivated tomato (Solanum lycopersicum L.) has narrow genetic diversity that makes it difficult to identify polymorphisms between elite germplasm. We explored array-based single feature polymorphism (SFP) discovery as a high-throughput approach for marker development in cultivated tomato.

RESULTS

Three varieties, FL7600 (fresh-market), OH9242 (processing), and PI114490 (cherry) were used as a source of genomic DNA for hybridization to oligonucleotide arrays. Identification of SFPs was based on outlier detection using regression analysis of normalized hybridization data within a probe set for each gene. A subset of 189 putative SFPs was sequenced for validation. The rate of validation depended on the desired level of significance (alpha) used to define the confidence interval (CI), and ranged from 76% for polymorphisms identified at alpha <or= 10-6 to 60% for those identified at alpha <or= 10-2. Validation percentage reached a plateau between alpha <or= 10-4 and alpha <or= 10-7, but failure to identify known SFPs (Type II error) increased dramatically at alpha <or= 10-6. Trough sequence validation, we identified 279 SNPs and 27 InDels in 111 loci. Sixty loci contained >or= 2 SNPs per locus. We used a subset of validated SNPs for genetic diversity analysis of 92 tomato varieties and accessions. Pairwise estimation of theta (Fst) suggested significant differentiation between collections of fresh-market, processing, vintage, Latin American (landrace), and S. pimpinellifolium accessions. The fresh-market and processing groups displayed high genetic diversity relative to vintage and landrace groups. Furthermore, the patterns of SNP variation indicated that domestication and early breeding practices have led to progressive genetic bottlenecks while modern breeding practices have reintroduced genetic variation into the crop from wild species. Finally, we examined the ratio of non-synonymous (Ka) to synonymous substitutions (Ks) for 20 loci with multiple SNPs (>or= 4 per locus). Six of 20 loci showed ratios of Ka/Ks >or= 0.9.

CONCLUSION

Array-based SFP discovery was an efficient method to identify a large number of molecular markers for genetics and breeding in elite tomato germplasm. Patterns of sequence variation across five major tomato groups provided insight into to the effect of human selection on genetic variation.

Development of SSR markers and analysis of diversity in Turkish populations of Brachypodium distachyon

BACKGROUND

Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to facilitate research aimed at improving grass crops for grain, forage and energy production. To characterize the natural diversity of Brachypodium and provide a valuable new tool to the growing list of resources available to Brachypodium researchers, we created and characterized a large, diverse collection of inbred lines.

RESULTS

We developed 84 inbred lines from eight locations in Turkey. To enable genotypic characterization of this collection, we created 398 SSR markers from BAC end and EST sequences. An analysis of 187 diploid lines from 56 locations with 43 SSR markers showed considerable genotypic diversity. There was some correlation between SSR genotypes and broad geographic regions, but there was also a high level of genotypic diversity at individual locations. Phenotypic analysis of this new germplasm resource revealed considerable variation in flowering time, seed size, and plant architecture. The inbreeding nature of Brachypodium was confirmed by an extremely high level of homozygosity in wild plants and a lack of cross-pollination under laboratory conditions.

CONCLUSION

Taken together, the inbreeding nature and genotypic diversity observed at individual locations suggest a significant amount of long-distance seed dispersal. The resources developed in this study are freely available to the research community and will facilitate experimental applications based on natural diversity.

Seq4SNPs: new software for retrieval of multiple, accurately annotated DNA sequences, ready formatted for SNP assay design

Background

In moderate-throughput SNP genotyping there was a gap in the workflow, between choosing a set of SNPs and submitting their sequences to proprietary assay design software, which was not met by existing software. Retrieval and formatting of sequences flanking each SNP, prior to assay design, becomes rate-limiting for more than about ten SNPs, especially if annotated for repetitive regions and adjacent variations. We routinely process up to 50 SNPs at once.

Implementation

We created Seq4SNPs, a web-based, walk-away software that can process one to several hundred SNPs given rs numbers as input. It outputs a file of fully annotated sequences formatted for one of three proprietary design softwares: TaqMan's Primer-By-Design FileBuilder, Sequenom's iPLEX or SNPstream's Autoprimer, as well as unannotated fasta sequences. We found genotyping assays to be inhibited by repetitive sequences or the presence of additional variations flanking the SNP under test, and in multiplexes, repetitive sequence flanking one SNP adversely affects multiple assays. Assay design software programs avoid such regions if the input sequences are appropriately annotated, so we used Seq4SNPs to provide suitably annotated input sequences, and improved our genotyping success rate. Adjacent SNPs can also be avoided, by annotating sequences used as input for primer design.

Conclusion

The accuracy of annotation by Seq4SNPs is significantly better than manual annotation (P < 1e-5).

Using Seq4SNPs to incorporate all annotation for additional SNPs and repetitive elements into sequences, for genotyping assay designer software, minimizes assay failure at the design stage, reducing the cost of genotyping. Seq4SNPs provides a rapid route for replacement of poor test SNP sequences. We routinely use this software for assay sequence preparation.

Seq4SNPs is available as a service at and , currently for human SNPs, but easily extended to include any species in dbSNP.

Defining global neuroendocrine gene expression patterns associated with reproductive seasonality in fish

BACKGROUND

Many vertebrates, including the goldfish, exhibit seasonal reproductive rhythms, which are a result of interactions between external environmental stimuli and internal endocrine systems in the hypothalamo-pituitary-gonadal axis. While it is long believed that differential expression of neuroendocrine genes contributes to establishing seasonal reproductive rhythms, no systems-level investigation has yet been conducted.

METHODOLOGY/PRINCIPAL FINDINGS

In the present study, by analyzing multiple female goldfish brain microarray datasets, we have characterized global gene expression patterns for a seasonal cycle. A core set of genes (873 genes) in the hypothalamus were identified to be differentially expressed between May, August and December, which correspond to physiologically distinct stages that are sexually mature (prespawning), sexual regression, and early gonadal redevelopment, respectively. Expression changes of these genes are also shared by another brain region, the telencephalon, as revealed by multivariate analysis. More importantly, by examining one dataset obtained from fish in October who were kept under long-daylength photoperiod (16 h) typical of the springtime breeding season (May), we observed that the expression of identified genes appears regulated by photoperiod, a major factor controlling vertebrate reproductive cyclicity. Gene ontology analysis revealed that hormone genes and genes functionally involved in G-protein coupled receptor signaling pathway and transmission of nerve impulses are significantly enriched in an expression pattern, whose transition is located between prespawning and sexually regressed stages. The existence of seasonal expression patterns was verified for several genes including isotocin, ependymin II, GABA(A) gamma2 receptor, calmodulin, and aromatase b by independent samplings of goldfish brains from six seasonal time points and real-time PCR assays.

CONCLUSIONS/SIGNIFICANCE

Using both theoretical and experimental strategies, we report for the first time global gene expression patterns throughout a breeding season which may account for dynamic neuroendocrine regulation of seasonal reproductive development.

Measuring differential gene expression by short read sequencing: quantitative comparison to 2-channel gene expression microarrays

Background

High-throughput cDNA synthesis and sequencing of poly(A)-enriched RNA is rapidly emerging as a technology competing to replace microarrays as a quantitative platform for measuring gene expression.

Results

Consequently, we compared full length cDNA sequencing to 2-channel gene expression microarrays in the context of measuring differential gene expression. Because of its comparable cost to a gene expression microarray, our study focused on the data obtainable from a single lane of an Illumina 1 G sequencer. We compared sequencing data to a highly replicated microarray experiment profiling two divergent strains of S. cerevisiae.

Conclusion

Using a large number of quantitative PCR (qPCR) assays, more than previous studies, we found that neither technology is decisively better at measuring differential gene expression. Further, we report sequencing results from a diploid hybrid of two strains of S. cerevisiae that indicate full length cDNA sequencing can discover heterozygosity and measure quantitative allele-specific expression simultaneously.

QuantPrime--a flexible tool for reliable high-throughput primer design for quantitative PCR

Background

Medium- to large-scale expression profiling using quantitative polymerase chain reaction (qPCR) assays are becoming increasingly important in genomics research. A major bottleneck in experiment preparation is the design of specific primer pairs, where researchers have to make several informed choices, often outside their area of expertise. Using currently available primer design tools, several interactive decisions have to be made, resulting in lengthy design processes with varying qualities of the assays.

Results

Here we present QuantPrime, an intuitive and user-friendly, fully automated tool for primer pair design in small- to large-scale qPCR analyses. QuantPrime can be used online through the internet or on a local computer after download; it offers design and specificity checking with highly customizable parameters and is ready to use with many publicly available transcriptomes of important higher eukaryotic model organisms and plant crops (currently 295 species in total), while benefiting from exon-intron border and alternative splice variant information in available genome annotations. Experimental results with the model plant Arabidopsis thaliana, the crop Hordeum vulgare and the model green alga Chlamydomonas reinhardtii show success rates of designed primer pairs exceeding 96%.

Conclusion

QuantPrime constitutes a flexible, fully automated web application for reliable primer design for use in larger qPCR experiments, as proven by experimental data. The flexible framework is also open for simple use in other quantification applications, such as hydrolyzation probe design for qPCR and oligonucleotide probe design for quantitative in situ hybridization. Future suggestions made by users can be easily implemented, thus allowing QuantPrime to be developed into a broad-range platform for the design of RNA expression assays.