Oligonucleotide microarray for the study of functional gene diversity in the nitrogen cycle in the environment

The analysis of functional diversity and its dynamics in the environment is essential for understanding the microbial ecology and biogeochemistry of aquatic systems. Here we describe the development and optimization of a DNA microarray method for the detection and quantification of functional genes in the environment and report on their preliminary application to the study of the denitrification gene nirS in the Choptank River-Chesapeake Bay system. Intergenic and intragenic resolution constraints were determined by an oligonucleotide (70-mer) microarray approach. Complete signal separation was achieved when comparing unrelated genes within the nitrogen cycle (amoA, nifH, nirK, and nirS) and detecting different variants of the same gene, nirK, corresponding to organisms with two different physiological modes, ammonia oxidizers and denitrifying halobenzoate degraders. The limits of intragenic resolution were investigated with a microarray containing 64 nirS sequences comprising 14 cultured organisms and 50 clones obtained from the Choptank River in Maryland. The nirS oligonucleotides covered a range of sequence identities from approximately 40 to 100%. The threshold values for specificity were determined to be 87% sequence identity and a target-to-probe perfect match-to-mismatch binding free-energy ratio of 0.56. The lower detection limit was 10 pg of DNA (equivalent to approximately 107 copies) per target per microarray. Hybridization patterns on the microarray differed between sediment samples from two stations in the Choptank River, implying important differences in the composition of the denitirifer community along an environmental gradient of salinity, inorganic nitrogen, and dissolved organic carbon. This work establishes a useful set of design constraints (independent of the target gene) for the implementation of functional gene microarrays for environmental applications.