Characterization of the genome and transcriptome of the blue tit Cyanistes caeruleus: polymorphisms, sex-biased expression and selection signals

Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

It is common that males and females display sexual dimorphisms, which usually result from sex-biased gene expression. Chinese hwamei (Garrulax canorus) is a good model for studying sex-biased gene expression because the song between the sexes is quite different. In this study, we analyze cerebrum and syrinx sex-biased gene expression and evolution using the de novo assembled Chinese hwamei transcriptome. In both the cerebrum and syrinx, our study revealed that most female-biased genes were actively expressed in females only, while most male-biased genes were actively expressed in both sexes. In addition, both male- and female-biased genes were enriched on the putative Z chromosome, suggesting the existence of sexually antagonistic genes and the insufficient dosage compensation of the Z-linked genes. We also identified a 9 Mb sex linkage region on the putative 4A chromosome which enriched more than 20% of female-biased genes. Resultantly, male-biased genes in both tissues had significantly higher Ka/Ks and effective number of codons (ENCs) than unbiased genes, and this suggested that male-biased genes which exhibit accelerated divergence may have resulted from positive selection. Taken together, our results initially revealed the reasons for the differences in singing behavior between males and females of Chinese hwamei.

Sequence Transpositions Restore Genes on the Highly Degenerated W Chromosomes of Songbirds

The female-specific W chromosomes of most Neognathae birds are highly degenerated and gene-poor. Previous studies have demonstrated that the gene repertoires of the Neognathae bird W chromosomes, despite being in small numbers, are conserved across bird species, likely due to purifying selection maintaining the regulatory and dosage-sensitive genes. Here we report the discovery of DNA-based sequence duplications from the Z to the W chromosome in birds-of-paradise (Paradisaeidae, Passeriformes), through sequence transposition. The original transposition involved nine genes, but only two of them (ANXA1 and ALDH1A1) survived on the W chromosomes. Both ANXA1 and ALDH1A1 are predicted to be dosage-sensitive, and the expression of ANXA1 is restricted to ovaries in all the investigated birds. These analyses suggest the newly transposed gene onto the W chromosomes can be favored for their role in restoring dosage imbalance or through female-specific selection. After examining seven additional songbird genomes, we further identified five other transposed genes on the W chromosomes of Darwin's finches and one in the great tit, expanding the observation of the Z-to-W transpositions to a larger range of bird species, but not all transposed genes exhibit dosage-sensitivity or ovary-biased expression We demonstrate a new mechanism by which the highly degenerated W chromosomes of songbirds can acquire genes from the homologous Z chromosomes, but further functional investigations are needed to validate the evolutionary forces underlying the transpositions.

Negative effects of individual heterozygosity on reproductive success in a wild bird population

The evolutionary consequences of individual genetic diversity are frequently studied by assessing heterozygosity-fitness correlations (HFCs). The prevalence of positive and negative HFCs and the predominance of general versus local effects in wild populations are far from understood, partly because comprehensive studies testing for both inbreeding and outbreeding depression are lacking. We studied a genetically diverse population of blue tits in southern Germany using a genome-wide set of 87 microsatellites to investigate the relationship between proxies of reproductive success and measures of multilocus and single-locus individual heterozygosity (MLH and SLH). We used complimentary measures of MLH and partitioned markers into functional categories according to their position in the blue tit genome. HFCs based on MLH were consistently negative for functional loci, whereas correlations were rather inconsistent for loci found in nonfunctional areas of the genome. Clutch size was the only reproductive variable showing a general effect. We found evidence for local effects for three measures of reproductive success: arrival date at the breeding site, the probability of breeding at the study site and male reproductive success. For these, we observed consistent, and relatively strong, negative effects at one functional locus. Remarkably, this marker had a similar effect in another blue tit population from Austria (~400 km to the east). We suggest that a genetic local effect on timing of arrival might be responsible for most negative HFCs detected, with carry-over effects on other reproductive traits. This effect could reflect individual differences in the distance between overwintering areas and breeding sites.

Regional epigenetic differentiation of the Z Chromosome between sexes in a female heterogametic system

In male heterogametic systems, the X Chromosome is epigenetically differentiated between males and females, to facilitate dosage compensation. For example, the X Chromosome in female mammals is largely inactivated. Relative to well-studied male heterogametic systems, the extent of epigenetic differentiation between male and female Z Chromosomes in female heterogametic species, which often lack complete dosage compensation, is poorly understood. Here, we examined the chromosomal DNA methylation landscapes of male and female Z Chromosomes in two distantly related avian species, namely chicken and white-throated sparrow. We show that, in contrast to the pattern in mammals, male and female Z Chromosomes in these species exhibit highly similar patterns of DNA methylation, which is consistent with weak or absent dosage compensation. We further demonstrate that the epigenetic differences between male and female chicken Z Chromosomes are localized to a few regions, including a previously identified male hypermethylated region 1 (MHM1; CGNC: 80601). We discovered a novel region with elevated male-to-female methylation ratios on the chicken Z Chromosome (male hypermethylated region 2 [MHM2]; CGNC: 80602). The MHM1 and MHM2, despite little sequence similarity between them, bear similar molecular features that are likely associated with their functions. We present evidence consistent with female hypomethylation of MHMs and up-regulation of nearby genes. Therefore, despite little methylation differentiation between sexes, extremely localized DNA methylation differences between male and female chicken Z Chromosomes have evolved and affect expression of nearby regions. Our findings offer new insights into epigenetic regulation of gene expression between sexes in female heterogametic systems.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

Ice-Age Climate Adaptations Trap the Alpine Marmot in a State of Low Genetic Diversity

Some species responded successfully to prehistoric changes in climate [1, 2], while others failed to adapt and became extinct [3]. The factors that determine successful climate adaptation remain poorly understood. We constructed a reference genome and studied physiological adaptations in the Alpine marmot (Marmota marmota), a large ground-dwelling squirrel exquisitely adapted to the "ice-age" climate of the Pleistocene steppe [4, 5]. Since the disappearance of this habitat, the rodent persists in large numbers in the high-altitude Alpine meadow [6, 7]. Genome and metabolome showed evidence of adaptation consistent with cold climate, affecting white adipose tissue. Conversely, however, we found that the Alpine marmot has levels of genetic variation that are among the lowest for mammals, such that deleterious mutations are less effectively purged. Our data rule out typical explanations for low diversity, such as high levels of consanguineous mating, or a very recent bottleneck. Instead, ancient demographic reconstruction revealed that genetic diversity was lost during the climate shifts of the Pleistocene and has not recovered, despite the current high population size. We attribute this slow recovery to the marmot's adaptive life history. The case of the Alpine marmot reveals a complicated relationship between climatic changes, genetic diversity, and conservation status. It shows that species of extremely low genetic diversity can be very successful and persist over thousands of years, but also that climate-adapted life history can trap a species in a persistent state of low genetic diversity.

Genomes of Three Closely Related Caribbean Amazons Provide Insight for Species History and Conservation

Islands have been used as model systems for studies of speciation and extinction since Darwin published his observations about finches found on the Galapagos. Amazon parrots inhabiting the Greater Antillean Islands represent a fascinating model of species diversification. Unfortunately, many of these birds are threatened as a result of human activity and some, like the Puerto Rican parrot, are now critically endangered. In this study we used a combination of de novo and reference-assisted assembly methods, integrating it with information obtained from related genomes to perform genome reconstruction of three amazon species. First, we used whole genome sequencing data to generate a new de novo genome assembly for the Puerto Rican parrot (Amazona vittata). We then improved the obtained assembly using transcriptome data from Amazona ventralis and used the resulting sequences as a reference to assemble the genomes Hispaniolan (A. ventralis) and Cuban (Amazona leucocephala) parrots. Finally, we, annotated genes and repetitive elements, estimated genome sizes and current levels of heterozygosity, built models of demographic history and provided interpretation of our findings in the context of parrot evolution in the Caribbean.

Life history shapes variation in egg composition in the blue tit Cyanistes caeruleus

Maternal investment directly shapes early developmental conditions and therefore has long-term fitness consequences for the offspring. In oviparous species prenatal maternal investment is fixed at the time of laying. To ensure the best survival chances for most of their offspring, females must equip their eggs with the resources required to perform well under various circumstances, yet the actual mechanisms remain unknown. Here we describe the blue tit egg albumen and yolk proteomes and evaluate their potential to mediate maternal effects. We show that variation in egg composition (proteins, lipids, carotenoids) primarily depends on laying order and female age. Egg proteomic profiles are mainly driven by laying order, and investment in the egg proteome is functionally biased among eggs. Our results suggest that maternal effects on egg composition result from both passive and active (partly compensatory) mechanisms, and that variation in egg composition creates diverse biochemical environments for embryonic development.

Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

Characterization of a de novo assembled transcriptome of the Common Blackbird (Turdus merula)

BACKGROUND

In recent years, next generation high throughput sequencing technologies have proven to be useful tools for investigations concerning the genomics or transcriptomics also of non-model species. Consequently, ornithologists have adopted these technologies and the respective bioinformatics tools to survey the genomes and transcriptomes of a few avian non-model species. The Common Blackbird is one of the most common bird species living in European cities, which has successfully colonized urban areas and for which no reference genome or transcriptome is publicly available. However, to target questions like genome wide gene expression analysis, a reference genome or transcriptome is needed.

METHODS

Therefore, in this study two Common Blackbirds were sacrificed, their mRNA was isolated and analyzed by RNA-Seq to de novo assemble a transcriptome and characterize it. Illumina reads (125 bp paired-end) and a Velvet/Oases pipeline led to 162,158 transcripts. For the annotation (using Blast+), an unfiltered protein database was used. SNPs were identified using SAMtools and BCFtools. Furthermore, mRNA from three single tissues (brain, heart and liver) of the same two Common Blackbirds were sequenced by Illumina (75 bp single-end reads). The draft transcriptome and the three single tissues were compared by their BLAST hits with the package VennDiagram in R.

RESULTS

Following the annotation against protein databases, we found evidence for 15,580 genes in the transcriptome (all well characterized hits after annotation). On 18% of the assembled transcripts, 144,742 SNPs were identified which are, consequently, 0.09% of all nucleotides in the assembled transcriptome. In the transcriptome and in the single tissues (brain, heart and liver), 10,182 shared genes were found.

DISCUSSION

Using a next-generation technology and bioinformatics tools, we made a first step towards the genomic investigation of the Common Blackbird. The de novo assembled transcriptome is usable for downstream analyses such as differential gene expression analysis and SNP identification. This study shows the importance of the approach to sequence single tissues to understand functions of tissues, proteins and the phenotype.

Assembly and RNA-free annotation of highly heterozygous genomes: The case of the thick-billed murre (Uria lomvia)

Thanks to a dramatic reduction in sequencing costs followed by a rapid development of bioinformatics tools, genome assembly and annotation have become accessible to many researchers in recent years. Among tetrapods, birds have genomes that display many features that facilitate their assembly and annotation, such as small genome size, low number of repeats and highly conserved genomic structure. However, we found that high genomic heterozygosity could have a great impact on the quality of the genome assembly of the thick-billed murre (Uria lomvia), an arctic colonial seabird. In this study, we tested the performance of three genome assemblers, ray/sscape, soapdenovo2 and platanus, in assembling the highly heterozygous genome of the thick-billed murre. Our results show that platanus, an assembler specifically designed for heterozygous genomes, outperforms the other two approaches and produces a highly contiguous (N50 = 15.8 Mb) and complete genome assembly (93% presence of genes from the Benchmarking Universal Single Copy Ortholog [BUSCO] gene set). Additionally, we annotated the thick-billed murre genome using a homology-based approach that takes advantage of the genomic resources available for birds and other taxa. Our study will be useful for those researchers who are approaching assembly and annotation of highly heterozygous genomes, or genomes of species of conservation concern, and/or who have limited financial resources.

rehh 2.0: a reimplementation of the R package rehh to detect positive selection from haplotype structure

Identifying genomic regions with unusually high local haplotype homozygosity represents a powerful strategy to characterize candidate genes responding to natural or artificial positive selection. To that end, statistics measuring the extent of haplotype homozygosity within (e.g. EHH, iHS) and between (Rsb or XP-EHH) populations have been proposed in the literature. The rehh package for r was previously developed to facilitate genome-wide scans of selection, based on the analysis of long-range haplotypes. However, its performance was not sufficient to cope with the growing size of available data sets. Here, we propose a major upgrade of the rehh package, which includes an improved processing of the input files, a faster algorithm to enumerate haplotypes, as well as multithreading. As illustrated with the analysis of large human haplotype data sets, these improvements decrease the computation time by more than one order of magnitude. This new version of rehh will thus allow performing iHS-, Rsb- or XP-EHH-based scans on large data sets. The package rehh 2.0 is available from the CRAN repository (http://cran.r-project.org/web/packages/rehh/index.html) together with help files and a detailed manual.

An ecological model organism flies into the genomics era

Despite the very rapid 'genomicization' of the field of Molecular Ecology in recent years, there have been relatively few annotated whole-genome assemblies of nonmodel organisms published. Instead, molecular ecologists have more frequently utilized next-generation sequencing technologies to develop genome-wide markers or to generate transcriptome data. Whole-genome assemblies are more expensive and require considerable computational resources and bioinformatic expertise. However, the availability of an annotated genome offers exciting opportunities to address fundamental questions in ecology and evolution that are difficult to address with moderate sets of markers or by transcriptome sequencing. Such questions include elucidating the roles of natural and sexual selection in shaping diversity, determining the roles of regulatory and protein-coding change in the evolution of traits, and determining the genomic architecture of sex-specific trait variation. Arguably, these questions are most tractable--and most interesting--in well-characterized species for which there is already some knowledge of natural and sexual selection, and of the traits that are most likely to link to fitness. In this issue, Mueller et al. (2016) present the assembly and annotation of the genome of the blue tit (Cyanistes caeruleus), a model ecological species. In addition, by sequencing the transcriptome of male and female blue tits, the authors identify and annotate sex-biased gene expression and conclude that noncoding RNA genes are likely to play a significant role in sex-biased expression. By making their assembly and annotation publically available and accessible via a genome browser, Mueller et al. (2016) offer exciting possibilities for further research into the genomic basis of adaptation, and investigation of the roles of natural and sexual selection, in this well-studied ecological model species.

Adaptation to nocturnality - learning from avian genomes

The recent availability of multiple avian genomes has laid the foundation for a huge variety of comparative genomics analyses including scans for changes and signatures of selection that arose from adaptions to new ecological niches. Nocturnal adaptation in birds, unlike in mammals, is comparatively recent, a fact that makes birds good candidates for identifying early genetic changes that support adaptation to dim-light environments. In this review, we give examples of comparative genomics analyses that could shed light on mechanisms of adaptation to nocturnality. We present advantages and disadvantages of both "data-driven" and "hypothesis-driven" approaches that lead to the discovery of candidate genes and genetic changes promoting nocturnality. We anticipate that the accessibility of multiple genomes from the Genome 10K Project will allow a better understanding of evolutionary mechanisms and adaptation in general.