Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome

Although there is growing interest to take genomics into the complex realms of natural populations, there is a general shortage of genomic resources and tools available for wild species. This applies not at least to birds, for which genomic approaches should be helpful to questions such as adaptation, speciation and population genetics. In this study, we describe a genome-wide reference set of conserved avian gene markers, broadly applicable across birds. By aligning protein-coding sequences from the recently assembled chicken genome with orthologous sequences in zebra finch, we identified particularly conserved exonic regions flanking introns of suitable size for subsequent amplification and sequencing. Primers were designed for 242 gene markers evenly distributed across the chicken genome, with a mean inter-marker interval of 4.2 Mb. Between 78% and 93% of the markers amplified a specific product in five species tested (chicken, peregrine falcon, collared flycatcher, great reed warbler and blue tit). Two hundred markers were sequenced in collared flycatcher, yielding a total of 122.41 kb of genomic DNA sequence (12096 bp coding sequence and 110 314 bp noncoding). Intron size of collared flycatcher and chicken was highly correlated, as was GC content. A polymorphism screening using these markers in a panel of 10 unrelated collared flycatchers identified 871 single nucleotide polymorphisms (pi = 0.0029) and 33 indels (mainly very short). Avian genome characteristics such as uniform genome size and low rate of syntenic rearrangements suggest that this marker set will find broad utility as a genome-wide reference resource for molecular ecological and population genomic analysis of birds. We envision that it will be particularly useful for obtaining large-scale orthologous targets in different species--important in, for instance, phylogenetics--and for large-scale identification of evenly distributed single nucleotide polymorphisms needed in linkage mapping or in studies of gene flow and hybridization.     

Eight PCR primers to amplify EST‐linked microsatellites in the Eastern oyster,Crassostrea virginica genome

Twenty‐four microsatellite repeat sequences were identified by screening a total of 4446 eastern oyster, Crassostrea virginica, expressed sequence tags. Polymerase chain reaction primers were designed to amplify 12 of these loci. After optimizing reaction parameters, eight loci showed high variability with consistent amplification that could be scored unambiguously. Ninety two C. virginica from the James River, VA, were genotyped at these loci. Number of alleles per locus ranged from 11 to 53, expected and observed heterozygosities ranged from 0.69 to 0.97, and from 0.30 to 0.99, respectively. Discrepancies between expected and observed heterozygosities were common and likely caused by null alleles.     

Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints

It is generally believed that the organization of avian genomes remains highly conserved in evolution as chromosome number is constant and comparative chromosome painting demonstrated there to be very few interchromosomal rearrangements. The recent sequencing of the zebra finch (Taeniopygia guttata) genome allowed an assessment of the number of intrachromosomal rearrangements between it and the chicken (Gallus gallus) genome, revealing a surprisingly high number of intrachromosomal rearrangements. With the publication of the turkey (Meleagris gallopavo) genome it has become possible to describe intrachromosomal rearrangements between these three important avian species, gain insight into the direction of evolutionary change and assess whether breakpoint regions are reused in birds. To this end, we aligned entire chromosomes between chicken, turkey and zebra finch, identifying syntenic blocks of at least 250 kb. Potential optimal pathways of rearrangements between each of the three genomes were determined, as was a potential Galliform ancestral organization. From this, our data suggest that around one-third of chromosomal breakpoint regions may recur during avian evolution, with 10% of breakpoints apparently recurring in different lineages. This agrees with our previous hypothesis that mechanisms of genome evolution are driven by hotspots of non-allelic homologous recombination.     

Variable visual habitats may influence the spread of colourful plumage across an avian hybrid zone

Several studies have shown that hybridization can be a creative process by acting as a conduit for the spread of adaptive traits between species, but few provide the mechanism that favours this spread. In the hybrid zone between the golden- (Manacus vitellinus) and white-collared (Manacus candei) manakins, sexual selection drives the introgression of golden/yellow plumage into the white species; however, the mechanism for the yellow male's mating advantage and the reasons why yellow plumage has not swept further into the white species remain mostly speculative. We quantified the colour properties of male plumage, the background and the ambient light at the hybrid zone, and allopatric golden and white populations. As measured by the perceived difference in colour between plumage and background, we found that yellow plumage appears more conspicuous than white plumage in the hybrid zone and allopatric golden-collar habitats, whereas white plumage appears more conspicuous than yellow plumage in the allopatric white-collared habitat. These results suggest a mechanism for the unidirectional spread of yellow plumage across the hybrid zone but slowed movement beyond it.     

Primers for 52 polymorphic regions in the Quercus rubra chloroplast, 47 of which amplify across 11 tracheophyte clades

Postglacial migration studies in Quercus rubra L. (northern red oak) are hampered by low levels of population differentiation in the widely used universal chloroplast primers. We sequenced the large single copy (LSC) regions of the Q. rubra and Quercus ellipsoidalis chloroplasts to enable us to query additional regions for future studies on migration and speciation. Using 454 sequencing of long-range PCR amplicons and Sanger sequencing for gap closure, we report 65 coding sequences from Q. rubra and 59 from Q. ellipsoidalis. Comparison of our de novo assembly of the LSC region sequence for Q. rubra to Q. rubra chloroplast sequence (NCBI Reference Sequence: NC_020152.1) from a different tree revealed 106 polymorphisms, all within intergenic regions, that can serve as tools for postglacial migration studies and taxonomic studies within the Lobatae. Sequence alignment for the 59 complete coding regions in common for theQ. rubrachloroplast reference sequence, our Q. rubra sequence and our Q. ellipsoidalis sequence revealed no sequence polymorphisms and no indels. We also report the 52 primer pairs we used for gap closure, including 53 new primer pairs not previously reported. We tested these 52 primer pairs against 11 species representing the Tracheophyta and detected 47 that produced amplicons in all 11 species. The new universal primers we have identified provide additional tools for resolving the taxonomic relationships among the congeneric taxa of forest trees in the temperate and subtropical forests of the Northern Hemisphere.     

A complete set of maize individual chromosome additions to the oat genome

All 10 chromosomes of maize (Zea mays, 2n = 2x = 20) were recovered as single additions to the haploid complement of oat (Avena sativa, 2n = 6x = 42) among F(1) plants generated from crosses involving three different lines of maize to eight different lines of oat. In vitro rescue culture of more than 4,300 immature F(1) embryos resulted in a germination frequency of 11% with recovery of 379 F(1) plantlets (8.7%) of moderately vigorous growth. Some F(1) plants were sectored with distinct chromosome constitutions among tillers of the same plant and also between root and shoot cells. Meiotic restitution facilitated development of un-reduced gametes in the F(1). Self-pollination of these partially fertile F(1) plants resulted in disomic additions (2n = 6x + 2 = 44) for maize chromosomes 1, 2, 3, 4, 6, 7, and 9. Maize chromosome 8 was recovered as a monosomic addition (2n = 6x + 1 = 43). Monosomic additions for maize chromosomes 5 and 10 to a haploid complement of oat (n = 3x + 1 = 22) were recovered several times among the F(1) plants. Although partially fertile, these chromosome 5 and 10 addition plants have not yet transmitted the added maize chromosome to F(2) offspring. We discuss the development and general utility of this set of oat-maize addition lines as a novel tool for maize genomics and genetics.     

A set of plastid DNA-specific universal primers for flowering plants

MatK and rbcL are recommended as the official barcode loci for higher plants but there remains a need for additional universal markers. We generated a series of 84 new universal primers targeting 42 plastid loci that all yielded single amplicons when applied to DNA templates from 19 diverse higher plant families. Marker utility ultimately depends on sequence variability, with rapidly evolving loci being useful for barcoding or biogeographic applications and more conserved loci being better suited to deep phylogeny reconstruction. Whereas excessive size variation is undesirable for many applications, modest size variability caused by indels and the sequence variation frequently associated with indels are highly desirable. We therefore performed a quick screen of the markers for size and sequence variation using pooled DNA templates from 96 taxonomically diverse species. All markers produced little or no size variation (consistent with the presence of minor indels). The seven regions exhibiting most size variation in pooled templates (rpl23&rpl2.1, 16S, 23S, 4.5S&5S, petB&D and rpl2, rpoC1 and trnK introns) were then amplified for all species individually, confirming the pooled template results. When the most variable loci (introns of trnK and rpoC1) were sequenced for all 96 species, a high level of sequence variation (nucleotide substitutions and indels) was observed among congeneric species groups for both loci. Both markers therefore have potential as supplementary barcode markers.     

A geographical spread of vaccine-resistance in avian influenza epidemics

Vaccination can be a useful tool for control of avian influenza outbreaks in poultry, but its use is reconsidered in most of the countries worldwide because of its negative effects on the disease control. One of the most important negative effects is the potential for emergence of vaccine-resistant viruses. Actually, in the vaccination program in China and Mexico, several vaccine-resistant strains were confirmed. Vaccine-resistant strains usually cause a loss of the protection effectiveness of vaccination. Therefore, a vaccination program that engenders the emergence of the resistant strain might promote the spread of the resistant strain and undermine the control of the infectious disease, even if the vaccination protects against the transmission of a vaccine-sensitive strain. We designed and analyzed a deterministic patch-structured model in heterogeneous areas (with or without vaccination) illustrating transmission of vaccine-sensitive and vaccine-resistant strains during a vaccination program. We found that the vaccination program can eradicate the vaccine-sensitive strain but lead to a prevalence of vaccine-resistant strain. Further, interestingly, the replacement of viral strain could occur in another area without vaccination through a migration of non-infectious individuals due to an illegal trade of poultry. It is also a novel result that only a complete eradication of both strains in vaccination area can achieve the complete eradication in another areas. Thus we can obtain deeper understanding of an effect of vaccination for better development of vaccination strategies to control avian influenza spread.     

A set of BAC clones spanning the human genome

Using the human bacterial artificial chromosome (BAC) fingerprint-based physical map, genome sequence assembly and BAC end sequences, we have generated a fingerprint-validated set of 32855 BAC clones spanning the human genome. The clone set provides coverage for at least 98% of the human fingerprint map, 99% of the current assembled sequence and has an effective resolving power of 79 kb. We have made the clone set publicly available, anticipating that it will generally facilitate FISH or array-CGH-based identification and characterization of chromosomal alterations relevant to disease.     

Development of primers to amplify mitochondrial DNA control region of Old World porcupines (subgenus Hystrix)

Eight primers were developed for the amplification of mitochondrial DNA control region of Old world porcupines (subgenus Hystrix). Successful amplifications of low‐quality DNA extracted from old (12 years old) and recent quills were performed, thus facilitating field sampling. Successful cross‐species amplifications were obtained for Hystrix africaeaustralis, H. cristata and H. indica. Length and structure of mitochondrial DNA control region were analysed and its usefulness as genetic marker for interspecific and population investigation was discussed.     

Identifying genomic regions for fine-mapping using genome scan meta-analysis (GSMA) to identify the minimum regions of maximum significance (MRMS) across populations

In order to detect linkage of the simulated complex disease Kofendrerd Personality Disorder across studies from multiple populations, we performed a genome scan meta-analysis (GSMA). Using the 7-cM microsatellite map, nonparametric multipoint linkage analyses were performed separately on each of the four simulated populations independently to determine p-values. The genome of each population was divided into 20-cM bin regions, and each bin was rank-ordered based on the most significant linkage p-value for that population in that region. The bin ranks were then averaged across all four studies to determine the most significant 20-cM regions over all studies. Statistical significance of the averaged bin ranks was determined from a normal distribution of randomly assigned rank averages. To narrow the region of interest for fine-mapping, the meta-analysis was repeated two additional times, with each of the 20-cM bins offset by 7 cM and 13 cM, respectively, creating regions of overlap with the original method. The 6-7 cM shared regions, where the highest averaged 20-cM bins from each of the three offsets overlap, designated the minimum region of maximum significance (MRMS). Application of the GSMA-MRMS method revealed genome wide significance (p-values refer to the average rank assigned to the bin) at regions including or adjacent to all of the simulated disease loci: chromosome 1 (p < 0.0001 for 160-167 cM, including D1), chromosome 3 (p-value < 0.0000001 for 287-294 cM, including D2), chromosome 5 (p-value < 0.001 for 0-7 cM, including D3), and chromosome 9 (p-value < 0.05 for 7-14 cM, the region adjacent to D4). This GSMA analysis approach demonstrates the power of linkage meta-analysis to detect multiple genes simultaneously for a complex disorder. The MRMS method enhances this powerful tool to focus on more localized regions of linkage.     

Design and evaluation of PCR primers to amplify bacterial 16S ribosomal DNA fragments used for community fingerprinting

Denaturing gradient gel electrophoresis of PCR-amplified 16S ribosomal DNA (rDNA) fragments has frequently been applied to the fingerprinting of natural bacterial populations (PCR/DGGE). In this study, sequences of bacterial universal primers frequently used in PCR/DGGE were compared with 16S rDNA sequences that represent recently proposed divisions in the domain Bacteria. We found mismatches in 16S rDNA sequences from some groups of bacteria. Inosine residues were then introduced into the bacterial universal primers to reduce amplification biases caused by these mismatches. Using the improved primers, phylotypes affiliated with Verrucomicrobia and candidate division OP11, were detected in DGGE fingerprints of groundwater populations, which have not been detected by PCR/DGGE with conventional universal primers.     

A genome walking strategy for the identification of nucleotide sequences adjacent to known regions

To identify the transposon insertion sites in a soil actinomycete, Saccharopolyspora spinosa, a genome walking approach, termed SPTA-PCR, was developed. In SPTA-PCR, a simple procedure consisting of TA cloning and a high stringency PCR, following the single primer-mediated, randomly-primed PCR, can eliminate non-target DNA fragments and obtain target fragments specifically. Using SPTA-PCR, the DNA sequence adjacent to the highly conserved region of lectin coding gene in onion plant, Allium chinense, was also cloned.     

Primaclade--a flexible tool to find conserved PCR primers across multiple species

Primaclade is a web-based application that accepts a multiple species nucleotide alignment file as input and identifies a set of polymerase chain reaction (PCR) primers that will bind across the alignment. Primaclade iteratively runs the Primer3 application for each alignment sequence and collates the results. Primaclade creates an HTML results page that recaps the original alignment, provides a consensus sequence and lists primers for each alignment area, with primers color-coded to reflect the level of degeneracy in the primer.     

Utility of chicken-specific microsatellite primers for mapping the turkey genome

As part of the University of Minnesota's initiative to map the turkey genome, we are currently evaluating chicken microsatellite loci for use in mapping the turkey genome. To date, 141 primer pairs have been tested for amplification at six different combinations of temperature and MgCl2 concentration. Microsatellite primer pairs from the Chicken Comprehensive Mapping Kit #2, and additional unpublished chromosome 1 and 2 primers were screened. Analyzable PCR products were produced from 78 of the 141 (55%) primer combinations. In the majority of cases (68%), PCR fragments obtained from the turkey were similar in size to respective chicken loci. The presence of dinucleotide repeats (CA/TG repeats) was determined by Southern hybridization with a (TG)15, oligonucleotide probe. Five of 12 (41.63%) turkey fragments hybridized under low stringency conditions. The length of the dinucleotide repeats in the turkey, relative to the chicken sequences, were found to correspond directly with hybridization intensity. Amplification of homologous loci was confirmed by direct sequencing and subsequent alignment of the turkey and chicken sequences. The results of this study indicate that the use of chicken-specific microsatellite primers will rapidly and significantly enhance construction of a genetic map for the turkey.