Comparative BAC-based mapping in the white-throated sparrow, a novel behavioral genomics model, using interspecies overgo hybridization

BACKGROUND The genomics era has produced an arsenal of resources from sequenced organisms allowing researchers to target species that do not have comparable mapping and sequence information. These new "non-model" organisms offer unique opportunities to examine environmental effects on genomic patterns and processes. Here we use comparative mapping as a first step in characterizing the genome organization of a novel animal model, the white-throated sparrow (Zonotrichia albicollis), which occurs as white or tan morphs that exhibit alternative behaviors and physiology. Morph is determined by the presence or absence of a complex chromosomal rearrangement. This species is an ideal model for behavioral genomics because the association between genotype and phenotype is absolute, making it possible to identify the genomic bases of phenotypic variation. FINDINGS We initiated a genomic study in this species by characterizing the white-throated sparrow BAC library via filter hybridization with overgo probes designed for the chicken, turkey, and zebra finch. Cross-species hybridization resulted in 640 positive sparrow BACs assigned to 77 chicken loci across almost all macro-and microchromosomes, with a focus on the chromosomes associated with morph. Out of 216 overgos, 36% of the probes hybridized successfully, with an average number of 3.0 positive sparrow BACs per overgo. CONCLUSIONS These data will be utilized for determining chromosomal architecture and for fine-scale mapping of candidate genes associated with phenotypic differences. Our research confirms the utility of interspecies hybridization for developing comparative maps in other non-model organisms

Molecular evolution of genes in avian genomes

Nam K, Mugal C, Nabholz B, et al. Molecular evolution of genes in avian genomes. Genome Biology. 2010;11(6): R68.Background: Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species. Results: We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an overrepresentation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton. Conclusions: Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes

Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor

Background The availability of multiple avian genome sequence assemblies greatly improves our ability to define overall genome organization and reconstruct evolutionary changes. In birds, this has previously been impeded by a near intractable karyotype and relied almost exclusively on comparative molecular cytogenetics of only the largest chromosomes. Here, novel whole genome sequence information from 21 avian genome sequences (most newly assembled) made available on an interactive browser (Evolution Highway) was analyzed. Results Focusing on the six best-assembled genomes allowed us to assemble a putative karyotype of the dinosaur ancestor for each chromosome. Reconstructing evolutionary events that led to each species' genome organization, we determined that the fastest rate of change occurred in the zebra finch and budgerigar, consistent with rapid speciation events in the Passeriformes and Psittaciformes. Intra- and interchromosomal changes were explained most parsimoniously by a series of inversions and translocations respectively, with breakpoint reuse being commonplace. Analyzing chicken and zebra finch, we found little evidence to support the hypothesis of an association of evolutionary breakpoint regions with recombination hotspots but some evidence to support the hypothesis that microchromosomes largely represent conserved blocks of synteny in the majority of the 21 species analyzed. All but one species showed the expected number of microchromosomal rearrangements predicted by the haploid chromosome count. Ostrich, however, appeared to retain an overall karyotype structure of 2n = 80 despite undergoing a large number (26) of hitherto un-described interchromosomal changes. Conclusions Results suggest that mechanisms exist to preserve a static overall avian karyotype/genomic structure, including the microchromosomes, with widespread interchromosomal change occurring rarely (e.g. in ostrich and budgerigar lineages). Of the species analyzed, the chicken lineage appeared to have undergone the fewest changes compared to the dinosaur ancestor

A multilocus assay reveals high nucleotide diversity and limited differentiation among Scandinavian willow grouse (Lagopus lagopus)

<p>Abstract</p> <p>Background</p> <p>There is so far very little data on autosomal nucleotide diversity in birds, except for data from the domesticated chicken and some passerines species. Estimates of nucleotide diversity reported so far in birds have been high (~10^-3) and a likely explanation for this is the generally higher effective population sizes compared to mammals. In this study, the level of nucleotide diversity has been examined in the willow grouse, a non-domesticated bird species from the order Galliformes, which also holds the chicken. The willow grouse (<it>Lagopus lagopus</it>) has an almost circumpolar distribution but is absent from Greenland and the north Atlantic islands. It primarily inhabits tundra, forest edge habitats and sub-alpine vegetation. Willow grouse are hunted throughout its range, and regionally it is a game bird of great cultural and economical importance.</p> <p>Results</p> <p>We sequenced 18 autosomal protein coding loci from approximately 15–18 individuals per population. We found a total of 127 SNP's, which corresponds to 1 SNP every 51 bp. 26 SNP's were amino acid replacement substitutions. Total nucleotide diversity (<it>π</it>_<it>t</it>) was between 1.30 × 10^-4and 7.66 × 10^-3(average <it>π</it>_{<it>t </it>}= 2.72 × 10^-3± 2.06 × 10^-3) and silent nucleotide diversity varied between 4.20 × 10^-4and 2.76 × 10^-2(average <it>π</it>_{<it>S </it>}= 9.22 × 10^-3± 7.43 × 10^-4). The synonymous diversity is approximately 20 times higher than in humans and two times higher than in chicken. Non-synonymous diversity was on average 18 times lower than the synonymous diversity and varied between 0 and 4.90 × 10^-3(average <it>π</it>_{<it>a </it>}= 5.08 × 10^-4± 7.43 × 10³), which suggest that purifying selection is strong in these genes. <it>F</it>_STvalues based on synonymous SNP's varied between -5.60 × 10^-4and 0.20 among loci and revealed low levels of differentiation among the four localities, with an overall value of <it>F</it>_ST= 0.03 (95% CI: 0.006 – 0.057) over 60 unlinked loci. Non-synonymous SNP's gave similar results. Low levels of linkage disequilibrium were observed within genes, with an average r²= 0.084 ± 0.110, which is expected for a large outbred population with no population differentiation. The mean per site per generation recombination parameter (ρ) was comparably high (0.028 ± 0.018), indicating high recombination rates in these genes.</p> <p>Conclusion</p> <p>We found unusually high levels of nucleotide diversity in the Scandinavian willow grouse as well as very little population structure among localities with up to 1647 km distance. There are also low levels of linkage disequilibrium within the genes and the population recombination rate is high, which is indicative of an old panmictic population, where recombination has had time to break up any haplotype blocks. The non-synonymous nucleotide diversity is low compared with the silent, which is in agreement with effective purifying selection, possibly due to the large effective population size.</p

Recombination dynamics of a human Y-chromosomal palindrome:rapid GC-biased gene conversion, multi-kilobase conversion tracts, and rare inversions

The male-specific region of the human Y chromosome (MSY) includes eight large inverted repeats (palindromes) in which arm-to-arm similarity exceeds 99.9%, due to gene conversion activity. Here, we studied one of these palindromes, P6, in order to illuminate the dynamics of the gene conversion process. We genotyped ten paralogous sequence variants (PSVs) within the arms of P6 in 378 Y chromosomes whose evolutionary relationships within the SNP-defined Y phylogeny are known. This allowed the identification of 146 historical gene conversion events involving individual PSVs, occurring at a rate of 2.9-8.4×10(-4) events per generation. A consideration of the nature of nucleotide change and the ancestral state of each PSV showed that the conversion process was significantly biased towards the fixation of G or C nucleotides (GC-biased), and also towards the ancestral state. Determination of haplotypes by long-PCR allowed likely co-conversion of PSVs to be identified, and suggested that conversion tract lengths are large, with a mean of 2068 bp, and a maximum in excess of 9 kb. Despite the frequent formation of recombination intermediates implied by the rapid observed gene conversion activity, resolution via crossover is rare: only three inversions within P6 were detected in the sample. An analysis of chimpanzee and gorilla P6 orthologs showed that the ancestral state bias has existed in all three species, and comparison of human and chimpanzee sequences with the gorilla outgroup confirmed that GC bias of the conversion process has apparently been active in both the human and chimpanzee lineages

Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor.

BACKGROUND: The availability of multiple avian genome sequence assemblies greatly improves our ability to define overall genome organization and reconstruct evolutionary changes. In birds, this has previously been impeded by a near intractable karyotype and relied almost exclusively on comparative molecular cytogenetics of only the largest chromosomes. Here, novel whole genome sequence information from 21 avian genome sequences (most newly assembled) made available on an interactive browser (Evolution Highway) was analyzed. RESULTS: Focusing on the six best-assembled genomes allowed us to assemble a putative karyotype of the dinosaur ancestor for each chromosome. Reconstructing evolutionary events that led to each species' genome organization, we determined that the fastest rate of change occurred in the zebra finch and budgerigar, consistent with rapid speciation events in the Passeriformes and Psittaciformes. Intra- and interchromosomal changes were explained most parsimoniously by a series of inversions and translocations respectively, with breakpoint reuse being commonplace. Analyzing chicken and zebra finch, we found little evidence to support the hypothesis of an association of evolutionary breakpoint regions with recombination hotspots but some evidence to support the hypothesis that microchromosomes largely represent conserved blocks of synteny in the majority of the 21 species analyzed. All but one species showed the expected number of microchromosomal rearrangements predicted by the haploid chromosome count. Ostrich, however, appeared to retain an overall karyotype structure of 2n=80 despite undergoing a large number (26) of hitherto un-described interchromosomal changes. CONCLUSIONS: Results suggest that mechanisms exist to preserve a static overall avian karyotype/genomic structure, including the microchromosomes, with widespread interchromosomal change occurring rarely (e.g., in ostrich and budgerigar lineages). Of the species analyzed, the chicken lineage appeared to have undergone the fewest changes compared to the dinosaur ancestor.The authors would like to thank Alain Vignal and Thomas Faraut of INRA Toulouse (France) for access to the duck chromosome assembly data. This research was funded in part by PL-Grid Infrastructure (DML), Biotechnology and Biological Sciences Research Council BB/K008161 (DML, DKG), BB/K008226/1 (DML), BB/J010170/1 (DML, MF) and a knowledge transfer partnership award (DKG and Cytocell Ltd). The authors are grateful to Malcolm Ferguson-Smith’s lab (Cambridge, UK) for producing the flow-sorted chicken microchromosome paints. We also thank Cytocell Ltd (Cambridge, UK) for technical support in FISH technologies.This article was originally published in BMC Genomics 2014, 15:1060 doi:10.1186/1471-2164-15-106

Autosomal and Z-linked microsatellite markers enhanced for cross-species utility and assessed in a range of birds, including species of conservation concern

Microsatellite markers were designed to be of utility for genotyping multiple species of birds, including those of conservation concern, hence saving resources and enabling species/genome comparisons. We used the proven approach of Dawson et al. (Mol Ecol Resour 10:475–494, 2010) and assessed markers in multiple species, including nine species of conservation interest. We ensured both primer sequences matched multiple species (13 loci) or designed primer sets from expressed sequence tags (2 loci). Eleven primer sets were 100 % identical to the zebra finch (Taeniopygia guttata) and a second passerine species and/or the chicken (Gallus gallus). All 15 loci were polymorphic when assessed in a non-source species (Gouldian finch, Erythrura gouldiae) suggesting utility in multiple species. Four of the five Z-linked loci were assessed in at least nine additional species each (including ratites). All were variable in multiple species, demonstrating cross-species utility and potential for identifying Z chromosome rearrangements

The genome of a songbird

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chickenthe only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat- based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour. © 2010 Macmillan Publishers Limited. All rights reserved

Evidence for Mito-Nuclear and Sex-Linked Reproductive Barriers between the Hybrid Italian Sparrow and Its Parent Species

Studies of reproductive isolation between homoploid hybrid species and their parent species have rarely been carried out. Here we investigate reproductive barriers between a recently recognized hybrid bird species, the Italian sparrow Passer italiae and its parent species, the house sparrow P. domesticus and Spanish sparrow P. hispaniolensis. Reproductive barriers can be difficult to study in hybrid species due to lack of geographical contact between taxa. However, the Italian sparrow lives parapatrically with the house sparrow and both sympatrically and parapatrically with the Spanish sparrow. Through whole-transcriptome sequencing of six individuals of each of the two parent species we identified a set of putatively parent species-diagnostic single nucleotide polymorphism (SNP) markers. After filtering for coverage, genotyping success (>97%) and multiple SNPs per gene, we retained 86 species-informative, genic, nuclear and mitochondrial SNP markers from 84 genes for analysis of 612 male individuals. We show that a disproportionately large number of sex-linked genes, as well as the mitochondria and nuclear genes with mitochondrial function, exhibit sharp clines at the boundaries between the hybrid and the parent species, suggesting a role for mito-nuclear and sex-linked incompatibilities in forming reproductive barriers. We suggest that genomic conflict via interactions between mitochondria and sex-linked genes with mitochondrial function ("mother's curse") at one boundary and centromeric drive at the other may best explain our findings. Hybrid speciation in the Italian sparrow may therefore be influenced by mechanisms similar to those involved in non-hybrid speciation, but with the formation of two geographically separated species boundaries instead of one. Spanish sparrow alleles at some loci have spread north to form reproductive barriers with house sparrows, while house sparrow alleles at different loci, including some on the same chromosome, have spread in the opposite direction to form barriers against Spanish sparrows

The Dynamics of Incomplete Lineage Sorting across the Ancient Adaptive Radiation of Neoavian Birds

The diversification of neoavian birds is one of the most rapid adaptive radiations of extant organisms. Recent whole-genome sequence analyses have much improved the resolution of the neoavian radiation and suggest concurrence with the Cretaceous-Paleogene (K-Pg) boundary, yet the causes of the remaining genome-level irresolvabilities appear unclear. Here we show that genome-level analyses of 2,118 retrotransposon presence/absence markers converge at a largely consistent Neoaves phylogeny and detect a highly differential temporal prevalence of incomplete lineage sorting (ILS), i.e., the persistence of ancestral genetic variation as polymorphisms during speciation events. We found that ILS-derived incongruences are spread over the genome and involve 35% and 34% of the analyzed loci on the autosomes and the Z chromosome, respectively. Surprisingly, Neoaves diversification comprises three adaptive radiations, an initial near-K-Pg super-radiation with highly discordant phylogenetic signals from near-simultaneous speciation events, followed by two post-K-Pg radiations of core landbirds and core waterbirds with much less pronounced ILS. We provide evidence that, given the extreme level of up to 100% ILS per branch in super-radiations, particularly rapid speciation events may neither resemble a fully bifurcating tree nor are they resolvable as such. As a consequence, their complex demographic history is more accurately represented as local networks within a species tree