Identification of development and tissue-specific gene expression in the fathead minnow Pimephales promelas, Rafinesque using computational and DNA microarray methods

The fathead minnow Pimephales promelas serves as a model organism for assessing the effects of environmental contaminants on early life stage growth and development. Yet, the utilization of genomic tools has been hindered by the lack of genome sequence and genomic information known from this model species. Utilizing published cDNA library sequences, the authors used sequence similarity to compare 4105 cDNAs isolated from fathead minnow fry (<14 days old) with over 250 000 adult cDNA sequences derived from whole body and various tissue types. The objectives of the computational subtraction were to (1) assess the extent of sequence similarity between developing and adult cDNA libraries and (2) predict which cDNA clones are expressed only in developing organisms. The results of the computational predictions were assessed through the construction of a development‐specific DNA microarray targeting all 4105 sequences in the fry cDNA library as well as 56 known mRNAs in P. promelas. Gene expression was determined by comparing total RNA isolated from fry with total RNA isolated from adult samples (whole animal, kidney, liver, brain, ovary and testes). The results showed that 1381 of the targeted fry cDNA sequences (34%) displayed expression across all sample comparisons, and of these, only 166 genes were found to harbour fry‐specific expression (i.e. no expression in adult samples). Of note, 69% of the genes computationally predicted to be fry specific were found across all experimental results; yet, only 27% of the computationally predicted fry‐specific sequences were experimentally confirmed to be fry specific. An important result was the identification of many novel mRNA sequences specific to the developing minnow, which lack homology with any other known sequence. In addition, the study results included tissue‐specific expression in adult samples. These results demonstrate the capabilities and limitations of inter‐library sequence comparisons as a predictor of gene activity in non‐sequenced organisms and tissues, as well as DNA microarray gene expression studies in non‐sequenced organisms.     

Diversity arrays: a solid state technology for sequence information independent genotyping

Here we present the successful application of the microarray technology platform to the analysis of DNA polymorphisms. Using the rice genome as a model, we demonstrate the potential of a high-throughput genome analysis method called Diversity Array Technology, DArT‘. In the format presented here the technology is assaying for the presence (or amount) of a specific DNA fragment in a representation derived from the total genomic DNA of an organism or a population of organisms. Two different approaches are presented: the first involves contrasting two representations on a single array while the second involves contrasting a representation with a reference DNA fragment common to all elements of the array. The Diversity Panels created using this method allow genetic fingerprinting of any organism or group of organisms belonging to the gene pool from which the panel was developed. Diversity Arrays enable rapid and economical application of a highly parallel, solid-state genotyping technology to any genome or complex genomic mixtures.     

Phylogeographic genomics of mitochondrial DNA: Highly-resolved patterns of intraspecific evolution and a multi-species, microarray-based DNA sequencing strategy for biodiversity studies

Phylogeographic genomics, based on multiple complete mtDNA genome sequences from within individual vertebrate species, provides highly-resolved intraspecific trees for the detailed study of evolutionary biology. We describe new biogeographic and historical insights from our studies of the genomes of codfish, wolffish, and harp seal populations in the Northwest Atlantic, and from the descendants of the founding human population of Newfoundland. Population genomics by conventional sequencing methods remains laborious. A new biotechnology, iterative DNA “re-sequencing”, uses a DNA microarray to recover 30–300 kb of contiguous DNA sequence in a single experiment. Experiments with a single-species mtDNA microarray show that the method is accurate and efficient, and sufficiently species-specific to discriminate mtDNA genomes of moderately-divergent taxa. Experiments with a multi-species DNA microarray (the “ArkChip”) show that simultaneous sequencing of species in different orders and classes detects SNPs within each taxon with equal accuracy as single-species-specific experiments. Iterative DNA sequencing offers a practical method for high-throughput biodiversity genomics that will enable standardized, coordinated investigation of multiple species of interest to Species at Risk and conservation biologists.     

CMG-Biotools,a Free Workbench for Basic Comparative Microbial Genomics

Background

Today, there are more than a hundred times as many sequenced prokaryotic genomes than were present in the year 2000. The economical sequencing of genomic DNA has facilitated a whole new approach to microbial genomics. The real power of genomics is manifested through comparative genomics that can reveal strain specific characteristics, diversity within species and many other aspects. However, comparative genomics is a field not easily entered into by scientists with few computational skills. The CMG-biotools package is designed for microbiologists with limited knowledge of computational analysis and can be used to perform a number of analyses and comparisons of genomic data.

Results

The CMG-biotools system presents a stand-alone interface for comparative microbial genomics. The package is a customized operating system, based on Xubuntu 10.10, available through the open source Ubuntu project. The system can be installed on a virtual computer, allowing the user to run the system alongside any other operating system. Source codes for all programs are provided under GNU license, which makes it possible to transfer the programs to other systems if so desired. We here demonstrate the package by comparing and analyzing the diversity within the class Negativicutes, represented by 31 genomes including 10 genera. The analyses include 16S rRNA phylogeny, basic DNA and codon statistics, proteome comparisons using BLAST and graphical analyses of DNA structures.

Conclusion

This paper shows the strength and diverse use of the CMG-biotools system. The system can be installed on a vide range of host operating systems and utilizes as much of the host computer as desired. It allows the user to compare multiple genomes, from various sources using standardized data formats and intuitive visualizations of results. The examples presented here clearly shows that users with limited computational experience can perform complicated analysis without much training.     

Selective Whole Genome Amplification for Resequencing Target Microbial Species from Complex Natural Samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host–microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10⁵-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2–9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.     

Mercury‐associated DNA hypomethylation in polar bear brains via the LUminometric Methylation Assay: a sensitive method to study epigenetics in wildlife

In this paper we describe a novel approach that may shed light on the genomic DNA methylation of organisms with non‐resolved genomes. The LUminometric Methylation Assay (LUMA) is permissive for genomic DNA methylation studies of any genome as it relies on the use of methyl‐sensitive and ‐insensitive restriction enzymes followed by polymerase extension via Pyrosequencing technology. Here, LUMA was used to characterize genomic DNA methylation in the lower brain stem region from 47 polar bears subsistence hunted in central East Greenland between 1999 and 2001. In these samples, average genomic DNA methylation was 57.9% ± 6.69 (SD; range was 42.0 to 72.4%). When genomic DNA methylation was related to brain mercury (Hg) exposure levels, an inverse association was seen between these two variables for the entire study population (P for trend = 0.17). After dichotomizing animals by gender and controlling for age, a negative trend was seen amongst male animals (P for trend = 0.07) but no associations were found in female bears. Such sexually dimorphic responses have been found in other toxicological studies. Our results show that genomic DNA methylation can be quantitatively studied in a highly reproducible manner in tissue samples from a wild organism with a non‐resolved genome. As such, LUMA holds great promise as a novel method to explore consequential questions across the ecological sciences that may require an epigenetic understanding.     

Applications of DNA and protein microarrays in comparative physiology

DNA microarrays have revolutionized gene expression studies and made large-scale parallel measurement of whole genome expression a feasible technique in model species where genomes are well characterized. Such studies are perfectly suited to unraveling the complex regulation and/or interaction of both genes and proteins likely involved in most physiological processes. Gene expression profiles are currently being used to identify genes underlying a range of physiological responses. Characterization of these genes will help to elucidate the pathways and processes regulating physiological processes. Expanding the use of DNA microarrays to non-model species that have been critical in elucidating certain physiological pathways will be valuable in determining the genes associated with these processes. Approaches that do not require complete genome information have recently been applied to "non-model" organisms. As whole genomes are sequenced for non-model organisms, the application of DNA microarrays to comparative physiology will expand even further. The recent development of protein microarrays will be critical in understanding the regulation of physiological processes not accounted for at the genomic level. Together, DNA and protein microarrays provide the most thorough and efficient method of understanding the molecular basis of physiological processes to date. In turn, classical physiological approaches will be vital in characterizing and verifying the function of the novel genes identified by microarray experiments. Ultimately, DNA and protein microarray expression profiles may be used to predict physiological responses.     

BAC-HAPPY Mapping (BAP Mapping): A New and Efficient Protocol for Physical Mapping

Physical and linkage mapping underpin efforts to sequence and characterize the genomes of eukaryotic organisms by providing a skeleton framework for whole genome assembly. Hitherto, linkage and physical “contig” maps were generated independently prior to merging. Here, we develop a new and easy method, BAC HAPPY MAPPING (BAP mapping), that utilizes BAC library pools as a HAPPY mapping panel together with an Mbp-sized DNA panel to integrate the linkage and physical mapping efforts into one pipeline. Using Arabidopsis thaliana as an exemplar, a set of 40 Sequence Tagged Site (STS) markers spanning ∼10% of chromosome 4 were simultaneously assembled onto a BAP map compiled using both a series of BAC pools each comprising 0.7x genome coverage and dilute (0.7x genome) samples of sheared genomic DNA. The resultant BAP map overcomes the need for polymorphic loci to separate genetic loci by recombination and allows physical mapping in segments of suppressed recombination that are difficult to analyze using traditional mapping techniques. Even virtual “BAC-HAPPY-mapping” to convert BAC landing data into BAC linkage contigs is possible.     

Improved Multiple Displacement Amplification (iMDA) and Ultraclean Reagents

Background

Next-generation sequencing sample preparation requires nanogram to microgram quantities of DNA; however, many relevant samples are comprised of only a few cells. Genomic analysis of these samples requires a whole genome amplification method that is unbiased and free of exogenous DNA contamination. To address these challenges we have developed protocols for the production of DNA-free consumables including reagents and have improved upon multiple displacement amplification (iMDA).

Results

A specialized ethylene oxide treatment was developed that renders free DNA and DNA present within Gram positive bacterial cells undetectable by qPCR. To reduce DNA contamination in amplification reagents, a combination of ion exchange chromatography, filtration, and lot testing protocols were developed. Our multiple displacement amplification protocol employs a second strand-displacing DNA polymerase, improved buffers, improved reaction conditions and DNA free reagents. The iMDA protocol, when used in combination with DNA-free laboratory consumables and reagents, significantly improved efficiency and accuracy of amplification and sequencing of specimens with moderate to low levels of DNA. The sensitivity and specificity of sequencing of amplified DNA prepared using iMDA was compared to that of DNA obtained with two commercial whole genome amplification kits using 10 fg (~1-2 bacterial cells worth) of bacterial genomic DNA as a template. Analysis showed >99% of the iMDA reads mapped to the template organism whereas only 0.02% of the reads from the commercial kits mapped to the template. To assess the ability of iMDA to achieve balanced genomic coverage, a non-stochastic amount of bacterial genomic DNA (1 pg) was amplified and sequenced, and data obtained were compared to sequencing data obtained directly from genomic DNA. The iMDA DNA and genomic DNA sequencing had comparable coverage 99.98% of the reference genome at ≥1X coverage and 99.9% at ≥5X coverage while maintaining both balance and representation of the genome.

Conclusions

The iMDA protocol in combination with DNA-free laboratory consumables, significantly improved the ability to sequence specimens with low levels of DNA. iMDA has broad utility in metagenomics, diagnostics, ancient DNA analysis, pre-implantation embryo screening, single-cell genomics, whole genome sequencing of unculturable organisms, and forensic applications for both human and microbial targets.     

Human-blind probes and primers for dengue virus identification

Reliable detection and identification of pathogens in complex biological samples, in the presence of contaminating DNA from a variety of sources, is an important and challenging diagnostic problem for the development of field tests. The problem is compounded by the difficulty of finding a single, unique genomic sequence that is present simultaneously in all genomes of a species of closely related pathogens and absent in the genomes of the host or the organisms that contribute to the sample background. Here we describe 'host-blind probe design'- a novel strategy of designing probes based on highly frequent genomic signatures found in the pathogen genomes of interest but absent from the host genome. Upon hybridization, an array of such informative probes will produce a unique pattern that is a genetic fingerprint for each pathogen strain. This multiprobe approach was applied to 83 dengue virus genome sequences, available in public databases, to design and perform in silico microarray experiments. The resulting patterns allow one to unequivocally distinguish the four major serotypes, and within each serotype to identify the most similar strain among those that have been completely sequenced. In an environment where dengue is indigenous, this would allow investigators to determine if a particular isolate belongs to an ongoing outbreak or is a previously circulating version. Using our probe set, the probability that misdiagnosis at the serotype level would occur is approximately 1 : 10(150).     

Vom Phänotyp zum Genotyp

Forward Genomics – a comparative genomics approach to link phenotype to genotype Despite availability of several sequenced genomes, we know very little about the specific changes in the DNA that underlie phenotypic differences between species. The main reason is that species differ by both numerous genomic and phenotypic changes. A new comparative genomics method addresses this question by for phenotypes with independent evolutionary losses by searching for genomic regions that exhibit an elevated number of mutations in exactly these phenotype‐loss species. The near future sequencing of thousands of novel genomes will make it possible to use comparative genomics to systematically search for such DNA changes that are associated with phenotypic differences.     

Two-dimensional DNA displays for comparisons of bacterial genomes

We have developed two whole genome-scanning techniques to aid in the discovery of polymorphisms as well as horizontally acquired genes in prokaryotic organisms.First, two-dimensional bacterial genomic display (2DBGD) was developed using restriction enzyme fragmentation to separate genomic DNA based on size, and then employing denaturing gradient gel electrophoresis (DGGE) in the second dimension to exploit differences in sequence composition. This technique was used to generate high-resolution displays that enable the direct comparison of >800 genomic fragments simultaneously and can be adapted for the high-throughput comparison of bacterial genomes. 2DBGDs are capable of detecting acquired and altered DNA, however, only in very closely related strains. If used to compare more distantly related strains (e.g. different species within a genus) numerous small changes (i.e. small deletions and point mutations) unrelated to the interesting phenotype, would encumber the comparison of 2DBGDs. For this reason asecond method, bacterial comparative genomic hybridization (BCGH), was developed to directly compare bacterial genomes to identify gain or loss of genomic DNA. BCGH relies on performing 2DBGD on a pooled sample of genomic DNA from 2 strains to be compared and subsequently hybridizing the resulting 2DBGD blot separately with DNA from each individual strain. Unique spots (hybridization signals) represent foreign DNA. The identification of novel DNA is easily achieved by excising the DNA from a dried gel followed by subsequent cloning and sequencing. 2DBGD and BCGH thus represent novel high resolution genome scanning techniques for directly identifying altered and/or acquired DNA. Published: June 15, 2003     

A PCR-based amplification method retaining the quantitative difference between two complex genomes

With the increasing emergence of genome-wide analysis technologies (including comparative genomic hybridization (CGH), expression profiling on microarrays, differential display (DD), subtractive hybridization, and representational difference analysis (RDA)), there is frequently a need to amplify entire genomes or cDNAs by PCR to obtain enough material for comparisons among target and control samples. A major problem with PCR is that amplification occurs in a nonlinear manner and reproducibility is influenced by stray impurities. As a result, when two complex DNA populations are amplified separately, the quantitative relationship between two genes after amplification is generally not the same as their relation before amplification. Here we describe balanced PCR, a procedure that faithfully retains the difference among corresponding amplified genes by using a simple principle. Two distinct genomic DNA samples are tagged with oligonucleotides containing both a common and a unique DNA sequence. The genomic DNA samples are pooled and amplified in a single PCR tube using the common DNA tag. By mixing the two genomes, PCR loses the ability to discriminate among the different alleles and the influence of impurities is eliminated. The PCR-amplified pooled samples can be separated using the DNA tag unique to each individual genomic DNA sample. The principle of this method has been validated with synthetic DNA, genomic DNA, and cDNA applied on microarrays. By removing the bias of PCR, this method allows a balanced amplification of allelic fragments from two complex DNAs even after three sequential rounds of PCR. This balanced PCR approach should allow genetic analysis in minute laser-microdissected tissues, paraffin-embedded archived material, or single cells.     

A new rice repetitive DNA shows sequence homology to both 5S RNA and tRNA.

Moderately repetitive DNA sequences are found in the genomes of all eucaryotes that have been examined. We now report the discovery of a novel, transcribed, moderately repetitive DNA sequence in a higher plant which is different from any of the known repetitive DNA sequences from any organism. We isolated a rice cDNA clone which hybridizes to multiple bands on genomic blot analysis. The sequence of this 352 bp cDNA contains four regions of homology to the wheat phenylalanine tRNA, including the polymerase III-type promoter. Unexpectedly, two regions of the same 352 bp sequence also show homology to the wheat 5S RNA sequence. Using the cDNA as a probe, we have isolated six genomic clones which contain long tandem repeats of 355 bp sequence, and have sequenced nine repeat units. Our findings suggest that the rice repetitive sequence may be an amplified pseudogene with sequence homology to both 5S RNA and tRNA, but organized as long tandem repeats resembling 5S RNA genes. This is the first example showing homology between the sequences of a moderately repetitive DNA with unknown function and 5S RNA.     

Plant genomics: an overview

Recent technological advancements have substantially expanded our ability to analyze and understand plant genomes and to reduce the gap existing between genotype and phenotype. The fast evolving field of genomics allows scientists to analyze thousand of genes in parallel, to understand the genetic architecture of plant genomes and also to isolate the genes responsible for mutations. Furthermore, whole genomes can now be sequenced. This review addresses these issues and also discusses ways to extract biological meaning from DNA data. Although genomic issuesare addressed from a plant perspective, this review provides insights into the genomic analyses of other organisms.     

RADSeq: next-generation population genetics

Next-generation sequencing technologies are making a substantial impact on many areas of biology, including the analysis of genetic diversity in populations. However, genome-scale population genetic studies have been accessible only to well-funded model systems. Restriction-site associated DNA sequencing, a method that samples at reduced complexity across target genomes, promises to deliver high resolution population genomic data-thousands of sequenced markers across many individuals-for any organism at reasonable costs. It has found application in wild populations and non-traditional study species, and promises to become an important technology for ecological population genomics.     

Genomic Treasure Troves: Complete Genome Sequencing of Herbarium and Insect Museum Specimens

Unlocking the vast genomic diversity stored in natural history collections would create unprecedented opportunities for genome-scale evolutionary, phylogenetic, domestication and population genomic studies. Many researchers have been discouraged from using historical specimens in molecular studies because of both generally limited success of DNA extraction and the challenges associated with PCR-amplifying highly degraded DNA. In today''s next-generation sequencing (NGS) world, opportunities and prospects for historical DNA have changed dramatically, as most NGS methods are actually designed for taking short fragmented DNA molecules as templates. Here we show that using a standard multiplex and paired-end Illumina sequencing approach, genome-scale sequence data can be generated reliably from dry-preserved plant, fungal and insect specimens collected up to 115 years ago, and with minimal destructive sampling. Using a reference-based assembly approach, we were able to produce the entire nuclear genome of a 43-year-old Arabidopsis thaliana (Brassicaceae) herbarium specimen with high and uniform sequence coverage. Nuclear genome sequences of three fungal specimens of 22–82 years of age (Agaricus bisporus, Laccaria bicolor, Pleurotus ostreatus) were generated with 81.4–97.9% exome coverage. Complete organellar genome sequences were assembled for all specimens. Using de novo assembly we retrieved between 16.2–71.0% of coding sequence regions, and hence remain somewhat cautious about prospects for de novo genome assembly from historical specimens. Non-target sequence contaminations were observed in 2 of our insect museum specimens. We anticipate that future museum genomics projects will perhaps not generate entire genome sequences in all cases (our specimens contained relatively small and low-complexity genomes), but at least generating vital comparative genomic data for testing (phylo)genetic, demographic and genetic hypotheses, that become increasingly more horizontal. Furthermore, NGS of historical DNA enables recovering crucial genetic information from old type specimens that to date have remained mostly unutilized and, thus, opens up a new frontier for taxonomic research as well.     

Restriction site tagged (RST) microarrays: a novel technique to study the species composition of complex microbial systems

We have developed a new type of microarray, restriction site tagged (RST), for example NotI, microarrays. In this approach only sequences surrounding specific restriction sites (i.e. NotI linking clones) were used for generating microarrays. DNA was labeled using a new procedure, NotI representation, where only sequences surrounding NotI sites were labeled. Due to these modifications, the sensitivity of RST microarrays increases several hundred-fold compared to that of ordinary genomic microarrays. In a pilot experiment we have produced NotI microarrays from Gram-positive and Gram-negative bacteria and have shown that even closely related Escherichia coli strains can be easily discriminated using this technique. For example, two E.coli strains, K12 and R2, differ by less than 0.1% in their 16S rRNA sequences and thus the 16S rRNA sequence would not easily discriminate between these strains. However, these strains showed distinctly different hybridization patterns with NotI microarrays. The same technique can be adapted to other restriction enzymes as well. This type of microarray opens the possibility not only for studies of the normal flora of the gut but also for any problem where quantitative and qualitative analysis of microbial (or large viral) genomes is needed.