De novo transcriptome sequencing in a songbird, the dark-eyed junco (Junco hyemalis): genomic tools for an ecological model system

Effects of short-term experimental manipulation of captive social environment on uropygial gland microbiome and preen oil volatile composition

Introduction

Avian preen oil, secreted by the uropygial gland, is an important source of volatile compounds that convey information about the sender’s identity and quality, making preen oil useful for the recognition and assessment of potential mates and rivals. Although intrinsic factors such as hormone levels, genetic background, and diet can affect preen oil volatile compound composition, many of these compounds are not the products of the animal’s own metabolic processes, but rather those of odor-producing symbiotic microbes. Social behavior affects the composition of uropygial microbial communities, as physical contact results in microbe sharing. We experimentally manipulated social interactions in captive dark-eyed juncos (Junco hyemalis) to assess the relative influence of social interactions, subspecies, and sex on uropygial gland microbial composition and the resulting preen oil odor profiles.

Methods

We captured 24 birds at Mountain Lake Biological Station in Virginia, USA, including birds from two seasonally sympatric subspecies – one resident, one migratory. We housed them in an outdoor aviary in three phases of social configurations: first in same-sex, same-subspecies flocks, then in male-female pairs, and finally in the original flocks. Using samples taken every four days of the experiment, we characterized their uropygial gland microbiome through 16S rRNA gene sequencing and their preen oil volatile compounds via GC-MS.

Results

We predicted that if social environment was the primary driver of uropygial gland microbiome composition, and if microbiome composition in turn affected preen oil volatile profiles, then birds housed together would become more similar over time. Our results did not support this hypothesis, instead showing that sex and subspecies were stronger predictors of microbiome composition. We observed changes in volatile compounds after the birds had been housed in pairs, which disappeared after they were moved back into flocks, suggesting that hormonal changes related to breeding condition were the most important factor in these patterns.

Discussion

Although early life social environment of nestlings and long-term social relationships have been shown to be important in shaping uropygial gland microbial communities, our study suggests that shorter-term changes in social environment do not have a strong effect on uropygial microbiomes and the resulting preen oil volatile compounds.

A High-Quality Genome Assembly of the North American Song Sparrow, Melospiza melodia

The song sparrow, Melospiza melodia, is one of the most widely distributed species of songbirds found in North America. It has been used in a wide range of behavioral and ecological studies. This species' pronounced morphological and behavioral diversity across populations makes it a favorable candidate in several areas of biomedical research. We have generated a high-quality de novo genome assembly of M. melodia using Illumina short read sequences from genomic and in vitro proximity-ligation libraries. The assembled genome is 978.3 Mb, with a physical coverage of 24.9×, N50 scaffold size of 5.6 Mb and N50 contig size of 31.7 Kb. Our genome assembly is highly complete, with 87.5% full-length genes present out of a set of 4,915 universal single-copy orthologs present in most avian genomes. We annotated our genome assembly and constructed 15,086 gene models, a majority of which have high homology to related birds, Taeniopygia guttata and Junco hyemalis In total, 83% of the annotated genes are assigned with putative functions. Furthermore, only ∼7% of the genome is found to be repetitive; these regions and other non-coding functional regions are also identified. The high-quality M. melodia genome assembly and annotations we report will serve as a valuable resource for facilitating studies on genome structure and evolution that can contribute to biomedical research and serve as a reference in population genomic and comparative genomic studies of closely related species.

Genomes to space stations: the need for the integrative study of migration for avian conservation

Ongoing changes to global weather patterns and human modifications of the environment have altered the breeding and non-breeding ranges of migratory species, the timing of their migrations, and even whether they continue to migrate at all. Animal movements are arguably one of the most difficult behaviours to study, particularly in smaller birds that migrate tens to thousands of kilometres seasonally, often moving hundreds of kilometres each day. The recent miniaturization of tracking and logging devices has led to a radical transformation in our understanding of avian migratory behaviour and migratory connectivity. While advances in technology have altered the way researchers study migratory behaviour in the field, advances in techniques related to the study of physiological and genetic mechanisms underlying migratory behaviour have rarely been integrated into field studies of tracking. To predict the capacity of migrants to adjust to a changing planet, it is essential that we combine avian migration data with physiological and genetic measurements taken at key time points prior to, during and after migration.

De novo annotation of the transcriptome of the Northern Wheatear (Oenanthe oenanthe)

We have sequenced a partial transcriptome of the Northern Wheatear (Oenanthe oenanthe), a species with one of the longest migrations on Earth. The transcriptome was constructed de novo using RNA-Seq sequence data from the pooled mRNA of six different tissues: brain, muscle, intestine, liver, adipose tissue and skin. The samples came from nine captive-bred wheatears collected at three different stages of the endogenous autumn migratory period: (1) lean birds prior the onset of migration, (2) during the fattening stage and (3) individuals at their migratory body mass plateau, when they have almost doubled their lean body mass. The sample structure used to build up the transcriptome of the Northern Wheatears concerning tissue composition and time guarantees the future survey of the regulatory genes involved in the development of the migratory phenotype. Through the pre-migratory period, birds accomplish outstanding physical and behavioural changes that involve all organ systems. Nevertheless, the molecular mechanisms through which birds synchronize and control hyperphagia, fattening, restlessness increase, immunity boosting and tuning the muscles for such endurance flight are still largely unknown. The use of RNA-Seq has emerged as a powerful tool to analyse complex traits on a broad scale, and we believe it can help to characterize the migratory phenotype of wheatears at an unprecedented level. The primary challenge to conduct quantitative transcriptomic studies in non-model species is the availability of a reference transcriptome, which we have constructed and described in this paper. The cDNA was sequenced by pyrosequencing using the Genome Sequencer Roche GS FLX System; with single paired-end reads of about 400 bp. We estimate the total number of genes at 15,640, of which  67% could be annotated using Turkey and Zebra Finch genomes, or protein sequence information from SwissProt and NCBI databases. With our study, we have made a first step towards understanding the migratory phenotype regarding gene expression of a species that has become a model to study birds long-distance migrations.

The effect of chronic and acute stressors, and their interaction, on testes function: an experimental test during testicular recrudescence

Organisms are expected to invest less in reproduction in response to a stressor, but theory predicts that this effect should depend on the frequency and duration of stressors in the environment. Here, we investigated how an acute stressor affected testes function in a songbird, and how chronic stressors influenced the acute stress response. We exposed male dark-eyed juncos (Junco hyemalis) either to chronic or minimal (control) disturbance during testicular recrudescence, after which we measured baseline testosterone, testosterone after an acute handling stressor, and capacity to produce testosterone after hormonal stimulation. In a 2×2 design, we then killed males from the two chronic treatment groups either immediately or after an acute stressor to investigate the effect of long- and short-term stressors on the testicular transcriptome. We found that chronically disturbed birds had marginally lower baseline testosterone. The acute stressor suppressed testosterone in control birds, but not in the chronic disturbance group. The ability to elevate testosterone did not differ between the chronic treatments. Surprisingly, chronic disturbance had a weak effect on the testicular transcriptome, and did not affect the transcriptomic response to the acute stressor. The acute stressor, on the other hand, upregulated the cellular stress response and affected expression of genes associated with hormonal stress response. Overall, we show that testicular function is sensitive to acute stressors but surprisingly robust to long-term stressors, and that chronic disturbance attenuates the decrease in testosterone in response to an acute stressor.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

De novo transcriptome analysis shows differential expression of genes in salivary glands of edible bird's nest producing swiftlets

BACKGROUND

Edible bird's nest (EBN), produced from solidified saliva secretions of specific swiftlet species during the breeding season, is one of the most valuable animal by-products in the world. The composition and medicinal benefits of EBN have been extensively studied, however, genomic and transcriptomic studies of the salivary glands of these birds have not been conducted.

RESULTS

The study described the transcriptomes of salivary glands from three swiftlet species (28 samples) generated by RNASeq. A total of 14,835 annotated genes and 428 unmapped genes were cataloged. The current study investigated the genes and pathways that are associated with the development of salivary gland and EBN composition. Differential expression and pathway enrichment analysis indicated that the expression of CREB3L2 and several signaling pathways involved in salivary gland development, namely, the EGFR, BMP, and MAPK signaling pathways, were up-regulated in swiftlets producing white EBN (Aerodramus fuciphagus) and black EBN (Aerodramus maximus) compared with non-EBN-producing swiftlets (Apus affinis). Furthermore, MGAT, an essential gene for the biosynthesis of N-acetylneuraminic acid (sialic acid), was highly expressed in both white- and black-nest swiftlets compared to non-EBN-producing swiftlets. Interspecies comparison between Aerodramus fuciphagus and Aerodramus maximus indicated that the genes involved in N-acetylneuraminic and fatty acid synthesis were up-regulated in Aerodramus fuciphagus, while alanine and aspartate synthesis pathways were up-regulated in Aerodramus maximus. Furthermore, gender-based analysis revealed that N-glycan trimming pathway was significantly up-regulated in male Aerodramus fuciphagus from its natural habitat (cave) compared to their female counterpart.

CONCLUSIONS

Transcriptomic analysis of salivary glands of different swiftlet species reveal differential expressions of candidate genes that are involved in salivary gland development and in the biosynthesis of various bioactive compounds found in EBN.

Widespread patterns of sexually dimorphic gene expression in an avian hypothalamic-pituitary-gonadal (HPG) axis

The hypothalamic-pituitary-gonadal (HPG) axis is a key biological system required for reproduction and associated sexual behaviors to occur. In the avian reproductive model of the rock dove (Columba livia), we characterized the transcript community of each tissue of the HPG axis in both sexes, thereby significantly expanding our mechanistic insight into HPG activity. We report greater sex-biased differential expression in the pituitary as compared to the hypothalamus, with multiple genes more highly expressed in the male pituitary being related to secretory function, and multiple genes more highly expressed in the female pituitary being related to reproduction, growth, and development. We report tissue-specific and sex-biased expression in genes commonly investigated when studying reproduction, highlighting the need for sex parity in future studies. In addition, we uncover new targets of investigation in both sexes, which could potentially change our understanding of HPG function.

De novo transcriptome sequencing of a non-model polychaete species

Transcriptome sequencing is a useful method for studying gene-sequences data, especially in non-model organisms whose genomic sequences are yet to be determined. Indeed, even without any genome reference, transcripts can be assembled de novo to produce a genome-scale transcription map. Here we describe a de novo transcriptome assembly for the polychaete worm Ophryotrocha diadema, generated from 454-sequencing (Roche GS-FLX Titanium Sequencing). We describe the sequencing, the assembly and functional annotation of EST sequences, and the level of transcriptome coverage provided by our sequence data. The sequences we assembled represent a wide depiction of expressed genes from throughout the lifespan of O. diadema. We found gene sequences to be part of a broad range of functions (i.e. biological and metabolic processes, immune system, reproductive processes, development and growth, signalling and response to stimuli) and we therefore built the first transcriptomic reference for Ophryotrocha polychaete worms. These results encourage us to focus our investigation on many aspects of O. diadema and to extend our research to co-generic species.

Comparative genomics reveals molecular features unique to the songbird lineage

BACKGROUND

Songbirds (oscine Passeriformes) are among the most diverse and successful vertebrate groups, comprising almost half of all known bird species. Identifying the genomic innovations that might be associated with this success, as well as with characteristic songbird traits such as vocal learning and the brain circuits that underlie this behavior, has proven difficult, in part due to the small number of avian genomes available until recently. Here we performed a comparative analysis of 48 avian genomes to identify genomic features that are unique to songbirds, as well as an initial assessment of function by investigating their tissue distribution and predicted protein domain structure.

RESULTS

Using BLAT alignments and gene synteny analysis, we curated a large set of Ensembl gene models that were annotated as novel or duplicated in the most commonly studied songbird, the Zebra finch (Taeniopygia guttata), and then extended this analysis to 47 additional avian and 4 non-avian genomes. We identified 10 novel genes uniquely present in songbird genomes. A refined map of chromosomal synteny disruptions in the Zebra finch genome revealed that the majority of these novel genes localized to regions of genomic instability associated with apparent chromosomal breakpoints. Analyses of in situ hybridization and RNA-seq data revealed that a subset of songbird-unique genes is expressed in the brain and/or other tissues, and that 2 of these (YTHDC2L1 and TMRA) are highly differentially expressed in vocal learning-associated nuclei relative to the rest of the brain.

CONCLUSIONS

Our study reveals novel genes unique to songbirds, including some that may subserve their unique vocal control system, substantially improves the quality of Zebra finch genome annotations, and contributes to a better understanding of how genomic features may have evolved in conjunction with the emergence of the songbird lineage.

Behavioral effects of social challenges and genomic mechanisms of social priming: What's testosterone got to do with it?

Social challenges from rival conspecifics are common in the lives of animals, and changes in an animal's social environment can influence physiology and behavior in ways that appear to be adaptive in the face of continued social instability (i.e. social priming). Recently, it has become clear that testosterone, long thought to be the primary mediator of these effects, may not always change in response to social challenges, an observation that highlights gaps in our understanding of the proximate mechanisms by which animals respond to their social environment. Here, our goal is to address the degree to which testosterone mediates organismal responses to social cues. To this end, we review the behavioral and physiological consequences of social challenges, as well as their underlying hormonal and gene regulatory mechanisms. We also present a new case study from a wild songbird, the dark-eyed junco (Junco hyemalis), in which we find largely divergent genome-wide transcriptional changes induced by social challenges and testosterone, respectively, in muscle and liver tissue. Our review underscores the diversity of mechanisms that link the dynamic social environment with an organisms' genomic, hormonal, and behavioral state. This diversity among species, and even among tissues within an organism, reveals new insights into the pattern and process by which evolution may alter proximate mechanisms of social priming.

Opening the genetic toolbox of niche model organisms with high throughput techniques: novel proteins in regeneration as a case study

Understanding in vivo regeneration of complex structures offers a fascinating perspective for translation into medical applications. Unfortunately, mammals in general lack large-scale regenerative capacity, whereas planarians, newts or Hydra can regenerate complete body parts. Such organisms are, however, poorly annotated because of the lack of sequence information. This leads to limited access for molecular biological investigations. In the last decade, high throughput technologies and new methods enabling the effective generation of transgenic animals have rapidly evolved. These developments have allowed the extensive use of niche model organisms as part of a trend towards the accessibility of a greater panel of model species for scientific research. The case study that follows provides an insight into the impact of high throughput techniques on the landscape of models of regeneration. The cases presented here give evidence of alternative stem cell maintenance pathways, the identification of new protein families and new stem cell markers.

Profile of whole blood gene expression following immune stimulation in a wild passerine

Background

Immunoecology aims to explain variation among hosts in the strength and efficacy of immunological defences in natural populations. This requires development of biomarkers of the activation of the immune system so that they can be collected non-lethally and sampled from small amounts of easily obtainable tissue. We used transcriptome profiling in wild greenfinches (Carduelis chloris) to detect whole blood transcripts that most profoundly indicate upregulation of antimicrobial defences during acute phase response. The more general aim of this study was to obtain a functional annotation of a substantial portion of the greenfinch transcriptome that would enable to gain access to more specific genomic tools in subsequent studies. The birds received either bacterial lipopolysaccharide or saline injections and RNA-seq transcriptional profiling was performed 12 h after treatment to provide initial functional annotation of the transcriptome and assess whole blood response to immune stimulation.

Results

A total of 66,084 transcripts were obtained from de novo Trinty assembly, out of which 23,153 could be functionally annotated. Only 1,911 of these were significantly upregulated or downregulated. The manipulation caused marked upregulation of several transcripts related to immune activation. These included avian-specific antimicrobial agents avidin and gallinacin, but also some more general host response genes, such as serum amyloid A protein, lymphocyte antigen 75 and copper-transporting ATPase 1. However, links with avian immunity for most differentially regulated transcripts remained rather hypothetical, as a large set of differentially expressed transcripts lacked functional annotation.

Conclusions

This appears to be the first large scale transcriptional profiling of immune function in passerine birds. The transcriptomic data obtained suggest novel markers for the assessment of the immunological state of wild passerines. Characterizing the function of those possible novel infection markers would assist future vertebrate genome annotation. The extensive sequence information collected enables to identify possible target and housekeeping genes needed to gain access to more specific genomic tools in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-533) contains supplementary material, which is available to authorized users.

A house finch (Haemorhous mexicanus) spleen transcriptome reveals intra- and interspecific patterns of gene expression, alternative splicing and genetic diversity in passerines

Background

With its plumage color dimorphism and unique history in North America, including a recent population expansion and an epizootic of Mycoplasma gallisepticum (MG), the house finch (Haemorhous mexicanus) is a model species for studying sexual selection, plumage coloration and host-parasite interactions. As part of our ongoing efforts to make available genomic resources for this species, here we report a transcriptome assembly derived from genes expressed in spleen.

Results

We characterize transcriptomes from two populations with different histories of demography and disease exposure: a recently founded population in the eastern US that has been exposed to MG for over a decade and a native population from the western range that has never been exposed to MG. We utilize this resource to quantify conservation in gene expression in passerine birds over approximately 50 MY by comparing splenic expression profiles for 9,646 house finch transcripts and those from zebra finch and find that less than half of all genes expressed in spleen in either species are expressed in both species. Comparative gene annotations from several vertebrate species suggest that the house finch transcriptomes contain ~15 genes not yet found in previously sequenced vertebrate genomes. The house finch transcriptomes harbour ~85,000 SNPs, ~20,000 of which are non-synonymous. Although not yet validated by biological or technical replication, we identify a set of genes exhibiting differences between populations in gene expression (n = 182; 2% of all transcripts), allele frequencies (76 F_ST ouliers) and alternative splicing as well as genes with several fixed non-synonymous substitutions; this set includes genes with functions related to double-strand break repair and immune response.

Conclusions

The two house finch spleen transcriptome profiles will add to the increasing data on genome and transcriptome sequence information from natural populations. Differences in splenic expression between house finch and zebra finch imply either significant evolutionary turnover of splenic expression patterns or different physiological states of the individuals examined. The transcriptome resource will enhance the potential to annotate an eventual house finch genome, and the set of gene-based high-quality SNPs will help clarify the genetic underpinnings of host-pathogen interactions and sexual selection.

De novo sequencing, assembly, and analysis of the root transcriptome of Persea americana (Mill.) in response to Phytophthora cinnamomi and flooding

Avocado is a diploid angiosperm containing 24 chromosomes with a genome estimated to be around 920 Mb. It is an important fruit crop worldwide but is susceptible to a root rot caused by the ubiquitous oomycete Phytophthora cinnamomi. Phytophthora root rot (PRR) causes damage to the feeder roots of trees, causing necrosis. This leads to branch-dieback and eventual tree death, resulting in severe losses in production. Control strategies are limited and at present an integrated approach involving the use of phosphite, tolerant rootstocks, and proper nursery management has shown the best results. Disease progression of PRR is accelerated under high soil moisture or flooding conditions. In addition, avocado is highly susceptible to flooding, with even short periods of flooding causing significant losses. Despite the commercial importance of avocado, limited genomic resources are available. Next generation sequencing has provided the means to generate sequence data at a relatively low cost, making this an attractive option for non-model organisms such as avocado. The aims of this study were to generate sequence data for the avocado root transcriptome and identify stress-related genes. Tissue was isolated from avocado infected with P. cinnamomi, avocado exposed to flooding and avocado exposed to a combination of these two stresses. Three separate sequencing runs were performed on the Roche 454 platform and produced approximately 124 Mb of data. This was assembled into 7685 contigs, with 106 448 sequences remaining as singletons. Genes involved in defence pathways such as the salicylic acid and jasmonic acid pathways as well as genes associated with the response to low oxygen caused by flooding, were identified. This is the most comprehensive study of transcripts derived from root tissue of avocado to date and will provide a useful resource for future studies.

Examining sources of variation in HPG axis function among individuals and populations of the dark-eyed junco

Gonadal steroids are important mediators of traits relevant to fitness, and thus may be targets of selection. However, more knowledge is needed about sources of variation along the endocrine axes that may contribute to functional variation in steroid levels. In a controlled captive environment, we studied males of two closely related subspecies of the dark-eyed junco (Junco hyemalis) that differ in testosterone-related phenotype, asking whether they also differ in testosterone (T), and assessing the contribution of the sequential links of the hypothalamic-pituitary-gonadal axis. When males of both subspecies were challenged with gonadotropin-releasing hormone (GnRH), they were similar in circulating luteinizing hormone (LH) and T responses. When challenged with exogenous LH, they again produced levels of T similar to one another, and to the levels produced in response to GnRH. However, the smaller, less ornamented, and less aggressive subspecies had greater abundance of mRNA for LH receptor in the testes and for androgen receptor in the rostral hypothalamus, suggesting potential differences in regulatory feedback. We suggest that circulating hormone levels may be less prone to evolutionary change than the responsiveness of individual hormone targets. Among individuals, T titers were highly repeatable whether males were challenged with GnRH or with LH, but LH produced in response to GnRH did not covary with T produced in response to LH. Testis mass, but not LH receptor transcript abundance, predicted individual variation in T responses. These data implicate the gonad, but not the pituitary, as an important source of individual variation in T production.

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.

Sources of variation in HPG axis reactivity and individually consistent elevation of sex steroids in a female songbird

Understanding sources of individual differences in steroid hormone production has important implications for the evolution of reproductive and social behaviors. In females in particular, little is known about the mechanistic sources of these individual differences, despite established linkages between sex steroids and a variety of fitness-related traits. Using captive female dark-eyed juncos (Junco hyemalis) from two subspecies, we asked how variation in different components of the hypothalamo-pituitary-gonadal (HPG) axis related to variation in testosterone production among females, and we compared females to males in multiple components of the HPG axis. We demonstrated consistent individual differences in testosterone elevation in response to challenges with luteinizing hormone (LH) and gonadotropin-releasing hormone (GnRH). These hormone challenges led to more LH production but less testosterone production in females than males, and the sexes differed in some but not all measures of sensitivity to hormones along the HPG axis. Similar to findings in males, variation in testosterone production among females was not related to variation in LH production, gonadal LH-receptor mRNA abundance, or hypothalamic abundance of androgen receptor mRNA or aromatase mRNA. Rather, the primary source of individual variation in circulating steroids appears to the gonad, a conclusion further supported by positive correlations between testosterone and estradiol production. Unlike males, females did not differ by subspecies in any of the endocrine parameters that we assessed, suggesting some degree of independent evolution between the two sexes. Our results highlight the sources of physiological variation that may underlie the evolution of hormone-mediated phenotypes in females.

Potential for sexual conflict assessed via testosterone-mediated transcriptional changes in liver and muscle of a songbird

Males and females can be highly dimorphic in metabolism and physiology despite sharing nearly identical genomes, and both sexes respond phenotypically to elevated testosterone, a steroid hormone that alters gene expression. Only recently has it become possible to learn how a hormone such as testosterone affects global gene expression in non-model systems, and whether it affects the same genes in males and females. To investigate the transcriptional mechanisms by which testosterone exerts its metabolic and physiological effects on the periphery, we compared gene expression by sex and in response to experimentally elevated testosterone in a well-studied bird species, the dark-eyed junco (Junco hyemalis). We identified 291 genes in the liver and 658 in the pectoralis muscle that were differentially expressed between males and females. In addition, we identified 1727 genes that were differentially expressed between testosterone-treated and control individuals in at least one tissue and sex. Testosterone treatment altered the expression of only 128 genes in both males and females in the same tissue, and 847 genes were affected significantly differently by testosterone treatment in the two sexes. These substantial differences in transcriptional response to testosterone suggest that males and females may employ different pathways when responding to elevated testosterone, despite the fact that many phenotypic effects of experimentally elevated testosterone are similar in both sexes. In contrast, of the 121 genes that were affected by testosterone treatment in both sexes, 78% were regulated in the same direction (e.g. either higher or lower in testosterone-treated than control individuals) in both males and females. Thus, it appears that testosterone acts through both unique and shared transcriptional pathways in males and females, suggesting multiple mechanisms by which sexual conflict can be mediated.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.

Brain transcriptome of the violet-eared waxbill Uraeginthus granatina and recent evolution in the songbird genome

Songbirds are important models for the study of social behaviour and communication. To complement the recent genome sequencing of the domesticated zebra finch, we sequenced the brain transcriptome of a closely related songbird species, the violet-eared waxbill (Uraeginthus granatina). Both the zebra finch and violet-eared waxbill are members of the family Estrildidae, but differ markedly in their social behaviour. Using Roche 454 RNA sequencing, we generated an assembly and annotation of 11 084 waxbill orthologues of 17 475 zebra finch genes (64%), with an average transcript length of 1555 bp. We also identified 5985 single nucleotide polymorphisms (SNPs) of potential utility for future population genomic studies. Comparing the two species, we found evidence for rapid protein evolution (ω) and low polymorphism of the avian Z sex chromosome, consistent with prior studies of more divergent avian species. An intriguing outlier was putative chromosome 4A, which showed a high density of SNPs and low evolutionary rate relative to other chromosomes. Genome-wide ω was identical in zebra finch and violet-eared waxbill lineages, suggesting a similar demographic history with efficient purifying natural selection. Further comparisons of these and other estrildid finches may provide insights into the evolutionary neurogenomics of social behaviour.

Variation in candidate genes CLOCK and ADCYAP1 does not consistently predict differences in migratory behavior in the songbird genus Junco

Recent studies exploring the molecular genetic basis for migratory variation in animals have identified polymorphisms in two genes ( CLOCK and ADCYAP1) that are linked to circadian rhythms and correlate with migratory propensity and phenology among individuals and populations. Results from these initial studies are mixed, however, and additional data are needed to assess the generality and diversity of the molecular mechanisms that regulate the biology of migration. We sequenced CLOCK and ADCYAP1 in 15 populations across the two species of the avian genus Junco, a North American lineage in which multiple recently diverged subspecies and populations range from sedentary to long-distance migrants. We found no consistent associations between allele length and migratory status across the genus for either CLOCK or ADCYAP1. However, within two subspecies groups, populations that migrate longer distances have longer CLOCK alleles on average. Additionally, there was a positive relationship between ADCYAP1 allele length and migratory restlessness (zugunruhe) among individuals within one of two captive populations studied—a result similar to those reported previously within captive blackcaps ( Sylvia atricapilla). We conclude that, while both ADCYAP1 and CLOCK may correlate with migratory propensity within or among certain populations or species, previously identified relationships between migratory behavior and sequence variants cannot be easily generalized across taxa.

Testosterone affects neural gene expression differently in male and female juncos: a role for hormones in mediating sexual dimorphism and conflict

Despite sharing much of their genomes, males and females are often highly dimorphic, reflecting at least in part the resolution of sexual conflict in response to sexually antagonistic selection. Sexual dimorphism arises owing to sex differences in gene expression, and steroid hormones are often invoked as a proximate cause of sexual dimorphism. Experimental elevation of androgens can modify behavior, physiology, and gene expression, but knowledge of the role of hormones remains incomplete, including how the sexes differ in gene expression in response to hormones. We addressed these questions in a bird species with a long history of behavioral endocrinological and ecological study, the dark-eyed junco (Junco hyemalis), using a custom microarray. Focusing on two brain regions involved in sexually dimorphic behavior and regulation of hormone secretion, we identified 651 genes that differed in expression by sex in medial amygdala and 611 in hypothalamus. Additionally, we treated individuals of each sex with testosterone implants and identified many genes that may be related to previously identified phenotypic effects of testosterone treatment. Some of these genes relate to previously identified effects of testosterone-treatment and suggest that the multiple effects of testosterone may be mediated by modifying the expression of a small number of genes. Notably, testosterone-treatment tended to alter expression of different genes in each sex: only 4 of the 527 genes identified as significant in one sex or the other were significantly differentially expressed in both sexes. Hormonally regulated gene expression is a key mechanism underlying sexual dimorphism, and our study identifies specific genes that may mediate some of these processes.

Neural steroid sensitivity and aggression: comparing individuals of two songbird subspecies

Hormones coordinate the expression of complex phenotypes and thus may play important roles in evolutionary processes. When populations diverge in hormone-mediated phenotypes, differences may arise via changes in circulating hormones, sensitivity to hormones or both. Determining the relative importance of signal and sensitivity requires consideration of both inter- and intrapopulation variation in hormone levels, hormone sensitivity and phenotype, but such studies are rare, particularly among closely related taxa. We compared males of two subspecies of the dark-eyed junco (Junco hyemalis) for territorial aggression and associations among behaviour, circulating testosterone (T), and gene expression of androgen receptor (AR), aromatase (AROM) and oestrogen receptor α in three behaviourally relevant brain regions. Thus, we examined the degree to which evolution may shape behaviour via changes in plasma T as compared with key sex steroid binding/converting molecules. We found that the white-winged junco (J. h. aikeni) was more aggressive than the smaller, less ornamented Carolina junco (J. h. carolinensis). The subspecies did not differ in circulating testosterone, but did differ significantly in the abundance of AR and AROM mRNA in key areas of the brain. Within populations, both gene expression and circulating T co-varied significantly with individual differences in aggression. Notably, the differences identified between populations were opposite to those predicted by the patterns among individuals within populations. These findings suggest that hormone-phenotype relationships may evolve via multiple pathways, and that changes that have occurred over evolutionary time do not necessarily reflect standing physiological variation on which current evolutionary processes may act.

Why so many unknown genes? Partitioning orphans from a representative transcriptome of the lone star tick Amblyomma americanum

Background

Genomic resources within the phylum Arthropoda are largely limited to the true insects but are beginning to include unexplored subphyla, such as the Crustacea and Chelicerata. Investigations of these understudied taxa uncover high frequencies of orphan genes, which lack detectable sequence homology to genes in pre-existing databases. The ticks (Acari: Chelicerata) are one such understudied taxon for which genomic resources are urgently needed. Ticks are obligate blood-feeders that vector major diseases of humans, domesticated animals, and wildlife. In analyzing a transcriptome of the lone star tick Amblyomma americanum, one of the most abundant disease vectors in the United States, we find a high representation of unannotated sequences. We apply a general framework for quantifying the origin and true representation of unannotated sequences in a dataset and for evaluating the biological significance of orphan genes.

Results

Expressed sequence tags (ESTs) were derived from different life stages and populations of A. americanum and combined with ESTs available from GenBank to produce 14,310 ESTs, over twice the number previously available. The vast majority (71%) has no sequence homology to proteins archived in UniProtKB. We show that poor sequence or assembly quality is not a major contributor to this high representation by orphan genes. Moreover, most unannotated sequences are functional: a microarray experiment demonstrates that 59% of functional ESTs are unannotated. Lastly, we attempt to further annotate our EST dataset using genomic datasets from other members of the Acari, including Ixodes scapularis, four other tick species and the mite Tetranychus urticae. We find low homology with these species, consistent with significant divergence within this subclass.

Conclusions

We conclude that the abundance of orphan genes in A. americanum likely results from 1) taxonomic isolation stemming from divergence within the tick lineage and limited genomic resources for ticks and 2) lineage-specific genes needing functional genomic studies to evaluate their association with the unique biology of ticks. The EST sequences described here will contribute substantially to the development of tick genomics. Moreover, the framework provided for the evaluation of orphan genes can guide analyses of future transcriptome sequencing projects.

Neural sensitivity to sex steroids predicts individual differences in aggression: implications for behavioural evolution

Testosterone (T) regulates many traits related to fitness, including aggression. However, individual variation in aggressiveness does not always relate to circulating T, suggesting that behavioural variation may be more closely related to neural sensitivity to steroids, though this issue remains unresolved. To assess the relative importance of circulating T and neural steroid sensitivity in predicting behaviour, we measured aggressiveness during staged intrusions in free-living male and female dark-eyed juncos (Junco hyemalis). We compared aggressiveness to plasma T levels and to the abundance of androgen receptor (AR), aromatase (AROM) and oestrogen receptor alpha (ORα) mRNA in behaviourally relevant brain areas (avian medial amygdala, hypothalamus and song control regions). We also asked whether patterns of covariation among behaviour and endocrine parameters differed in males and females, anticipating that circulating T may be a better predictor of behaviour in males than in females. We found that circulating T related to aggressiveness only in males, but that gene expression for ORα, AR and AROM covaried with individual differences in aggressiveness in both sexes. These findings are among the first to show that individual variation in neural gene expression for three major sex steroid-processing molecules predicts individual variation in aggressiveness in both sexes in nature. The results have broad implications for our understanding of the mechanisms by which aggressive behaviour may evolve.