The zebra finch is an important model organism in several fields 1,2 with unique relevance to hum... more The zebra finch is an important model organism in several fields 1,2 with unique relevance to human neuroscience . Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken 5 -the only bird with a sequenced genome until now 6 . Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes 7 . We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.
Background: Avian β-defensins (AvBDs) represent a group of innate immune genes with broad antimic... more Background: Avian β-defensins (AvBDs) represent a group of innate immune genes with broad antimicrobial activity. Within the chicken genome, previous work identified 14 AvBDs in a cluster on chromosome three. The release of a second bird genome, the zebra finch, allows us to study the comparative evolutionary history of these gene clusters between from two species that shared a common ancestor about 100 million years ago. Results: A phylogenetic analysis of the β-defensin gene clusters in the chicken and the zebra finch identified several cases of gene duplication and gene loss along their ancestral lines. In the zebra finch genome a cluster of 22 AvBD genes were identified, all located within 125 Kbp on chromosome three. Ten of the 22 genes were found to be highly conserved with orthologous genes in the chicken genome. The remaining 12 genes were all located within a cluster of 58 Kbp and are suggested to be a result of recent gene duplication events that occurred after the galliformes-passeriformes split (G-P split). Within the chicken genome, AvBD6 was found to be a duplication of AvBD7, whereas the gene AvDB14 seems to have been lost along the ancestral line of the zebra finch. The duplicated β-defensin genes have had a significantly higher accumulation of non-synonymous over synonymous substitutions compared to the genes that have not undergone duplication since the G-P split. The expression patterns of avian β-defensin genes seem to be well conserved between chicken and zebra finch.
The MHC (Major Histocompatibility Complex) molecules are responsible for presenting pathogens to ... more The MHC (Major Histocompatibility Complex) molecules are responsible for presenting pathogens to the immune system of the body. The aim of this study was to investigate the MHC of the lekking Great Snipe. Because of its mating system, this species could easily be used to monitor mate choice in the wild. Parts of the second exon of the class II gene was amplified with PCR using degenerate primers and sequenced after cloning of the individual PCR products. From five individuals seven different sequence types were found, all of them were more similar to each other than to the B-LBII sequence from chicken. The fact that three different sequences were found in one of the individuals indicates that more than one locus was amplified. No frameshift substitutions or amino acid shifts in traditionally conserved regions were found suggesting that the sequences could come from functional genes. The rates of synonymous substitutions were smaller than the rates non-synonymous substitutions when comparing the different sequence types to each other, this supports the view of balancing selection on the MHC.
Genome sequencing projects were long confined to biomedical model organisms and required the conc... more Genome sequencing projects were long confined to biomedical model organisms and required the concerted effort of large consortia. Rapid progress in high-throughput sequencing technology and the simultaneous development of bioinformatic tools have democratized the field. It is now within reach for individual research groups in the eco-evolutionary and conservation community to generate de novo draft genome sequences for any organism of choice. Because of the cost and considerable effort involved in such an endeavour, the important first step is to thoroughly consider whether a genome sequence is necessary for addressing the biological question at hand. Once this decision is taken, a genome project requires careful planning with respect to the organism involved and the intended quality of the genome draft. Here, we briefly review the state of the art within this field and provide a step-by-step introduction to the workflow involved in genome sequencing, assembly and annotation with pa...
Genome and transcriptome sequencing applications that rely on variation in sequence depth can be ... more Genome and transcriptome sequencing applications that rely on variation in sequence depth can be negatively affected if there are systematic biases in coverage. We have investigated patterns of local variation in sequencing coverage by utilising ultra-deep sequencing (>100,000X) of mtDNA obtained during sequencing of two vertebrate genomes, wolverine (Gulo gulo) and collared flycatcher (Ficedula albicollis). With such extreme depth, stochastic variation in coverage should be negligible, which allows us to provide a very detailed, fine-scale picture of sequence dependent coverage variation and sequencing error rates. Sequencing coverage showed up to six-fold variation across the complete mtDNA and this variation was highly repeatable in sequencing of multiple individuals of the same species. Moreover, coverage in orthologous regions was correlated between the two species and was negatively correlated with GC content. We also found a negative correlation between the site-specific s...
Parasites and diseases constitute major evolutionary forces in many natural populations, and thus... more Parasites and diseases constitute major evolutionary forces in many natural populations, and thus having an efficient immune defense to resist infections is crucial for many organisms. Properties of the immune response may also influence mate choice decisions in many animals. Theory predicts several advantages for females when choosing males with superior immune systems. These benefits can be both direct (e.g.
The house sparrow (Passer domesticus) is an important model species in ecology and evolution. How... more The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.
By next generation transcriptome sequencing, it is possible to obtain data on both nucleotide seq... more By next generation transcriptome sequencing, it is possible to obtain data on both nucleotide sequence variation and gene expression. We have used this approach (RNA-Seq) to investigate the genetic basis for differences in plumage coloration and mating strategies in a non-model bird species, the ruff (Philomachus pugnax). Ruff males show enormous variation in the coloration of ornamental feathers, used for individual recognition. This polymorphism is linked to reproductive strategies, with dark males (Independents) defending territories on leks against other Independents, whereas white morphs (Satellites) co-occupy Independent's courts without agonistic interactions. Previous work found a strong genetic component for mating strategy, but the genes involved were not identified. We present feather transcriptome data of more than 6,000 de-novo sequenced ruff genes (although with limited coverage for many of them). None of the identified genes showed significant expression divergence between males, but many genetic markers showed nucleotide differentiation between different color morphs and mating strategies. These include several feather keratin genes, splicing factors, and the Xg blood-group gene. Many of the genes with significant genetic structure between mating strategies have not yet been annotated and their functions remain to be elucidated. We also conducted in-depth investigations of 28 pre-identified coloration candidate genes. Two of these (EDNRB and TYR) were specifically expressed in black-and rust-colored males, respectively. We have demonstrated the utility of next generation transcriptome sequencing for identifying and genotyping large number of genetic markers in a non-model species without previous genomic resources, and highlight the potential of this approach for addressing the genetic basis of ecologically important variation. 2011). Strong genetic differences (both nucleotide divergence and variation in gene expression rates) were also found between normal and dwarf morphs of lake whitefish (Coregonus clupeaformis) .
The genomic organisation of the major histocompatibility complex (MHC) seems to vary considerably... more The genomic organisation of the major histocompatibility complex (MHC) seems to vary considerably between different bird species. In order to understand this variation it is important to gather information from different species. We have, for the first time, investigated MHC class II polymorphism in a wader species, the great snipe (Gallinago media). Eleven alleles were found in five sequenced individuals; these come from at least three different loci, but RFLP data suggest that a larger number of genes may be present. For MHC genes, amino acid substitutions followed the, for MHC genes, general pattern of high non-synonymous substitution rates in peptide-binding regions, suggesting that the sequenced alleles may be expressed. The number of genes, lengths of introns and exon sequences of the great snipe MHC seem to be intermediate between those of chicken and passerine birds.
The global loss of biodiversity continues at an alarming rate. Genomic approaches have been sugge... more The global loss of biodiversity continues at an alarming rate. Genomic approaches have been suggested as a promising tool for conservation practice as scaling up to genome-wide data can improve traditional conservation genetic inferences and provide qualitatively novel insights. However, the generation of genomic data and subsequent analyses and interpretations remain challenging and largely confined to academic research in ecology and evolution. This generates a gap between basic research and applicable solutions for conservation managers faced with multifaceted problems. Before the real-world conservation potential of genomic research can be realized, we suggest that current infrastructures need to be modified, methods must mature, analytical pipelines need to be developed, and successful case studies must be disseminated to practitioners.
The black grouse (Tetrao tetrix) is a galliform bird species that is important for both ecologica... more The black grouse (Tetrao tetrix) is a galliform bird species that is important for both ecological studies and conservation genetics. Here, we report the sequencing of the spleen transcriptome of black grouse using 454 GS FLX Titanium sequencing. We performed a large-scale gene discovery analysis with a focus on genes that might be related to fitness in this species and also identified a large set of microsatellites. In total, we obtained 182 179 quality-filtered sequencing reads that we assembled into 9035 contigs. Using these contigs and 15 794 length-filtered (greater than 200 bp) singletons, we identified 7762 transcripts that appear to be homologues of chicken genes. A specific BLAST search with an emphasis on immune genes found 308 homologous chicken genes that have immune function, including ten major histocompatibility complex-related genes located on chicken chromosome 16. We also identified 1300 expressed sequence tag microsatellites and were able to design suitable flanki...
Indirect mate choice is any behavior that restricts the individual&am... more Indirect mate choice is any behavior that restricts the individual's set of potential mates without discrimination of mate attributes directly, for example, by having preferences about where to mate. We analyzed a 14-year data set from great snipe (Gallinago media) leks for evidence of indirect mate choice based on relative and absolute position of lek territories. We found little or no effect of the centrality of territories on mating and no between-year consistency in the spatial distribution of matings within leks. Instead, the probability of matings occurring at a particular site increased if the current territory owner had mated the previous year. Furthermore, individual females returned in later seasons to mate with the same male as previously rather than at the same site. Previous work found that male interactions and dominance do not control matings and that females are very choosy about which territory they mate in. Here we show that this is because of the male occupying the territory rather than its position. We therefore conclude that direct female mate choice is the main behavioral process affecting variation in mating success among great snipe males, unlike in some lekking mammals where male competition and/or indirect mate choice appears more important.
This article is protected by copyright. All rights reserved. evolution. However, until recently g... more This article is protected by copyright. All rights reserved. evolution. However, until recently genomic resources for molecular ecological projects have been lacking in this species. Here we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximise the yield of immune-related gene functions. After de-novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8,756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3,177 single nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.
The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen rec... more The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and triggering of the adaptive immune response, and are therefore likely to be under selection from parasites. These selection regimes may vary in space and time. Here we report a strong geographical structure in MHC class II B genes of a migrating bird, the great snipe ( Gallinago media ). Genetic differentiation in the MHC between two ecologically distinct distributional regions (Scandinavian mountain populations vs. East European lowland populations) was still present after statistically controlling for the effect of selectively neutral variation (microsatellites) using partial Mantel tests. This suggests a role for selection in generating this spatial structure and that it represents local adaptation to different environments. Differentiation between populations within the two regions was negligible. Overall, we found a high number of MHC alleles (50, from 175 individuals). This, together with a tendency for a higher rate of nonsynonymous than synonymous substitutions in the peptide binding sites, and high Tajima's D in certain regions of the gene, suggests a history of balancing selection. MHC variation is often thought to be maintained by some form of balancing selection, but the nature of this selection remains unclear. Our results support the hypothesis that spatial variation in selection regimes contributes to the high polymorphism.
The zebra finch is an important model organism in several fields 1,2 with unique relevance to hum... more The zebra finch is an important model organism in several fields 1,2 with unique relevance to human neuroscience . Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken 5 -the only bird with a sequenced genome until now 6 . Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes 7 . We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.
Background: Avian β-defensins (AvBDs) represent a group of innate immune genes with broad antimic... more Background: Avian β-defensins (AvBDs) represent a group of innate immune genes with broad antimicrobial activity. Within the chicken genome, previous work identified 14 AvBDs in a cluster on chromosome three. The release of a second bird genome, the zebra finch, allows us to study the comparative evolutionary history of these gene clusters between from two species that shared a common ancestor about 100 million years ago. Results: A phylogenetic analysis of the β-defensin gene clusters in the chicken and the zebra finch identified several cases of gene duplication and gene loss along their ancestral lines. In the zebra finch genome a cluster of 22 AvBD genes were identified, all located within 125 Kbp on chromosome three. Ten of the 22 genes were found to be highly conserved with orthologous genes in the chicken genome. The remaining 12 genes were all located within a cluster of 58 Kbp and are suggested to be a result of recent gene duplication events that occurred after the galliformes-passeriformes split (G-P split). Within the chicken genome, AvBD6 was found to be a duplication of AvBD7, whereas the gene AvDB14 seems to have been lost along the ancestral line of the zebra finch. The duplicated β-defensin genes have had a significantly higher accumulation of non-synonymous over synonymous substitutions compared to the genes that have not undergone duplication since the G-P split. The expression patterns of avian β-defensin genes seem to be well conserved between chicken and zebra finch.
The MHC (Major Histocompatibility Complex) molecules are responsible for presenting pathogens to ... more The MHC (Major Histocompatibility Complex) molecules are responsible for presenting pathogens to the immune system of the body. The aim of this study was to investigate the MHC of the lekking Great Snipe. Because of its mating system, this species could easily be used to monitor mate choice in the wild. Parts of the second exon of the class II gene was amplified with PCR using degenerate primers and sequenced after cloning of the individual PCR products. From five individuals seven different sequence types were found, all of them were more similar to each other than to the B-LBII sequence from chicken. The fact that three different sequences were found in one of the individuals indicates that more than one locus was amplified. No frameshift substitutions or amino acid shifts in traditionally conserved regions were found suggesting that the sequences could come from functional genes. The rates of synonymous substitutions were smaller than the rates non-synonymous substitutions when comparing the different sequence types to each other, this supports the view of balancing selection on the MHC.
Genome sequencing projects were long confined to biomedical model organisms and required the conc... more Genome sequencing projects were long confined to biomedical model organisms and required the concerted effort of large consortia. Rapid progress in high-throughput sequencing technology and the simultaneous development of bioinformatic tools have democratized the field. It is now within reach for individual research groups in the eco-evolutionary and conservation community to generate de novo draft genome sequences for any organism of choice. Because of the cost and considerable effort involved in such an endeavour, the important first step is to thoroughly consider whether a genome sequence is necessary for addressing the biological question at hand. Once this decision is taken, a genome project requires careful planning with respect to the organism involved and the intended quality of the genome draft. Here, we briefly review the state of the art within this field and provide a step-by-step introduction to the workflow involved in genome sequencing, assembly and annotation with pa...
Genome and transcriptome sequencing applications that rely on variation in sequence depth can be ... more Genome and transcriptome sequencing applications that rely on variation in sequence depth can be negatively affected if there are systematic biases in coverage. We have investigated patterns of local variation in sequencing coverage by utilising ultra-deep sequencing (>100,000X) of mtDNA obtained during sequencing of two vertebrate genomes, wolverine (Gulo gulo) and collared flycatcher (Ficedula albicollis). With such extreme depth, stochastic variation in coverage should be negligible, which allows us to provide a very detailed, fine-scale picture of sequence dependent coverage variation and sequencing error rates. Sequencing coverage showed up to six-fold variation across the complete mtDNA and this variation was highly repeatable in sequencing of multiple individuals of the same species. Moreover, coverage in orthologous regions was correlated between the two species and was negatively correlated with GC content. We also found a negative correlation between the site-specific s...
Parasites and diseases constitute major evolutionary forces in many natural populations, and thus... more Parasites and diseases constitute major evolutionary forces in many natural populations, and thus having an efficient immune defense to resist infections is crucial for many organisms. Properties of the immune response may also influence mate choice decisions in many animals. Theory predicts several advantages for females when choosing males with superior immune systems. These benefits can be both direct (e.g.
The house sparrow (Passer domesticus) is an important model species in ecology and evolution. How... more The house sparrow (Passer domesticus) is an important model species in ecology and evolution. However, until recently, genomic resources for molecular ecological projects have been lacking in this species. Here, we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximize the yield of immune-related gene functions. After de novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3177 single-nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.
By next generation transcriptome sequencing, it is possible to obtain data on both nucleotide seq... more By next generation transcriptome sequencing, it is possible to obtain data on both nucleotide sequence variation and gene expression. We have used this approach (RNA-Seq) to investigate the genetic basis for differences in plumage coloration and mating strategies in a non-model bird species, the ruff (Philomachus pugnax). Ruff males show enormous variation in the coloration of ornamental feathers, used for individual recognition. This polymorphism is linked to reproductive strategies, with dark males (Independents) defending territories on leks against other Independents, whereas white morphs (Satellites) co-occupy Independent's courts without agonistic interactions. Previous work found a strong genetic component for mating strategy, but the genes involved were not identified. We present feather transcriptome data of more than 6,000 de-novo sequenced ruff genes (although with limited coverage for many of them). None of the identified genes showed significant expression divergence between males, but many genetic markers showed nucleotide differentiation between different color morphs and mating strategies. These include several feather keratin genes, splicing factors, and the Xg blood-group gene. Many of the genes with significant genetic structure between mating strategies have not yet been annotated and their functions remain to be elucidated. We also conducted in-depth investigations of 28 pre-identified coloration candidate genes. Two of these (EDNRB and TYR) were specifically expressed in black-and rust-colored males, respectively. We have demonstrated the utility of next generation transcriptome sequencing for identifying and genotyping large number of genetic markers in a non-model species without previous genomic resources, and highlight the potential of this approach for addressing the genetic basis of ecologically important variation. 2011). Strong genetic differences (both nucleotide divergence and variation in gene expression rates) were also found between normal and dwarf morphs of lake whitefish (Coregonus clupeaformis) .
The genomic organisation of the major histocompatibility complex (MHC) seems to vary considerably... more The genomic organisation of the major histocompatibility complex (MHC) seems to vary considerably between different bird species. In order to understand this variation it is important to gather information from different species. We have, for the first time, investigated MHC class II polymorphism in a wader species, the great snipe (Gallinago media). Eleven alleles were found in five sequenced individuals; these come from at least three different loci, but RFLP data suggest that a larger number of genes may be present. For MHC genes, amino acid substitutions followed the, for MHC genes, general pattern of high non-synonymous substitution rates in peptide-binding regions, suggesting that the sequenced alleles may be expressed. The number of genes, lengths of introns and exon sequences of the great snipe MHC seem to be intermediate between those of chicken and passerine birds.
The global loss of biodiversity continues at an alarming rate. Genomic approaches have been sugge... more The global loss of biodiversity continues at an alarming rate. Genomic approaches have been suggested as a promising tool for conservation practice as scaling up to genome-wide data can improve traditional conservation genetic inferences and provide qualitatively novel insights. However, the generation of genomic data and subsequent analyses and interpretations remain challenging and largely confined to academic research in ecology and evolution. This generates a gap between basic research and applicable solutions for conservation managers faced with multifaceted problems. Before the real-world conservation potential of genomic research can be realized, we suggest that current infrastructures need to be modified, methods must mature, analytical pipelines need to be developed, and successful case studies must be disseminated to practitioners.
The black grouse (Tetrao tetrix) is a galliform bird species that is important for both ecologica... more The black grouse (Tetrao tetrix) is a galliform bird species that is important for both ecological studies and conservation genetics. Here, we report the sequencing of the spleen transcriptome of black grouse using 454 GS FLX Titanium sequencing. We performed a large-scale gene discovery analysis with a focus on genes that might be related to fitness in this species and also identified a large set of microsatellites. In total, we obtained 182 179 quality-filtered sequencing reads that we assembled into 9035 contigs. Using these contigs and 15 794 length-filtered (greater than 200 bp) singletons, we identified 7762 transcripts that appear to be homologues of chicken genes. A specific BLAST search with an emphasis on immune genes found 308 homologous chicken genes that have immune function, including ten major histocompatibility complex-related genes located on chicken chromosome 16. We also identified 1300 expressed sequence tag microsatellites and were able to design suitable flanki...
Indirect mate choice is any behavior that restricts the individual&am... more Indirect mate choice is any behavior that restricts the individual's set of potential mates without discrimination of mate attributes directly, for example, by having preferences about where to mate. We analyzed a 14-year data set from great snipe (Gallinago media) leks for evidence of indirect mate choice based on relative and absolute position of lek territories. We found little or no effect of the centrality of territories on mating and no between-year consistency in the spatial distribution of matings within leks. Instead, the probability of matings occurring at a particular site increased if the current territory owner had mated the previous year. Furthermore, individual females returned in later seasons to mate with the same male as previously rather than at the same site. Previous work found that male interactions and dominance do not control matings and that females are very choosy about which territory they mate in. Here we show that this is because of the male occupying the territory rather than its position. We therefore conclude that direct female mate choice is the main behavioral process affecting variation in mating success among great snipe males, unlike in some lekking mammals where male competition and/or indirect mate choice appears more important.
This article is protected by copyright. All rights reserved. evolution. However, until recently g... more This article is protected by copyright. All rights reserved. evolution. However, until recently genomic resources for molecular ecological projects have been lacking in this species. Here we present transcriptome sequencing data (RNA-Seq) from three different house sparrow tissues (spleen, blood and bursa). These tissues were specifically chosen to obtain a diverse representation of expressed genes and to maximise the yield of immune-related gene functions. After de-novo assembly, 15,250 contigs were identified, representing sequence data from a total of 8,756 known avian genes (as inferred from the closely related zebra finch). The transcriptome assembly contain sequence data from nine manually annotated MHC genes, including an almost complete MHC class I coding sequence. There were 407, 303 and 68 genes overexpressed in spleen, blood and bursa, respectively. Gene ontology terms related to ribosomal function were associated with overexpression in spleen and oxygen transport functions with overexpression in blood. In addition to the transcript sequences, we provide 327 gene-linked microsatellites (SSRs) with sufficient flanking sequences for primer design, and 3,177 single nucleotide polymorphisms (SNPs) within genes, that can be used in follow-up molecular ecology studies of this ecological well-studied species.
The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen rec... more The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and triggering of the adaptive immune response, and are therefore likely to be under selection from parasites. These selection regimes may vary in space and time. Here we report a strong geographical structure in MHC class II B genes of a migrating bird, the great snipe ( Gallinago media ). Genetic differentiation in the MHC between two ecologically distinct distributional regions (Scandinavian mountain populations vs. East European lowland populations) was still present after statistically controlling for the effect of selectively neutral variation (microsatellites) using partial Mantel tests. This suggests a role for selection in generating this spatial structure and that it represents local adaptation to different environments. Differentiation between populations within the two regions was negligible. Overall, we found a high number of MHC alleles (50, from 175 individuals). This, together with a tendency for a higher rate of nonsynonymous than synonymous substitutions in the peptide binding sites, and high Tajima's D in certain regions of the gene, suggests a history of balancing selection. MHC variation is often thought to be maintained by some form of balancing selection, but the nature of this selection remains unclear. Our results support the hypothesis that spatial variation in selection regimes contributes to the high polymorphism.
Uploads
Papers by Robert Ekblom