Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

This article is part of the supplement: Proceedings of the Avian Genomics Conference and Gene Ontology Annotation Workshop

Open Access Highly Accessed Proceedings

The value of avian genomics to the conservation of wildlife

Michael N Romanov18*, Elaina M Tuttle2, Marlys L Houck1, William S Modi1, Leona G Chemnick1, Marisa L Korody2, Emily M Stremel Mork1, Christie A Otten1, Tanya Renner1, Kenneth C Jones3, Sugandha Dandekar4, Jeanette C Papp4, Yang Da5, NISC Comparative Sequencing Program6, Eric D Green6, Vincent Magrini7, Matthew T Hickenbotham7, Jarret Glasscock7, Sean McGrath7, Elaine R Mardis7 and Oliver A Ryder1

Author Affiliations

1 Genetics Division, San Diego Zoo's Institute for Conservation Research, Zoological Society of San Diego, Arnold and Mabel Beckman Center for Conservation Research, 15600 San Pasqual Valley Rd., Escondido, CA 92027, USA

2 Dept. of Biology, Indiana State University, Science 285E, 403-25 North 6th St., Terre Haute, IN 47809, USA

3 Genetic Identification Services, 9552 Topanga Canyon Blvd., Chatsworth, CA 91311, USA

4 UCLA Sequencing and Genotyping Core, Dept. of Human Genetics, 5309 Gonda, 695 Charles E. Young Dr. South, David Geffen School of Medicine at UCLA, Los Angeles, CA 90095, USA

5 Dept. of Animal Science, 265D Haecker Hall, University of Minnesota, 1364 Eckles Ave., St. Paul, MN 55108, USA

6 Genome Technology Branch and NIH Intramural Sequencing Center (NISC), National Human Genome Research Institute, National Institutes of Health, 50 South Dr., Bldg. 50, Rm. 5222, Bethesda, MD 20892, USA

7 Washington University Genome Sequencing Center, Washington University School of Medicine, Campus Box 8501, 4444 Forest Park Ave., St. Louis, MO 63108, USA

8 BACPAC Resources, Children's Hospital Oakland Research Institute, 747 52nd St., Oakland, CA 94609, USA

For all author emails, please log on.

BMC Genomics 2009, 10(Suppl 2):S10  doi:10.1186/1471-2164-10-S2-S10

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/10/S2/S10


Published:14 July 2009

© 2009 Romanov et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Genomic studies in non-domestic avian models, such as the California condor and white-throated sparrow, can lead to more comprehensive conservation plans and provide clues for understanding mechanisms affecting genetic variation, adaptation and evolution.

Developing genomic tools and resources including genomic libraries and a genetic map of the California condor is a prerequisite for identification of candidate loci for a heritable embryonic lethal condition. The white-throated sparrow exhibits a stable genetic polymorphism (i.e. chromosomal rearrangements) associated with variation in morphology, physiology, and behavior (e.g., aggression, social behavior, sexual behavior, parental care).

In this paper we outline the utility of these species as well as report on recent advances in the study of their genomes.

Results

Genotyping of the condor resource population at 17 microsatellite loci provided a better assessment of the current population's genetic variation. Specific New World vulture repeats were found in the condor genome. Using condor BAC library and clones, chicken-condor comparative maps were generated. A condor fibroblast cell line transcriptome was characterized using the 454 sequencing technology.

Our karyotypic analyses of the sparrow in combination with other studies indicate that the rearrangements in both chromosomes 2m and 3a are complex and likely involve multiple inversions, interchromosomal linkage, and pleiotropy. At least a portion of the rearrangement in chromosome 2m existed in the common ancestor of the four North American species of Zonotrichia, but not in the one South American species, and that the 2m form, originally thought to be the derived condition, might actually be the ancestral one.

Conclusion

Mining and characterization of candidate loci in the California condor using molecular genetic and genomic techniques as well as linkage and comparative genomic mapping will eventually enable the identification of carriers of the chondrodystrophy allele, resulting in improved genetic management of this disease.

In the white-throated sparrow, genomic studies, combined with ecological data, will help elucidate the basis of genic selection in a natural population. Morphs of the sparrow provide us with a unique opportunity to study intraspecific genomic differences, which have resulted from two separate yet linked evolutionary trajectories. Such results can transform our understanding of evolutionary and conservation biology.

Background

Introduction

Genomic studies in a variety of mammalian taxa have contributed to the development of more comprehensive plans for their conservation, as well as to our understanding of the generation and maintenance of genetic diversity in general (e.g., [1-3]). With the advances of genomic resources for other species, it is now feasible to expand investigation to additional taxonomic realms.

Non-domesticated species of birds are emerging as new animal models with great potential to advance comparative avian genomics and contribute to conservation efforts for threatened species. Birds show incredible diversity in morphology, physiology, and behavior – all of which are analogous to phenotypic variation encountered in other species. With the completion of the chicken genome, as well as significant advances in turkey and zebra finch genomic resources, it is now possible to examine the genetic bases of complex traits using "tools" borrowed from these species. There is a special need for genomic studies in endangered birds to secure their recovery in natural habitats and to increase their resistance to potential threats (e.g., disease, anthropogenic effects). In addition, because many avian species are relatively easy to observe in the wild, they have been intensely studied for decades and, as a result, several long-term studies of their behavior, ecology, and evolution now exist (e.g., the pied flycatcher, Ficedula hypoleuca, the great tit, Parus major, and the blue tit, Parus caeruleus; see for review [4,5]), making it possible to examine adaptive genetic variation. Finally, many avian species are the proverbial "canary in the coal mine", providing us with indicators of the health of an ecosystem and sentinels for the impacts of environmental changes.

In this paper, we review the recent genomics advances in two model avian species, the California condor (Gymnogyps californianus) and the white-throated sparrow (Zonotrichia albicollis). The critically endangered condor has been affected by habitat destruction and an extreme bottleneck. The white-throated sparrow is an indicator species for the Northern Boreal forest, an essential habitat for avian biodiversity and a key carbon sink [6]. Through studies of their genomes, we will develop a deeper understanding of the genetic mechanisms affecting variation, adaptation, evolution, and behavior in these species, thereby obtaining new information relevant to conservation management of wild populations.

California condor

The management of wild avian species, including those that are critically threatened and/or endangered, has greatly benefited from the recent advances in genetics and genomics. Among the most striking examples of genetic application to conservation of an endangered species is the California condor (Gymnogyps californianus). The condor is an avian icon of California that belongs to the New World vulture family (Cathartidae). Over the last two centuries, there has been a rapid decline in the condor population as their natural habitat dramatically shrank. In the 1950s, wild condors were thought to be restricted to Los Padres National Forest in central California. By the early 1980's, a California condor rescue and restoration program was set in motion. Captive propagation was quite successful, and condors have been reintroduced to central and northern California (including Los Padres National Forest and Monterey County), the Grand Canyon area in Arizona, and Baja California del Norte, Mexico [7].

Captive breeding necessitated the initiation of genetic studies to assess the diversity of the remaining population. Early genetic studies of the condor were done at the San Diego Zoo's Institute for Conservation Research, and included the genotyping of the original population with RFLPs [8] and mtDNA [9] analyses (Figure 1). Within the founder population, three unrelated groups (clans) of closely related individuals were identified based on the analysis of multi-locus RFLP fingerprinting [8]. Using mtDNA markers, four distinct maternal lineages were then identified [9]. Finally, karyotype and molecular techniques were developed for sexing condors, involving first an RFLP approach and later utilizing conserved avian sex chromosome sequences for PCR primers [9]. Together, these findings were used to recreate the hypothetical founding generation of California condors, an essential step in developing an effective population management strategy [10].

thumbnailFigure 1. California condor pedigree of 69 individuals used for genetic studies. A pedigree includes founders and early generation offspring used for the DNA fingerprint [8], mtDNA [9] and microsatellite (msat; current study) analyses. All individuals are designated using their studbook numbers.

Genetically small populations are at increased risk of impacts of deleterious mutations as a result of founder effect and genetic drift. A lethal mutation was identified in the expanding pedigree of California condors that negatively affected limb development [11]. The mutant phenotype, called chondrodystrophy, is similar to genetic disorders observed in the chicken (nanomelia) [12], turkey [13], Japanese quail [14], mouse [15,16], and human [17]. Ralls et al. [11] found that the embryonic lethal condition of chondrodystrophy segregates in the affected pedigree as a Mendelian character in a manner consistent with an autosomal recessive mode of inheritance.

With the advent of the genomic era, a detailed comparative molecular cytogenetic analysis of condor chromosomes (GCA) was carried out using chicken macrochromosome (GGA1-GGA9, GGAZ, and GGAW) paints, revealing that the condor has 80 chromosomes [18], comparable to the chicken and turkey karyotypes. This study showed a great homology of chicken and condor macrochromosomes, except for GGA4 that corresponds to GCA4 and GCA9, and suggested an incomplete differentiation of the condor sex chromosomes in the evolution [18].

These studies laid the groundwork for application of genomic technologies in support of California condor conservation efforts. Construction of a genetic map of the California condor and identification of linked markers for the chondrodystrophy mutation were primary goals. The availability of the chicken genome sequence [19] and a large-insert BAC library that was generated at the BACPAC Resources Center, Children's Hospital Oakland Research Institute (CHORI) [20] were crucial components of these efforts. The BAC library consists of 89,665 clones and provides approximately 14-fold coverage of the condor genome. Based upon 172 probes, mostly available from the chicken genome project, we screened the library and produced a first-generation BAC-based chicken-condor comparative map that contained 93 loci.

Among the genes identified in the condor BAC library, several are candidates for chondrodystrophy because they are involved in bone and tissue formation. One such candidate is aggrecan (ACAN), which is an integral part of the extracellular matrix [21]. The aggrecan molecule consists of three globular domains, two chondroitin sulfate domains, and one keratin region. Aggrecan molecules attach to a long hyaluronan molecule and form aggregates. Moreover, this functional candidate gene was found to cause skeletal dysplasia in model organisms [12,13,15-17]. We detected two BAC clones positive for the condor ACAN gene, one of which was sequenced. As expected, the condor gene sequence showed a great similarity to the chicken homolog [20]. Building on these data, we have continued to investigate the genetic basis of disease in this species. The development and application of practical approaches for management of deleterious genes identified in small populations, including diagnosis of risk factors, identification of unaffected carriers, and breeding strategies to retain genetic variation in such populations with segregating lethal mutations will find precedent in these initial studies.

White-throated sparrow

The white-throated sparrow (Zonotrichia albicollis; Family Emberizidae) is a socially monogamous passerine that breeds almost exclusively (>85% of the breeding population) in the North American boreal forest. With large declines in species abundances, increasing evidence for global warming and acid rain, and with the realization that the amount of boreal forest currently protected will not sustain migratory bird populations [6], conservation efforts have begun to focus on such key areas before it is too late. An essential component of a valid conservation plan will include the study of the genomes of boreal species, such as the white-throated sparrow, for the amount and kind of adaptive genetic variation they actually possess, and its distribution across varying habitats.

The white-throated sparrow has several important characteristics that make it an ideal model for the study of conservation genomics. The species is polymorphic, and both sexes occur as either tan (T) or white (W) morphs based on the color of their crown stripes (Figure 2). In addition to morphology, behavioral and life-history characteristics differ between the morphs as well [22-25] (Table 1). In general, white birds are more aggressive than tan birds, yet white are not dominant to tan. Males of the two morphs employ alternative reproductive strategies based on the tradeoffs between current and future reproduction. White males invest heavily in securing matings through song, territorial intrusion, bigamy, and the pursuit of promiscuity (i.e. extra-pair copulations) at the expense of mate-guarding and paternal care [22,23,26-28]. White males have higher levels of circulating testosterone than tan males early in the breeding season [29], most likely because they tend to establish territories in areas of high density where encounters with conspecifics are frequent [24]. By contrast, tan males invest in monogamy and high levels of parental care, and tend to establish territories in more isolated areas. White and tan females seem to exhibit similar trade-offs between investment in parental effort and investment in reproductive effort.

Table 1. Characteristics of morph/sex classes in the white-throated sparrow. Characteristics of morph/sex classes; all data from [22-24] except where noted: (1), [76-80]; (2), [26-28,81]; (3), [22,82]. "?" indicates where current knowledge is lacking.

thumbnailFigure 2. Morphs of the white-throated sparrow. Photographs of the morphs of the white-throated sparrow; A) white morphs are heterozygous for a chromosomal rearrangement on autosome 2, while B) tan morphs are homozygous for the non-inverted form of autosome 2.

Interestingly, morph was absolutely correlated with a chromosomal polymorphism resulting from what was believed to be a pericentric inversion of the 2nd largest autosome [30,31]. White birds are heterozygous for the rearrangement (i.e. 2m/2, where 2m represents the inverted chromosome 2), and tan birds are homozygous non-carriers (i.e. 2/2). The autosomal inversion segregates in a Mendelian fashion [31]. White-throated sparrows mate disassortatively with respect to the polymorphism [22,23,31], which maintains both morphs in the population [31] and results in parental "phenotypes" that differ in the amount and type of care provided. These parental types respond differently to environmental conditions (e.g., predation, disease, and long-term climatic factors) by altering the condition, sex, and morph of their own offspring through mechanisms that persist across generations.

Results gathered from detailed cytogenetic and molecular analyses in the white-throated sparrow can be combined with long-term data on behavior, physiology, and ecology to reveal the bases of adaptive variation in a natural population. Morphs of the sparrow are somewhat analogous to natural "inbred" lines, where selective factors affecting fitness can be linked to specific genes. In this species, we have the unique opportunity to transcend from genotype to ecosystem. Such results can transform our understanding of evolutionary and conservation biology. In addition, results can be generalized to endangered and at-risk species, providing conservation biologists with a strong foundation with which to build their captive and reintroduction plans.

Methods

DNA samples

A condor female named Molloko, studbook (SB) #45, the first ever condor chick hatched in captivity, was selected as the DNA source for constructing the condor microsatellite and BAC libraries. Sequences of microsatellite and BAC clones from Molloko have been deposited in the GenBank and Trace Archive.

For genetic diversity and linkage studies, a set of 121 condor DNA samples including that of Molloko was chosen. These samples were previously isolated for fingerprinting studies and/or gender determination and represent a pedigree of related individuals from four generations and seven families (Additional file 1, created using a computer program Pedigraph, Version 2.4 [32]). This pedigree serves as a condor resource population for further microsatellite and linkage analyses, and also involves three chondrodystrophic individuals, SB #60 and eggs #1405 and #2537.

Additional file 1. Pedigree of a condor resource population chosen for microsatellite and linkage analyses, with 121 individuals.

Format: PNG Size: 62KB Download fileOpen Data

White-throated sparrow blood samples were collected from individuals that were either part of the long-term study population at the Cranberry Lake Biological Station, New York, or from a captive population housed in research aviaries at Indiana State University. Additional blood samples were also collected from junco (Junco hyemalis hyemalis), white-crowned sparrow (Zonotrichia leucophrys), song sparrow (Melospiza melodia), swamp sparrow (M. georgiana), and purple finch (Carpodacus purpureus). All blood was separated via centrifugation, and the red blood cells were stored in a lysis buffer at 4°C until needed for further analyses.

Tissue samples from the golden-crowned sparrow (Z. atricapilla), Harris's sparrow (Z. querula), as well as the majority of white-crowned sparrows (Z. leucophrys) and juncos (J. h. oreganus, J. h. caniceps), were provided by the Museum of Vertebrate Zoology (University of California, Berkeley) and the University of Alaska Museum of the North. Finally, blood samples from the rufous-collared sparrow (Z. capensis) were kindly provided by Ignacio T. Moore, Department of Biological Sciences, Virginia Polytechnic Institute and State University.

Microsatellite-based analyses

A California condor microsatellite-enriched library was created at Genetic Identification Services, following the company's protocol, and includes eight sublibraries, A to H, with enrichments for repeat motifs CA, ATG, TACA, TAGA, AAG, AAT, AAAT, and CATC, respectively. For marker design, approximately 180 random clones were sequenced. Dye-labeled primers for a set of 17 informative microsatellite markers identified by Genetic Identification Services (Table 2) were used to amplify condor DNA. These included 13 tetra-, two tri- and two dinucleotide repeat loci. Three more avian microsatellite markers, FhU2, HrU2, and HrU6 (Accession Numbers X84361, X84087, and X84091) [33,34], were tested using Schuelke's [35] procedure for the fluorescent labeling of PCR fragments with the M13(-21) universal primer.

Table 2. Summary of genetic population statistics for 17 condor microsatellite loci.

The PCR products were run on an ABI 3100, ABI 3130 or ABI 3730 Genetic Analyzer instruments (Applied Biosystems), the fragments or alleles generated were verified (Figure 3), and the data analyzed for population genetic parameters and linkage relationships. The following computer programs were used for estimating population genetic statistics: Microsatellite Tools for Excel, Version 3.1 [36], and POPGENE, Version 1.31 [37]. A preliminary linkage analysis was done using Locusmap, Version 1.1 [38].

thumbnailFigure 3. An example of microsatellite genotyping profile at the loci C5 and D6. The DNA from a condor male (SB #162) was amplified with the C5 and D5 primers labeled with two different fluorescent dyes. The bird is heterozygous for C5 (green alleles, 179 and 191 bp) and homozygous for D5 (blue alleles, 256 bp both). Red peaks are size standards. Values along the y-axis represent the intensity of fluorescence (a.u.), values along the x-axis mark the size of the DNA fragments (a function of time).

Additionally, over 1,900 random clones from the same microsatellite-enriched library were bi-directionally sequenced at High-Throughput Sequencing Solutions (University of Washington, Seattle). These were then used to develop additional polymorphic microsatellite markers using the Primer3 online program [39] for primer design. All sequences have been deposited in the GenBank (Accession Numbers DQ471953, DQ483109DQ484036, EF108178EF108181, EF116886EF116903) and analyzed for the presence of repetitive elements using a package of NCBI BLAST programs [40] and RepeatMasker http://www.repeatmasker.org/ webcite[41].

BAC library screen

A second library screening was carried out using the CHORI protocol http://bacpac.chori.org/overgohyb.htm webcite and a four-dimensional filter hybridization procedure [42], with modifications [20]. The latter, in brief, included the simultaneous 32P-labeling of 216 overgo probes (Additional file 2) and their hybridization to BAC filters, with variable 36 overgo pools per filter per dimension. This strategy allowed the completion of a total of 24 hybridizations in a relatively short amount of time. The overgo probes were available from the previous chicken and zebra finch genome projects [43,44] or newly designed on the basis of available condor and chicken sequences (Additional file 2). Other screening procedure details and overgo sequences can be provided by the authors upon request.

Additional file 2. List of overgo probes used for the second screening of the condor BAC library.

Format: XLS Size: 102KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Thirteen condor BAC clones were subsequently sequenced at NIH Intramural Sequencing Center (NISC) following protocols described in [45], and their sequences are publicly available from GenBank (Accession Numbers AC171379, AC171743, AC172166, AC183448AC183451, AC183845, AC183846, AC188502AC188504, and AC191192).

Cytogenetic techniques

Avian fibroblast cell lines from feather pulp were established at the San Diego Zoo's Institute for Conservation Research as described in [46]; these cell lines were used for karyotyping and fluorescence in situ hybridization (FISH).

Metaphase chromosomes were obtained by arresting cell division with colcemid (Gibco BRL), hypotonic treatment with 0.075 M KCl, and fixation in 3:1 methanol:acetic acid. To assess chromosome number, we stained chromosome preparations with 5% Giemsa (Sigma) in 0.07 M phosphate buffer (pH 6.8). To examine banding patterns, some slide preparations were trypsin G-banded following [47].

BAC DNA was used as FISH probes and extracted following standard procedures. After extraction from a BAC clone, DNA was purified and labeled with a fluorochrome according to the manufacturer's guidelines (Vysis, Abbott Laboratories).

Fragment and sequence analysis

Fragments from coding regions of the ACAN candidate gene were PCR amplified and sequenced using an ABI 3100 or ABI 3130 Genetic Analyzer. Sequencher (Gene Codes) was employed for post-sequencing analysis.

A condor fibroblast cell line was utilized for mRNA extraction and consequent generation of cDNA sequences at Washington University Genome Sequencing Center using the 454 sequencing technology [48]. The total RNA was used to prepare a polyT-primed cDNA library by an in-house version of the Clontech SMART protocol without normalization [49]. The 454 sequence data set has been deposited in Trace Archive.

For analyses of the white-throated sparrow morph, we modified the Michopoulos et al. [50] protocol to amplify a portion of the vasoactive intestinal peptide (VIP) gene located on chromosome 2 that either has a DraI restriction site (the "2m" state) or lacks it (the "2" state). Our modifications included fluorescently labeling both the forward and reverse primers so that all products could be visualized on an ABI PRISM 310 Genetic Analyzer, and reducing PCR reactions to a final volume of 10 μl. Figure 4 shows the genetic analyzer output for white and tan morphs. To examine evolutionary history, we used this same "morph" assay on DNA samples collected from several related passerine species.

thumbnailFigure 4. Molecular assay for morph determination. Molecular assay for morph determination in the white-throated sparrow, modified in the Tuttle laboratory for use on automated sequencers from primers described in [50]. Panel A shows a tan morph with a single peak at 284 bp amplified from chromosome 2. Panel B shows a white morph with a band for chromosome 2 (at 284 bp) and two additional bands at 88 bp and 189 bp, which occur when PCR product amplified from chromosome 2m is digested with restriction endonuclease DraI. If they occur in nature, we assume that birds homozygous for 2m ("super whites") would produce a pattern showing only the 88 bp and 189 bp bands. This technique has been verified in adults of known plumage. Values along the y-axis represent the intensity of fluorescence (a.u.), values along the x-axis mark the size of the DNA fragments (a function of time).

Results

California condor studies

Genetic diversity and linkage analysis

For genetic diversity assessment, linkage analysis and construction of the condor genetic map, we chose 121 related individuals that formed a condor resource population. Their genotyping at 17 microsatellite loci resulted in the estimates of within-population genetic variation shown in the Table 2. The average number of alleles was 2.41 ± 0.80 per locus, the effective number of alleles 1.85 ± 0.50, and the Shannon's information index 0.66 ± 0.26. The average heterozygosity for all loci was 0.45 ± 0.01, and the genetic diversity 0.42 ± 0.04. For eight of the 17 loci, there were deviations from Hardy-Weinberg equilibrium.

A preliminary linkage analysis suggested that among the 17 microsatellite loci we may expect a linkage between loci D10 and D6 (LOD score = 21.67), as well as between A20 and D9 (LOD score = 5.12) and between B7 and H238 (LOD score = 5.12).

To expand the microsatellite mapping subproject, we tested three more avian microsatellite markers, FhU2, HrU2, and HrU6 [34,35] that previously showed satisfactory amplification results in a wide range of other birds, including some New World vultures. We were able to amplify the condor FhU2 and HrU2 PCR fragments, while HrU6 did not work on condor DNA in our hands.

Additionally, we obtained over 900 clone sequences from the microsatellite-enriched library and derived from them around 300 new markers for genotyping the resource population at the next project stage.

In silico mapping

In total, we obtained approximately 1000 short genomic sequences from the microsatellite-enriched library. In the course of in silico mapping, we found that 30% of these are homologous to chicken sequences across almost all chicken chromosomes (GGA1-GGA15, GGA17, GGA19-GGA21, GGA23, GGA24, GGA26-GGA28, GGAZ, and UN) and some other avian sequences. Many of these sequences contain microsatellites and other repetitive elements (LINEs, retroviral LTRs, and satellites) known in chicken. We also found satellite sequences previously detected by Keyser et al. [51] that only occur in New World vultures.

BAC library screening

The second screen of a 2.8× subset of the condor BAC library was done using 216 overgo probes. The latter were mostly available from the zebra finch and chicken genome projects [20] or derived from the condor microsatellite library clone sequences (Additional file 2). As a result, we added 100 more loci onto the BAC-based chicken-condor comparative map, and it currently contains 199 loci (Additional file 3). At present, the total number of condor BAC-gene or other sequence assignments is almost 670 (Additional file 4).

Additional file 3. Chicken chromosome map (in Mb) showing location of the 199 loci aligned with the California condor BAC clones (numbers in parentheses). Created with MapChart [83] and GenomePixelizer [84] software. With the colored symbols, chicken loci are given that were anchored to condor BACs using California condor (red), New World vulture (gold) and zebra finch (green) sequence derived OVERGOs or Uprobes (blue) followed up by direct sequencing. GGAUn_random, an unknown virtual chicken chromosome sequence generated by the International Chicken Genome Sequencing Consortium [19] that is yet to be assigned to a named chromosome. LOC770630 (multiple GGAUn_random?), the probe for this locus may also align with the chicken sequence at additional unresolved locations. *, the CRES0014 and GAPR1 sequences align with several chicken chromosomes, GGA1 and GGA2 being the best matches, respectively.

Format: PPT Size: 296KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Additional file 4. Database of condor BACs positive for known avian genes, markers and other sequences

Format: XLS Size: 539KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

FISH mapping

To verify the efficiency of cross-species hybridization of the condor BAC library and to build the condor cytogenetic map, we performed FISH using more than 50 condor BAC clones that were positive for targeted condor genes (Figure 5). Most BACs were found to be orthologous to the appropriate chicken genes and chromosomes (Additional file 5). We identified one intrachromosomal rearrangement on chromosome 4 (the FGF2 locus), with more intrachromosomal rearrangements on the Z chromosome.

Additional file 5. Chicken-condor comparative cytogenetic map showing a considerable degree of conserved synteny.

GGA, Gallus gallus (chicken) chromosomes; GCA, Gymnogyps californianus (California condor) chromosomes.

Format: PPT Size: 147KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

thumbnailFigure 5. California condor FISH mapping. An example of FISH using metaphase chromosomes of the California condor and a biotinylated BAC clone containing the ACVR2B gene located on GCA2.

Candidate gene hunting

Among the positive clones for almost 200 genes identified in the condor BAC library, we selected and sequenced one that harbors the ACAN gene (CH262-21P20, Acc. No. AC171743). Based on the sequence information, we designed PCR primers to re-sequence the entire gene in affected chicks, their parents, as well as unaffected non-relatives. After sequencing coding regions of the gene, we have yet to detect SNPs that could be associated with the chondrodystrophy.

Comparative human-condor mapping

To date, NISC has generated sequence of 13 condor BACs, each 100–200 kb in size. These sequences are orthologous to human chromosome 7 and six chicken chromosomes. Based on the sequences, a condor-human comparative physical map for a region on HSA7 was designed and is available online through the NISC web site (Additional file 6). After comparing against the chicken whole genome sequence using BLAST, the 13 sequenced BACs were also added onto the BAC-based chicken-condor comparative map (Additional file 3).

Additional file 6. A condor-human comparative physical map for HSA7q31

Format: PPT Size: 53KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

454 sequencing

At the Washington University Genome Sequencing Center, almost 440,000 cDNA sequences were generated from a fibroblast cell line using the 454 technology. The total length of these short transcripts (about 100 bp each) is 43,379,069 bp. We found that 78% transcripts were homologous to chicken genes, while being distributed across almost all chicken chromosomes (Figure 6). After assembly and clustering, we were left with around 15,000 contigs that contained about 190,000 reads.

thumbnailFigure 6. Distribution of the condor 454 transcripts across the chicken chromosomes. The distribution of the condor transcripts based on their orthology to the chicken chromosomes.

Surprisingly, for a number of genes expressed in this cell line, we had unusually high numbers of reads, with the most extreme transcript having over 55,000 reads (Additional file 7). Amongst the most abundant transcripts (Additional file 8), there were several proteins that are expected to be expressed in fibroblast cell lines, like proteoglycans, a keratin associated protein, fibroblast growth factor receptor 1 and alpha-actin, and other housekeeping proteins. However, one protein, the mitotic arrest deficient-like 1 (MAD1L1), had especially higher levels of expression in this cell line.

Additional file 7. Distribution of the condor 454 transcripts homologous to chicken genes. The histogram illustrates how the condor reads fall into the chicken gene bins where the '0' bin represents all the genes with 0 to 9 members, the '10' bin genes with 10 to 19 reads and so on. There were a lot of rare transcripts (~45,000), with 0 to 20 members, vs. small numbers of other abundant and extremely abundant transcripts, with up to 55,000+ members, in this fibroblast cell line.

Format: PPT Size: 62KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Additional file 8. Most abundant transcripts in a transformed condor fibroblast cell line. Proteins that are expected to be expressed in a fibroblast cell line are highlighted with gray, while one protein (MAD1L1) that could be involved in abnormal cell features and functions is given in bold.

Format: DOC Size: 49KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

White-throated sparrow studies

Chromosomal studies

To confirm previous studies and as a necessary first step in genomic studies, we established eight white-throated sparrow cell lines (six white, and two tan; KB15792 – KB15799) accessioned in the Frozen Zoo® at the Zoological Society of San Diego. Using standard Giemsa staining, we confirmed that there is a total diploid chromosome number of 82, including the Z and W sex chromosomes (Figure 7). G-banding comparisons indicate that the rearrangements distinguishing 2m from 2 and 3a from 3 are complex. The banding differences in each pair involve the centromere and cannot be resolved by inverting only one segment of the chromosome and thus are likely the result of more than one inversion on each chromosome (Figure 8).

thumbnailFigure 7. Complete karyotype of white female. Complete karyotype of white (2m/2) female (Z/W) KB15793, showing 82 macro- and microchromosomes, including 40 autosome pairs and two sex chromosomes (Z and W); data generated at the San Diego Zoo's Institute for Conservation Research, Zoological Society of San Diego. In addition to the rearrangement of 2, chromosome 3 also sometimes exhibits an alternate form (3a; as seen here).

thumbnailFigure 8. G-banded chromosomes 1–9 and Z from the sparrow. G-banded chromosomes 1–9 and Z from a 2m/2 3a/3 white-throated sparrow (KB15794; San Diego Zoo's Institute for Conservation Research, Zoological Society of San Diego).

Karyotypic analyses of 46 randomly selected white-throated sparrows (Table 3) indicate that with regards to chromosome 2 and 3, the most common genotypes are 2/2 3a/3a (tan; type C) and 2m/2 3a/3a (white; type F). Homozygous genotypes 2m/2m and 3/3 are rare [31] or nonexistent (this study). Most of the chicks from tan male – white female pairs (T × W) were the type F genotype (2m/2 3a/3a); few of their chicks (approximately 15%) were heterozygous for chromosome 3 (3a/3). By contrast, offspring from white male – tan female pairs (W × T) were fairly evenly split between types B (2/2 3a/3), C (2/2 3a/3a), and F (2m/2 3a/3a); and, over 33% of their chicks were 3a/3 heterozygotes.

Table 3. The percentage of white-throated sparrows sampled having various arrangements of chromosomes 2 and 3. Table adapted from [31]. Types A-C are all tan, whereas types D-G are all white.

Comparative analyses

The VIP assay for morph is based on the fact that there is a SNP on chromosome 2m that forms a restriction site for DraI, whereas the VIP fragment on chromosome 2 does not contain that site [50]. In an assay of 546 white-throated sparrows (Table 4), 50.9% were 2m/2 (white) and 49.1% were tan (2/2); no individuals were 2m/2m ("superwhites"). When we examined this site in the four other species of Zonotrichia (Z. atricapilla, Z. querula, Z. leucophrys, Z. capensis), only the North American species of Zonotrichia show polymorphism for the VIP fragment suggesting that the common ancestor of all four had the polymorphism (Table 4; Figure 9). All 10 rufous-collared sparrows (Z. capensis) tested were monomorphic for the DraI restriction site. In addition, other passerine birds tested (three subspecies of junco, including Junco hyemalis hyemalis, J. h. oreganos, and J. h. caniceps, swamp sparrow (Melospiza georgiana), song sparrow (M. melodia), and purple finch (Carpodacus purpureus)) were also monomorphic for the DraI restriction site (Table 4). Preliminary sequence analysis of the VIP intron indicates that in the rufous-collared sparrow, the junco, the song sparrow, and the swamp sparrow, the SNP is the same as that found in chromosome 2m of the white-throated sparrow (Figure 10).

Table 4. Results of the white-throated sparrow "morph" assay using the DraI restriction site located in the VIP gene [50]. Species heterozygous for the DraI restriction site show 3 fragments (F1-F3); species homozygous for the DraI restriction site show 2 bands, one smaller band at approximately 72–89 bp and one larger band at 185–190 bp (F1 and F2); finally, homozygotes that do not carry the DraI restriction site show only a single band at 284 bp (F3).

thumbnailFigure 9. VIP restriction site mapped on Zonotrichia tree. Phylogenetic tree for Zonotrichia from [75] showing branch lengths (% nucleotide differentiation). Mapped on the tree are presence (+), absence (-), or both presence and absence (+/-) of the DraI restriction site in the VIP fragment. Data for the restriction site were derived from Table 1. Only the four North American Zonotrichia species show polymorphism for the VIP fragment, suggesting that the polymorphism arose in the common ancestor of those four species.

thumbnailFigure 10. VIP sequence for several passerine species. Sequence alignment for an intron of the vasoactive intestinal peptide (VIP) gene (primers described in [50]). Bases highlighted in gray show differences between the white-throated sparrow (Zonotrichia albicollis; WTSP of both morphs), Harris's sparrow (Z. querula; HASP), the white-crowned sparrow (Z. leucophrys, WCSP), the rufous-collared sparrow (Z. capensis, RUFS), the dark-eyed junco (Junco hyemalis; DEJU), the song sparrow (Melospiza melodia; SOSP), and the swamp sparrow (M. georgiana; SWSP). Yellow highlighting indicates the DraI restriction site, which is found in dark-eyed juncos and song sparrows, as well as white morphs (W) white-throated sparrows; green highlighting shows the lack of a DraI restriction site found in tan (T) morph white-throated sparrows.

Discussion

California condor genetics and genomics

Genetic management of California condors aims to maximize the retention of genetic diversity by minimizing mean kinship within the population. We used a series of mostly tetranucleotide repeat markers to evaluate genetic diversity among individuals obtained from the wild and ascertain linkage relationships between microsatellite loci. The observed level of heterozygosity (0.45) found in the resource population could be considered as a good sign of a successful captive propagation program and population recovery. On the other hand, observation of Hardy-Weinberg equilibrium deviations at eight out of 17 loci might suggest that the current condor population is impacted by inbreeding. These results can be used for validating genetic management recommendations currently used for this endangered species.

Taking advantage of recent progress in genome studies of the domestic chicken, we are developing a genetic map for the California condor. Among the 17 microsatellite loci already used, there are three potential linkages established. Approximately 300 more markers have been designed and, after validation, will be used for genotyping the resource population. This would be the first attempt to make a genetic linkage map for condors, which will be later supplemented with SNP markers to locate and characterize candidate loci for the chondrodystrophy mutation, to identify carriers of the chondrodystrophy allele and to improve genetic management of this and other heritable disease risk-factor loci.

The conventional re-sequencing of ACAN, a likely chondrodystrophy candidate gene, was complicated by the compound structure of this gene, including the large number of exons (18 in the chicken) and a complex repetitive structure of its longest exon that encodes for chondroitin sulfate attachment domain. The latter contains in the chicken a variable repeat region that includes 19 repeats, each of 60 nucleotides, with the reported nanomelic stop codon [12]. Although we have not yet been able to identify a causative SNP in coding regions of the condor ACAN gene, we continue this investigation in both exonic and intronic sequences as well as other regulatory regions that could also be associated with changes in the ACAN regulation. We plan to use a high throughput 454 sequencing technology that may identify or eliminate this locus as the locus of the chondrodystrophy mutation. Construction of the condor linkage map will bring us closer to the mapping of the chondrodystrophy locus and it will enable confirmation of a role for the ACAN gene or other (marker) loci linked to the disease. In either scenario, we will be able to develop a molecular assay for heterozygous carriers of the chondrodystrophy mutation.

The 454 cDNA sequence data reported here provides a first insight into the condor transcriptome. Interestingly, ACAN was found to be expressed in this cell line. We expect that establishment of fibroblast cell lines for normal and affected chicks will enable direct comparison of full-length cDNA sequences to significantly assist in the identification of the potential chondrodystrophy mutation in ACAN or other candidate genes. Such data will provide one more functional genomics approach to pinpoint the disorder mutation.

MAD1L1 was shown to be highly overexpressed in the condor fibroblast cell line analyzed. MAD1L1 is a component of the mitotic spindle-assembly checkpoint, and its mutations are suggested to play a role in the pathogenesis of various types of human cancer (e.g., [52]). Abnormal overexpression of MAL1L1 gene may be linked to a particular phenotype observed in this condor cell line. The cell line has a continuous long-term proliferation and heteroploid features, an indication that the line is transformed possibly due to a mutation affecting MAD1L1 expression. As MAD1L1 is known to be involved in neoplasia, further investigation of this cell line can provide a source of comparative studies and shed light on de-regulation of cell proliferation in avian species.

We were able to obtain dense coverage of the Z chromosome (Additional file 3) with the second BAC library screen. We chose to cytogenetically assign the Z-linked BACs to condor chromosomes using FISH because preliminary indications suggested that there could be intra- and interchromosomal rearrangements in the condor relative to the chicken Z chromosomes. Furthermore, we consider it useful for the future studies to compare Z-linked gene expression in condor males and females to see if there is any difference or if we have a dosage compensation effect similar to what was shown for the zebra finch and chicken [53,54].

Overall, we detected a few changes in gene location and order between chicken and condor chromosomes, indicating a high degree of conserved synteny. Our data support a higher similarity of the cathartid genome to the avian ancestral karyotype and a more basal position of Cathartidae as compared to Accipitridae and Falconidae [55], as well as the general stability of avian genomes over the course of evolution [56]. Another piece of evidence in support of those ideas comes from the fact that the chicken and condor genomes demonstrate a high degree of sequence homology and share many common avian repetitive elements including LINEs, retroviral LTRs, and satellites. However, the presence of satellite sequences specific to New World vultures [51] suggests that some distinctive genomic features have emerged in the cathartids over almost the 100-million years of evolution in this avian lineage. We hypothesize that the specific satellite sequences inherent to the Cathartidae may represent unique genomic signatures that could help resolve the long disputed taxonomic position of this avian family (e.g., [55,57,58]).

The 454 cDNA sequence data, as well as condor microsatellite library clone sequences, have been added to a new avian genomic database Gallus GBrowse that includes the whole genome sequence of the chicken as a reference http://birdbase.net/cgi-bin/gbrowse/gallus08/ webcite[59]. Along with sequence information for turkey and zebra finch, this database serves as a new powerful tool for avian comparative genomics. All the sequence information available for the condors is also being deposited in GenBank and Trace Archive, adding the condor to the ever-expanding list of avian genomes (Table 5). Finally, there is a web site we have developed for the condor genomics project: https://msu.edu/~romanoff/index2.htm webcite.

Table 5. Representation of avian species in the NCBI databases (as of June 8, 2009).

White-throated sparrow as a behavioral model

Advances in white-throated sparrow genomics have begun at an obvious point – the characterization of the chromosomal rearrangement affecting morphic differences. Thorneycroft [30,31] first karyotyped the species, reporting a diploid chromosome number of 82 or 84, and based on centromere position, presumed pericentric inversions in both chromosome 2 and 3. Here, our karyotypic analyses confirm a diploid number of 82, with the typical avian pattern of approximately 30 pairs of microchromosomes (Figure 7; see [60,61] for reviews of avian microchromosomes). G-banding analyses suggest that the rearrangements in chromosomes 2 and 3 are not simple, but instead consist of multiple inversions involving the centromere. These results are consistent with those reported in Thomas et al. [62], who used comparative cytogenetic mapping to show that 2m and 2 differ by a pair of pericentric inversions spanning at least 98 Mb or greater than 86% of the chromosome. Chromosomal inversions are often adaptive because they can result in co-adapted gene complexes [63], therefore further investigation of gene structure in this species is warranted. However, to be valuable to evolutionary and conservation biology, such studies must be done within the context of the ecological and social environments that selected for the maintenance of polymorphism.

In the white-throated sparrow, recombination is suppressed in heterozygotes for chromosome 2m. Thorneycroft [31] found that chiasma did form between the p-arm of chromosome 2 and one arm of chromosome 2m, but they did not form between the q-arm of chromosome 2 and chromosome 2 m. In addition, chiasma formed between the p-arm of chromosome 3 and the p-arm of chromosome 3a, but not between the q-arms of these chromosomes. Using nine loci within the 2m rearrangement, Thomas et al. [62] confirmed that recombination only occurred in the telomere area between the p-arm of chromosome 2 and chromosome 2m, thus restricting gene flow between the two. In addition, they found that chromosome 2 had five times the nucleotide diversity of chromosome 2m. Based on these results, Thorneycroft's [31] meiotic data, and our data on population genotype frequencies (Table 3), we predict that chromosome 3a has also evolved under limited recombination and reduced gene flow.

Chromosomal inversions may be disadvantageous if they foster the accumulation of suites of deleterious alleles within the area of low recombination. Since the inversion in the white-throated sparrow is relatively large and covers at least 86% of chromosome 2m, it is likely that some negative fitness affects are associated with homozygosity; the same is likely true for chromosome 3a. We did not find 2m/2m or 3a/3a individuals in our karyotypic analyses, nor did we find any 2m/2m birds in over 546 individuals sampled in the VIP assay. Thorneycroft [31] found that parents with a single 2m chromosome passed this on to half of their offspring however, parents with a single chromosome 3 passed this on to only a quarter of their offspring. Given his results, we would expect that white birds (2m/2) pairing with other white birds (2m/2) would result in up to 25% lethal or semi-lethal homozygotes. It is surprising then that in our study, the frequencies of particular genotypes differed in the sexes and between the two disassortative pair types (Table 3). Certain combinations of chromosomes 2 and 3 also appeared more viable than others. Together, these factors suggest non-Mendelian transmission, interchromosomal linkage and pleiotropy, and a strong interaction between autosomal genotype and sex. Our laboratory is continuing to investigate the relationship between chromosomes 2 and 3 using a combination of genomics techniques and population data.

Comparative analyses between the white-throated sparrow and its relatives can be quite revealing about genomic evolution. In the rufous-collared sparrow (Zonotrichia capensis), Rocha et al. [64] reported polymorphisms in chromosomes 3 and 5 – also presumably due to pericentric inversions. The closely related junco exhibits polymorphisms for pericentric inversions in chromosomes 2 and 5 [65] and the white-crowned sparrow (Z. leucophrys) has centric rearrangements in chromosomes 3, 5, and 12 [66]. Thus, chromosomal inversions in Zonotrichia and its congeners seem to be relatively common [65,66]. Our comparative analysis of the VIP intron [50] revealed that the DraI polymorphism likely existed in the common ancestor of all four North American species of Zonotrichia, but not in the South American species, Z. capensis. In accordance with these results, sequence data from VIP (Figures 9 and 10) and other areas of chromosome 2m suggest that the first inversion predated the divergence of white-throated sparrow from both white-crowned sparrow and Harris's sparrow, originating approximately 2.2 ± 0.3 MYA [62]. Surprisingly, our comparative data on the VIP intron suggests that, contrary to Thorneycroft's [31] expectations, chromosome 2m (or at least a portion of it) might actually be "ancestral". We are undertaking further comparative mapping to resolve these alternatives.

In the past few years, genomic resources for avian species have advanced by leaps and bounds, making it now possible to use comparative genomics to determine the genes responsible for differences in white and tan morphology, behavior, and physiology. Chicken chromosome paints revealed that white-throated sparrow chromosome 2 was analogous to chicken chromosome 3 [62], now allowing us to pinpoint candidate genes associated with phenotypic differences. In addition, major advances in the zebra finch genome [67] as well as a white-throated sparrow BAC library (CHORI-264; http://bacpac.chori.org/library.php?id=469 webcite), mean that detailed mapping is now possible. Over 1000 microsatellite markers, developed for parentage and population genetics analyses in songbirds [68,69], have been used to construct a predicted passerine genome map based on sequence similarity to the chicken genome. Using microsatellite markers that cross-amplify in the Zonotrichia with the detailed pedigrees of over 350 families (data not shown) we will be able to generate linkage maps for this species. Finally, in order to advance conservation and evolutionary genomics in this species, it will be important to identify adaptive genetic variation and the role that the environment has had in selecting for the evolution of alternative genotypes. We have also established a website for the sparrow project: http://www.whitethroatedsparrow.org webcite.

Conclusion

The retention and capacity for generation of novel adaptive genetic variation in a population is key to conservation efforts. Without maintaining this diversity, given the increased anthropogenic impacts of populations and their habitats, many species will face increased vulnerability to extinction in the near future [70,71]. Conservation genetics has traditionally relied on assays of neutral genetic variation, assuming that it accurately reflects adaptive variation [72]. However, the correlation between variation in adaptive and non-adaptive genetic elements can sometimes be low [72-74], and so a true understanding of the interplay between genes and biodiversity will benefit from the additional data that a broader genomics approach can provide. Here we show the utility of applying such a "genomics approach" to study of two key avian species.

Conservation efforts in the endangered California condor will benefit from the characterization and deeper understanding of the segregating lethal chondrodystrophy mutation. Management of the species' genetic diversity under circumstances of maintaining a viable, self-sustaining, managed population that provides birds for reintroduction, and genetic and demographic augmentation of newly established wild California condor populations in the United States and Mexico is a conservation paradigm that merits the fullest support of scientific prowess in genomics, the health sciences, and behavioral ecology, as well as through public policy, public education and law enforcement. Results gathered from the condor studies will help identify recovery options for other at-risk species as their populations continue to decline and experience genetic bottlenecks.

Genomic studies in the white-throated sparrow are still in their infancy although the potential benefits of this avenue of study are great. Morphs of the white-throated sparrow provide a novel intraspecific comparison of selective pressures, where factors affecting fitness can be linked to specific genes – in other words, we have the opportunity to directly assess adaptive variation. In addition, genomic studies in the sparrow will contribute to our understanding of chromosomal inversions, where suites of genes are often inherited as co-adapted gene complexes and the "success" of a phenotype is highly dependent on orchestrated gene cascades and environmental effects.

Birds are highly influential model organisms in conservation and evolutionary research, and the foundation of work on many avian species is already strong. With the advances of genomic resources for chicken, turkey, zebra finch, and now, the California condor and white-throated sparrow, it is now feasible to expand genomics studies to other avian groups.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MNR, as a principal researcher for the California condor genome project, performed BAC library screen, PCR amplification and sequencing of the ACAN gene and other avian microsatellites, analyzed the results and co-authored the manuscript writing. EMT is principal investigator of the white-throated sparrow project, coordinated the laboratory and field research, and co-authored the manuscript writing. MLH provided karyotyping and G-banding for the white-throated sparrow project, avian cell culture preparation and maintenance. WSM carried out FISH mapping and cell culture preparation for the condor project. LGC was involved in DNA extraction and curation, molecular sexing, and genotype work. MLK assisted in PCR amplification and sequencing for the white-throated sparrow project, and sample collection in the field. EMSM and CAO performed DNA extraction for the condor project and microsatellite genotyping. TR contributed to the re-sequencing of the ACAN gene (exon 12), microsatellite genotyping, and molecular sexing. KCJ was responsible for microsatellite library construction supervision, and primer design and validation. SD and JCP directed microsatellite genotyping for the condor project. YD provided resource population construction for the condor project and linkage analysis. EDG and NISC Comparative Sequencing Program supported the sequencing of condor BACs. VM, MTH, JG, SM, ERM carried out 454 sequencing and sequence analyses for the condor project. OAR managed overall condor genome project development and supervision, and contributed to writing the manuscript.

Acknowledgements

We would like to thank the many institutions and people because without their support these projects would be still at the conceptual level. California condor project funding mostly comes from the private Seaver Institute as well as the Zoological Society of San Diego. Funding for the white-throated sparrow work was provided by the Department of Life Sciences at Indiana State University (ISU), the Center for Public Service and Community Engagement (ISU), the Indiana Academy of Sciences, the Indiana Space Grant Consortium, the San Diego Zoo, the Lilly Foundation, and the National Science Foundation (DEB-0217463 to E.M. Tuttle). This work was supported in part by the Intramural Program of the National Human Genome Research Institute, National Institutes of Health; we are grateful to the members of the NISC Comparative Sequencing Program for their efforts in condor BAC sequencing.

We are also grateful to the following collaborators: Maxim Koriabine, Mikhail Nefedov and Pieter J. de Jong from CHORI; Jerry Dodgson and Bill Payne from Michigan State University; Suellen Charter, Julie Fronczek, Andrea Johnson, and Michael Mace, Zoological Soicety of San Diego; Kent Reed from University of Minnesota; Michael Dorschner from University of Washington and High-Throughput Sequencing Solutions; Tim Harkins and Jane Hutchinson from Roche; Eduard Akhunov from Kansas State University; Carl Schmidt from University of Delaware; Lincoln Stein from Cold Spring Harbor Laboratory; Denise Gonzalez, currently a UCLA undergrad student; Jonathan D. Ballou and Katherine Ralls from Smithsonian's National Zoological Park; Elizabeth A. Thompson and Adele Mitchell from University of Washington; Rusty A. Gonser and Gary Stuart from Indiana State University; Ignacio T. Moore from Virginia Polytechnic Institute and State University; Teri Lear from the University of Kentucky; and the Cranberry Lake Biological Station. We thank two anonymous reviewers for careful reading and useful comments on the manuscript.

This article has been published as part of BMC Genomics Volume 10 Supplement 2, 2009: Proceedings of the Avian Genomics Conference and Gene Ontology Annotation Workshop. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/10?issue=S2

References

  1. O'Brien SJ, Weinberg J, Lyons LA: Comparative genomics: lessons from cats.

    Trends Genet 1997, 13:393-399. PubMed Abstract | Publisher Full Text OpenURL

  2. Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ 3rd, Zody MC, et al.: Genomic sequence, comparative analysis and haplotype structure of the domestic dog.

    Nature 2005, 438:803-819. PubMed Abstract | Publisher Full Text OpenURL

  3. Ryder OA: Conservation genomics: applying whole genome studies to species conservation efforts.

    Cytogenet Genome Res 2005, 108:6-15. PubMed Abstract | Publisher Full Text OpenURL

  4. Lundberg A, Alatalo RV: The Pied Flycatcher. London, T. & A. D. Poyser Ltd; 1992. OpenURL

  5. Otter KA, (ed): The Ecology and Behavior of Chickadees and Titmice: An Integrated Approach. New York, Oxford University Press; 2007. OpenURL

  6. Blancher P, Wells J: The Boreal Forest Region: North America's Bird Nursery. Boreal Songbird Initiative and the Canadian Boreal Initiative; 2005. OpenURL

  7. Snyder N, Snyder H: The California Condor: a Saga of Natural History and Conservation. San Diego, Academic Press; 2000. OpenURL

  8. Geyer CJ, Ryder OA, Chemnick LG, Thompson EA: Analysis of relatedness in the California condors: from DNA fingerprints.

    Mol Biol Evol 1993, 10:571-589. OpenURL

  9. Chemnick LG, Kumamoto AT, Ryder OA: Genetic analyses in support of conservation efforts for the Californian condor.

    Int Zoo Yb 2000, 37:330-339. Publisher Full Text OpenURL

  10. Ralls K, Ballou JD: Genetic status and management of California condors.

    Condor 2004, 106:215-228. Publisher Full Text OpenURL

  11. Ralls K, Ballou JD, Rideout BA, Frankham R: Genetic management of chondrodystrophy in California condors.

    Anim Conserv 2000, 3:145-153. Publisher Full Text OpenURL

  12. Li H, Schwartz NB, Vertel BM: cDNA cloning of chick cartilage chondroitin sulfate (aggrecan) core protein and identification of a stop codon in the aggrecan gene associated with the chondrodystrophy, nanomelia.

    J Biol Chem 1993, 268:23504-23511. PubMed Abstract | Publisher Full Text OpenURL

  13. Dannenberg A, Buss EG, Goetinck PF: Proteoglycans in the cartilage of the turkey chondrodystrophy mutant.

    Prog Clin Biol Res 1982, 110(Pt B):85-89. PubMed Abstract OpenURL

  14. Hermes JC, Abbott UK, Johnston E, Owens M: A new chondrodystrophy mutation in Japanese quail (Coturnix japonica).

    J Hered 1990, 81:222-224. PubMed Abstract | Publisher Full Text OpenURL

  15. Watanabe H, Kimata K, Line S, Strong D, Gao LY, Kozak CA, Yamada Y: Mouse cartilage matrix deficiency (cmd) caused by a 7 bp deletion in the aggrecan gene.

    Nat Genet 1994, 7:154-157. PubMed Abstract | Publisher Full Text OpenURL

  16. Krueger RC Jr, Kurima K, Schwartz NB: Completion of the mouse aggrecan gene structure and identification of the defect in the cmd-Bc mouse as a near complete deletion of the murine aggrecan gene.

    Mamm Genome 1999, 10:1119-1125. PubMed Abstract | Publisher Full Text OpenURL

  17. Gleghorn L, Ramesar R, Beighton P, Wallis G: A mutation in the variable repeat region of the aggrecan gene (AGC1) causes a form of spondyloepiphyseal dysplasia associated with severe, premature osteoarthritis.

    Am J Hum Genet 2005, 77:484-490. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Raudsepp T, Houck ML, O'Brien PC, Ferguson-Smith MA, Ryder OA, Chowdhary BP: Cytogenetic analysis of California condor (Gymnogyps californianus) chromosomes: comparison with chicken (Gallus gallus) macrochromosomes.

    Cytogenet Genome Res 2002, 98:54-60. PubMed Abstract | Publisher Full Text OpenURL

  19. International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution.

    Nature 2004, 432:695-716. PubMed Abstract | Publisher Full Text OpenURL

  20. Romanov MN, Koriabine M, Nefedov M, de Jong PJ, Ryder OA: Construction of a California condor BAC library and first-generation chicken-condor comparative physical map as an endangered species conservation genomics resource.

    Genomics 2006, 88:711-718. PubMed Abstract | Publisher Full Text OpenURL

  21. Kiani C, Chen L, Wu YJ, Yee AJ, Yang BB: Structure and function of aggrecan.

    Cell Res 2002, 12:19-32. PubMed Abstract | Publisher Full Text OpenURL

  22. Tuttle EM: Mate choice and the maintenance of stable polymorphism in the white-throated sparrow. PhD Dissertation, State University of New York at Albany; 1993. OpenURL

  23. Tuttle EM: Alternative reproductive strategies in the polymorphic white-throated sparrow: behavioral and genetic evidence.

    Behav Ecol 2003, 14:425-432. Publisher Full Text OpenURL

  24. Formica VA, Gonser RA, Ramsay SM, Tuttle EM: Spatial dynamics of alternative reproductive strategies: the role of neighbors.

    Ecology 2004, 85:1125-1136. Publisher Full Text OpenURL

  25. Tuttle EM, Jensen RR, Formica VA, Gonser RA: Using remote sensing image texture to study habitat use patterns: A case study using the polymorphic white-throated sparrow (Zonotrichia albicollis).

    Glob Ecol Biogeogr 2006, 15:349-357. Publisher Full Text OpenURL

  26. Knapton RW, Falls JB: Differences in parental contribution among pair types in the polymorphic white-throated sparrow.

    Can J Zool 1983, 61:1288-1292. OpenURL

  27. Kopachena JG, Falls JB: Aggressive performance as a behavioral correlate of plumage polymorphism in the white-throated sparrow (Zonotrichia albicollis).

    Behaviour 1993, 124:249-266. Publisher Full Text OpenURL

  28. Kopachena JG, Falls JB: Postfledgling parental care in the white-throated sparrow (Zonotrichia albicollis).

    Can J Zool 1993, 71:227-232. Publisher Full Text OpenURL

  29. Spinney LH, Bentley GE, Hau M: Endocrine correlates of alternative phenotypes in the white-throated sparrow (Zonotrichia albicollis).

    Horm Behav 2006, 50:762-771. PubMed Abstract | Publisher Full Text OpenURL

  30. Thorneycroft HB: Chromosomal polymorphism in the white-throated sparrow, Zonotrichia albicollis.

    Science 1966, 154:1571-1572. PubMed Abstract | Publisher Full Text OpenURL

  31. Thorneycroft HB: A cytogenetic study of the white-throated sparrow, Zonotrichia albicollis (Gmelin).

    Evolution 1975, 29:611-621. Publisher Full Text OpenURL

  32. Ellegren H: Polymerase-chain reaction (PCR) analysis of microsatellites – a new approach to studies of genetic relationships in birds.

    Auk 1992, 109:886-895. OpenURL

  33. Garbe JR, Da Y: Locusmap User Manual. Version 1.1. St. Paul, Department of Animal Science, University of Minnesota; 2003. OpenURL

  34. Primmer CR, Møller AP, Ellegren H: Resolving genetic relationships with microsatellite markers: a parentage testing system for the swallow Hirundo rustica.

    Mol Ecol 1995, 4:493-498. PubMed Abstract | Publisher Full Text OpenURL

  35. Schuelke M: An economic method for the fluorescent labelling of PCR fragments.

    Nat Biotechnol 2000, 18:233-234. PubMed Abstract | Publisher Full Text OpenURL

  36. Park SDE: Trypanotolerance in West African cattle and the population genetic effects of selection. PhD Thesis, University of Dublin; 2001. OpenURL

  37. Yeh FC, Yang R-C, Boyle TBJ, Ye Z-H, Mao JX: POPGENE, the User-Friendly Shareware for Population Genetic Analysis. Version 1.31. Edmonton, Molecular Biology and Biotechnology Centre, University of Alberta; 1999. OpenURL

  38. Garbe JR, Da Y: Pedigraph User Manual. Version 2.4. St. Paul, Department of Animal Science, University of Minnesota; 2008. OpenURL

  39. Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by Krawetz S, Misener S. Totowa, Humana Press; 2000:365-386. OpenURL

  40. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract | Publisher Full Text OpenURL

  41. Smit AFA, Hubley R, Green P: [http://www.repeatmasker.org/] webcite

    RepeatMasker Open-3.0. Seattle, Institute for Systems Biology; 1996. OpenURL

  42. Romanov MN, Price JA, Dodgson JB: Integration of animal linkage and BAC contig maps using overgo hybridization.

    Cytogenet Genome Res 2003, 102:277-281. PubMed Abstract | Publisher Full Text OpenURL

  43. Wallis JW, Aerts J, Groenen M, Crooijmans R, Layman D, Graves T, Scheer D, Kremitzki C, Higgenbotham J, Gaige T, et al.: A physical map of the chicken genome.

    Nature 2004, 432:761-764. PubMed Abstract | Publisher Full Text OpenURL

  44. Romanov MN, Dodgson JB: Cross-species overgo hybridization and comparative physical mapping within avian genomes.

    Anim Genet 2006, 37:397-399. PubMed Abstract | Publisher Full Text OpenURL

  45. Blakesley RW, Hansen NF, Mullikin JC, Thomas PJ, McDowell JC, Maskeri B, Young AC, Benjamin B, Brooks SY, Coleman BI, et al.: An intermediate grade of finished genomic sequence suitable for comparative analyses.

    Genome Res 2004, 14:2235-2244. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Sasaki M, Ikechi T, Makino S: A feather pulp culture technique for avian chromosomes, with notes on the chromosomes of the peafowl and the ostrich.

    Experientia 1968, 24:1292-1293. PubMed Abstract | Publisher Full Text OpenURL

  47. Seabright M: The use of proteolytic enzymes for the mapping of structural rearrangements in the chromosomes of man.

    Chromosoma 1972, 36:204-210. PubMed Abstract | Publisher Full Text OpenURL

  48. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, et al.: Genome sequencing in microfabricated high-density picolitre reactors.

    Nature 2005, 437:376-380. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Warren WC, Hillier LW, Marshall Graves JA, Birney E, Ponting CP, Grützner F, Belov K, Miller W, Clarke L, Chinwalla AT, et al.: Genome analysis of the platypus reveals unique signatures of evolution.

    Nature 2008, 453:175-183. PubMed Abstract | Publisher Full Text OpenURL

  50. Michopoulos V, Maney DL, Morehouse CB, Thomas JW: A genotyping assay to determine plumage morph in the white-throated sparrow (Zonotrichia albicollis).

    Auk 2007, 124:1330-1335. Publisher Full Text OpenURL

  51. Keyser C, Montagnon D, Schlee M, Ludes B, Pfitzinger H, Mangin P: First isolation of tandemly repeated DNA sequences in New World vultures and phylogenetic implications.

    Genome 1996, 39:31-39. PubMed Abstract | Publisher Full Text OpenURL

  52. Tsukasaki K, Miller CW, Greenspun E, Eshaghian S, Kawabata H, Fujimoto T, Tomonaga M, Sawyers C, Said JW, Koeffler HP: Mutations in the mitotic check point gene, MAD1L1, in human cancers.

    Oncogene 2001, 20:3301-3305. PubMed Abstract | Publisher Full Text OpenURL

  53. Itoh Y, Melamed E, Yang X, Kampf K, Wang S, Yehya N, Van Nas A, Replogle K, Band MR, Clayton DF, et al.: Dosage compensation is less effective in birds than in mammals.

    J Biol 2007, 6:2. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  54. Arnold AP, Itoh Y, Melamed E: A bird's-eye view of sex chromosome dosage compensation.

    Annu Rev Genomics Hum Genet 2008, 9:109-127. PubMed Abstract | Publisher Full Text OpenURL

  55. Tagliarini MM, Pieczarka JC, Nagamachi CY, Rissino J, de Oliveira EH: Chromosomal analysis in Cathartidae: distribution of heterochromatic blocks and rDNA, and phylogenetic considerations.

    Genetica 2008, in press. PubMed Abstract | Publisher Full Text OpenURL

  56. Griffin DK, Robertson LB, Tempest HG, Vignal A, Fillon V, Crooijmans RP, Groenen MA, Deryusheva S, Gaginskaya E, Carré W, et al.: Whole genome comparative studies between chicken and turkey and their implications for avian genome evolution.

    BMC Genomics 2008, 9:168. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  57. Ericson PG, Anderson CL, Britton T, Elzanowski A, Johansson US, Kallersjo M, Ohlson JI, Parsons TJ, Zuccon D, Mayr G: Diversification of Neoaves: integration of molecular sequence data and fossils.

    Biol Lett 2006, 2:543-547. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Brown JW, Rest JS, García-Moreno J, Sorenson MD, Mindell DP: Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages.

    BMC Biol 2008, 6:6. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Schmidt C, Romanov M, Ryder O, Magrini V, Hickenotham M, Glasscock J, Mardis E, McGrath S, Stein L: Gallus Genome GBrowse: a unified genomic database for the chicken.

    Nucleic Acids Res 2008, 36(Suppl 1):D719-D723. OpenURL

  60. Burt DW: Origin and evolution of avian microchromosomes.

    Cytogenet Genome Res 2002, 96:97-112. PubMed Abstract | Publisher Full Text OpenURL

  61. Axelsson E, Webster MT, Smith NGC, Burt DW, Ellegren H: Comparison of the chicken and turkey genomes reveals a higher rate of nucleotide divergence on microchromosomes than macrochromosomes.

    Genome Res 2005, 15:120-125. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Thomas JW, Cáceres M, Lowman JJ, Morehouse CB, Short ME, Baldwin EL, Maney DL, Martin CL: The chromosomal polymorphism linked to variation in social behavior in the white-throated sparrow (Zonotrichia albicollis) is a complex rearrangement and suppressor of recombination.

    Genetics 2008, 179:1455-1468. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Hoffman AA, Sgro CM, Weeks AR: Chromosomal inversion polymorphisms and adaptation.

    Trends Ecol Evol 2004, 19:482-488. PubMed Abstract | Publisher Full Text OpenURL

  64. Rocha GT, DeLucca EJ, DeSouza EB: Chromosome polymorphism due to pericentric inversion in Zonotrichia capensis (Emberizidae-Passeriformes-Aves).

    Genetica 1990, 80:201-207. Publisher Full Text OpenURL

  65. Shields GF: Chromosomal polymorphism common to several species on Junco (Aves).

    Can J Genet Cytol 1973, 15:461-471. OpenURL

  66. Shields GF: Comparative avian cytogenetics: A review.

    Condor 1982, 84:45-58. Publisher Full Text OpenURL

  67. Replogle K, Arnold AP, Ball GF, Band M, Bensch S, Brenowitz EA, Dong S, Drnevich J, Ferris M, George JM, et al.: The Songbird Neurogenomics (SoNG) Initiative: Community-based tools and strategies for study of brain gene function and evolution.

    BMC Genomics 2008, 9:131. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  68. Dawson DA, Burke T, Hansson B, Pandhal J, Hale MC, Hinten GN, Slate J: A predicted microsatellite map based on chicken-passerine sequence similarity.

    Mol Ecol 2006, 15:1299-1320. PubMed Abstract | Publisher Full Text OpenURL

  69. Dawson DA, Akesson M, Burke T, Pemberton JM, Slate J, Hansson B: Gene order and recombination rate in homologous chromosome regions of the chicken and a passerine bird.

    Mol Biol Evol 2007, 24:1537-1552. PubMed Abstract | Publisher Full Text OpenURL

  70. Lande R, Shannon S: The role of genetic variation in adaptive and population persistence in a changing environment.

    Evol 1996, 50:434-437. Publisher Full Text OpenURL

  71. Spielman D, Brook BW, Frankham R: Most species are not driven to extinction before genetic factors impact them.

    Proc Nat Acad Sci USA 2004, 101:15261-15264. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Kohn MH, Murphy WJ, Ostrander EA, Wayne RK: Genomics and conservation genetics.

    Trends Ecol Evol 2006, 21:629-637. PubMed Abstract | Publisher Full Text OpenURL

  73. Hedrick PW: Conservation genetics: where are we now?

    Trends Ecol Evol 2001, 16:629-636. Publisher Full Text OpenURL

  74. Reed DH, Frankham R: Correlation between fitness and genetic diversity.

    Conserv Biol 2003, 17:230-237. Publisher Full Text OpenURL

  75. Zink RM, Dittmann DL, Rootes WL: Mitochondrial DNA variation and the phylogeny of Zonotrichia.

    Auk 1991, 108:578-584. OpenURL

  76. Hailman JP: Analysis of aggression in white-throated sparrow types of different proportions.

    Bird-Banding 1975, 46:236-240. OpenURL

  77. Ficken RW, Ficken MS, Hailman JP: Differential aggression in genetically different morphs of the white-throated sparrow (Zonotrichia albicollis).

    Z Tierpsychol 1978, 46:43-57. PubMed Abstract OpenURL

  78. Watt DJ, Ralph CJ, Atkinson CT: The role of plumage polymorphism in dominance relationships of the white-throated sparrow.

    Auk 1984, 101:110-120. OpenURL

  79. Piper WH, Wiley RH: Distinguishing morphs of the white-throated sparrow in basic plumage.

    J Field Ornithol 1989, 60:73-83. OpenURL

  80. Piper WH, Wiley RH: Correlates of dominance in wintering white-throated sparrows: age, sex, and location.

    Anim Behav 1989, 37:298-310. Publisher Full Text OpenURL

  81. Whillans KV, Falls JB: Effects of male removal on parental care of female white-throated sparrows, Zonotrichia albicollis.

    Anim Behav 1990, 39:869-878. Publisher Full Text OpenURL

  82. Houtman AM, Falls JB: Negative assortative mating in the white-throated sparrow: the effect of mate choice and intrasexual competition.

    Anim Behav 1994, 48:377-383. Publisher Full Text OpenURL

  83. Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTLs.

    J Hered 2002, 93:77-78. PubMed Abstract | Publisher Full Text OpenURL

  84. Kozik A, Kochetkova E, Michelmore R: GenomePixelizer – a visualization program for comparative genomics within and between species.

    Bioinformatics 2002, 18:335-336. PubMed Abstract | Publisher Full Text OpenURL