Search

Heat waves are predicted to be detrimental for organismal physiology with costs for survival that could be reflected in markers of biological state such as telomeres. Changes in early life telomere dynamics driven by thermal stress are of particular interest during the early post-natal stages of altricial birds because nestlings quickly shift from being ectothermic to endothermic after hatching. Telomeres of ectothermic and endothermic organisms respond differently to environmental temperature, but few investigations within species that transition from ectothermy to endothermy are available. Also, ambient temperature influences parental brooding behaviour, which will alter the temperature experienced by offspring and thereby, potentially, their telomeres. We exposed zebra finch nestlings to experimental heat waves and compared their telomere dynamics to that of a control group at 5, 12 and 80 days of age that encapsulate the transition from the ectothermic to the endothermic thermoregulatory stage; we also recorded parental brooding, offspring sex, mass, growth rates, brood size and hatch order. Nestling mass showed an inverse relationship with telomere length, and nestlings exposed to heat waves showed lower telomere attrition during their first 12 days of life (ectothermic stage) compared to controls. Additionally, parents of heated broods reduced the time they spent brooding offspring (at 5 days old) compared to controls. Our results indicate that the effect of heat waves on telomere dynamics likely varies depending on age and thermoregulatory stage of the offspring in combination with parental brooding behaviour during growth., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Telomere length and DNA methylation (DNAm) are two promising biomarkers of biological age. Environmental factors and life history traits are known to affect variation in both these biomarkers, especially during early life, yet surprisingly little is known about their reciprocal association, especially in natural populations. Here, we explore how variation in DNAm, growth rate, and early-life conditions are associated with telomere length changes during development. We tested these associations by collecting data from wild, nestling zebra finches in the Australian desert. We found that increases in the level of DNAm were negatively correlated with telomere length changes across early life. We also confirm previously documented effects of post hatch growth rate and clutch size on telomere length in a natural ecological context for a species that has been extensively studied in the laboratory. However, we did not detect any effect of ambient temperature during developmental on telomere length dynamics. We also found that the absolute telomere length of wild zebra finches, measured using the in-gel TRF method, was similar to that of captive birds. Our findings highlight exciting new opportunities to link and disentangle potential relationships between DNA based biomarkers of ageing, and of physiological reactions to environmental change., (© 2021 John Wiley & Sons Ltd.)

Emerging pathogens can have devastating effects on naïve hosts, but disease outcomes often vary among host species. Comparing the cellular response of different hosts to infection can provide insight into mechanisms of host defence. Here, we used RNA-seq to characterize the transcriptomic response of Darwin's finches to avian poxvirus, a disease of concern in the Galápagos Islands. We tested whether gene expression differs between infected and uninfected birds, and whether transcriptomic differences were related either to known antiviral mechanisms and/or the co-option of the host cellular environment by the virus. We compared two species, the medium ground finch (Geospiza fortis) and the vegetarian finch (Platyspiza crassirostris), to determine whether endemic Galápagos species differ in their response to pox. We found that medium ground finches had a strong transcriptomic response to infection, upregulating genes involved in the innate immune response including interferon production, inflammation, and other immune signalling pathways. In contrast, vegetarian finches had a more limited response, and some changes in this species were consistent with viral manipulation of the host's cellular function and metabolism. Many of the transcriptomic changes mirrored responses documented in model and in vitro studies of poxviruses. Our results thus indicate that many pathways of host defence against poxviruses are conserved among vertebrates and present even in hosts without a long evolutionary history with the virus. At the same time, the differences we observed between closely related species suggests that some endemic species of Galápagos finch could be more susceptible to avian pox than others., (© 2022 John Wiley & Sons Ltd.)

In a rapidly warming world, exposure to high temperatures may impact fitness, but the gene regulatory mechanisms that link sublethal heat to sexually selected traits are not well understood, particularly in endothermic animals. Our experiment used zebra finches (Taeniopygia guttata), songbirds that experience extreme temperature fluctuations in their native Australia. We exposed captive males to an acute thermal challenge (43°C) compared with thermoneutral (35°C) and lower (27°C) temperatures. We found significantly more heat dissipation behaviours at 43°C, a temperature previously shown to reduce song production and fertility, and more heat retention behaviours at 27°C. Next, we characterized transcriptomic responses in tissues important for mating effort-the posterior telencephalon, for its role in song production, and the testis, for its role in fertility and hormone production. Differential expression of hundreds of genes in the testes, but few in the brain, suggests the brain is less responsive to extreme temperatures. Nevertheless, gene network analyses revealed that expression related to dopaminergic signalling in the brain covaried with heat dissipation behaviours, providing a mechanism by which temporary thermal challenges may alter motivational circuits for song production. In both brain and testis, we observed correlations between thermally sensitive gene networks and individual differences in thermoregulatory behaviour. Although we cannot directly relate these gene regulatory changes to mating success, our results suggest that individual variation in response to thermal challenges could impact sexually selected traits in a warming world., (© 2022 John Wiley & Sons Ltd.)

Genomic diversity is the evolutionary foundation for adaptation to environmental change and thus is essential to consider in conservation planning. Island species are ideal for investigating the evolutionary drivers of genomic diversity, in part because of the potential for biological replicates. Here, we use genome data from 180 individuals spread among 27 island populations from 17 avian species to study the effects of island area, body size, demographic history and conservation status on contemporary genomic diversity. Our study expands earlier work on a small number of neutral loci to the entire genome and from a few species to many. We find significant positive correlation between island size and genomic diversity, a significant negative correlation between body size and genomic diversity, and that historical population declines significantly reduced contemporary genomic diversity. Our study shows that island size is the key factor in determining genomic diversity, indicating that habitat conservation is key to maintaining adaptive potential in the face of global environmental change. We found that threatened species generally had a significantly smaller values of Watterson's theta (θ W  = 4N e μ) compared to nonthreatened species, suggesting that θ W may be useful as a conservation indicator for at-risk species. Overall, these findings (a) provide biological insights into how genomic diversity scales with ecological, morphological and demographic factors; and (b) illustrate how population genomic data can be leveraged to better inform conservation efforts., (© 2019 John Wiley & Sons Ltd.)

Competitive environments promote high testosterone levels, produce oxidative stress and, consequently, impair cellular homeostasis. The regulation of genes involved in the synthesis of the pigment pheomelanin in melanocytes seems to help to maintain homeostasis against environmental oxidative stress. Here, we experimentally increased social interactions in some zebra finch (Taeniopygia guttata) males by keeping them in groups of six birds during feather growth, while others were kept alone, to test if melanocytes show epigenetic lability under a competitive social environment. As these changes may depend on the oxidative status, we administrated buthionine sulfoximine (BSO) to decrease the antioxidant capacity of some birds. The competitive environment downregulated a gene involved in pheomelanin synthesis (Slc7a11) by changing the level of DNA methylation in feather melanocytes. In other genes involved in pheomelanin synthesis (Slc45a2, MC1R and AGRP), DNA methylation was also affected, but no changes in expression were detected. Exposure to the competitive environment did not affect systemic oxidative stress and damage, indicating that a protective epigenetic mechanism that changes the expression of Slc7a11 may have been activated. However, no changes to the pigmentation phenotype of birds were found, probably due to the short duration or low intensity of the competitive environment. BSO treatment did not affect the epigenetic mechanism, suggesting that the antioxidant capacity of birds was high enough to deal with the competitive environment. An epigenetic mechanism limiting pheomelanin synthesis therefore becomes activated under exposure to a competitive environment in male zebra finches, which may help to avoid damage caused by competitive interactions., (© 2019 John Wiley & Sons Ltd.)

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier-beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild., (© 2017 John Wiley & Sons Ltd.)

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., (© 2017 John Wiley & Sons Ltd.)

Pheomelanin is a sulphur-containing yellow-to-reddish pigment whose synthesis consumes the main intracellular antioxidant (glutathione; GSH) and its precursor cysteine. Cysteine used for pheomelanogenesis cannot be used for antioxidant protection. We tested whether the expression of Slc7a11, the gene regulating the transport of cysteine to melanocytes for pheomelanogenesis, is environmentally influenced when cysteine/GSH are most required for antioxidant protection. We found that zebra finches Taeniopygia guttata developing pheomelanin-pigmented feathers during a 12-day exposure to the pro-oxidant diquat dibromide downregulated the expression of Slc7a11 in feather melanocytes, but not the expression of other genes that affect pheomelanogenesis by mechanisms different from cysteine transport such as MC1R and Slc45a2. Accordingly, diquat-treated birds did not suffer increased oxidative stress. This indicates that some animals have evolved an adaptive epigenetic lability that avoids damage derived from pheomelanogenesis. This mechanism should be explored in human Slc7a11 to help combat some cancer types related to cysteine consumption., (© 2016 John Wiley & Sons Ltd.)

Adaptive radiation unfolds as selection acts on the genetic variation underlying functional traits. The nature of this variation can be revealed by studying the tips of an ongoing adaptive radiation. We studied genomic variation at the tips of the Darwin's finch radiation; specifically focusing on polymorphism within, and variation among, three sympatric species of the genus Geospiza. Using restriction site-associated DNA (RAD-seq), we characterized 32 569 single-nucleotide polymorphisms (SNPs), from which 11 outlier SNPs for beak and body size were uncovered by a genomewide association study (GWAS). Principal component analysis revealed that these 11 SNPs formed four statistically linked groups. Stepwise regression then revealed that the first PC score, which included 6 of the 11 top SNPs, explained over 80% of the variation in beak size, suggesting that selection on these traits influences multiple correlated loci. The two SNPs most strongly associated with beak size were near genes associated with beak morphology across deeper branches of the radiation: delta-like 1 homologue (DLK1) and high-mobility group AT-hook 2 (HMGA2). Our results suggest that (i) key adaptive traits are associated with a small fraction of the genome (11 of 32 569 SNPs), (ii) SNPs linked to the candidate genes are dispersed throughout the genome (on several chromosomes), and (iii) micro- and macro-evolutionary variation (roots and tips of the radiation) involve some shared and some unique genomic regions., (© 2016 John Wiley & Sons Ltd.)

The two parental alleles at a specific locus are usually inherited with equal probability to the offspring. However, at least three processes can lead to an apparent departure from fair segregation: early viability selection, biased gene conversion and various kinds of segregation distortion. Here, we conduct a genome-wide scan for transmission distortion in a captive population of zebra finches (Taeniopygia guttata) using 1302 single-nucleotide polymorphisms (SNPs) followed by confirmatory analyses on independent samples from the same population. In the initial genome-wide scan, we found significant distortion at three linked loci on chromosome Tgu2 and we were able to replicate this finding in each of two follow-up data sets [overall transmission ratio = 0.567 (95% CI = 0.536-0.600), based on 1101 informative meioses]. Although the driving allele was preferentially transmitted by both heterozygous females [ratio = 0.560 (95% CI = 0.519-0.603)] and heterozygous males [ratio = 0.575 (95% CI = 0.531-0.623)], we could rule out postzygotic viability selection and biased gene conversion as possible mechanisms. Early postzygotic viability selection is unlikely, because it would result in eggs with no visible embryo and hence no opportunity for genotyping, and we confirmed that both females and males heterozygous for the driving allele did not produce a larger proportion of such eggs than homozygous birds. Biased gene conversion is expected to be rather localized, while we could trace transmission distortion in haplotypes of several megabases in a recombination desert. Thus, we here report the rare case of a prezygotically active transmission distorter operating equally effectively in female and male meioses., (© 2015 John Wiley & Sons Ltd.)

Many endangered species suffer from the loss of genetic diversity, but some populations may be able to thrive even if genetically depleted. To investigate the underlying genetic processes of population bottlenecks, we apply an innovative approach for assessing genetic diversity in the last known population of the endangered Pale-headed Brushfinch (Atlapetes pallidiceps) in Ecuador. First, we measure genetic diversity at eleven neutral microsatellite loci and adaptive SNP variation in five Toll-like receptor (TLR) immune system genes. Bottleneck tests confirm genetic drift as the main force shaping genetic diversity in this species and indicate a 99 % reduction in population size dating back several hundred years. Second, we compare contemporary microsatellite diversity with historic museum samples of A. pallidiceps, finding no change in genetic diversity. Third, we compare genetic diversity in the Pale-headed Brushfinch with two co-occurring-related brushfinch species (Atlapetes latinuchus, Buarremon torquatus), finding a reduction of up to 91% diversity in the immune system genes but not in microsatellites. High TLR diversity is linked to decreased survival probabilities in A. pallidiceps. Low TLR diversity is thus probably an adaptation to the specific selection regime within its currently very restricted distribution (approximately 200 ha), but could severely restrict the adaptive potential of the species in the long run. Our study illustrates the importance of investigating both neutral and adaptive markers to assess the effect of population bottlenecks and for recommending specific management plans in endangered species., (© 2014 John Wiley & Sons Ltd.)

Transcriptome Shotgun Sequencing (RNA-seq) has been readily embraced by geneticists and molecular ecologists alike. As with all high-throughput technologies, it is critical to understand which analytic strategies are best suited and which parameters may bias the interpretation of the data. Here we use a comprehensive simulation approach to explore how various features of the transcriptome (complexity, degree of polymorphism π, alternative splicing), technological processing (sequencing error ε, library normalization) and bioinformatic workflow (de novo vs. mapping assembly, reference genome quality) impact transcriptome quality and inference of differential gene expression (DE). We find that transcriptome assembly and gene expression profiling (EdgeR vs. BaySeq software) works well even in the absence of a reference genome and is robust across a broad range of parameters. We advise against library normalization and in most situations advocate mapping assemblies to an annotated genome of a divergent sister clade, which generally outperformed de novo assembly (Trans-Abyss, Trinity, Soapdenovo-Trans). Transcriptome complexity (size, paralogs, alternative splicing isoforms) negatively affected the assembly and DE profiling, whereas the effects of sequencing error and polymorphism were almost negligible. Finally, we highlight the challenge of gene name assignment for de novo assemblies, the importance of mapping strategies and raise awareness of challenges associated with the quality of reference genomes. Overall, our results have significant practical and methodological implications and can provide guidance in the design and analysis of RNA-seq experiments, particularly for organisms where genomic background information is lacking., (© 2012 Blackwell Publishing Ltd.)

Protective immunity is expected to evolve when the costs of mounting an immune response are less than those of harbouring pathogens. Estimating the costs of immunity vs. pathogenesis in natural systems is challenging, however, because they are typically closely linked. Here we attempt to disentangle the relative cost of each using experimental infections in a natural host-parasite system in which hosts (house finches, Carpodacus mexicanus) differ in resistance to a bacterium (Mycoplasma gallisepticum, MG), depending on whether they originate from co-evolved or unexposed populations. Experimental infection with a 2007-strain of MG caused finches from co-evolved populations to lose significantly more mass relative to controls, than those from unexposed populations. In addition, infected co-evolved finches that lost the most mass harboured the least amounts of MG, whereas the reverse was true in finches from unexposed populations. Finally, within co-evolved populations, individuals that displayed transcriptional evidence of higher protective immune activity, as indicated by changes in the expression of candidate immune and immune-related genes in a direction consistent with increased resistance to MG, showed greater mass loss and lower MG load. Thus, mass loss appeared to reflect the costs of immunity vs. pathogenesis in co-evolved and unexposed populations, respectively. Our results suggest that resistance can evolve even when the short-term energetic costs of protective immunity exceed those of pathogenesis, providing the longer-term fitness costs of infection are sufficiently high., (© 2012 Blackwell Publishing Ltd.)

The intra- and interspecific diversity of avian beak morphologies is one of the most compelling examples for the power of natural selection acting on a morphological trait. The development and diversification of the beak have also become a textbook example for evolutionary developmental biology, and variation in expression levels of several genes is known to causally affect beak shape. However, until now, no genomic polymorphisms have been identified, which are related to beak morphology in birds. QTL mapping does reveal the location of causal polymorphisms, albeit with poor spatial resolution. Here, we estimate heritability and genetic correlations for beak length, depth and width and perform a QTL linkage analysis for these traits based on 1404 informative single-nucleotide polymorphisms genotyped in a four-generation pedigree of 992 captive zebra finches (Taeniopygia guttata). Beak size, relative to body size, was sexually dimorphic (larger in males). Heritability estimates ranged from 0.47 for beak length to 0.74 for beak width. QTL mapping revealed four to five regions of significant or suggestive genome-wide linkage for each of the three beak dimensions (nine different regions in total). Eight out of 11 genes known to influence beak morphology are located in these nine peak regions. Five QTL do not cover known candidates demonstrating that yet unknown genes or regulatory elements may influence beak morphology in the zebra finch., (© 2012 Blackwell Publishing Ltd.)

Numerous studies have reported associations between heterozygosity in microsatellite markers and fitness-related traits (heterozygosity-fitness correlations, HFCs). However, it has often been questioned whether HFCs reflect general inbreeding depression, because a small panel of microsatellite markers does not reflect very well an individual's inbreeding coefficient (F) as calculated from a pedigree. Here, we challenge this prevailing view. Because of chance events during Mendelian segregation, an individual's realized proportion of the genome that is identical by descent (IBD) may substantially deviate from the pedigree-based expectation (i.e. F). This Mendelian noise may result in a weak correlation between F and multi-locus heterozygosity, but this does not imply that multi-locus heterozygosity is a bad estimator of realized IBD. We examined correlations between 11 fitness-related traits measured in up to 1192 captive zebra finches and three measures of inbreeding: (i) heterozygosity across 11 microsatellite markers, (ii) heterozygosity across 1359 single-nucleotide polymorphism (SNP) markers and (iii) F, based on a 5th-generation pedigree. All 11 phenotypic traits showed positive relationships with measures of heterozygosity, especially traits that are most closely related to fitness. Remarkably, the small panel of microsatellite markers produced equally strong HFCs as the large panel of SNP markers. Both marker-based approaches produced stronger correlations with phenotypes than the pedigree-based F, and this did not seem to result from the shortness of our pedigree. We argue that a small panel of microsatellites with high allelic richness may better reflect an individual's realized IBD than previously appreciated, especially in species like the zebra finch, where much of the genome is inherited in large blocks that rarely experience cross-over during meiosis., (© 2012 Blackwell Publishing Ltd.)

Innate immunity is expected to play a primary role in conferring resistance to novel infectious diseases, but few studies have attempted to examine its role in the evolution of resistance to emerging pathogens in wild vertebrate populations. Here, we used experimental infections and cDNA microarrays to examine whether changes in the innate and/or acquired immune responses likely accompanied the emergence of resistance in house finches (Carpodacus mexicanus) in the eastern United States subject to a recent outbreak of conjunctivitis-causing bacterium (Mycoplasma gallisepticum-MG). Three days following experimental infection with MG, we observed differences in the splenic transcriptional responses between house finches from eastern U.S. populations, with a 12-year history of MG exposure, versus western U.S. populations, with no history of exposure to MG. In particular, western birds down-regulated gene expression, while eastern finches showed no expression change relative to controls. Studies involving poultry have shown that MG can manipulate host immunity, and our observations suggest that pathogen manipulation occurred only in finches from the western populations, outside the range of MG. Fourteen days after infection, eastern finches, but not western finches, up-regulated genes associated with acquired immunity (cell-mediated immunity) relative to controls. These observations suggest population differences in the temporal course of the response to infection with MG and imply that innate immune processes were targets of selection in response to MG in the eastern U.S. population., (© 2012 Blackwell Publishing Ltd.)

The divergence-with-gene-flow model of speciation has a strong theoretical basis with a growing number of plausible examples in nature, but remains hotly debated. Darwin's finches of the Galápagos Archipelago have played an important role in our understanding of speciation processes. Recent studies suggest that this group may also provide insights into speciation via divergence with gene flow. On the island of Santa Cruz, recent studies found evidence for adaptive divergence in Darwin's small ground finch, Geospiza fuliginosa, between ecologically contrasting arid and humid zones. Despite the short geographical distance between these zones, strong disruptive selection during low rainfall periods is expected to generate and maintain adaptive divergence. Conversely, during high rainfall periods, when disruptive selection is predicted to be weakened, population divergence in adaptive traits is expected to break down. Because periods of low and high rainfall irregularly alternate, the geographical pattern of adaptive divergence can be assumed to break down and, importantly, regenerate in situ. Here, we use microsatellite allele frequency data to assess the genetic population structure of G. fuliginosa on Santa Cruz. We sample 21 sites and four ecological zones across the island. We reject hypotheses of population substructure linked to ecological and geographical differences among sites in favour of a single panmictic population. Panmixia implies high levels of gene flow within Santa Cruz, which favours selection over genetic drift as a valid process generating phenotypic divergence in G. fuliginosa on Santa Cruz. We discuss how our findings may support classic adaptation, phenotypic plasticity, matching habitat choice or any combination of these three processes., (© 2012 Blackwell Publishing Ltd.)

Avian wing length is an important trait that covaries with the ecology and migratory behaviour of a species and tends to change rapidly when the conditions are altered. Long-distance migrants typically have longer wings than short-distance migrants and sedentary species, and long-winged species also tend to be more dispersive. Although the substantial heritability of avian wing length is well established, the identification of causal genes has remained elusive. Based on large-scale genotyping of 1404 informative single nucleotide polymorphisms (SNP) in a captive population of 1067 zebra finches, we here show that the within-population variation of relative wing length (h(2) = 0.74 ± 0.05) is associated with standing genetic variation in at least six genomic regions (one genome-wide significant and five suggestive). The variance explained by these six quantitative trait loci (QTL) sums to 36.8% of the phenotypic variance (half of the additive genetic variance), although this likely is an overestimate attributable to the Beavis effect. As avian wing length is primarily determined by the length of the primary feathers, we then searched for candidate genes that are related to feather growth. Interestingly, all of the QTL signals co-locate with Wnt growth factors and closely interacting genes (Wnt3a, Wnt5a, Wnt6, Wnt7a, Wnt9a, RhoU and RhoV). Our findings therefore suggest that standing genetic variation in the Wnt genes might be linked to avian wing morphology, although there are many other genes that also fall within the confidence regions., (© 2011 Blackwell Publishing Ltd.)

The relationship between genetic diversity and fitness has important implications in evolutionary and conservation biology. This relationship has been widely investigated at the individual level in studies of heterozygosity-fitness correlations (HFC). General effects caused by inbreeding and/or local effects at single loci have been used as explanations of HFC, but the debate about the causes of HFC in open, natural populations is still ongoing. Study designs that control for variation in the inbreeding level of the individuals, and knowledge on the function and location of the markers used to measure heterozygosity, are fundamental to understand the causes of HFC. Here we investigated correlations between individual heterozygosity and estimates of survival at different life-history stages in an open population of blue tits (Cyanistes caeruleus). For survival at the embryo, nestling and fledgling stage, we used a full-sibling approach, i.e. we controlled for the level of inbreeding. We genotyped 1496 individuals with 79 microsatellites mapped across 25 chromosomes in the zebra finch (Taeniopygia guttata) that were classified either as potentially functional (58 loci) or as neutral (21 loci). We found different effects of standardized multilocus heterozygosity (SH): SH(functional) had a negative effect on the probability of hatching and local recruitment of females, whereas SH(neutral) had a positive effect on adult survival. The negative effects of functional loci are better explained by local effects, whereas the positive effects of neutral markers could reflect inbreeding effects in the population. Our results highlight the importance of considering the characteristics of the markers used in HFC studies and confirm the mixed effects of heterozygosity in different contexts (e.g. sex and life-history stage)., (© 2011 Blackwell Publishing Ltd.)

In recent years there has been a dramatic increase in the availability of high density genetic marker data for both model and non-model organisms. A potential application of these data is to infer relatedness in the absence of a complete pedigree. Using a marker panel of 771 SNPs genotyped in three generations of an extensive zebra finch pedigree, correlations between pedigree relatedness and seven marker-based estimates of relatedness were examined, as was the relationship between heterozygosity and inbreeding. Although marker-based and pedigree relatedness were highly correlated, the variance in estimated relatedness was high. Further, the correlation between heterozygosity and inbreeding was weak, even though mean inbreeding coefficient is typical of that seen in wild vertebrate pedigrees; the weak relationship was in part due to the small variance in inbreeding in the pedigree. Our data suggest that using marker information to reconstruct the pedigree, and then calculating relatedness from the pedigree, is likely to give more accurate relatedness estimates than using marker-based estimators directly.

Since the days of allozyme analysis, we have been enamored with the idea that if we just had enough polymorphic mendelian loci, we could gauge the inbreeding level of individuals by measuring heterozygosity and simultaneously measure the degree of genetic relatedness between pairs of individuals. Given Mendel's Laws, we have always known that we would need numerous independently segregating loci to achieve any reasonable degree of accuracy. Santure et al. (2010, this issue) use a 771 marker SNP panel to assess heterozygosity levels and to assess pairwise relatedness, and compare both with theoretical expectations obtained from a carefully recorded pedigree of a zebra finch breeding colony, as a function of increasing numbers of SNP markers. They also compare the SNP results with those from a 20-locus microsatellite panel, showing that adding SNPs to a fairly large microsatellite panel improves accuracy, but given an existing panel of 125 SNPs, little is to be gained by adding microsatellites. They show that the accuracy available for estimating individual levels of inbreeding is somewhat limited. They also show that the average pairwise relatedness measures bracket pedigree relationship very nicely, but the variances for individual pairs remain substantial, even with a very large panel.

Identifying genes that underlie ecological traits will open exiting possibilities to study gene-environment interactions in shaping phenotypes and in measuring natural selection on genes. Evolutionary ecology has been pursuing these objectives for decades, but they come into reach now that next generation sequencing technologies have dramatically lowered the costs to obtain the genomic sequence information that is currently lacking for most ecologically important species. Here we describe how we generated over 2 billion basepairs of novel sequence information for an ecological model species, the great tit Parus major. We used over 16 million short sequence reads for the de novo assembly of a reference sequence consisting of 550 000 contigs, covering 2.5% of the genome of the great tit. This reference sequence was used as the scaffold for mapping of the sequence reads, which allowed for the detection of over 20 000 novel single nucleotide polymorphisms. Contigs harbouring 4272 of the single nucleotide polymorphisms could be mapped to a unique location on the recently sequenced zebra finch genome. Of all the great tit contigs, significantly more were mapped to the microchromosomes than to the intermediate and the macrochromosomes of the zebra finch, indicating a higher overall level of sequence conservation on the microchromosomes than on the other types of chromosomes. The large number of great tit contigs that can be aligned to the zebra finch genome shows that this genome provides a valuable framework for large scale genetics, e.g. QTL mapping or whole genome association studies, in passerines.

Sex-chromosomes are thought to play an important role in speciation, but few studies of non-model organisms have investigated the relative influence of multiple sex-linked markers on reproductive isolation. We collected 222 individuals along a geographical transect spanning the hybrid zone between Passerina amoena and P. cyanea (Aves: Cardinalidae). Using maximum-likelihood cline fitting methods, we estimated locus-specific introgression rates for 10 z-linked markers. Although the cline width estimates ranged from 2.8 to 584 km, eight of 10 loci had cline widths between 224 and 271 km. We also used coalescent-based estimates of locus-specific divergence times between P. amoena and P. cyanea to test a recently proposed hypothesis of an inverse relationship between divergence time and cline width but did not find a significant association. The narrow width (2.8 km) of the cline estimated from the VLDLR9 locus indicates strong selection retarding introgression of alleles at this locus across the hybrid zone. Interestingly, a mutation in the very low density lipoprotein receptor (VLDLR) gene, in which VLDLR9 is an intron, is known to reduce the egg-laying ability of some chickens, suggesting a possible link between this gene region and reproductive isolation between P. amoena and P. cyanea. These results underscore the importance of sampling multiple loci to investigate introgression patterns across a chromosome or genome and support previous findings of the importance of sex-linked genes in speciation.

The zebra finch (Taeniopygia guttata) is a small Australian grassland songbird that has been domesticated over the past two centuries. Because it is easy to breed in captivity, it has become a widely used study organism, especially in behavioural research. Most work has been conducted on domesticated populations maintained at numerous laboratories in Europe and North America. However, little is known about the extent to which, during the process of domestication, captive populations have gone through bottlenecks in population size, leading to inbred and potentially genetically differentiated study populations. This is an important issue, because (i) behavioural studies on captive populations might suffer from artefacts arising from high levels of inbreeding or lack of genetic variation in such populations, and (ii) it may hamper the comparability of research findings. To address this issue, we genotyped 1000 zebra finches from 18 captive and two wild populations at 10 highly variable microsatellite loci. We found that all captive populations have lost some of the genetic variability present in the wild, but there is no evidence that they have gone through a severe bottleneck, as the average captive population still showed a mean of 11.7 alleles per locus, compared to a mean of 19.3 alleles/locus for wild zebra finches. We found significant differentiation between the captive populations (F(ST) = 0.062). Patterns of genetic similarity closely match geographical relationships, so the most pronounced differences occur between the three continents: Australia, North America, and Europe. By providing a tree of the genetic similarity of the different captive populations, we hope to contribute to a better understanding of variation in research findings obtained by different laboratories.

In 1994, the bacterial parasite Mycoplasma gallisepticum expanded its host range and swept through populations of a novel host--eastern US populations of the house finch (Carpodacus mexicanus). This epizootic caused a dramatic decline in finch population numbers, has been shown to have caused strong selection on house finch morphology, and presumably caused evolutionary change at the molecular level as finches evolved enhanced resistance. As a first step toward identifying finch genes that respond to infection by Mycoplasma and which may have experienced natural selection by this parasite, we used suppression subtractive hybridization (SSH) and cDNA macroarray approaches to identify differentially expressed genes regulated by the Mycoplasma parasite. Two subtractive cDNA libraries consisting of 16,512 clones were developed from spleen using an experimentally uninfected bird as the 'tester' and an infected bird as 'driver', and vice versa. Two hundred and twenty cDNA clones corresponding 34 genes with known vertebrate homologues and a large number of novel transcripts were found to be qualitatively up- or down-regulated genes by high-density filter hybridization. These gene expression changes were further confirmed by a high throughout reverse Northern blot approach and in specific cases by targeted Northern analysis. blast searches show that heat shock protein (HSP) 90, MHC II-associated invariant chain (CD74), T-cell immunoglobulin mucin 1 (TIM1), as well as numerous novel expressed genes not found in the databases were up- or down-regulated by the host in response to this parasite. Our results and macroarray resources provide a foundation for molecular co-evolutionary studies of the Mycoplasma parasite and its recently colonized avian host.

The impact of founder events on levels of genetic variation in natural populations remains a topic of significant interest. Well-documented introductions provide a valuable opportunity to examine how founder events influence genetic diversity in invasive species. House finches (Carpodacus mexicanus) are passerine birds native to western North America, with the large eastern North American population derived from a small number of captive individuals released in the 1940s. Previous comparisons using amplified fragment length polymorphism (AFLP) markers found equivalent levels of diversity in eastern and western populations, suggesting that any genetic effects of the founder event were ameliorated by the rapid growth of the newly established population. We used an alternative marker system, 10 highly polymorphic microsatellites, to compare levels of genetic diversity between four native and five introduced house finch populations. In contrast to the AFLP comparisons, we found significantly lower allelic richness and heterozygosity in introduced populations across all loci. Three out of five introduced populations showed significant reductions in the ratio of the number of alleles to the allele size range, a within-population characteristic of recent bottlenecks. Finally, native and introduced populations showed significant pairwise differences in allele frequencies in every case, with stronger isolation by distance within the introduced than native range. Overall, our results provide compelling molecular evidence for a founder effect during the introduction of eastern house finches that reduced diversity levels at polymorphic microsatellite loci and may have contributed to the emergence of the Mycoplasma epidemic which recently swept the eastern range of this species.

We use genetic divergence at 16 microsatellite loci to investigate how geographical features of the Galápagos landscape structure island populations of Darwin's finches. We compare the three most genetically divergent groups of Darwin's finches comprising morphologically and ecologically similar allopatric populations: the cactus finches (Geospiza scandens and Geospiza conirostris), the sharp-beaked ground finches (Geospiza difficilis) and the warbler finches (Certhidea olivacea and Certhidea fusca). Evidence of reduced genetic diversity due to drift was limited to warbler finches on small, peripheral islands. Evidence of low levels of recent interisland migration was widespread throughout all three groups. The hypothesis of distance-limited dispersal received the strongest support in cactus and sharp-beaked ground finches as evidenced by patterns of isolation by distance, while warbler finches showed a weaker relationship. Support for the hypothesis that gene flow constrains morphological divergence was only found in one of eight comparisons within these groups. Among warbler finches, genetic divergence was relatively high while phenotypic divergence was low, implicating stabilizing selection rather than constraint due to gene flow. We conclude that the adaptive radiation of Darwin's finches has occurred in the presence of ongoing but low levels of gene flow caused by distance-dependent interisland dispersal. Gene flow does not constrain phenotypic divergence, but may augment genetic variation and facilitate evolution due to natural selection. Both microsatellites and mtDNA agree in that subsets of peripheral populations of two older groups are genetically more similar to other species that underwent dramatic morphological change. The apparent decoupling of morphological and molecular evolution may be accounted for by a modification of Lack's two-stage model of speciation: relative ecological stasis in allopatry followed by secondary contact, ecological interactions and asymmetric phenotypic divergence.

Most phylogeographical studies of postglacial colonization focus on high latitude locations in the Northern Hemisphere. Here, we studied the phylogeographical structure of the red-necked snow finch Pyrgilauda ruficollis, an endemic species of the Tibetan plateau. We analysed 879 base pairs (bp) of the mitochondrial cytochrome b gene and 529 bp of the control region in 41 birds from four regional groups separated by mountain ranges. We detected 34 haplotypes, 31 of which occurred in a single individual and only three of which were shared among sampling sites within regional groups or among regional groups. Haplotype diversity was high (h = 0.94); nucleotide diversity was low (eth = 0.00415) and genetic differentiation was virtually non-existent. Analyses of mismatch distributions and geographically nested clades yielded results consistent with contiguous range expansion, and the expansion times were estimated as 0.07-0.19 million years ago (Ma). Our results suggest that P. ruficollis colonized the Tibetan plateau after the extensive glacial period (0.5-0.175 Ma), expanding from the eastern margin towards the inner plateau. Thus, in contrast to many of the post-glacial phylogeographical structures known at high latitudes, this colonization occurred without matrilineal population structuring. This might be due to the short glacial cycles typical of the Tibetan plateau, adaptation of P. ruficollis to cold conditions, or refugia and colonized habitat being semicontinuous and thus promoting population mixing.