Journal of Heredity 2008:99(2):187–192 Ó The American Genetic Association. 2008. All rights reserved. doi:10.1093/jhered/esm107 For permissions, please email: [email protected]. Advance Access publication January 24, 2008 Sex Identification of (Family Strigidae) Using Oligonucleotide Microarrays

LIH CHIANN WANG,LUCIA LIU SEVERINGHAUS,CHI TSONG CHEN,LU YUAN LIU,CHU HSIANG PAN, DEAN HUANG,HSIAO YUAN LEE,JIHN TSAIR LIR,SHIH CHIEN CHIN,CHANG EN PU, AND CHING HO WANG

From the Graduate Institute of Veterinary Medicine, National University, 1 Sec. 4, Roosevelt Road, Taipei 106, Taiwan (L. C. Wang and C. H. Wang); Taipei Zoo, Taipei, Taiwan (L. C. Wang, Chin, Huang, and Lee); the Research Center for Biodiversity, Academia Sinica, Taipei, Taiwan (Severinghaus); the Scientific and Technical Research Center, Ministry Justice Investigation Bureau, Taipei, Taiwan (Chen and Pu); the Graduate Institute of Plant Science, National Pintung University of Science and Technology, Taipei, Taiwan (Liu); the Department of Hog Cholera, Health Research Institute, Council of Agriculture, Taipei, Taiwan (Pan); and the Chung-Shan Institute of Science and Technology, Taipei, Taiwan (Lir).

Address correspondence to C. H. Wang at the address above, or e-mail: [email protected]. Downloaded from

Molecular sexing of the diversified avian family Strigidae is Z sexual chromosomes has made sexing in many nonratite jhered.oxfordjournals.org difficult. Sex identification using the intron length difference possible (Griffiths et al. 1998; Jensen et al. 2003). between W and Z chromosomal CHD1 genes, as visualized Because female birds are ZW, whereas males are ZZ, by agarose gel electrophoreses, often produces ambiguous polymerase chain reaction (PCR) products from males show results. Here we describe a simple method for sexing a single band, whereas those from females show 2 bands on a variety of Strigidae using oligonucleotide micro- a running gel (Russello and Amato 2001). The CHD1 gene by guest on July 5, 2011 arrays, on which several sex-specific probes operated has several introns (Griffiths and Korn 1997). Two introns complementarily or in concert. The sex of 8 species have been utilized extensively in different species, for was identified clearly on the microarrays through sequence example, the upstream intron flanked by primers 2550F/ recognition. This sequence-directed method can be easily 2718R (Fridolfsson and Ellegren 1999), and the downstream applied to a wider range of Strigidae species. intron flanked by primers 1237L/1272H (Kahn et al. 1998) or P2/P8 (Griffiths et al. 1998). However, there are still a few species that cannot be identified because of genetic diversity. DNA microarrays are miniaturized microsystems based on the Some efforts were made, such as using restriction enzymes ability of DNA to spontaneously find and bind to its (Bermudez-Humaran et al. 2002) and designing new primers complementary sequence by hybridization. Labeled DNA (Ito et al. 2003). The genetic diversity of Strigidae has also molecules in a sample are analyzed using DNA probes tethered made sexing results doubtful. The length polymorphism of at distinct sites on a solid support (Cuzin 2001; Vernet 2002). the 1237L/1272H-amplified products of some species was There are 2 main types of DNA microarrays: cDNA microarrays poor, making identification difficult on agarose gels (Kahn and oligonucleotide microarrays. The probe size of the former is et al. 1998). Although the 2550F/2718R-amplified products usually between 500 and 2000 bp and that of the later is around carried better length polymorphism (Fridolfsson and 25 bp (Kim and Watkinson 2002). DNA microarrays have been Ellegren 1999), the sex of some Strigidae species could not referred to as ‘‘fishing expeditions’’ (Lockhart and Winzeler be determined. None of these methods have been found 2000; Maughan et al. 2001) that allow widespread application: generally applicable to all Strigidae members thus far. gene expression, genotyping, infectious, and genetic disease This study describes sexing a variety of Strigidae species diagnoses, etc (Petrik 2000; Heller 2002). using oligonucleotide microarrays. The sex of 8 Strigidae Strigidae are sexually monomorphic, making sex identi- species was distinctly identified on microarrays, compared with fication difficult via their appearance. Sex identification is the ambiguous results obtained on running gel. The comple- especially important for their breeding and conservation as mentary binding on microarrays between the bird DNA and most of them are endangered species. the sex-specific probes made the sexing results more reliable. The intron length difference of the CHD1 (chromo- The sex of a wider range of Strigidae species was successfully helicase–DNA-binding protein 1) gene between the W and determined using this sequence-directed approach.

187 Journal of Heredity 2008:99(2)

Materials and Methods Owl samples for sex identification included muscle and blood. All the birds used in this study were of known sex, as determined by autopsy and behavior. Muscle tissues were obtained from carcasses during necropsy. Blood was collected during health checkups by veterinarians. Genomic DNA was extracted by incubating 0.1 g of muscle or 100 ll of blood with 20 ll of proteinase K (Amresco, Solon, OH) and 500 ll of digestion buffer (10 mM Tris–HCL, 2 mM ethylenediaminetetraacetic acid, 10 mM NaCl, 1% sodium dodecyl sulfate, 10 mg/ml dithiothreitol; Amresco) at 56 °C overnight, followed by phenol/chloroform isolation 3 times. The aqueous phase was precipitated in 100% ethanol at 4 °C Figure 1. Sex identification based on intronic length for 30 min. After centrifugation, the sample was washed polymorphism on agarose gel. Multiplex PCR was performed using 70% ethanol, resuspended in 100 llof1Â tris- with 2 primer pairs, 1237L/1272H and 2550F/2718R, flanking ethylenediaminetetraacetic acid buffer (Amresco), and passed 2 respective introns of the CHD1 gene. The known-sex through a purifying tube (Centricon-100, Amicon, Bedford, individuals are indicated as M (Male) and F (Female). Length England). Two pairs of 5# end-biotinylated primers, 1237L/ polymorphism appears in the Eurasian- and Elegant 1272H (Kahn et al. 1998) and 2550F/2718R (Fridolfsson and Owl in both 1237L/1272H- and 2550F/2718R-amplified Ellegren 1999), were used in this study. Multiplex PCR product regions, where males have 1 band and females have 2 amplification was carried out in a reaction volume of 25 ll bands. However, the polymorphism is only shown in the 2550F/2718R region of the Brown-wood Owl, Tawny-fish

containing 5 ll of each primer (1 lM), 0.5 llofTaq DNA Downloaded from polymerase (2.5 U/ll), 2.5 llof10Â PCR buffer (500 mM Owl, and Collared-scops Owl. Sex is not identified here for the Eurasian-eagle Owl, Mountain-scops Owl, and Short-eared KCl, 100 mM Tris–HCl pH 8.3, 15 mM MgCl2, 1% Triton X-100, 1600 lg/ml bovine serum albumin and 2 mM each Owl because of the lack of length polymorphism in both deoxynucleoside triphosphate, pH 8.2–8.4), and 2 llof 1237L/1272H and 2550F/2718R regions. jhered.oxfordjournals.org template DNA. The thermal profile for amplification was 94 °C for 11 min, 30Â (94 °Cfor30s,50°C for 30 s, and product was denatured at 95 °C for 10 min and cooled in an 72 °C for 60 s), 72 °C for 5 min. PCR products were ice bath for 2 min. To the microarray chamber was added separated in 3% agarose gels (Gibco, Grand Island, NY), run 200 ll of hybridization buffer (containing the 5# end- in 0.5Â tris-acetate-ethylenediaminetetraacetic acid buffer biotinylated oligonucleotide complementary to the sequence by guest on July 5, 2011 with 0.5 lg/ml of ethidium bromide (Gibco) at 120 V for of positive control probe) and 15 ll of denatured PCR 1.5 h, and visualized under UV light. product, incubated at 50 °C with vibration for 50 min, and The corresponding W and Z bands amplified using washed twice with wash buffer. The blocking reaction was primers 2550F and 2718R in Tawny fish Owls (Ketupa then performed by mixing 0.2 ll of Strep-AP (Streptavidin flavipes) and Elegant Owls (Otus elegans) on gels were cut, conjugate alkaline phosphatase) and 200 ll of blocking eluted, and sequenced using an automated sequencer (ABI reagent at room temperature for 30 min and washing twice 3100, Applied Biosystems, Foster, CA). The PCR products with wash buffer. The colorimetric reaction was then from Eurasian-eagle Owls (Bubo bubo) were separated in 3% implemented by adding 4 ll of nitro blue tetrazolium/ agarose gels at 80 V for 6 h, and the corresponding W and Z 5-bromo-4-chloro-3-indolyl phosphate indolyl phosphate bands amplified using primers 1237L and 1272H were cut, and 196 ll of detection buffer in the chamber, developing in eluted, and sequenced. These sequences were aligned further the dark at room temperature for 5 min, and washing twice for probe designing. with distilled water. The hybridization result was indicated A tail composed of 19 T bases was added to each 5# end as the developed pattern on the microarray, which was read of the oligonucleotide probe, including the positive control directly with the naked eyes. probe (an oligonucleotide from capsid protein VP1 of human enterovirus 71 gene, 5#-ATGAAGCATGTCAGG- # GCTTGGATACCTCG-3 ). Ten lM of each probe was Results then spotted to each specific position on the microarray polymer substrate using an automatic spotting machine The sexing results based on intronic length polymorphism (DR. Easy spotter, Miao-Li, Taiwan) and immobilized using are shown in Figure 1. Fragments amplified using primers a UV Crosslinker (Vilber Lourmat BLX-254, ECC, Marne, 1237L/1272H were located between 200 and 300 bp, and France) with 1.2 J for 5 min. those amplified using primers 2550F/2718R were distrib- The hybridization reaction between each DNA template uted between 600 and 1100 bp. This method was successful and probes was carried out with DR. Chip DIY Kit (DR. in 2 species of Strigidae, the Eurasian-scops Owl (Otus scops) Chip Biotech, Miao-Li, Taiwan). The procedures followed and Elegant Owl, using either of these 2 primer pairs. the manual and are briefly described below. The PCR However, the sex of the Brown-wood Owl (

188 Brief Communications

Figure 2. Alignment of amplified DNA using primers 1237L/1272H. The position of primers and the sequences of the probes are indicated. The consensus sequence (Con) is shown above the sequence alignment with nucleotide identities depicted as dots (.) and sequence deletions as dashes (-). The numbers refer to the nucleotide positions of each amplified fragment. Eu refers to the Eurasian-eagle Owl, whereas W and Z refer to the W and Z chromosome, respectively. leptogrammica), Tawny-fish Owl, and Collared-scops Owl Z chromosome. Sex identification with CHD1 must (Otus bakkamoena) was not identified using primers 1237L/ distinguish between these 2 duplicates. A given CHD1 1272H. Furthermore, the sex of the Eurasian-eagle Owl, intron, whose length differs between the W and Z Mountain-scops Owl (Otus spilocephalus), and Short-eared chromosomes, could be utilized for sex identification. There

Owl (Asio flammeus) was not identified using either primers are several introns on the CHD1 gene. The 1237L/1272H- Downloaded from 1237L/1272H or 2550F/2718R in our sexing procedures. amplified intron, whose length polymorphism between the DNA samples obtained from Eurasian-eagle Owls, W and Z chromosomes is generally inadequate in Strigidae Tawny-fish Owls, and Elegant Owls were sequenced, compared with other avian families, often leads to aligned, and analyzed further (Figures 2 and 3). Five misjudging a female (ZW) as a male (ZZ). Although the jhered.oxfordjournals.org nucleotide probes were designed in order to detect sex- sex of some Strigidae species could be determined through specific fragments of diversified Strigidae species. Probe high-resolution electrophoresis, for example, Tawny-fish W was employed to recognize female individuals because Owl and Collared-scops Owl, this method is time con- the W chromosome exists only in females. In addition, suming and cannot be applied to all species (data not shown). Probe ZW was utilized to detect both females and Some researchers suggested separating the PCR products by guest on July 5, 2011 males, which could be taken as a positive control. When both through polyacrylamide gels (Kahn et al. 1998) or single- Probe W and Probe ZW dots emerged on a microarray after strand conformation polymorphism gels (Cortes et al. 1999). hybridization, it indicated that the tested individual was However, both procedures were harder to perform than a female. If only Probe ZW dots presented, it indicated that agarose gels. Although the 2550F/2718R-amplified intron the specimen was a male. Several sex-specific probes operated has better length polymorphism between the W and Z complementarily or performed as multiple confirmation of chromosomes, this intronic polymorphism does not exist in the sex, making the obtained result more reliable. all Strigidae species. The sex of some species, moreover, The sexing results of 8 Strigidae species using oligonu- cannot be identified on agarose gel, regardless of which cleotide microarrays are shown in Figure 4. All females primer pairs are employed (Figure 1). exhibited at least one of the 3 Probe W dots as well as the Two primer pairs flanking 2 respective CHD1 introns ZW dots. All males displayed only the Probe ZW dots. were employed in this study. Sex identification became more These findings indicate that the sex of Strigidae can be efficient using multiplex PCR. The PCR products of identified easily and reliably using oligonucleotide micro- 3 Strigidae species were further sequenced and analyzed. arrays. It also demonstrates that W probes designed from Five probes were designed to generally detect CHD1W and a few species can be successfully performed in other species, CHD1Z of various Strigidae species. The sex of 8 Strigidae as a result of the sequence homology of CHD1W among species in this study was clearly identified through different Strigidae species. A more accurate method for sex hybridization between DNA templates and probes on identification is therefore achieved using this sequence- oligonucleotide microarrays. The hybridization signals on directed approach, which can be applicable to a wider range microarrays were easily read with the naked eye, requiring of Strigidae species. no laser scanning or imaging systems. The entire microarray manipulation time, including the data readout, was less than 2 h. This was nearly equal to the time needed for gel loading, Discussion electrophoresis, ethidium bromide staining, and visualizing of intron lengths under UV light. The CHD1 gene has been used for sex identification in The characteristic of sequence-directed hybridization on numerous avian species. It exists on the W and also on the microarrays, not the length difference on gels, allowed sex

189 Journal of Heredity 2008:99(2) Downloaded from jhered.oxfordjournals.org by guest on July 5, 2011

Figure 3. Sequence alignment of amplified DNA using primers 2550F/2718R. The position of primers and the sequences of the probes are indicated. Ta and El represent the Tawny-fish Owl and Elegant Owl, whereas W and Z symbolize the W and Z chromosome, respectively. Other illustrations are the same as Figure 2.

190 Brief Communications Downloaded from

Figure 4. Sex identification of 8 Strigidae species using oligonucleotide microarrays. All samples came from known-sex jhered.oxfordjournals.org birds. (A) Microarray map. Each dot indicates the spotted position of each probe. P: Positive control, 1: Probe ZW1, 2: Probe ZW2, 3: Probe W1, 4: Probe W2, 5: Probe W3. (B) The sexing results on the microarrays. M: Male, F: Female. 1: Brown-wood Owl, 2: Tawny-fish Owl, 3: Collared-scops Owl, 4: Eurasian-scops Owl, 5: Elegant owl, 6: Eurasian-eagle Owl, 7: Mountain-scops Owl, 8: Short-eared Owl, 9: Negative control. by guest on July 5, 2011 determination for the species whose CHD1W and CHD1Z Funding introns are equal or very close in length. The detection sensitivity of the oligonucleotide microarray was also higher This work was funded by the National Science Council than that for agarose gel, making fragments invisible on the grant 94N-1205. gel, visible on the array (Figures 1 and 4, 2250F/2718R- amplified product of Short-eared Owl). Sex identification of a wider Strigidae species range was thus achieved using Acknowledgments oligonucleotide microarrays. This approach has also suc- We thank the technicians and keepers at the Taipei Zoo for their kind help cessfully expanded its application to other avian families. with sampling. The sex of various avian families, including Strigidae, could be determined on 1 oligonucleotide microarray. Probes designed from only a few species were capable of References determining the sex of other species on microarrays. This was probably due to the sequence homology among Bermudez-Humaran LG, Chavez-Zamarripa P, Guzman-Velasco A, various species. Some efforts are still needed for those Leal-Garza CH, Montes-de-Oca-Luna R. 2002. Loss of restriction site species whose corresponding sequences are much diverged DdeI, used for avian molecular sexing in Oreophasis derbianus. Reprod from the general. We recommend that sequencing genes Domest Anim. 37:321–323. from those species be used to design additional probes. The Cortes O, Barroso A, Dunner S. 1999. Avian sexing: an optimized protocol microarray can efficiently expand its efficacy by recruiting using polymerase chain reaction single-strand conformation polymorphism. J Vet Diagn Invest. 11:297–299. new probes, which would make this device more practicable. Cuzin M. 2001. DNA chips: a new tool for genetic analysis and diagnostics. Transfus Clin Biol. 8:291–296. Supplementary Material Fridolfsson AK, Ellegren H. 1999. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 30:116–121. Supplementary material can be found at http://www.jhered. Griffiths R, Double MC, Orr K, Dawson RJG. 1998. A DNA test to sex oxfordjournals.org/. most birds. Mol Ecol. 7:1071–1075.

191 Journal of Heredity 2008:99(2)

Griffiths R, Korn RM. 1997. A CHD1 gene is Z chromosome linked in the Lockhart DJ, Winzeler EA. 2000. Genomics, gene expression and DNA chicken Gallus domesticus. Gene. 197:225–229. arrays. Nature. 405:827–836. Heller MJ. 2002. DNA microarray technology: devices, system, and Maughan NJ, Lewis FA, Smith V. 2001. An introduction to arrays. J Pathol. applications. Annu Rev Biomed Eng. 4:129–153. 195:3–6. Ito H, Sudo-Yamaji A, Abe M, Murase T, Tsubota T. 2003. Sex Petrik J. 2000. Micro array technology: the future of blood testing. Vox identification by alternative polymerase chain reaction methods in Sang. 80:1–11. Falconiformes. Zool Sci. 20:339–344. Russello MA, Amato G. 2001. Application of a noninvasive, PCR-based Jensen T, Pernasetti FM, Durrant B. 2003. Conditions for rapid sex test for sex identification in an endangered parrot, Amazona guildingii. Zoo determination in 47 avian species by PCR of genomic DNA from Biol. 20:41–45. blood, shell-membrane blood vessels, and feathers. Zoo Biol. 22: Vernet G. 2002. DNA-chip technology and infectious diseases. Virus Res. 561–571. 82:65–71. Kahn NW, Judith SJ, Quinn TW. 1998. Chromosome-specific intron size differences in the avian CHD gene provide a simple and efficient method Received December 8, 2006 for sex identification in birds. AUK. 115:1074–1078. Accepted October 12, 2007 Kim DS, Watkinson JC. 2002. Gene chip expression analysis in head and neck cancer. Clin Otolaryngol. 27:296–303. Corresponding Editor: Oliver Ryder Downloaded from jhered.oxfordjournals.org by guest on July 5, 2011

192