Search

The clouded apollo (Parnassius mnemosyne) used to have a wide distribution in Fennoscandia. Recent population declines have, however, led to regional extinctions and in Sweden it is currently one of the most endangered butterflies, confined to three geographically separated metapopulations: Blekinge, Roslagen and Västernorrland. Especially the Blekinge population has declined dramatically and few imagines have been observed during recent census efforts (< 10 in some localities). The clouded apollo is subject to a species action plan which includes both habitat restorations and captive breeding to produce individuals for release and reintroductions. Here, we apply whole-genome resequencing of clouded apollo individuals collected in the three natural populations and the captive population in Sweden and apply population genomic approaches to get a better understanding of the genetic structure and levels of genetic diversity in the species. We find that the clouded apollo populations in the different geographic regions have similar, but comparatively low levels of genetic diversity and we find evidence for significant genetic differentiation between the northernmost population and the populations in southern Sweden. Additional analysis, including previously available mitochondrial data, unveil that a bi-directional re-colonization of Fennoscandia after the latest glacial maximum most likely is the explanation for the considerable differentiation between some Swedish populations. Finally, we find evidence for population sub-structure in one of the Swedish populations. The results provide insights into the genetic consequences of population size declines and fragmentation in general and provide important information for direct conservation actions for the clouded apollo in Sweden in particular.

Migration is typically associated with risk and uncertainty at the population level, but little is known about its cost¿benefit trade-offs at the species level. Migratory insects in particular often exhibit strong demographic fluctuations due to local bottlenecks and outbreaks. Here, we use genomic data to investigate levels of heterozygosity and long-term population size dynamics in migratory insects, as an alternative to classical local and short-term approaches such as regional field monitoring. We analyse whole-genome sequences from 97 Lepidoptera species and show that individuals of migratory species have significantly higher levels of genome-wide heterozygosity, a proxy for effective population size, than do nonmigratory species. Also, we contribute whole-genome data for one of the most emblematic insect migratory species, the painted lady butterfly (Vanessa cardui), sampled across its worldwide distributional range. This species exhibits one of the highest levels of genomic heterozygosity described in Lepidoptera (2.95 ± 0.15%). Coalescent modelling (PSMC) shows historical demographic stability in V. cardui, and high effective population size estimates of 2 -20 million individuals 10,000 years ago. The study reveals that the high risks associated with migration and local environmental fluctuations do not seem to decrease overall genetic diversity and demographic stability in migratory Lepidoptera. We propose a "compensatory" demographic model for migratory r-strategist organisms in which local bottlenecks are counterbalanced by reproductive success elsewhere within their typically large distributional ranges. Our findings highlight that the boundaries of populations are substantially different for sedentary and migratory insects, and that, in the latter, local and even regional field monitoring results may not reflect whole population dynamics. Genomic diversity patterns may elucidate key aspects of an insect's migratory nature and population dynamics at large spatiotemporal scales., This work was funded by the National Geographic Society (grant WW1-300R-18) and by the grant PID2020-117739GA-I00 from MCIN/AEI/10.13039/501100011033 to G.T., by fellowship FPU19/01593 to A.G.-B., by project LINKA20399 from the CSIC i-LINK-2021 programme to G.T., N.E.P., R.V. and N.B., by grants from the Putnam Expeditionary Fund of the Museum of Comparative Zoology to G.T. and N.E.P., by projects PID2019-107078GB-I00/MCIN/AEI/10.13039/501100011033 and 2017-SGR-991 (Generalitat de Catalunya) to R.V. and G.T., by the University of Malaya (grant H50001-A-000027) to K.G.C., by the Swedish Collegium for Advanced Science (SCAS; Natural Sciences Programme, Knut and Alice Wallenberg Foundation, Postdoc funding) to D.S., by NSF grant DEB-1541560 to N.E.P., and by the Swedish Research Council FORMAS (grant 2019-00670) to N.B., G.T. and R.V., 1 INTRODUCTION2 MATERIALS AND METHODS2.1 Heterozygosity estimates2.2 Behavioural scoring2.3 Comparative analyses2.4 Vanessa cardui sampling and sequencing2.5 Vanessa cardui reference-based assemblies2.6 Vanessa cardui demographic inference3 RESULTS3.1 Heterozygosity rates in Lepidoptera3.2 Vanessa cardui demographic history4 DISCUSSION4.1 Genomic heterozygosity, effective population size and migration4.2 Short-term migratory cost but long-term demographic stability4.3 Delimiting populations in migratory insect speciesAUTHOR CONTRIBUTIONSACKNOWLEDGEMENTSCONFLICT OF INTEREST

Seasonal environmental fluctuations provide formidable challenges for living organisms, especially small ectotherms such as butterflies. A common strategy to cope with harsh environments is to enter diapause, but some species avoid the unsuitable conditions by migrating. Despite a growing understanding of migration in the life cycles of some butterfly species, it is unknown how individuals register and store environmental cues to determine whether and where to migrate. Here, we explored how competition and host plant availability during larval development affect patterns of DNA methylation in the migratory painted lady (Vanessa cardui) butterfly. We identify a set of potentially functional methylome shifts associated with differences in the environment, indicating that DNA methylation is involved in the response to different conditions during larval development. By analysing the transcriptome for the same samples used for methylation profiling, we also uncovered a non-linear relationship between gene body methylation and gene expression. Our results provide a starting point for understanding the interplay between DNA methylation and gene expression in butterflies in general and how differences in environmental conditions during development can trigger unique epigenetic marks that might be important for behavioural decisions in the adult stage.

Genetic recombination is a key molecular mechanism that has profound implications on both micro- and macro-evolutionary processes. However, the determinants of recombination rate variation in holocentric organisms are poorly understood, in particular in Lepidoptera (moths and butterflies). The wood white butterfly (Leptidea sinapis) shows considerable intraspecific variation in chromosome numbers and is a suitable system for studying regional recombination rate variation and its potential molecular underpinnings. Here, we developed a large wholegenome resequencing data set from a population of wood whites to obtain high-resolution recombination maps using linkage disequilibrium information. The analyses revealed that larger chromosomes had a bimodal recombination landscape, potentially due to interference between simultaneous chiasmata. The recombination rate was significantly lower in subtelomeric regions, with exceptions associated with segregating chromosome rearrangements, showing that fissions and fusions can have considerable effects on the recombination landscape. There was no association between the inferred recombination rate and base composition, supporting a negligible influence of GC-biased gene conversion in butterflies. We found significant but variable associations between the recombination rate and the density of different classes of transposable elements (TEs), most notably a significant enrichment of SINEs in genomic regions with higher recombination rate. Finally, the analyses unveiled significant enrichment of genes involved in farnesyltranstransferase activity in recombination cold-spots, potentially indicating that expression of transferases can inhibit formation of chiasmata during meiotic division. Our results provide novel information about recombination rate variation in holocentric organisms and has particular implications for forthcoming research in population genetics, molecular/genome evolution and speciation.

Migration is typically associated with risk and uncertainty at the population level, but little is known about its cost and benefit trade-offs at the species level. Migratory insects often exhibit strong demographic fluctuations due to local bottlenecks and outbreaks. Here, we use genomic data to investigate levels of heterozygosity and long-term population size dynamics in migratory insects, as an alternative to classical local and short-term approaches such as regional field monitoring. We analyze whole-genome sequences from 97 Lepidoptera species and show that migratory species have significantly higher levels of genome-wide heterozygosity, a proxy for effective population size, than non-migratory species. Also, we contribute whole-genome data for one of the most emblematic insect migratory species, the painted lady butterfly (Vanessa cardui), sampled across its worldwide distribution range. This species exhibits one of the highest levels of heterozygosity described in Lepidoptera (2.95% ± 0.15). Coalescent modeling (PSMC) shows historical demographic stability in V. cardui, and high effective population size estimates of 2 to 20 million individuals 10,000 years ago. The study reveals that the high risks associated with migration and local environmental fluctuations do not seem to decrease overall genetic diversity and demographic stability in migratory Lepidoptera. We propose a “compensatory” demographic model for migratory r-strategist organisms in which local bottlenecks are counterbalanced by reproductive success elsewhere within their typically large distributional ranges. Our findings highlight that the delimitations of populations are substantially different for sedentary and migratory insects, and that, in the latter, local and even regional field monitoring results may not reflect whole population dynamics. Genomic diversity patterns may elucidate key aspects of an insect’s migratory nature and population dynamics at large spatiotemporal scales.

The clouded apollo Clouded Apollo (Parnassius mnemosyne) is formerly used to have a widely distributed, Eurasian butterfly species that shows declining population numbers in many different regions. Within Sweden, both thewider distribution range and the population size has decreased dramatically over the latest decades and the species is now only found in Fenno Scandia but is in Sweden at present one of the most endangered butterflies confined to three geographically isolated regions. separated regional metapopulations: Blekinge, Roslagen and Västernorrland. Especially the Blekinge population have in recent years declined and show low numbers of adults (N

Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements., Financial support for this project was provided by FORMAS (Research grant 2019-00670 to N.B.) and The Swedish Collegium for Advanced Science (Natural Sciences Programme, Knut and Alice Wallenberg Foundation, Postdoc funding for D.S.). R.V. was supported by the grant PID2019-107078GB-I00 funded by MCIN/AEI/10.13039/501100011033. G.T. was supported by the grant PID2020-117739GA-I00 funded by MCIN/AEI/10.13039/501100011033 and by “La Caixa” Foundation (ID 100010434) through the grant LCF/BQ/PR19/11700004., 1. Introduction 2. Results 2.1. Linkage map and genome annotation 2.2. Synteny 2.3. Gene family evolution 2.4. Patterns of recombination rate variation 2.4.1. Global and chromosome specific recombination rates 2.4.2. Intra-chromosomal variation in recombination rate and associations with genomic elements 3. Discussion 3.1. The genome of the painted lady butterfly 3.2. Sex chromosomes 3.3. Gene family analysis 3.4. Patterns of recombination rate variation 4. Conclusions 5. Methods 5.1. Linkage map 5.1.1. Sampling and DNA-extraction 5.1.2. Building the linkage map 5.2. Genome annotation and whole genome statistics 5.2.1. Genome assembly statistics 5.2.2. Gene and repeat annotation 5.3. Manual curation 5.4. Gene family evolution 5.5. Gene ontology enrichment 5.6. Recombination rate analysis 5.6.1. Chromosome level analysis 5.6.2. Window-based analysis Data access CRediT authorship contribution statement Declaration of Competing Interest Acknowledgements Appendix A. Supplementary data Data availability References

Recombination reshuffles the alleles of a population through crossover and gene conversion. These mechanisms have considerable consequences on the evolution and maintenance of genetic diversity. Crossover, for example, can increase genetic diversity by breaking the linkage between selected and nearby neutral variants. Bias in favor of G or C alleles during gene conversion may instead promote the fixation of one allele over the other, thus decreasing diversity. Mutation bias from G or C to A and T opposes GC-biased gene conversion (gBGC). Less recognized is that these two processes may –when balanced– promote genetic diversity. Here we investigate how gBGC and mutation bias shape genetic diversity patterns in wood white butterflies (Leptideasp.). This constitutes the first in-depth investigation of gBGC in butterflies. Using 60 re-sequenced genomes from six populations of three species, we find substantial variation in the strength of gBGC across lineages. When modeling the balance of gBGC and mutation bias and comparing analytical results with empirical data, we reject gBGC as the main determinant of genetic diversity in these butterfly species. As alternatives, we consider linked selection and GC content. We find evidence that high values of both reduce diversity. We also show that the joint effects of gBGC and mutation bias can give rise to a diversity pattern which resembles the signature of linked selection. Consequently, gBGC should be considered when interpreting the effects of linked selection on levels of genetic diversity.

In temperate latitudes, many insects enter diapause (dormancy) during the cold season, a period during which developmental processes come to a standstill. The wood white (Leptidea sinapis) is a butterfly species distributed across western Eurasia that shows photoperiod-induced diapause with variation in critical day-length across populations at different latitudes. We assembled transcriptomes and estimated gene expression levels at different developmental stages in experimentally induced directly developing and diapausing cohorts of a single Swedish population of L. sinapis to investigate the regulatory mechanisms underpinning diapause initiation. Different day lengths resulted in expression changes of developmental genes and affected the rate of accumulation of signal molecules, suggesting that diapause induction might be controlled by increased activity of monoamine neurotransmitters in larvae reared under short-day light conditions. Expression differences between light treatment groups of two monoamine regulator genes (DDC and ST) were observed already in instar III larvae. Once developmental pathways were irreversibly set at instar V, a handful of genes related to dopamine production were differentially expressed leading to a significant decrease in expression of global metabolic genes and increase in expression of genes related to fatty acid synthesis and sequestration. This is in line with a time-dependent (hour-glass) model of diapause regulation where a gradual shift in the concentration of monoamine neurotransmitters and their metabolites during development of larvae under short-day conditions leads to increased storage of fat, decreased energy expenditures, and ultimately developmental stasis at the pupal stage.

Changes in interacting cis- and trans-regulatory elements are important candidates for Dobzhansky-Muller hybrid incompatibilities and may contribute to hybrid dysfunction by giving rise to misexpression in hybrids. To gain insight into the molecular mechanisms and determinants of gene expression evolution in natural populations, we analyzed the transcriptome from multiple tissues of two recently diverged Ficedula flycatcher species and their naturally occurring F(1) hybrids. Differential gene expression analysis revealed that the extent of differentiation between species and the set of differentially expressed genes varied across tissues. Common to all tissues, a higher proportion of Z-linked genes than autosomal genes showed differential expression, providing evidence for a fast-Z effect. We further found clear signatures of hybrid misexpression in brain, heart, kidney, and liver. However, while testis showed the highest divergence of gene expression among tissues, it showed no clear signature of misexpression in F(1) hybrids, even though these hybrids were found to be sterile. It is therefore unlikely that incompatibilities between cis-trans regulatory changes explain the observed sterility. Instead, we found evidence that cis-regulatory changes play a significant role in the evolution of gene expression in testis, which illustrates the tissue-specific nature of cis-regulatory evolution bypassing constraints associated with pleiotropic effects of genes.

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large-scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence—whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole-genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)—a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (F) regions generally had low genetic diversity (θ), but increased absolute divergence (D) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression., This work was supported by a junior research grant from the Swedish Research Council (VR) to NB. The authors acknowledge support from the National Genomics Infrastructure in Stockholm and Uppsala funded by the Science for Life Laboratory, the Knut and Alice Wallenberg Foundation and the Swedish Research Council, and SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science for assistance with massively parallel sequencing and access to the UPPMAX computational infrastructure. RV was supported by project CGL2016‐76322‐P (AEI/FEDER, UE).

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock‐like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between‐species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression PST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20‐kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated.

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM/Mb and correlated closely with chromosome size, from 2 cM/Mb for chromosomes >100 Mb to >10 cM/Mb for chromosomes

With the access to draft genome sequence assemblies and whole-genome resequencing data from population samples, molecular ecology studies will be able to take truly genome-wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single-nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10-fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F1 hybrids but no later-generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as

The complexity of eukaryote genomes makes assembly errors inevitable in the process of constructing reference genomes. Next-generation sequencing (NGS) could provide an efficient way to validate previously assembled genomes. Here, we exploited NGS data to interrogate the chicken reference genome and identified 35 pairs of nearly identical regions with >99.5 % sequence similarity and a median size of 109 kb. Several lines of evidence, including read depth, the composition of junction sequences, and sequence similarity, suggest that these regions present genome assembly errors and should be excluded from forthcoming genomic studies.

Characterization of the genetic basis of fitness traits in natural populations is important for understanding how organisms adapt to the changing environment and to novel events, such as epizootics. However, candidate fitness‐influencing loci, such as regulatory regions, are usually unavailable in nonmodel species. Here, we analyze sequence data from targeted resequencing of the cis‐regulatory regions of three candidate genes for disease resistance (CD74, HSP90α, and LCP1) in populations of the house finch (Carpodacus mexicanus) historically exposed (Alabama) and naïve (Arizona) to Mycoplasma gallisepticum. Our study, the first to quantify variation in regulatory regions in wild birds, reveals that the upstream regions of CD74 and HSP90α are GC‐rich, with the former exhibiting unusually low sequence variation for this species. We identified two SNPs, located in a GC‐rich region immediately upstream of an inferred promoter site in the gene HSP90α, that were significantly associated with Mycoplasma pathogen load in the two populations. The SNPs are closely linked and situated in potential regulatory sequences: one in a binding site for the transcription factor nuclear NFYα and the other in a dinucleotide microsatellite ((GC)6). The genotype associated with pathogen load in the putative NFYα binding site was significantly overrepresented in the Alabama birds. However, we did not see strong effects of selection at this SNP, perhaps because selection has acted on standing genetic variation over an extremely short time in a highly recombining region. Our study is a useful starting point to explore functional relationships between sequence polymorphisms, gene expression, and phenotypic traits, such as pathogen resistance that affect fitness in the wild.

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Identifying causal genetic variants underlying heritable phenotypic variation is a long-standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole-genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild-caught zebra finches. Then, we validated the most promising SNP-phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype-phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.

Anolis carolinensis is an emerging model species and the sole member of its genus native to the United States. Considerable morphological and physiological variation has been described in the species, and the recent sequencing of its genome makes it an attractive system for studies of genome variation. To inform future studies of molecular and phenotypic variation within A. carolinensis, a rigorous account of intraspecific population structure and relatedness is needed. Here, we present the most extensive phylogeographic study of this species to date. Phylogenetic analyses of mitochondrial DNA sequence data support the previous hypothesis of a western Cuban origin of the species. We found five well-supported, geographically distinct mitochondrial haplotype clades throughout the southeastern United States. Most Florida populations fall into one of three divergent clades, whereas the vast majority of populations outside Florida belong to a single, shallowly diverged clade. Genetic boundaries do not correspond to major rivers, but may reflect effects of Pleistocene glaciation events and the Appalachian Mountains on migration and expansion of the species. Phylogeographic signal should be examined using nuclear loci to complement these findings.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

Theoretical and empirical data suggest that genes located on sex chromosomes may play an important role both for sexually selected traits and for traits involved in the build-up of hybrid incompatibilities. We investigated patterns of genetic variation in 73 genes located on the Z chromosomes of two species of the flycatcher genus Ficedula, the pied flycatcher and the collared flycatcher. Sequence data were evaluated for signs of selection potentially related to genomic differentiation in these young sister species, which hybridize despite reduced fitness of hybrids. Seven loci were significantly more divergent between the two species than expected under neutrality and they also displayed reduced nucleotide diversity, consistent with having been influenced by directional selection. Two of the detected candidate regions contain genes that are associated with plumage coloration in birds. Plumage characteristics play an important role in species recognition in these flycatchers suggesting that the detected genes may have been involved in the evolution of sexual isolation between the species.

Chromosomal rearrangements and copy number variants (CNVs) play key roles in genome evolution and genetic disease; however, the molecular mechanisms underlying these types of structural genomic variation are not fully understood. The availability of complete genome sequences for two bird species, the chicken and the zebra finch, provides, for the first time, an ideal opportunity to analyze the relationship between structural genomic variation (chromosomal and CNV) and recombination on a genome-wide level. The aims of this study were therefore threefold: (1) to combine bioinformatics, physical mapping to produce comprehensive comparative maps of the genomes of chicken and zebra finch. In so doing, this allowed the identification of evolutionary chromosomal rearrangements distinguishing them. The previously reported interchromosomal conservation of synteny was confirmed, but a larger than expected number of intrachromosomal rearrangements were reported; (2) to hybridize zebra finch genomic DNA to a chicken tiling path microarray and identify CNVs in the zebra finch genome relative to chicken; 32 interspecific CNVs were identified; and (3) to test the hypothesis that there is an association between CNV, chromosomal rearrangements, and recombination by correlating data from (1) and (2) with recombination rate data from a high-resolution genetic linkage map of the zebra finch. We found a highly significant association of both chromosomal rearrangements and CNVs with elevated recombination rates. The results thus provide support for the notion of recombination-based processes playing a major role in avian genome evolution. © 2010 by Cold Spring Harbor Laboratory Press.

Next-generation sequencing technology provides an attractive means to obtain large-scale sequence data necessary for comparative genomic analysis. To analyse the patterns of mutation rate variation and selection intensity across the avian genome, we performed brain transcriptome sequencing using Roche 454 technology of 10 different non-model avian species. Contigs from de novo assemblies were aligned to the two available avian reference genomes, chicken and zebra finch. In total, we identified 6499 different genes across all 10 species, with approximately 1000 genes found in each full run per species. We found evidence for a higher mutation rate of the Z chromosome than of autosomes (male-biased mutation) and a negative correlation between the neutral substitution rate (d(S)) and chromosome size. Analyses of the mean d(N)/d(S) ratio (omega) of genes across chromosomes supported the Hill-Robertson effect (the effect of selection at linked loci) and point at stochastic problems with omega as an independent measure of selection. Overall, this study demonstrates the usefulness of next-generation sequencing for obtaining genomic resources for comparative genomic analysis of non-model organisms.

Abstract We have used phylogeographic methods to investigate the genetic structure and population history of the endangered Himalayan snowcock (Tetraogallus himalayensis) in northwestern China. The mitochondrial cytochrome b gene was sequenced of 102 individuals sampled throughout the distribution range. In total, we found 26 different haplotypes defined by 28 polymorphic sites. Phylogenetic analyses indicated that the samples were divided into two major haplogroups corresponding to one western and one eastern clade. The divergence time between these major clades was estimated to be approximately one million years. An analysis of molecular variance showed that 40% of the total genetic variability was found within local populations, 12% among populations within regional groups and 48% among groups. An analysis of the demographic history of the populations suggested that major expansions have occurred in the Himalayan snowcock populations and these correlate mainly with the first and the second largest glaciations during the Pleistocene. In addition, the data indicate that there was a population expansion of the Tianshan population during the uplift of the Qinghai-Tibet Plateau, approximately 2 million years ago.

By taking advantage of a recently developed reference marker set for avian genome analysis we have constructed a gene-based genetic map of the collared flycatcher, an important “ecological model” for studies of life-history evolution, sexual selection, speciation, and quantitative genetics. A pedigree of 322 birds from a natural population was genotyped for 384 single nucleotide polymorphisms (SNPs) from 170 protein-coding genes and 71 microsatellites. Altogether, 147 gene markers and 64 microsatellites form 33 linkage groups with a total genetic distance of 1787 cM. Male recombination rates are, on average, 22% higher than female rates (total distance 1982 vs. 1627 cM). The ability to anchor the collared flycatcher map with the chicken genome via the gene-based SNPs revealed an extraordinary degree of both synteny and gene-order conservation during avian evolution. The great majority of chicken chromosomes correspond to a single linkage group in collared flycatchers, with only a few cases of inter- and intrachromosomal rearrangements. The rate of chromosomal diversification, fissions/fusions, and inversions combined is thus considerably lower in birds (0.05/MY) than in mammals (0.6–2.0/MY). A dearth of repeat elements, known to promote chromosomal breakage, in avian genomes may contribute to their stability. The degree of genome stability is likely to have important consequences for general evolutionary patterns and may explain, for example, the comparatively slow rate by which genetic incompatibility among lineages of birds evolves.

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.

The HINTW gene on the female-specific W chromosome of chicken and other birds is amplified and present in numerous copies. Moreover, as HINTW is distinctly different from its homolog on the Z chromosome (HINTZ), is a candidate gene in avian sex determination, and evolves rapidly under positive selection, it shows several common features to ampliconic and testis-specific genes on the mammalian Y chromosome. A phylogenetic analysis within galliform birds (chicken, turkey, quail, and pheasant) shows that individual HINTW copies within each species are more similar to each other than to gene copies of related species. Such convergent evolution is most easily explained by recurrent events of gene conversion, the rate of which we estimated at 10(-6)-10(-5) per site and generation. A significantly higher GC content of HINTW than of other W-linked genes is consistent with biased gene conversion increasing the fixation probability of mutations involving G and C nucleotides. Furthermore, and as a likely consequence, the neutral substitution rate is almost twice as high in HINTW as in other W-linked genes. The region on W encompassing the HINTW gene cluster is not covered in the initial assembly of the chicken genome, but analysis of raw sequence reads indicates that gene copy number is significantly higher than a previous estimate of 40. While sexual selection is one of several factors that potentially affect the evolution of ampliconic, male-specific genes on the mammalian Y chromosome, data from HINTW provide evidence that gene amplification followed by gene conversion can evolve in female-specific chromosomes in the absence of sexual selection. The presence of multiple and highly similar copies of HINTW may be related to protein function, but, more generally, amplification and conversion offers a means to the avoidance of accumulation of deleterious mutations in nonrecombining chromosomes.

There is increasing evidence that dosage compensation is not a ubiquitous feature following sex chromosome evolution, especially not in organisms where females are the heterogametic sex, like in birds. Even when it occurs, compensation can be incomplete and limited to dosage-sensitive genes. However, previous work has mainly studied transcriptional regulation of sex-linked genes, which may not reflect expression at the protein level. Here, we used liquid chromatography‐ tandem mass spectrometry to detect and quantify expressed levels of more than 2,400 proteins in ten different tissues of male and female chicken embryos. For comparison, transcriptome sequencing was performed in the same individuals, five of each sex. The proteomic analysis revealed that dosage compensation was incomplete, with a mean male-to-female (M:F) expression ratio of Z-linked genes of 1.32 across tissues, similar to that at the RNA level (1.29). The mean Z chromosome-to-autosome expression ratio was close to 1 in males and lower than 1 in females, consistent with partly reduced Z chromosome expression in females. Although our results exclude a general mechanism for chromosome-wide dosage compensation at translation, 30% of all proteins encoded from Z-linked genes showed a significant change in the M:F ratio compared with the corresponding ratio at the RNA level. This resulted in a pattern where some genes showed balanced expression between sexes and some close to 2-fold higher expression in males. This suggests that proteomic analyses will be necessary to reveal a more complete picture of gene regulation and sex chromosome evolution.

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-1060

Adaptiveevolutionistheprocesswherebymutationsthat provide the carrier with a selective advantage increase in frequency in a population via the process of natural selection. Passerines are widespread, common and longterm targets for field study and they demonstrate a copious diversity in physiological and morphological adaptations to varying habitats, for example, beak size and wing shape, and they are, therefore, an important study system to understand adaptive evolution. Recent technological advancements have made it easier to investigate the mechanistic and evolutionary underpinnings of adaptive evolution by allowing genome sequence data to be generated in almost any species of interest. However, it is important to assess the contribution of neutral forces like demographic events and GCbiased gene conversion before concluding that selection has shaped the patterns observed in genomic data. Initial analyses in passerines have identified candidate genes that might be involved in, for example, song learning, beak morphology, disease resistance, high-altitude adaptation and exploratory behavior, but functional verifications are needed to establish a causative relationship between the identified genes and the traits.

With its plumage color dimorphism and unique history in North America, including a recent population expansion and an epizootic of Mycoplasma gallisepticum (MG), the house finch (Haemorhous mexicanus) is a model species for studying sexual selection, plumage coloration and host-parasite interactions. As part of our ongoing efforts to make available genomic resources for this species, here we report a transcriptome assembly derived from genes expressed in spleen. We characterize transcriptomes from two populations with different histories of demography and disease exposure: a recently founded population in the eastern US that has been exposed to MG for over a decade and a native population from the western range that has never been exposed to MG. We utilize this resource to quantify conservation in gene expression in passerine birds over approximately 50 MY by comparing splenic expression profiles for 9,646 house finch transcripts and those from zebra finch and find that less than half of all genes expressed in spleen in either species are expressed in both species. Comparative gene annotations from several vertebrate species suggest that the house finch transcriptomes contain ~15 genes not yet found in previously sequenced vertebrate genomes. The house finch transcriptomes harbour ~85,000 SNPs, ~20,000 of which are non-synonymous. Although not yet validated by biological or technical replication, we identify a set of genes exhibiting differences between populations in gene expression (n = 182; 2% of all transcripts), allele frequencies (76 F ST ouliers) and alternative splicing as well as genes with several fixed non-synonymous substitutions; this set includes genes with functions related to double-strand break repair and immune response. The two house finch spleen transcriptome profiles will add to the increasing data on genome and transcriptome sequence information from natural populations. Differences in splenic expression between house finch and zebra finch imply either significant evolutionary turnover of splenic expression patterns or different physiological states of the individuals examined. The transcriptome resource will enhance the potential to annotate an eventual house finch genome, and the set of gene-based high-quality SNPs will help clarify the genetic underpinnings of host-pathogen interactions and sexual selection.

Identifying genes influenced by natural selection can provide information about lineage-specific adaptations, and transcriptomes generated by next-generation sequencing are a useful resource for identifying such genes. Here, we utilize a spleen transcriptome for the house finch (Haemorhous mexicanus), an emerging model for sexual selection and disease ecology, together with previously sequenced avian genomes (chicken, turkey, and zebra finch), to investigate lineage-specific adaptations within birds. An analysis of 4,398 orthologous genes revealed a significantly higher ratio of nonsynonymous to synonymous substitutions and significantly higher GC content in passerines than in galliforms, an observation deviating from strictly neutral expectations but consistent with an effect of biased gene conversion on the evolutionary rate in passerines. These data also showed that genes exhibiting signs of positive selection and fast evolution in passerines have functional roles related to fat metabolism, neurodevelopment, and ion binding.

Background Inference of population and species histories and population stratification using genetic data is important for discriminating between different speciation scenarios and for correct interpretation of genome scans for signs of adaptive evolution and trait association. Here we use data from 24 intronic loci re-sequenced in population samples of two closely related species, the pied flycatcher and the collared flycatcher. Results We applied Isolation-Migration models, assignment analyses and estimated the genetic differentiation and diversity between species and between populations within species. The data indicate a divergence time between the species of

To test the validity of subspecies status within the common pheasant (Phasianus colchicus), we used genetic data to assess differences among five subspecies distributed across the Loess Plateau in central China. Samples were collected from 191 individuals representing subspecies P. c. strauchi, P. c. satscheuensis, P. c. kiangsuensis, P. c. suehschanensis, and P. c. edzinensis. Analysis of mitochondrial genetic variation suggested that there were no phylogenetically distinct subspecies clades, and we could not exclude that hybridization has occurred among subspecies. A potential explanation of the observed pattern is that an increase in common pheasant rearing farms may have led to enforcement of hybridization between different subspecies. Consistent with phenotypic data, our analysis suggests that a previously identified common pheasant subspecies (P. c. strauchi) represents an intermediate form between the subspecies P. c. satscheuensis and P. c. suehschanensis.

Background Recent theoretical and empirical work points toward a significant role for sex-chromosome linked genes in the evolution of traits that induce reproductive isolation and for traits that evolve under influence of sexual selection. Empirical studies including recently diverged (Pleistocene), short-lived avian species pairs with short generation times have found that introgression occurs on the autosomes but not on the Z-chromosome. Here we study genetic differentiation and gene flow in the long-lived greater spotted eagle (Aquila clanga) and lesser spotted eagle (A. pomarina), two species with comparatively long generation times. Results Our data suggest that there is a directional bias in migration rates between hybridizing spotted eagles in eastern Europe. We find that a model including post divergence gene flow fits our data best for both autosomal and Z-chromosome linked loci but, for the Z-chromosome, the rate is reduced in the direction from A. pomarina to A. clanga. Conclusions The fact that some introgression still occurs on the Z-chromosome between these species suggests that the differentiation process is in a more premature phase in our study system than in previously studied avian species pairs and that could be explained by a shorter divergence time and/or a longer average generation time in the spotted eagles. The results are in agreement with field observations and provide further insight into the role of sex-linked loci for the build-up of barriers to gene flow among diverging populations and species.

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post-zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z-chromosome. Could rearrangements of the Z-chromosome between the species explain their reproductive isolation? We developed high coverage Z-chromosome linkage maps for both species, using gene-based markers and large-scale SNP genotyping. Best order maps contained 57-62 gene markers with an estimated average density of one every 1-1.5 Mb. We estimated the recombination rates in flycatcher Z-chromosomes to 1.1-1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co-linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post-zygotic reproductive isolation between these two Ficedula species.

The view of species as entities subjected to natural selection and amenable to change put forth by Charles Darwin and Alfred Wallace laid the conceptual foundation for understanding speciation. Initially marred by a rudimental understanding of hereditary principles, evolutionists gained appreciation of the mechanistic underpinnings of speciation following the merger of Mendelian genetic principles with Darwinian evolution. Only recently have we entered an era where deciphering the molecular basis of speciation is within reach. Much focus has been devoted to the genetic basis of intrinsic postzygotic isolation in model organisms and several hybrid incompatibility genes have been successfully identified. However, concomitant with the recent technological advancements in genome analysis and a newfound interest in the role of ecology in the differentiation process, speciation genetic research is becoming increasingly open to non-model organisms. This development will expand speciation research beyond the traditional boundaries and unveil the genetic basis of speciation from manifold perspectives and at various stages of the splitting process. This review aims at providing an extensive overview of speciation genetics. Starting from key historical developments and core concepts of speciation genetics, we focus much of our attention on evolving approaches and introduce promising methodological approaches for future research venues.

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 over-representation 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 omega 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.

Understanding the causes and consequences of variation in the rate of recombination is essential since this parameter is considered to affect levels of genetic variability, the efficacy of selection, and the design of association and linkage mapping studies. However, there is limited knowledge about the factors governing recombination rate variation. We genotyped 1920 single nucleotide polymorphisms in a multigeneration pedigree of more than 1000 zebra finches (Taeniopygia guttata) to develop a genetic linkage map, and then we used these map data together with the recently available draft genome sequence of the zebra finch to estimate recombination rates in 1 Mb intervals across the genome. The average zebra finch recombination rate (1.5 cM/Mb) is higher than in humans, but significantly lower than in chicken. The local rates of recombination in chicken and zebra finch were only weakly correlated, demonstrating evolutionary turnover of the recombination landscape in birds. The distribution of recombination events was heavily biased toward ends of chromosomes, with a stronger telomere effect than so far seen in any organism. In fact, the recombination rate was as low as 0.1 cM/Mb in intervals up to 100 Mb long in the middle of the larger chromosomes. We found a positive correlation between recombination rate and GC content, as well as GC-rich sequence motifs. Levels of linkage disequilibrium (LD) were significantly higher in regions of low recombination, showing that heterogeneity in recombination rates have left a footprint on the genomic landscape of LD in zebra finch populations.

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.

V. Dinca and P.D.N. Hebert. Biodiversity Institute of Ontario, University of Guelph, 50 Stone Road East, ON N1G 2W1, Canada. N. Backstrom. Department of Evolutionary Biology, Uppsala University, Norbyvagen 18D, 75236 Uppsala, Sweden. L. Dapporto. Department of Biological and Medical Sciences, Oxford Brookes University, Headington, Oxford, OX3 0BP, UK. M. Friberg. Department of Ecology and Genetics, Plant Ecology and Evolution, Uppsala University, Norbyvagen 18D, 75236 Uppsala, Sweden. E. Garcia-Barros, J. Hernandez-Roldan, and M.L. Munguira. Department of Biology, Universidad Autonoma de Madrid, Campus Cantoblanco 28049, Madrid, Spain. E. Hornett. Department of Zoology, University of Cambridge, Cambridge, CB2 3EJ, UK. V. Lukhtanov. Department of Karyosystematics, Zoological Institute of Russian Academy of Sciences, Universitetskaya nab. 1, 199034 St. Petersburg, Russia. F. Marec and J. Sichova Faculty of Science, University of South Bohemia, Ceske Budějovice, 370 05, Czech Republic. S. Montagud. Institut Cavanilles de Biodiversitat i Biologia Evolutiva (ICBiBE) Universitat de Valencia, Carrer Catedratic Jose Beltran 2, 46980, Paterna, Spain. M. Olofsson and C. Wiklund. Department of Zoology, Stockholm University, Svante Arrhenius vag 18B, S-106 91, Stockholm, Sweden. G. Talavera. Department of Organismic and Evolutionary Biology, Museum of Comparative Zoology, Harvard University, 26 Oxford St., 02138 Cambridge, MA, USA. J.C. Vicente-Arranz. Asociacion Espanola para la Proteccion de las Mariposas y su Medio (Zerynthia), Logrono (La Rioja), Spain. R. Vila. Institut de Biologia Evolutiva (CSIC-Universitat Pompeu Fabra), Passeig Maritim de la Barceloneta 37, 08003, Barcelona, Spain. Corresponding author: Vlad Dinca (e-mail: vdinca@uoguelph.ca). 391

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (alpha(m)) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of alpha(m) and c. Alpha(m) is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of alpha(m) based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G15 = 37.6, p = 0.0010) as well as within W chromosome introns (G15 = 44.0, p = 0.0001). This heterogeneity also affected the estimates of alpha(m), which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F(13,14) = 2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of alpha(m).