Digital gene expression analysis of the zebra finch genome

Peripheral inflammation-induced changes in songbird brain gene expression: 3' mRNA transcriptomic approach

Species-specific neural inflammation can be induced by profound immune signalling from periphery to brain. Recent advances in transcriptomics offer cost-effective approaches to study this regulation. In a population of captive zebra finch (Taeniopygia guttata), we compare the differential gene expression patterns in lipopolysaccharide (LPS)-triggered peripheral inflammation revealed by RNA-seq and QuantSeq. The RNA-seq approach identified more differentially expressed genes but failed to detect any inflammatory markers. In contrast, QuantSeq results identified specific expression changes in the genes regulating inflammation. Next, we adopted QuantSeq to relate peripheral and brain transcriptomes. We identified subtle changes in the brain gene expression during the peripheral inflammation (e.g. up-regulation in AVD-like and ACOD1 expression) and detected co-structure between the peripheral and brain inflammation. Our results suggest benefits of the 3'end transcriptomics for association studies between peripheral and neural inflammation in genetically heterogeneous models and identify potential targets for the future brain research in birds.

Varied impacts of social relationships on neuroendocrine state

Social relationships, affiliative social attachments, are important for many species. The best studied types of relationships are monogamous pair bonds. However, it remains unclear how generalizable models of pair bonding are across types of social attachments. Zebra finches are a fascinating system to explore the neurobiology of social relationships because they form various adult bonds with both same- and opposite-sex partners. To test whether different bonds are supported by a single brain network, we quantified individuals' neuroendocrine state after either 24 h or 2 weeks of co-housing with a novel same- or opposite-sex partner. We defined neuroendocrine state by the expression of 22 genes related to 4 major signaling pathways (dopamine, steroid, nonapeptide, and opioid) in six brain regions associated with affiliation or communication [nucleus accumbens (NAc), nucleus taeniae of the amygdala (TnA), medial preoptic area (POM), and periaqueductal gray (PAG), ventral tegmental area, and auditory cortex]. Overall, we found dissociable effects of social contexts (same- or opposite-sex partnerships) and duration of co-housing. Social bonding impacted the neuroendocrine state of four regions in males (NAc, TnA, POM, and PAG) and three regions in females (NAc, TnA, and POM). Monogamous pair bonding specifically appeared to impact male NAc. However, the patterns of gene expression in zebra finches were different than has previously been reported in mammals. Together, our results support the view that there are numerous mechanisms regulating social relationships and highlight the need to further our understanding of how social interactions shape social bonds.

High-throughput sequencing analysis identified microRNAs associated with egg production in ducks ovaries

BACKGROUND

MicroRNAs (miRNAs) exist widely and are involved in multiple biological processes in ducks, whereas the regulatory mechanism of miRNAs in egg laying of ducks has remained unclear. This study aims to reveal key miRNAs involved in the regulation of egg production in duck ovaries.

METHODS

High-throughput sequencing was performed on four egg-type duck ovaries and four egg-meat-type duck ovaries at the start of the egg-laying stage. Quantitative reverse transcription PCR (qRT-PCR) validation was performed on differentially expressed miRNAs (DE miRNAs). Gene network of DEmiRNA-mRNA-pathway was constructed by Cytoscape.

RESULTS

A total of 251 know miRNAs and 1,972 novel miRNAs were obtained from whole clean reads. Among the known miRNAs, we identified 21 DEmiRNAs, including eight down-regulated and 13 up-regulated miRNAs in egg-type ducks compared with egg-meat-type ducks. Among the novel miRNAs, we identified 70 DEmiRNAs, including 58 down-regulated and 12 up-regulated in egg-type ducks compared with egg-meat-type ducks. The expression patterns of four miRNAs were verified by qRT-PCR. The DEmiRNAs were involved in the function of response to folic acid and the pathway of valine, leucine and isoleucine degradation. Specific target genes of DEmiRNAs enrichment was found in some egg-laying regulation pathways, such as dopaminergic synapse, ovarian steroidogenesis and oocyte meiosis. The DEmiRNA-mRNA-pathway network including three DEmiRNAs, nine mRNAs and 11 pathways. apl-miR-194-5p and apl-miR-215-5p may be potential key miRNAs in regulating egg laying.

CONCLUSIONS

This study provided miRNAs profiles in ducks about egg laying and establish a theoretical basis for subsequent selection or modification of duck phenotypes at the molecular level.

Tissue-specific expression profiles and positive selection analysis in the tree swallow (Tachycineta bicolor) using a de novo transcriptome assembly

Tree swallows (Tachycineta bicolor) are one of the most commonly studied wild birds in North America. They have advanced numerous research areas, including life history, physiology, and organismal responses to global change; however, transcriptomic resources are scarce. To further advance the utility of this system for biologists across disciplines, we generated a transcriptome for the tree swallow using six tissues (brain, blood, ovary, spleen, liver, and muscle) collected from breeding females. We de novo assembled 207,739 transcripts, which we aligned to 14,717 high confidence protein-coding genes. We then characterized each tissue with regard to its unique genes and processes and applied this transcriptome to two fundamental questions in evolutionary biology and endocrinology. First, we analyzed 3,015 single-copy orthologs and identified 46 genes under positive selection in the tree swallow lineage, including those with putative links to adaptations in this species. Second, we analyzed tissue-specific expression patterns of genes involved in sex steroidogenesis and processing. Enzymes capable of synthesizing these behaviorally relevant hormones were largely limited to the ovary, whereas steroid binding genes were found in nearly all other tissues, highlighting the potential for local regulation of sex steroid-mediated traits. These analyses provide new insights into potential sources of phenotypic variation in a free-living female bird and advance our understanding of fundamental questions in evolutionary and organismal biology.

Research Resource: Hormones, Genes, and Athleticism: Effect of Androgens on the Avian Muscular Transcriptome

Male vertebrate social displays vary from physically simple to complex, with the latter involving exquisite motor command of the body and appendages. Studies of these displays have, in turn, provided substantial insight into neuromotor mechanisms. The neotropical golden-collared manakin (Manacus vitellinus) has been used previously as a model to investigate intricate motor skills because adult males of this species perform an acrobatic and androgen-dependent courtship display. To support this behavior, these birds express elevated levels of androgen receptors (AR) in their skeletal muscles. Here we use RNA sequencing to explore how testosterone (T) modulates the muscular transcriptome to support male manakin courtship displays. In addition, we explore how androgens influence gene expression in the muscles of the zebra finch (Taenopygia guttata), a model passerine bird with a limited courtship display and minimal muscle AR. We identify androgen-dependent, muscle-specific gene regulation in both species. In addition, we identify manakin-specific effects that are linked to muscle use during the manakin display, including androgenic regulation of genes associated with muscle fiber contractility, cellular homeostasis, and energetic efficiency. Overall, our results point to numerous genes and gene networks impacted by androgens in male birds, including some that underlie optimal muscle function necessary for performing acrobatic display routines. Manakins are excellent models to explore gene regulation promoting athletic ability.

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.

Highly variable recombinational landscape modulates efficacy of natural selection in birds

Determining the rate of protein evolution and identifying the causes of its variation across the genome are powerful ways to understand forces that are important for genome evolution. By using a multitissue transcriptome data set from great tit (Parus major), we analyzed patterns of molecular evolution between two passerine birds, great tit and zebra finch (Taeniopygia guttata), using the chicken genome (Gallus gallus) as an outgroup. We investigated whether a special feature of avian genomes, the highly variable recombinational landscape, modulates the efficacy of natural selection through the effects of Hill-Robertson interference, which predicts that selection should be more effective in removing deleterious mutations and incorporating beneficial mutations in high-recombination regions than in low-recombination regions. In agreement with these predictions, genes located in low-recombination regions tend to have a high proportion of neutrally evolving sites and relaxed selective constraint on sites subject to purifying selection, whereas genes that show strong support for past episodes of positive selection appear disproportionally in high-recombination regions. There is also evidence that genes located in high-recombination regions tend to have higher gene expression specificity than those located in low-recombination regions. Furthermore, more compact genes (i.e., those with fewer/shorter introns or shorter proteins) evolve faster than less compact ones. In sum, our results demonstrate that transcriptome sequencing is a powerful method to answer fundamental questions about genome evolution in nonmodel organisms.

Using genomics to characterize evolutionary potential for conservation of wild populations

Genomics promises exciting advances towards the important conservation goal of maximizing evolutionary potential, notwithstanding associated challenges. Here, we explore some of the complexity of adaptation genetics and discuss the strengths and limitations of genomics as a tool for characterizing evolutionary potential in the context of conservation management. Many traits are polygenic and can be strongly influenced by minor differences in regulatory networks and by epigenetic variation not visible in DNA sequence. Much of this critical complexity is difficult to detect using methods commonly used to identify adaptive variation, and this needs appropriate consideration when planning genomic screens, and when basing management decisions on genomic data. When the genomic basis of adaptation and future threats are well understood, it may be appropriate to focus management on particular adaptive traits. For more typical conservations scenarios, we argue that screening genome-wide variation should be a sensible approach that may provide a generalized measure of evolutionary potential that accounts for the contributions of small-effect loci and cryptic variation and is robust to uncertainty about future change and required adaptive response(s). The best conservation outcomes should be achieved when genomic estimates of evolutionary potential are used within an adaptive management framework.

Identification of differentially expressed genes in breast muscle and skin fat of postnatal Pekin duck

Lean-type Pekin duck is a commercial breed that has been obtained through long-term selection. Investigation of the differentially expressed genes in breast muscle and skin fat at different developmental stages will contribute to a comprehensive understanding of the potential mechanisms underlying the lean-type Pekin duck phenotype. In the present study, RNA-seq was performed on breast muscle and skin fat at 2-, 4- and 6-weeks of age. More than 89% of the annotated duck genes were covered by our RNA-seq dataset. Thousands of differentially expressed genes, including many important genes involved in the regulation of muscle development and fat deposition, were detected through comparison of the expression levels in the muscle and skin fat of the same time point, or the same tissue at different time points. KEGG pathway analysis showed that the differentially expressed genes clustered significantly in many muscle development and fat deposition related pathways such as MAPK signaling pathway, PPAR signaling pathway, Calcium signaling pathway, Fat digestion and absorption, and TGF-beta signaling pathway. The results presented here could provide a basis for further investigation of the mechanisms involved in muscle development and fat deposition in Pekin duck.

Brain transcriptome sequencing and assembly of three songbird model systems for the study of social behavior

Emberizid sparrows (emberizidae) have played a prominent role in the study of avian vocal communication and social behavior. We present here brain transcriptomes for three emberizid model systems, song sparrow Melospiza melodia, white-throated sparrow Zonotrichia albicollis, and Gambel's white-crowned sparrow Zonotrichia leucophrys gambelii. Each of the assemblies covered fully or in part, over 89% of the previously annotated protein coding genes in the zebra finch Taeniopygia guttata, with 16,846, 15,805, and 16,646 unique BLAST hits in song, white-throated and white-crowned sparrows, respectively. As in previous studies, we find tissue of origin (auditory forebrain versus hypothalamus and whole brain) as an important determinant of overall expression profile. We also demonstrate the successful isolation of RNA and RNA-sequencing from post-mortem samples from building strikes and suggest that such an approach could be useful when traditional sampling opportunities are limited. These transcriptomes will be an important resource for the study of social behavior in birds and for data driven annotation of forthcoming whole genome sequences for these and other bird species.

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.

EST sequencing and gene expression profiling in Scutellaria baicalensis

Scutellaria baicalensis is an important medicinal plant, but few genomic resources are available for this species, as well as for other non-model plants. One of the major new directions in genome research is to discover the full spectrum of genes transcribed from the whole genome. Here, we report extensive transcriptome data of the early growth stage of S. baicalensis. This transcriptome consensus sequence was constructed by de novo assembly of shotgun sequencing data, obtained using multiple next-generation DNA sequencing (NGS) platforms (Roche/454 GS_FLX+ and Illumina/Solexa HiSeq2000). We show that this new approach to obtain extensive mRNA is an efficient strategy for genome-wide transcriptome analysis. We obtained 1,226,938 and 161,417,646 reads using the GS_FLX and the Illumina/Solexa HiSeq2000, respectively. De novo assembly of the high-quality GS_FLX and Illumina reads (95 % and 75 %) resulted in more than 82 Mb of mRNA consensus sequence, which we assembled into 51,188 contigs, with at least 500 bp per contig. Of these contigs, 39,581 contained known genes, as determined by BLASTX searches against non-redundant NCBI database. Of these, 20,498 different genes were expressed during the early growth stage of S. baicalensis. We have made the expressed sequences available on a public database. Our results demonstrate the utility of combining NGS technologies as a basis for the development of genomic tools in non-model, medicinal plant species. Knowledge of all described genes and quantitation of the expressed genes, including the transcription factors involved, will be useful in studies of the biology of S. baicalensis gene regulation.

Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection

Hemileia vastatrix is the causal agent of coffee leaf rust, the most important disease of coffee Arabica. In this work, a 454-pyrosequencing transcriptome analysis of H. vastatrix germinating urediniospores (gU) and appressoria (Ap) was performed and compared to previously published in planta haustoria-rich (H) data. A total of 9234 transcripts were identified and annotated. Ca. 50% of these transcripts showed no significant homology to international databases. Only 784 sequences were shared by the three conditions, and 75% were exclusive of either gU (2146), Ap (1479) or H (3270). Relative transcript abundance and RT-qPCR analyses for a selection of genes indicated a particularly active metabolism, translational activity and production of new structures in the appressoria and intense signaling, transport, secretory activity and cellular multiplication in the germinating urediniospores, suggesting the onset of a plant-fungus dialogue as early as at the germ tube stage. Gene expression related to the production of carbohydrate-active enzymes and accumulation of glycerol in germinating urediniospores and appressoria suggests that combined lytic and physical mechanisms are involved in appressoria-mediated penetration. Besides contributing to the characterization of molecular processes leading to appressoria-mediated infection by rust fungi, these results point toward the identification of new H. vastatrix candidate virulence factors, with 516 genes predicted to encode secreted proteins.

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.

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.

Expression profiles for genes in the turkey major histocompatibility complex B-locus

The major histocompatibility complex (MHC) is a highly polymorphic region of the genome essential to immune responses and animal health. In galliforms, the MHC is divided into 2 genetically unlinked regions (MHC-B and MHC-Y). Many MHC-B genes are involved in adaptive or innate immunity, yet others have nonimmune or unknown functions. The sequenced MHC-B region of the turkey (Meleagris gallopavo) contains 40 genes, the majority of which are predicted transcripts based on comparison with the chicken or quail, without direct evidence for expression. This study was designed to test for the presence of MHC-B gene transcripts in a panel of immune and nonimmune system tissues from domestic turkeys. This analysis provides the first locus-wide examination of MHC-B gene expression in any avian species. Most MHC-B genes were broadly expressed across tissues. Expression of all predicted genes was verified by reverse-transcription PCR, including B-butyrophilin 2 (BTN2), a predicted gene with no previous evidence for expression in any species. Previously undescribed splice variants were also detected and sequenced from 3 genes. Characterization of MHC-B expression patterns helps elucidate unknown gene functions and potential gene coregulation. Determining turkey MHC-B expression profiles increases our overall understanding of the avian MHC and provides a necessary resource for future research on the immunological response of these genes.

Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus

BACKGROUND

Animal migration requires adaptations in morphological, physiological and behavioural traits. Several of these traits have been shown to possess a strong heritable component in birds, but little is known about their genetic architecture. Here we used 454 sequencing of brain-derived transcriptomes from two differentially migrating subspecies of the willow warbler Phylloscopus trochilus to detect genes potentially underlying traits associated with migration.

RESULTS

The transcriptome sequencing resulted in 1.8 million reads following filtering steps. Most of the reads (84%) were successfully mapped to the genome of the zebra finch Taeniopygia gutatta. The mapped reads were situated within at least 12,101 predicted zebra finch genes, with the greatest sequencing depth in exons. Reads that were mapped to intergenic regions were generally located close to predicted genes and possibly located in uncharacterized untranslated regions (UTRs). Out of 85,000 single nucleotide polymorphisms (SNPs) with a minimum sequencing depth of eight reads from each of two subspecies-specific pools, only 55 showed high differentiation, confirming previous studies showing that most of the genetic variation is shared between the subspecies. Validation of a subset of the most highly differentiated SNPs using Sanger sequencing demonstrated that several of them also were differentiated between an independent set of individuals of each subspecies. These SNPs were clustered in two chromosome regions that are likely to be influenced by divergent selection between the subspecies and that could potentially be associated with adaptations to their different migratory strategies.

CONCLUSIONS

Our study represents the first large-scale sequencing analysis aiming at detecting genes underlying migratory phenotypes in birds and provides new candidates for genes potentially involved in migration.

Transcriptome of the dead: characterisation of immune genes and marker development from necropsy samples in a free-ranging marine mammal

Background

Transcriptomes are powerful resources, providing a window on the expressed portion of the genome that can be generated rapidly and at low cost for virtually any organism. However, because many genes have tissue-specific expression patterns, developing a complete transcriptome usually requires a 'discovery pool' of individuals to be sacrificed in order to harvest mRNA from as many different types of tissue as possible. This hinders transcriptome development in large, charismatic and endangered species, many of which stand the most to gain from such approaches. To circumvent this problem in a model pinniped species, we 454 sequenced cDNA from testis, heart, spleen, intestine, kidney and lung tissues obtained from nine adult male Antarctic fur seals (Arctocephalus gazella) that died of natural causes at Bird Island, South Georgia.

Results

After applying stringent quality control criteria based on length and annotation, we obtained 12,397 contigs which, in combination with 454 data previously obtained from skin, gave a total of 23,096 unique contigs. Homology was found to 77.0% of dog (Canis lupus familiaris) transcripts, suggesting that the combined assembly represents a substantial proportion of this species' transcriptome. Moreover, only 0.5% of transcripts revealed sequence similarity to bacteria, implying minimal contamination, and the percentage of transcripts involved in cell death was low at 2.6%. Transcripts with immune-related annotations were almost five-fold enriched relative to skin and represented 13.2% of all spleen-specific contigs. By reference to the dog, we also identified transcripts revealing homology to five class I, ten class II and three class III genes of the Major Histocompatibility Complex and derived the putative genomic distribution of 17,121 contigs, 2,119 in silico mined microsatellites and 9,382 single nucleotide polymorphisms.

Conclusions

Our findings suggest that transcriptome development based on samples collected post mortem may greatly facilitate genomic studies, not only of marine mammals but also more generally of species that are of conservation concern.

Conservation genomics of threatened animal species

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.

RNA-seq transcriptome analysis of male and female zebra finch cell lines

The derivation of stably cultured cell lines has been critical to the advance of molecular biology. We profiled gene expression in the first two generally available cell lines derived from the zebra finch. Using Illumina RNA-seq, we generated ~93 million reads and mapped the majority to the recently assembled zebra finch genome. Expression of most Ensembl-annotated genes was detected, but over half of the mapped reads aligned outside annotated genes. The male-derived G266 line expressed Z-linked genes at a higher level than did the female-derived ZFTMA line, indicating persistence in culture of the distinctive lack of avian sex chromosome dosage compensation. Although these cell lines were not derived from neural tissue, many neurobiologically relevant genes were expressed, although typically at lower levels than in a reference sample from auditory forebrain. These cell lines recapitulate fundamental songbird biology and will be useful for future studies of songbird gene regulation and function.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

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

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.

Transcriptome analysis of a North American songbird, Melospiza melodia

An effective way to understand the genomics of divergence in non-model organisms is to use the transcriptome to identify genes associated with divergence. We examine the transcriptome of the song sparrow (Melospiza melodia) and contrast it with the avian models zebra finch (Taeniopygia guttata) and chicken (Gallus gallus). We aimed to (i) obtain a functional annotation of a substantial portion of the song sparrow transcriptome; (ii) compare transcript divergence; (iii) efficiently characterize single nucleotide polymorphism/indel markers possibly fixed between song sparrow subspecies; and (iv) identify the most common set of transcripts in birds using the zebra finch as a reference. Using two individuals from each of three populations, whole-body mRNA was normalized and sequenced (110Mb total). The assembly yielded 38 539 contigs [N50 (the length–weighted median) = 482 bp]; 4574 were orthologous to both model genomes and 3680 are functionally annotated. This low-coverage scan of the song sparrow transcriptome revealed 29 982 SNPs/indels, 1402 fixed between populations and subspecies. Referencing zebra finch and chicken, we identified 43 and 5 fast-evolving genes, respectively. We also identified the most common set of transcripts present in birds with respect to zebra finch. This study provides new insight into songbird transcriptomes, and candidate markers identified here may help research in songbirds (oscine Passeriformes), a frequently studied group.

Digital gene expression analysis based on integrated de novo transcriptome assembly of sweet potato [Ipomoea batatas (L.) Lam]

BACKGROUND

Sweet potato (Ipomoea batatas L. [Lam.]) ranks among the top six most important food crops in the world. It is widely grown throughout the world with high and stable yield, strong adaptability, rich nutrient content, and multiple uses. However, little is known about the molecular biology of this important non-model organism due to lack of genomic resources. Hence, studies based on high-throughput sequencing technologies are needed to get a comprehensive and integrated genomic resource and better understanding of gene expression patterns in different tissues and at various developmental stages.

METHODOLOGY/PRINCIPAL FINDINGS

Illumina paired-end (PE) RNA-Sequencing was performed, and generated 48.7 million of 75 bp PE reads. These reads were de novo assembled into 128,052 transcripts (≥ 100 bp), which correspond to 41.1 million base pairs, by using a combined assembly strategy. Transcripts were annotated by Blast2GO and 51,763 transcripts got BLASTX hits, in which 39,677 transcripts have GO terms and 14,117 have ECs that are associated with 147 KEGG pathways. Furthermore, transcriptome differences of seven tissues were analyzed by using Illumina digital gene expression (DGE) tag profiling and numerous differentially and specifically expressed transcripts were identified. Moreover, the expression characteristics of genes involved in viral genomes, starch metabolism and potential stress tolerance and insect resistance were also identified.

CONCLUSIONS/SIGNIFICANCE

The combined de novo transcriptome assembly strategy can be applied to other organisms whose reference genomes are not available. The data provided here represent the most comprehensive and integrated genomic resources for cloning and identifying genes of interest in sweet potato. Characterization of sweet potato transcriptome provides an effective tool for better understanding the molecular mechanisms of cellular processes including development of leaves and storage roots, tissue-specific gene expression, potential biotic and abiotic stress response in sweet potato.

Transcriptome sequencing of black grouse (Tetrao tetrix) for immune gene discovery and microsatellite development

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 flanking primers for 526 of these. A preliminary test of the polymorphism of the microsatellites found 10 polymorphic microsatellites of the 102 tested. Genomic resources generated in this study should greatly benefit future ecological, evolutionary and conservation genetic studies on this species.

Survey of the Applications of NGS to Whole-Genome Sequencing and Expression Profiling

Recently, the technologies of DNA sequence variation and gene expression profiling have been used widely as approaches in the expertise of genome biology and genetics. The application to genome study has been particularly developed with the introduction of the next-generation DNA sequencer (NGS) Roche/454 and Illumina/Solexa systems, along with bioinformation analysis technologies of whole-genome de novo assembly, expression profiling, DNA variation discovery, and genotyping. Both massive whole-genome shotgun paired-end sequencing and mate paired-end sequencing data are important steps for constructing de novo assembly of novel genome sequencing data. It is necessary to have DNA sequence information from a multiplatform NGS with at least 2× and 30× depth sequence of genome coverage using Roche/454 and Illumina/Solexa, respectively, for effective an way of de novo assembly. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. Whole-genome expression profile data are useful to approach genome system biology with quantification of expressed RNAs from a whole-genome transcriptome, depending on the tissue samples. The hybrid mRNA sequences from Rohce/454 and Illumina/Solexa are more powerful to find novel genes through de novo assembly in any whole-genome sequenced species. The 20× and 50× coverage of the estimated transcriptome sequences using Roche/454 and Illumina/Solexa, respectively, is effective to create novel expressed reference sequences. However, only an average 30× coverage of a transcriptome with short read sequences of Illumina/Solexa is enough to check expression quantification, compared to the reference expressed sequence tag sequence.

Comparison between Normalised and Unnormalised 454-Sequencing Libraries for Small-Scale RNA-Seq Studies

Next-generation sequencing of transcriptomes (RNA-Seq) is being used increasingly in studies of nonmodel organisms. Here, we evaluate the effectiveness of normalising cDNA libraries prior to sequencing in a small-scale study of the zebra finch. We find that assemblies produced from normalised libraries had a larger number of contigs but used fewer reads compared to unnormalised libraries. Considerably more genes were also detected using the contigs produced from normalised cDNA, and microsatellite discovery was up to 73% more efficient in these. There was a positive correlation between the detected expression level of genes in normalised and unnormalised cDNA, and there was no difference in the number of genes identified as being differentially expressed between blood and spleen for the normalised and unnormalised libraries. We conclude that normalised cDNA libraries are preferable for many applications of RNA-Seq and that these can also be used in quantitative gene expression studies.

Genomic toolboxes for conservation biologists

Conservation genetics is expanding its research horizon with a genomic approach, by incorporating the modern techniques of next-generation sequencing (NGS). Application of NGS overcomes many limitations of conservation genetics. First, NGS allows for genome-wide screening of markers, which may lead to a more representative estimation of genetic variation within and between populations. Second, NGS allows for distinction between neutral and non-neutral markers. By screening populations on thousands of single nucleotide polymorphism markers, signals of selection can be found for some markers. Variation in these markers will give insight into functional rather than neutral genetic variation. Third, NGS facilitates the study of gene expression. Conservation genomics will increase our insight in how the environment and genes interact to affect phenotype and fitness. In addition, the NGS approach opens a way to study processes such as inbreeding depression and local adaptation mechanistically. Conservation genetics programs are directed to a fundamental understanding of the processes involved in conservation genetics and should preferably be started in species for which large databases on ecology, demography and genetics are available. Here, we describe and illustrate the connection between the application of NGS technologies and the research questions in conservation. The perspectives of conservation genomics programs are also discussed.

Overview of the transcriptome profiles identified in hagfish, shark, and bichir: current issues arising from some nonmodel vertebrate taxa

Because of their crucial phylogenetic positions, hagfishes, sharks, and bichirs are recognized as key taxa in our understanding of vertebrate evolution. The expression patterns of the regulatory genes involved in developmental patterning have been analyzed in the context of evolutionary developmental studies. However, in a survey of public sequence databases, we found that the large-scale sequence data for these taxa are still limited. To address this deficit, we used conventional Sanger DNA sequencing and a next-generation sequencing technology based on 454 GS FLX sequencing to obtain expressed sequence tags (ESTs) of the Japanese inshore hagfish (Eptatretus burgeri; 161,482 ESTs), cloudy catshark (Scyliorhinus torazame; 165,819 ESTs), and gray bichir (Polypterus senegalus; 34,336 ESTs). We deposited the ESTs in a newly constructed database, designated the "Vertebrate TimeCapsule." The ESTs include sequences from genes that can be effectively used in evolutionary developmental studies; for instance, several encode cartilaginous extracellular matrix proteins, which are central to an understanding of the ways in which evolutionary processes affected the skeletal elements, whereas others encode regulatory genes involved in craniofacial development and early embryogenesis. Here, we discuss how hagfishes, sharks, and bichirs contribute to our understanding of vertebrate evolution, we review the current status of the publicly available sequence data for these three taxa, and we introduce our EST projects and newly developed database.

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

Genetic mapping of the major histocompatibility complex in the zebra finch (Taeniopygia guttata)

Genes of the major histocompatibility complex (MHC) have received much attention in immunology, genetics, and ecology because they are highly polymorphic and play important roles in parasite resistance and mate choice. Until recently, the MHC of passerine birds was not well-described. However, the genome sequencing of the zebra finch (Taeniopygia guttata) has partially redressed this gap in our knowledge of avian MHC genes. Here, we contribute further to the understanding of the zebra finch MHC organization by mapping SNPs within or close to known MHC genes in the zebra finch genome. MHC class I and IIB genes were both mapped to zebra finch chromosome 16, and there was no evidence that MHC class I genes are located on chromosome 22 (as suggested by the genome assembly). We confirm the location in the MHC region on chromosome 16 for several other genes (BRD2, FLOT1, TRIM7.2, GNB2L1, and CSNK2B). Two of these (CSNK2B and FLOT1) have not previously been mapped in any other bird species. In line with previous results, we also find that orthologs to the immune-related genes B-NK and CLEC2D, which are part of the MHC region in chicken, are situated on zebra finch chromosome Z and not among other MHC genes in the zebra finch.

Next generation sequencing and analysis of a conserved transcriptome of New Zealand's kiwi

BACKGROUND

Kiwi is a highly distinctive, flightless and endangered ratite bird endemic to New Zealand. To understand the patterns of molecular evolution of the nuclear protein-coding genes in brown kiwi (Apteryx australis mantelli) and to determine the timescale of avian history we sequenced a transcriptome obtained from a kiwi embryo using next generation sequencing methods. We then assembled the conserved protein-coding regions using the chicken proteome as a scaffold.

RESULTS

Using 1,543 conserved protein coding genes we estimated the neutral evolutionary divergence between the kiwi and chicken to be ~45%, which is approximately equal to the divergence computed for the human-mouse pair using the same set of genes. A large fraction of genes was found to be under high selective constraint, as most of the expressed genes appeared to be involved in developmental gene regulation. Our study suggests a significant relationship between gene expression levels and protein evolution. Using sequences from over 700 nuclear genes we estimated the divergence between the two basal avian groups, Palaeognathae and Neognathae to be 132 million years, which is consistent with previous studies using mitochondrial genes.

CONCLUSIONS

The results of this investigation revealed patterns of mutation and purifying selection in conserved protein coding regions in birds. Furthermore this study suggests a relatively cost-effective way of obtaining a glimpse into the fundamental molecular evolutionary attributes of a genome, particularly when no closely related genomic sequence is available.

Evolutionary analysis and expression profiling of zebra finch immune genes

Genes of the immune system are generally considered to evolve rapidly due to host-parasite coevolution. They are therefore of great interest in evolutionary biology and molecular ecology. In this study, we manually annotated 144 avian immune genes from the zebra finch (Taeniopygia guttata) genome and conducted evolutionary analyses of these by comparing them with their orthologs in the chicken (Gallus gallus). Genes classified as immune receptors showed elevated d(N)/d(S) ratios compared with other classes of immune genes. Immune genes in general also appear to be evolving more rapidly than other genes, as inferred from a higher d(N)/d(S) ratio compared with the rest of the genome. Furthermore, ten genes (of 27) for which sequence data were available from at least three bird species showed evidence of positive selection acting on specific codons. From transcriptome data of eight different tissues, we found evidence for expression of 106 of the studied immune genes, with primary expression of most of these in bursa, blood, and spleen. These immune-related genes showed a more tissue-specific expression pattern than other genes in the zebra finch genome. Several of the avian immune genes investigated here provide strong candidates for in-depth studies of molecular adaptation in birds.

Ecological genetics and genomics of plant defenses: Evidence and approaches

Herbivores exert significant selection on plants, and plants have evolved a variety of constitutive and inducible defenses to resist and tolerate herbivory. Assessing the genetic mechanisms that influence defenses against herbivores will deepen our understanding of the evolution of essential phenotypic traits.Ecogenomics is a powerful interdisciplinary approach that can address fundamental questions about the ecology and evolutionary biology of species, such as: which evolutionary forces maintain variation within a population? and What is the genetic architecture of adaptation? This field seeks to identify gene regions that influence ecologically-important traits, assess the fitness consequences under natural conditions of alleles at key quantitative trait loci (QTLs), and test how the abiotic and biotic environment affects gene expression.Here, we review ecogenomics techniques and emphasize how this framework can address long-standing and emerging questions relating to anti-herbivore defenses in plants. For example, ecogenomics tools can be used to investigate: inducible vs. constitutive defenses; tradeoffs between resistance and tolerance; adaptation to the local herbivore community; selection on alleles that confer resistance and tolerance in natural populations; and whether different genes are activated in response to specialist vs. generalist herbivores and to different types of damage.Ecogenomic studies can be conducted with model species, such as Arabidopsis, or their relatives, in which case myriad molecular tools are already available. Burgeoning sequence data will also facilitate ecogenomic studies of non-model species. Throughout this paper, we highlight approaches that are particularly suitable for ecological studies of non-model organisms, discuss the benefits and disadvantages of specific techniques, and review bioinformatic tools for analyzing data.We focus on established and promising techniques, such as QTL mapping with pedigreed populations, genome wide association studies, transcription profiling strategies, population genomics, and transgenic methodologies. Many of these techniques are complementary and can be used jointly to investigate the genetic architecture of defense traits and selection on alleles in nature.

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

BACKGROUND

Due to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.

RESULTS

The zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.

CONCLUSION

The zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.