Gene expression divergence and nucleotide differentiation between males of different color morphs and mating strategies in the ruff

Gene Expression Modification by an Autosomal Inversion Associated With Three Male Mating Morphs

Chromosomal inversions are structural rearrangements that frequently provide genomic substrate for phenotypic diversity. In the ruff Philomachus pugnax, three distinct male reproductive morphs (Independents, Satellites and Faeders) are genetically determined by a 4.5 Mb autosomal inversion. Here we test how this stable inversion polymorphism affects gene expression in males during the lekking season. Gene expression may be altered through disruptions at the breakpoints and the accumulation of mutations due to suppressed recombination. We used quantitative PCR to measure expression of 11 candidate inversion genes across three different tissues (liver, adrenal glands and gonads) and tested for allelic imbalance in four inversion genes across 12 males of all three morphs (8 Independents, 2 Satellites, 2 Faeders). We quantified transcripts of CENPN, an essential gene disrupted by the inversion at the proximal breakpoint, at different exons distributed near and across the breakpoint region. Consistent with dosage dependent gene expression for the breakpoint gene CENPN, we found that expression in Independents was broadly similar for transcripts segments from inside and outside the inversion regions, whereas for Satellites and Faeders, transcript segments outside of the inversion showed at least twofold higher expression than those spanning over the breakpoint. Within the inversion, observed expression differences for inversion males across all four genes with allele-specific primers were consistent with allelic imbalance. We further analyzed gonadal expression of two inversion genes, HSD17B2 and SDR42E1, along with 12 non-inversion genes related to steroid metabolism and signaling in 25 males (13 Independents, 7 Satellites, 5 Faeders). Although we did not find clear morph differentiation for many individual genes, all three morphs could be separated based on gene expression differences when using linear discriminant analysis (LDA), regardless of genomic location (i.e., inside or outside of the inversion). This was robust to the removal of genes with the highest loadings. Pairwise correlations in the expression of genes showed significant correlations for 9–18 pairs of genes within morphs. However, between morphs, we only found a single association between genes SDR42E1 and AROM for Independents and Satellites. Our results suggest complex and wide-ranging changes in gene expression caused by structural variants.

Variation in pigmentation gene expression is associated with distinct aposematic color morphs in the poison frog Dendrobates auratus

Background

Color and pattern phenotypes have clear implications for survival and reproduction in many species. However, the mechanisms that produce this coloration are still poorly characterized, especially at the genomic level. Here we have taken a transcriptomics-based approach to elucidate the underlying genetic mechanisms affecting color and pattern in a highly polytypic poison frog. We sequenced RNA from the skin from four different color morphs during the final stage of metamorphosis and assembled a de novo transcriptome. We then investigated differential gene expression, with an emphasis on examining candidate color genes from other taxa.

Results

Overall, we found differential expression of a suite of genes that control melanogenesis, melanocyte differentiation, and melanocyte proliferation (e.g., tyrp1, lef1, leo1, and mitf) as well as several differentially expressed genes involved in purine synthesis and iridophore development (e.g., arfgap1, arfgap2, airc, and gart).

Conclusions

Our results provide evidence that several gene networks known to affect color and pattern in vertebrates play a role in color and pattern variation in this species of poison frog.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1410-7) contains supplementary material, which is available to authorized users.

De novo assembly of the dual transcriptomes of a polymorphic raptor species and its malarial parasite

BACKGROUND

Studies of non-model species are important for understanding the molecular processes underpinning phenotypic variation under natural ecological conditions. The common buzzard (Buteo buteo; Aves: Accipitriformes) is a widespread and common Eurasian raptor with three distinct plumage morphs that differ in several fitness-related traits, including parasite infestation. To provide a genomic resource for plumage polymorphic birds in general and to search for candidate genes relating to fitness, we generated a transcriptome from a single dead buzzard specimen plus easily accessible, minimally invasive samples from live chicks.

RESULTS

We not only de novo assembled a near-complete buzzard transcriptome, but also obtained a significant fraction of the transcriptome of its malaria-like parasite, Leucocytozoon buteonis. By identifying melanogenesis-related transcripts that are differentially expressed in light ventral and dark dorsal feathers, but which are also expressed in other regions of the body, we also identified a suite of candidate genes that could be associated with fitness differences among the morphs. These include several immune-related genes, providing a plausible link between melanisation and parasite load. qPCR analysis of a subset of these genes revealed significant differences between ventral and dorsal feathers and an additional effect of morph.

CONCLUSION

This new resource provides preliminary insights into genes that could be involved in fitness differences between the buzzard colour morphs, and should facilitate future studies of raptors and their malaria-like parasites.

A supergene determines highly divergent male reproductive morphs in the ruff

Three strikingly different alternative male mating morphs (aggressive “Independents”, semi-cooperative “Satellites” and female mimic “Faeders”) coexist as a balanced polymorphism in the ruff, Philomachus pugnax, a lek-breeding wading bird^1,2,3. Major differences in body size, ornamentation, and aggressive and mating behaviour are inherited as an autosomal polymorphism^4,5. We show that development into Satellites and Faeders is determined by a supergene^6,7,8 consisting of divergent alternative, dominant, non-recombining haplotypes of an inversion on chromosome 11, which contains 125 predicted genes. Independents are homozygous for the ancestral sequence. One breakpoint of the inversion disrupts the essential Centromere protein N (CENP-N) gene, and pedigree analysis confirms lethality of inversion homozygotes. We describe novel behavioural, testes size, and steroid metabolic differences among morphs, and identify polymorphic genes within the inversion that are likely to contribute to the differences among morphs in reproductive traits.

The Genome of the "Great Speciator" Provides Insights into Bird Diversification

Among birds, white-eyes (genus Zosterops) have diversified so extensively that Jared Diamond and Ernst Mayr referred to them as the "great speciator." The Zosterops lineage exhibits some of the fastest rates of species diversification among vertebrates, and its members are the most prolific passerine island colonizers. We present a high-quality genome assembly for the silvereye (Zosterops lateralis), a white-eye species consisting of several subspecies distributed across multiple islands. We investigate the genetic basis of rapid diversification in white-eyes by conducting genomic analyses at varying taxonomic levels. First, we compare the silvereye genome with those of birds from different families and searched for genomic features that may be unique to Zosterops. Second, we compare the genomes of different species of white-eyes from Lifou island (South Pacific), using whole genome resequencing and restriction site associated DNA. Third, we contrast the genomes of two subspecies of silvereye that differ in plumage color. In accordance with theory, we show that white-eyes have high rates of substitutions, gene duplication, and positive selection relative to other birds. Below genus level, we find that genomic differentiation accumulates rapidly and reveals contrasting demographic histories between sympatric species on Lifou, indicative of past interspecific interactions. Finally, we highlight genes possibly involved in color polymorphism between the subspecies of silvereye. By providing the first whole-genome sequence resources for white-eyes and by conducting analyses at different taxonomic levels, we provide genomic evidence underpinning this extraordinary bird radiation.

Major breeding plumage color differences of male ruffs (Philomachus pugnax) are not associated with coding sequence variation in the MC1R gene

Sequence variation in the melanocortin-1 receptor (MC1R) gene explains color morph variation in several species of birds and mammals. Ruffs (Philomachus pugnax) exhibit major dark/light color differences in melanin-based male breeding plumage which is closely associated with alternative reproductive behavior. A previous study identified a microsatellite marker (Ppu020) near the MC1R locus associated with the presence/absence of ornamental plumage. We investigated whether coding sequence variation in the MC1R gene explains major dark/light plumage color variation and/or the presence/absence of ornamental plumage in ruffs. Among 821bp of the MC1R coding region from 44 male ruffs we found 3 single nucleotide polymorphisms, representing 1 nonsynonymous and 2 synonymous amino acid substitutions. None were associated with major dark/light color differences or the presence/absence of ornamental plumage. At all amino acid sites known to be functionally important in other avian species with dark/light plumage color variation, ruffs were either monomorphic or the shared polymorphism did not coincide with color morph. Neither ornamental plumage color differences nor the presence/absence of ornamental plumage in ruffs are likely to be caused entirely by amino acid variation within the coding regions of the MC1R locus. Regulatory elements and structural variation at other loci may be involved in melanin expression and contribute to the extreme plumage polymorphism observed in this species.

Evolution of endothelin receptors in vertebrates

Endothelin receptors are G protein coupled receptors (GPCRs) of the β-group of rhodopsin receptors that bind to endothelin ligands, which are 21 amino acid long peptides derived from longer prepro-endothelin precursors. The most basal Ednr-like GPCR is found outside vertebrates in the cephalochordate amphioxus, but endothelin ligands are only present among vertebrates, including the lineages of jawless vertebrates (lampreys and hagfishes), cartilaginous vertebrates (sharks, rays, and chimaeras), and bony vertebrates (ray-finned fishes and lobe-finned vertebrates including tetrapods). A bona fide endothelin system is thus a vertebrate-specific innovation with important roles for regulating the cardiovascular system, renal and pulmonary processes, as well as for the development of the vertebrate-specific neural crest cell population and its derivatives. Expectedly, dysregulation of endothelin receptors and the endothelin system leads to a multitude of human diseases. Despite the importance of different types of endothelin receptors for vertebrate development and physiology, current knowledge on endothelin ligand-receptor interactions, on the expression of endothelin receptors and their ligands, and on the functional roles of the endothelin system for embryonic development and in adult vertebrates is very much biased towards amniote vertebrates. Recent analyses from a variety of vertebrate lineages, however, have shown that the endothelin system in lineages such as teleost fish and lampreys is more diverse and is divergent from the mammalian endothelin system. This diversity is mainly based on differential evolution of numerous endothelin system components among vertebrate lineages generated by two rounds of whole genome duplication (three in teleosts) during vertebrate evolution. Here we review current understanding of the evolutionary history of the endothelin receptor family in vertebrates supplemented with surveys on the endothelin receptor gene complement of newly available genome assemblies from phylogenetically informative taxa. Our assessment further highlights the diversity of the vertebrate endothelin system and calls for detailed functional and pharmacological analyses of the endothelin system beyond tetrapods.

A comparison of next generation sequencing technologies for transcriptome assembly and utility for RNA-Seq in a non-model bird

De novo assembled transcriptomes, in combination with RNA-Seq, are powerful tools to explore gene sequence and expression level in organisms without reference genomes. Investigators must first choose which high throughput sequencing platforms will provide data most suitable for their experimental goals. In this study, we explore the utility of 454 and Illumina sequences for de novo transcriptome assembly and downstream RNA-Seq applications in a reproductive gland from a non-model bird species, the Japanese quail (Coturnix japonica). Four transcriptomes composed of either pure 454 or Illumina reads or mixtures of read types were assembled and evaluated for the same cost. Illumina assemblies performed best for de novo transcriptome characterization in terms of contig length, transcriptome coverage, and complete assembly of gene transcripts. Improvements over the Hybrid assembly were marginal, with the exception that the addition of 454 data significantly increased the number of genes annotated. The Illumina assembly provided the best reference to align an independent set of RNA-Seq data as ∼84% of reads mapped to single genes in the transcriptome. Contigs constructed solely from 454 data may impose problems for RNA-Seq as our 454 transcriptome revealed a high number of indels and many ambiguously mapped reads. Correcting the 454 transcriptome with Illumina reads was an effective strategy to deal with indel and frameshift errors inherent to the 454 transcriptome, but at the cost of transcriptome coverage. In the absence of a reference genome, we find that Illumina reads alone produced a high quality transcriptome appropriate for RNA-Seq gene expression analyses.

Characterization of the house sparrow (Passer domesticus) transcriptome: a resource for molecular ecology and immunogenetics

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.

Pleiotropy in the melanocortin system: expression levels of this system are associated with melanogenesis and pigmentation in the tawny owl (Strix aluco)

The adaptive function of melanin-based coloration is a long-standing debate. A recent genetic model suggested that pleiotropy could account for covariations between pigmentation, behaviour, morphology, physiology and life history traits. We explored whether the expression levels of genes belonging to the melanocortin system (MC1R, POMC, PC1/3, PC2 and the antagonist ASIP), which have many pleiotropic effects, are associated with melanogenesis (through variation in the expression of the genes MITF, SLC7A11, TYR, TYRP1) and in turn melanin-based coloration. We considered the tawny owl (Strix aluco) because individuals vary continuously from light to dark reddish, and thus, colour variation is likely to stem from differences in the levels of gene expression. We measured gene expression in feather bases collected in nestlings at the time of melanin production. As expected, the melanocortin system was associated with the expression of melanogenic genes and pigmentation. Offspring of darker reddish fathers expressed PC1/3 to lower levels but tended to express PC2 to higher levels. The convertase enzyme PC1/3 cleaves the POMC prohormone to obtain ACTH, while the convertase enzyme PC2 cleaves ACTH to produce α-melanin-stimulating hormone (α-MSH). ACTH regulates glucocorticoids, hormones that modulate stress responses, while α-MSH induces eumelanogenesis. We therefore conclude that the melanocortin system, through the convertase enzymes PC1/3 and PC2, may account for part of the interindividual variation in melanin-based coloration in nestling tawny owls. Pleiotropy may thus account for the covariation between phenotypic traits involved in social interactions (here pigmentation) and life history, morphology, behaviour and physiology.

Genetics of colouration in birds

Establishing the links between phenotype and genotype is of great importance for resolving key questions about the evolution, maintenance and adaptive function of phenotypic variation. Bird colouration is one of the most studied systems to investigate the role of natural and sexual selection in the evolution of phenotypic diversity. Given the recent advances in molecular tools that allow discovering genetic polymorphisms and measuring gene and protein expression levels, it is timely to review the literature on the genetics of bird colouration. The present study shows that melanin-based colour phenotypes are often associated with mutations at melanogenic genes. Differences in melanin-based colouration are caused by switches of eumelanin to pheomelanin production or by changes in feather keratin structure, melanoblast migration and differentiation, as well as melanosome structure. Similar associations with other types of colourations are difficult to establish, because our knowledge about the molecular genetics of carotenoid-based and structural colouration is quasi inexistent. This discrepancy stems from the fact that only melanin-based colouration shows pronounced heritability estimates, i.e. the resemblance between related individuals is usually mainly explained by genetic factors. In contrast, the expression of carotenoid-based colouration is phenotypically plastic with a high sensitivity to variation in environmental conditions. It therefore appears that melanin-based colour traits are prime systems to understand the genetic basis of phenotypic variation. In this context, birds have a great potential to bring us to new frontiers where many exciting discoveries will be made on the genetics of phenotypic traits, such as colouration. In this context, a major goal of our review is to suggest a number of exciting future avenues.