Conserved syntenic clusters of protein coding genes are missing in birds

Curation of microarray oligonucleotides and corresponding ESTs/cDNAs used for gene expression analysis in zebra finches

OBJECTIVES

Zebra finches are a major model organism for investigating mechanisms of vocal learning, a trait that enables spoken language in humans. The development of cDNA collections with expressed sequence tags (ESTs) and microarrays has allowed for extensive molecular characterizations of circuitry underlying vocal learning and production. However, poor database curation can lead to errors in transcriptome and bioinformatics analyses, limiting the impact of these resources. Here we used genomic alignments and synteny analysis for orthology verification to curate and reannotate ~ 35% of the oligonucleotides and corresponding ESTs/cDNAs that make-up Agilent microarrays for gene expression analysis in finches.

DATA DESCRIPTION

We found that: (1) 5475 out of 43,084 oligos (a) failed to align to the zebra finch genome, (b) aligned to multiple loci, or (c) aligned to Chr_un only, and thus need to be flagged until a better genome assembly is available, or (d) reflect cloning artifacts; (2) Out of 9635 valid oligos examined further, 3120 were incorrectly named, including 1533 with no known orthologs; and (3) 2635 oligos required name update. The resulting curated dataset provides a reference for correcting gene identification errors in previous finch microarrays studies, and avoiding such errors in future studies.

Characterization of Chicken Tumor Necrosis Factor-α, a Long Missed Cytokine in Birds

Tumor necrosis factor-α (TNF-α) is a pleiotropic cytokine playing critical roles in host defense and acute and chronic inflammation. It has been described in fish, amphibians, and mammals but was considered to be absent in the avian genomes. Here, we report on the identification and functional characterization of the avian ortholog. The chicken TNF-α (chTNF-α) is encoded by a highly GC-rich gene, whose product shares with its mammalian counterpart 45% homology in the extracellular part displaying the characteristic TNF homology domain. Orthologs of chTNF-α were identified in the genomes of 12 additional avian species including Palaeognathae and Neognathae, and the synteny of the closely adjacent loci with mammalian TNF-α orthologs was demonstrated in the crow (Corvus cornix) genome. In addition to chTNF-α, we obtained full sequences for homologs of TNF-α receptors 1 and 2 (TNFR1, TNFR2). chTNF-α mRNA is strongly induced by lipopolysaccharide (LPS) stimulation of monocyte derived, splenic and bone marrow macrophages, and significantly upregulated in splenic tissue in response to i.v. LPS treatment. Activation of T-lymphocytes by TCR crosslinking induces chTNF-α expression in CD4⁺ but not in CD8⁺ cells. To gain insights into its biological activity, we generated recombinant chTNF-α in eukaryotic and prokaryotic expression systems. Both, the full-length cytokine and the extracellular domain rapidly induced an NFκB-luciferase reporter in stably transfected CEC-32 reporter cells. Collectively, these data provide strong evidence for the existence of a fully functional TNF-α/TNF-α receptor system in birds thus filling a gap in our understanding of the evolution of cytokine systems.

The constitutive differential transcriptome of a brain circuit for vocal learning

Background

The ability to imitate the vocalizations of other organisms, a trait known as vocal learning, is shared by only a few organisms, including humans, where it subserves the acquisition of speech and language, and 3 groups of birds. In songbirds, vocal learning requires the coordinated activity of a set of specialized brain nuclei referred to as the song control system. Recent efforts have revealed some of the genes that are expressed in these vocal nuclei, however a thorough characterization of the transcriptional specializations of this system is still missing. We conducted a rigorous and comprehensive analysis of microarrays, and conducted a separate analysis of 380 genes by in situ hybridizations in order to identify molecular specializations of the major nuclei of the song system of zebra finches (Taeniopygia guttata), a songbird species.

Results

Our efforts identified more than 3300 genes that are differentially regulated in one or more vocal nuclei of adult male birds compared to the adjacent brain regions. Bioinformatics analyses provided insights into the possible involvement of these genes in molecular pathways such as cellular morphogenesis, intrinsic cellular excitability, neurotransmission and neuromodulation, axonal guidance and cela-to-cell interactions, and cell survival, which are known to strongly influence the functional properties of the song system. Moreover, an in-depth analysis of specific gene families with known involvement in regulating the development and physiological properties of neuronal circuits provides further insights into possible modulators of the song system.

Conclusion

Our study represents one of the most comprehensive molecular characterizations of a brain circuit that evolved to facilitate a learned behavior in a vertebrate. The data provide novel insights into possible molecular determinants of the functional properties of the song control circuitry. It also provides lists of compelling targets for pharmacological and genetic manipulations to elucidate the molecular regulation of song behavior and vocal learning.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4578-0) contains supplementary material, which is available to authorized users.

The curious case of peroxiredoxin-5: what its absence in aves can tell us and how it can be used

Background

Peroxiredoxins are ubiquitous thiol-dependent peroxidases that represent a major antioxidant defense in both prokaryotic cells and eukaryotic organisms. Among the six vertebrate peroxiredoxin isoforms, peroxiredoxin-5 (PRDX5) appears to be a particular peroxiredoxin, displaying a different catalytic mechanism, as well as a wider substrate specificity and subcellular distribution. In addition, several evolutionary peculiarities, such as loss of subcellular targeting in certain species, have been reported for this enzyme.

Results

Western blotting analyses of 2-cys PRDXs (PRDX1–5) failed to identify the PRDX5 isoform in chicken tissue homogenates. Thereafter, via in silico analysis of PRDX5 orthologs, we went on to show that the PRDX5 gene is conserved in all branches of the amniotes clade, with the exception of aves. Further investigation of bird genomic sequences and expressed tag sequences confirmed the disappearance of the gene, though TRMT112, a gene located closely to the 5′ extremity of the PRDX5 gene, is conserved. Finally, using in ovo electroporation to overexpress the long and short forms of human PRDX5, we showed that, though the gene is lost in birds, subcellular targeting of human PRDX5 is conserved in the chick.

Conclusions

Further adding to the distinctiveness of this enzyme, this study reports converging evidence supporting loss of PRDX5 in aves. In-depth analysis revealed that this absence is proper to birds as PRDX5 appears to be conserved in non-avian amniotes. Finally, taking advantage of the in ovo electroporation technique, we validate the subcellular targeting of human PRDX5 in the chick embryo and bring forward this gain-of-function model as a potent way to study PRDX5 functions in vivo.

Inferring synteny between genome assemblies: a systematic evaluation

BACKGROUND

Genome assemblies across all domains of life are being produced routinely. Initial analysis of a new genome usually includes annotation and comparative genomics. Synteny provides a framework in which conservation of homologous genes and gene order is identified between genomes of different species. The availability of human and mouse genomes paved the way for algorithm development in large-scale synteny mapping, which eventually became an integral part of comparative genomics. Synteny analysis is regularly performed on assembled sequences that are fragmented, neglecting the fact that most methods were developed using complete genomes. It is unknown to what extent draft assemblies lead to errors in such analysis.

RESULTS

We fragmented genome assemblies of model nematodes to various extents and conducted synteny identification and downstream analysis. We first show that synteny between species can be underestimated up to 40% and find disagreements between popular tools that infer synteny blocks. This inconsistency and further demonstration of erroneous gene ontology enrichment tests raise questions about the robustness of previous synteny analysis when gold standard genome sequences remain limited. In addition, assembly scaffolding using a reference guided approach with a closely related species may result in chimeric scaffolds with inflated assembly metrics if a true evolutionary relationship was overlooked. Annotation quality, however, has minimal effect on synteny if the assembled genome is highly contiguous.

CONCLUSIONS

Our results show that a minimum N50 of 1 Mb is required for robust downstream synteny analysis, which emphasizes the importance of gold standard genomes to the science community, and should be achieved given the current progress in sequencing technology.

A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model

The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies.

Avian genomics lends insights into endocrine function in birds

The genomics era has brought along the completed sequencing of a large number of bird genomes that cover a broad range of the avian phylogenetic tree (>30 orders), leading to major novel insights into avian biology and evolution. Among recent findings, the discovery that birds lack a large number of protein coding genes that are organized in highly conserved syntenic clusters in other vertebrates is very intriguing, given the physiological importance of many of these genes. A considerable number of them play prominent endocrine roles, suggesting that birds evolved compensatory genetic or physiological mechanisms that allowed them to survive and thrive in spite of these losses. While further studies are needed to establish the exact extent of avian gene losses, these findings point to birds as potentially highly relevant model organisms for exploring the genetic basis and possible therapeutic approaches for a wide range of endocrine functions and disorders.

Chicken Is a Useful Model to Investigate the Role of Adipokines in Metabolic and Reproductive Diseases

Reproduction is a complex and essential physiological process required by all species to produce a new generation. This process involves strict hormonal regulation, depending on a connection between the hypothalamus-pituitary-gonadal axis and peripheral organs. Metabolic homeostasis influences the reproductive functions, and its alteration leads to disturbances in the reproductive functions of humans as well as animals. For a long time, adipose tissue has been recognised as an endocrine organ but its ability to secrete and release hormones called adipokines is now emerging. Adipokines have been found to play a major role in the regulation of metabolic and reproductive processes at both central and peripheral levels. Leptin was initially the first adipokine that has been described to be the most involved in the metabolism/reproduction interrelation in mammals. In avian species, the role of leptin is still under debate. Recently, three novel adipokines have been discovered: adiponectin (ADIPOQ, ACRP30), visfatin (NAMPT, PBEF), and chemerin (RARRES2, TIG2). However, their mode of action between mammalian and nonmammalian species is different due to the different reproductive and metabolic systems. Herein, we will provide an overview of the structure and function related to metabolic and reproductive mechanisms of the latter three adipokines with emphasis on avian species.

Androgen and estrogen sensitivity of bird song: a comparative view on gene regulatory levels

Singing of songbirds is sensitive to testosterone and its androgenic and estrogenic metabolites in a species-specific way. The hormonal effects on song pattern are likely mediated by androgen receptors (AR) and estrogen receptor alpha (ERα), ligand activated transcription factors that are expressed in neurons of various areas of the songbirds' vocal control circuit. The distribution of AR in this circuit is rather similar between species while that of ERα is species variant and concerns a key vocal control area, the HVC (proper name). We discuss the regulation of the expression of the cognate AR and ERα and putative splice variants. In particular, we suggest that transcription factor binding sites in the promoter of these receptors differ between bird species. Further, we suggest that AR- and ERα-dependent gene regulation in vocal areas differs between species due to species-specific DNA binding sites of putative target genes that are required for the transcriptional activity of the receptors. We suggest that species differences in the distribution of AR and ERα in vocal areas and in the genomic sensitivity to these receptors contribute to species-specific hormonal regulation of the song.

Bat Biology, Genomes, and the Bat1K Project: To Generate Chromosome-Level Genomes for All Living Bat Species

Bats are unique among mammals, possessing some of the rarest mammalian adaptations, including true self-powered flight, laryngeal echolocation, exceptional longevity, unique immunity, contracted genomes, and vocal learning. They provide key ecosystem services, pollinating tropical plants, dispersing seeds, and controlling insect pest populations, thus driving healthy ecosystems. They account for more than 20% of all living mammalian diversity, and their crown-group evolutionary history dates back to the Eocene. Despite their great numbers and diversity, many species are threatened and endangered. Here we announce Bat1K, an initiative to sequence the genomes of all living bat species (n∼1,300) to chromosome-level assembly. The Bat1K genome consortium unites bat biologists (>148 members as of writing), computational scientists, conservation organizations, genome technologists, and any interested individuals committed to a better understanding of the genetic and evolutionary mechanisms that underlie the unique adaptations of bats. Our aim is to catalog the unique genetic diversity present in all living bats to better understand the molecular basis of their unique adaptations; uncover their evolutionary history; link genotype with phenotype; and ultimately better understand, promote, and conserve bats. Here we review the unique adaptations of bats and highlight how chromosome-level genome assemblies can uncover the molecular basis of these traits. We present a novel sequencing and assembly strategy and review the striking societal and scientific benefits that will result from the Bat1K initiative.

LXRα and LXRβ Nuclear Receptors Evolved in the Common Ancestor of Gnathostomes

Nuclear receptors (NRs) regulate numerous aspects of the endocrine system. They mediate endogenous and exogenous cues, ensuring a homeostatic control of development and metabolism. Gene duplication, loss and mutation have shaped the repertoire and function of NRs in metazoans. Here, we examine the evolution of a pivotal orchestrator of cholesterol metabolism in vertebrates, the liver X receptors (LXRs). Previous studies suggested that LXRα and LXRβ genes emerged in the mammalian ancestor. However, we show through genome analysis and functional assay that bona fide LXRα and LXRβ orthologues are present in reptiles, coelacanth and chondrichthyans but not in cyclostomes. These findings show that LXR duplicated before gnathostome radiation, followed by asymmetric paralogue loss in some lineages. We suggest that a tighter control of cholesterol levels in vertebrates was achieved through the exploitation of a wider range of oxysterols, an ability contingent on ligand-binding pocket remodeling.

Identification, chromosomal arrangements and expression analyses of the evolutionarily conserved prmt1 gene in chicken in comparison with its vertebrate paralogue prmt8

Nine protein arginine methyltransferases (PRMTs) are conserved in mammals and fish. Among these, PRMT1 is the major type I PRMT for asymmetric dimethylarginine (ADMA) formation and is the most conserved and widely distributed one. Two chicken prmt1 splicing variants were assembled and confirmed by RT-PCR experiments. However, only two scaffolds containing single separate prmt1 exon with high GC contents are present in the current chicken genome assembly. Besides, prmt1 exons are scattered in separate small scaffolds in most avian species. Complete prmt1 gene has only been predicted from two falcon species with few neighboring genes. Crocodilians are considered close to the common ancestor shared by crocodilians and birds. The gene arrangements around prmt1 in American alligator are different from that in birds but are largely conserved in human. Orthologues of genes in a large segment of human chromosomal 19 around PRMT1 are missing or not assigned to the current chicken chromosomes. In comparison, prmt8, the prmt1 paralogue, is on chicken chromosome 1 with the gene arrangements downstream of prmt8 highly conserved in birds, crocodilians, and human. However, the ones upstream vary greatly in birds. Biochemically, we found that though prmt1 transcripts were detected, limited or none PRMT1 protein was present in chicken tissues. Moreover, a much higher level of PRMT8 protein was detected in chicken brain than in mouse brain. While PRMT8 is brain specific in other vertebrate species studied, low level of PRMT8 was present in chicken but not mouse liver and muscle. We also showed that the ADMA level in chicken was similar to that in mouse. This study provides the critical information of chicken PRMT1 and PRMT8 for future analyses of the function of protein arginine methyltransferases in birds.

Mapping of leptin and its syntenic genes to chicken chromosome 1p

Background

Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly.

Results

A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = − 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes.

Conclusion

Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0543-1) contains supplementary material, which is available to authorized users.

Correspondence on Lovell et al.: response to Bornelöv et al

While the analysis of Bornelöv et al. is informative, they provide evidence for the existence of only 3% of the reported avian missing genes set, and thus do not significantly challenge our main findings that specific groups of syntenic protein-coding genes are missing in birds.

This is a response to the Correspondence article: https://www.dx.doi.org/10.1186/s13059-017-1231-1

Correspondence on Lovell et al.: identification of chicken genes previously assumed to be evolutionarily lost

Through RNA-Seq analyses, we identified 137 genes that are missing in chicken, including the long-sought-after nephrin and tumor necrosis factor genes. These genes tended to cluster in GC-rich regions that have poor coverage in genome sequence databases. Hence, the occurrence of syntenic groups of vertebrate genes that have not been observed in Aves does not prove the evolutionary loss of such genes.Please see related Research article: http://dx.doi.org/10.1186/s13059-014-0565-1 and Please see response from Lovell et al: https://www.dx.doi.org/10.1186/s13059-017-1234-y.

Evolution of the functionally conserved DCC gene in birds

Understanding the loss of conserved genes is critical for determining how phenotypic diversity is generated. Here we focus on the evolution of DCC, a gene that encodes a highly conserved neural guidance receptor. Disruption of DCC in animal models and humans results in major neurodevelopmental defects including commissural axon defects. Here we examine DCC evolution in birds, which is of particular interest as a major model system in neurodevelopmental research. We found the DCC containing locus was disrupted several times during evolution, resulting in both gene losses and faster evolution rate of salvaged genes. These data suggest that DCC had been lost independently twice during bird evolution, including in chicken and zebra finch, whereas it was preserved in many other closely related bird species, including ducks. Strikingly, we observed that commissural axon trajectory appeared similar regardless of whether DCC could be detected or not. We conclude that the DCC locus is susceptible to genomic instability leading to independent disruptions in different branches of birds and a significant influence on evolution rate. Overall, the phenomenon of loss or molecular evolution of a highly conserved gene without apparent phenotype change is of conceptual importance for understanding molecular evolution of key biological processes.

Dynamics of genome size evolution in birds and mammals

Genome size in mammals and birds shows remarkably little interspecific variation compared with other taxa. However, genome sequencing has revealed that many mammal and bird lineages have experienced differential rates of transposable element (TE) accumulation, which would be predicted to cause substantial variation in genome size between species. Thus, we hypothesize that there has been covariation between the amount of DNA gained by transposition and lost by deletion during mammal and avian evolution, resulting in genome size equilibrium. To test this model, we develop computational methods to quantify the amount of DNA gained by TE expansion and lost by deletion over the last 100 My in the lineages of 10 species of eutherian mammals and 24 species of birds. The results reveal extensive variation in the amount of DNA gained via lineage-specific transposition, but that DNA loss counteracted this expansion to various extents across lineages. Our analysis of the rate and size spectrum of deletion events implies that DNA removal in both mammals and birds has proceeded mostly through large segmental deletions (>10 kb). These findings support a unified "accordion" model of genome size evolution in eukaryotes whereby DNA loss counteracting TE expansion is a major determinant of genome size. Furthermore, we propose that extensive DNA loss, and not necessarily a dearth of TE activity, has been the primary force maintaining the greater genomic compaction of flying birds and bats relative to their flightless relatives.

A New Chicken Genome Assembly Provides Insight into Avian Genome Structure

The importance of the Gallus gallus (chicken) as a model organism and agricultural animal merits a continuation of sequence assembly improvement efforts. We present a new version of the chicken genome assembly (Gallus_gallus-5.0; GCA_000002315.3), built from combined long single molecule sequencing technology, finished BACs, and improved physical maps. In overall assembled bases, we see a gain of 183 Mb, including 16.4 Mb in placed chromosomes with a corresponding gain in the percentage of intact repeat elements characterized. Of the 1.21 Gb genome, we include three previously missing autosomes, GGA30, 31, and 33, and improve sequence contig length 10-fold over the previous Gallus_gallus-4.0. Despite the significant base representation improvements made, 138 Mb of sequence is not yet located to chromosomes. When annotated for gene content, Gallus_gallus-5.0 shows an increase of 4679 annotated genes (2768 noncoding and 1911 protein-coding) over those in Gallus_gallus-4.0. We also revisited the question of what genes are missing in the avian lineage, as assessed by the highest quality avian genome assembly to date, and found that a large fraction of the original set of missing genes are still absent in sequenced bird species. Finally, our new data support a detailed map of MHC-B, encompassing two segments: one with a highly stable gene copy number and another in which the gene copy number is highly variable. The chicken model has been a critical resource for many other fields of study, and this new reference assembly will substantially further these efforts.

Novel Insights into Chromosome Evolution in Birds, Archosaurs, and Reptiles

Homologous synteny blocks (HSBs) and evolutionary breakpoint regions (EBRs) in mammalian chromosomes are enriched for distinct DNA features, contributing to distinct phenotypes. To reveal HSB and EBR roles in avian evolution, we performed a sequence-based comparison of 21 avian and 5 outgroup species using recently sequenced genomes across the avian family tree and a newly-developed algorithm. We identified EBRs and HSBs in ancestral bird, archosaurian (bird, crocodile, and dinosaur), and reptile chromosomes. Genes involved in the regulation of gene expression and biosynthetic processes were preferably located in HSBs, including for example, avian-specific HSBs enriched for genes involved in limb development. Within birds, some lineage-specific EBRs rearranged genes were related to distinct phenotypes, such as forebrain development in parrots. Our findings provide novel evolutionary insights into genome evolution in birds, particularly on how chromosome rearrangements likely contributed to the formation of novel phenotypes.

Deep landscape update of dispersed and tandem repeats in the genome model of the red jungle fowl, Gallus gallus, using a series of de novo investigating tools

BACKGROUND

The program RepeatMasker and the database Repbase-ISB are part of the most widely used strategy for annotating repeats in animal genomes. They have been used to show that avian genomes have a lower repeat content (8-12 %) than the sequenced genomes of many vertebrate species (30-55 %). However, the efficiency of such a library-based strategies is dependent on the quality and completeness of the sequences in the database that is used. An alternative to these library based methods are methods that identify repeats de novo. These alternative methods have existed for a least a decade and may be more powerful than the library based methods. We have used an annotation strategy involving several complementary de novo tools to determine the repeat content of the model genome galGal4 (1.04 Gbp), including identifying simple sequence repeats (SSRs), tandem repeats and transposable elements (TEs).

RESULTS

We annotated over one Gbp. of the galGal4 genome and showed that it is composed of approximately 19 % SSRs and TEs repeats. Furthermore, we estimate that the actual genome of the red jungle fowl contains about 31-35 % repeats. We find that library-based methods tend to overestimate TE diversity. These results have a major impact on the current understanding of repeats distributions throughout chromosomes in the red jungle fowl.

CONCLUSIONS

Our results are a proof of concept of the reliability of using de novo tools to annotate repeats in large animal genomes. They have also revealed issues that will need to be resolved in order to develop gold-standard methodologies for annotating repeats in eukaryote genomes.

Corticotropin-releasing hormone family evolution: five ancestral genes remain in some lineages

The evolution of the peptide family consisting of corticotropin-releasing hormone (CRH) and the three urocortins (UCN1-3) has been puzzling due to uneven evolutionary rates. Distinct gene duplication scenarios have been proposed in relation to the two basal rounds of vertebrate genome doubling (2R) and the teleost fish-specific genome doubling (3R). By analyses of sequences and chromosomal regions, including many neighboring gene families, we show here that the vertebrate progenitor had two peptide genes that served as the founders of separate subfamilies. Then, 2R resulted in a total of five members: one subfamily consists of CRH1, CRH2, and UCN1. The other subfamily contains UCN2 and UCN3. All five peptide genes are present in the slowly evolving genomes of the coelacanth Latimeria chalumnae (a lobe-finned fish), the spotted gar Lepisosteus oculatus (a basal ray-finned fish), and the elephant shark Callorhinchus milii (a cartilaginous fish). The CRH2 gene has been lost independently in placental mammals and in teleost fish, but is present in birds (except chicken), anole lizard, and the nonplacental mammals platypus and opossum. Teleost 3R resulted in an additional surviving duplicate only for crh1 in some teleosts including zebrafish (crh1a and crh1b). We have previously reported that the two vertebrate CRH/UCN receptors arose in 2R and that CRHR1 was duplicated in 3R. Thus, we can now conclude that this peptide-receptor system was quite complex in the ancestor of the jawed vertebrates with five CRH/UCN peptides and two receptors, and that crh and crhr1 were duplicated in the teleost fish tetraploidization.

Characterization of Chicken Thyroid Hormone Transporters

Thyroid hormone (TH) transmembrane transporters are key regulators of TH availability in target cells where correct TH signaling is essential for normal development. Although the chicken embryo is a valuable model for developmental studies, the only functionally characterized chicken TH transporter so far is the organic anion transporting polypeptide 1C1 (OATP1C1). We therefore cloned the chicken L-type amino acid transporter 1 (LAT1) and the monocarboxylate transporters 8 (MCT8) and 10 (MCT10), and functionally characterized them, together with OATP1C1, in JEG3, COS1, and DF-1 cells. In addition, we used in situ hybridization to study their mRNA expression pattern during development. MCT8 and OATP1C1 are both high affinity transporters for the prohormone T4, whereas receptor-active T3 is preferably transported by MCT8 and MCT10. The latter one shows lower affinity but has a high Vmax and seems to be especially good at T3 export. Also, LAT1 has a lower affinity for its preferred substrate 3,3'-diiodothyronine. Reverse T3 is transported by all 4 TH transporters and is a good export product for OATP1C1. TH transporters are strongly expressed in eye (LAT1, MCT8, MCT10), pancreas (LAT1, MCT10), kidney, and testis (MCT8). Their extensive expression in the central nervous system, especially at the brain barriers, indicates an important role in brain development. In conclusion, we show TH transport by chicken MCT8, MCT10, and LAT1. Together with OATP1C1, these transporters have functional characteristics similar to their mammalian orthologs and are interesting target genes to further elucidate the role of THs during embryonic development.

Adaptation to nocturnality - learning from avian genomes

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

Contraction of the type I IFN locus and unusual constitutive expression of IFN-α in bats

Bats harbor many emerging and reemerging viruses, several of which are highly pathogenic in other mammals but cause no clinical signs of disease in bats. To determine the role of interferons (IFNs) in the ability of bats to coexist with viruses, we sequenced the type I IFN locus of the Australian black flying fox, Pteropus alecto, providing what is, to our knowledge, the first gene map of the IFN region of any bat species. Our results reveal a highly contracted type I IFN family consisting of only 10 IFNs, including three functional IFN-α loci. Furthermore, the three IFN-α genes are constitutively expressed in unstimulated bat tissues and cells and their expression is unaffected by viral infection. Constitutively expressed IFN-α results in the induction of a subset of IFN-stimulated genes associated with antiviral activity and resistance to DNA damage, providing evidence for a unique IFN system that may be linked to the ability of bats to coexist with viruses.

Identification of the Long-Sought Leptin in Chicken and Duck: Expression Pattern of the Highly GC-Rich Avian leptin Fits an Autocrine/Paracrine Rather Than Endocrine Function

More than 20 years after characterization of the key regulator of mammalian energy balance, leptin, we identified the leptin (LEP) genes of chicken (Gallus gallus) and duck (Anas platyrhynchos). The extreme guanine-cytosine content (∼70%), the location in a genomic region with low-complexity repetitive and palindromic sequence elements, the relatively low sequence conservation, and low level of expression have hampered the identification of these genes until now. In vitro-expressed chicken and duck leptins specifically activated signaling through the chicken leptin receptor in cell culture. In situ hybridization demonstrated expression of LEP mRNA in granular and Purkinje cells of the cerebellum, anterior pituitary, and in embryonic limb buds, somites, and branchial arches, suggesting roles in adult brain control of energy balance and during embryonic development. The expression patterns of LEP and the leptin receptor (LEPR) were explored in chicken, duck, and quail (Coturnix japonica) using RNA-sequencing experiments available in the Short Read Archive and by quantitative RT-PCR. In adipose tissue, LEP and LEPR were scarcely transcribed, and the expression level was not correlated to adiposity. Our identification of the leptin genes in chicken and duck genomes resolves a long lasting controversy regarding the existence of leptin genes in these species. This identification was confirmed by sequence and structural similarity, conserved exon-intron boundaries, detection in numerous genomic, and transcriptomic datasets and characterization by PCR, quantitative RT-PCR, in situ hybridization, and bioassays. Our results point to an autocrine/paracrine mode of action for bird leptin instead of being a circulating hormone as in mammals.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Rapid genome reshaping by multiple-gene loss after whole-genome duplication in teleost fish suggested by mathematical modeling

Whole-genome duplication (WGD) is believed to be a significant source of major evolutionary innovation. Redundant genes resulting from WGD are thought to be lost or acquire new functions. However, the rates of gene loss and thus temporal process of genome reshaping after WGD remain unclear. The WGD shared by all teleost fish, one-half of all jawed vertebrates, was more recent than the two ancient WGDs that occurred before the origin of jawed vertebrates, and thus lends itself to analysis of gene loss and genome reshaping. Using a newly developed orthology identification pipeline, we inferred the post-teleost-specific WGD evolutionary histories of 6,892 protein-coding genes from nine phylogenetically representative teleost genomes on a time-calibrated tree. We found that rapid gene loss did occur in the first 60 My, with a loss of more than 70-80% of duplicated genes, and produced similar genomic gene arrangements within teleosts in that relatively short time. Mathematical modeling suggests that rapid gene loss occurred mainly by events involving simultaneous loss of multiple genes. We found that the subsequent 250 My were characterized by slow and steady loss of individual genes. Our pipeline also identified about 1,100 shared single-copy genes that are inferred to have become singletons before the divergence of clupeocephalan teleosts. Therefore, our comparative genome analysis suggests that rapid gene loss just after the WGD reshaped teleost genomes before the major divergence, and provides a useful set of marker genes for future phylogenetic analysis.

The Origin and Diversity of Cpt1 Genes in Vertebrate Species

The Carnitine palmitoyltransferase I (Cpt1) gene family plays a crucial role in energy homeostasis since it is required for the occurrence of fatty acid β-oxidation in the mitochondria. The exact gene repertoire in different vertebrate lineages is variable. Presently, four genes are documented: Cpt1a, also known as Cpt1a1, Cpt1a2; Cpt1b and Cpt1c. The later is considered a mammalian innovation resulting from a gene duplication event in the ancestor of mammals, after the divergence of sauropsids. In contrast, Cpt1a2 has been found exclusively in teleosts. Here, we reassess the overall evolutionary relationships of Cpt1 genes using a combination of approaches, including the survey of the gene repertoire in basal gnathostome lineages. Through molecular phylogenetics and synteny studies, we find that Cpt1c is most likely a rapidly evolving orthologue of Cpt1a2. Thus, Cpt1c is present in other lineages such as cartilaginous fish, reptiles, amphibians and the coelacanth. We show that genome duplications (2R) and variable rates of sequence evolution contribute to the history of Cpt1 genes in vertebrates. Finally, we propose that loss of Cpt1b is the likely cause for the unusual energy metabolism of elasmobranch.

Early Vertebrate Evolution of the Host Restriction Factor Tetherin

Tetherin is an interferon-inducible restriction factor targeting a broad range of enveloped viruses. Its antiviral activity depends on an unusual topology comprising an N-terminal transmembrane domain (TMD) followed by an extracellular coiled-coil region and a C-terminal glycosylphosphatidylinositol (GPI) anchor. One of the two membrane anchors is inserted into assembling virions, while the other remains in the plasma membrane of the infected cell. Thus, tetherin entraps budding viruses by physically bridging viral and cellular membranes. Although tetherin restricts the release of a large variety of diverse human and animal viruses, only mammalian orthologs have been described to date. Here, we examined the evolutionary origin of this protein and demonstrate that tetherin orthologs are also found in fish, reptiles, and birds. Notably, alligator tetherin efficiently blocks the release of retroviral particles. Thus, tetherin emerged early during vertebrate evolution and acquired its antiviral activity before the mammal/reptile divergence. Although there is only limited sequence homology, all orthologs share the typical topology. Two unrelated proteins of the slime mold Dictyostelium discoideum also adopt a tetherin-like configuration with an N-terminal TMD and a C-terminal GPI anchor. However, these proteins showed no evidence for convergent evolution and failed to inhibit virion release. In summary, our findings demonstrate that tetherin emerged at least 450 million years ago and is more widespread than previously anticipated. The early evolution of antiviral activity together with the high topology conservation but low sequence homology suggests that restriction of virus release is the primary function of tetherin.

IMPORTANCE

The continuous arms race with viruses has driven the evolution of a variety of cell-intrinsic immunity factors that inhibit different steps of the viral replication cycle. One of these restriction factors, tetherin, inhibits the release of newly formed progeny virions from infected cells. Although tetherin targets a broad range of enveloped viruses, including retro-, filo-, herpes-, and arenaviruses, the evolutionary origin of this restriction factor and its antiviral activity remained obscure. Here, we examined diverse vertebrate genomes for genes encoding cellular proteins that share with tetherin the highly unusual combination of an N-terminal transmembrane domain and a C-terminal glycosylphosphatidylinositol anchor. We show that tetherin orthologs are found in fish, reptiles, and birds and demonstrate that alligator tetherin efficiently inhibits the release of retroviral particles. Our findings identify tetherin as an evolutionarily ancient restriction factor and provide new important insights into the continuous arms race between viruses and their hosts.

Living without DAT: Loss and compensation of the dopamine transporter gene in sauropsids (birds and reptiles)

The dopamine transporter (DAT) is a major regulator of synaptic dopamine (DA) availability. It plays key roles in motor control and motor learning, memory formation, and reward-seeking behavior, is a major target of cocaine and methamphetamines, and has been assumed to be conserved among vertebrates. We have found, however, that birds, crocodiles, and lizards lack the DAT gene. We also found that the unprecedented loss of this important gene is compensated for by the expression of the noradrenaline transporter (NAT) gene, and not the serotonin transporter genes, in dopaminergic cells, which explains the peculiar pharmacology of the DA reuptake activity previously noted in bird striatum. This unexpected pattern contrasts with that of ancestral vertebrates (e.g. fish) and mammals, where the NAT gene is selectively expressed in noradrenergic cells. DA circuits in birds/reptiles and mammals thus operate with an analogous reuptake mechanism exerted by different genes, bringing new insights into gene expression regulation in dopaminergic cells and the evolution of a key molecular player in reward and addiction pathways.

Response to Hron et al

Hron et al. provide transcriptome evidence that three (1.1 %) of the 274 genes reported by Lovell et al. as missing in birds may actually be 'hidden' as a result of high GC content. Although this factor may explain some gene absences from genomic assemblies, we believe it is insufficient to account for the extensive syntenic losses described in Lovell et al. Please see related article: www.dx.doi.org/10.1186/s13059-015-0724-z.

Hidden genes in birds

We report that a subset of avian genes is characterized by very high GC content and long G/C stretches. These sequence characteristics correlate with the frequent absence of these genes from genomic databases. We provide several examples where genes in this subset are mistakenly reported as missing in birds.

www.dx.doi.org/10.1186/s13059-015-0725-y

The serine/threonine kinase 33 is present and expressed in palaeognath birds but has become a unitary pseudogene in neognaths about 100 million years ago

Background

Serine/threonine kinase 33 (STK33) has been shown to be conserved across all major vertebrate classes including reptiles, mammals, amphibians and fish, suggesting its importance within vertebrates. It has been shown to phosphorylate vimentin and might play a role in spermatogenesis and organ ontogenesis. In this study we analyzed the genomic locus and expression of stk33 in the class Aves, using a combination of large scale next generation sequencing data analysis and traditional PCR.

Results

Within the subclass Palaeognathae we analyzed the white-throated tinamou (Tinamus guttatus), the African ostrich (Struthio camelus) and the emu (Dromaius novaehollandiae). For the African ostrich we were able to generate a 62,778 bp long genomic contig and an mRNA sequence that encodes a protein showing highly significant similarity to STK33 proteins from other vertebrates. The emu has been shown to encode and transcribe a functional STK33 as well. For the white-throated tinamou we were able to identify 13 exons by sequence comparison encoding a protein similar to STK33 as well.

In contrast, in all 28 neognath birds analyzed, we could not find evidence for the existence of a functional copy of stk33 or its expression. In the genomes of these 28 bird species, we found only remnants of the stk33 locus carrying several large genomic deletions, leading to the loss of multiple exons. The remaining exons have acquired various indels and premature stop codons.

Conclusions

We were able to elucidate and describe the genomic structure and the transcription of a functional stk33 gene within the subclass Palaeognathae, but we could only find degenerate remnants of stk33 in all neognath birds analyzed. This led us to the conclusion that stk33 became a unitary pseudogene in the evolutionary history of the class Aves at the paleognath-neognath branch point during the late cretaceous period about 100 million years ago. We hypothesize that the pseudogenization of stk33 might have become fixed in neognaths due to either genetic redundancy or a non-orthologous gene displacement and present potential candidate genes for such an incident.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1769-9) contains supplementary material, which is available to authorized users.