Conserved syntenic clusters of protein coding genes are missing in birds

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents

BACKGROUND

Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches.

RESULTS

We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought.

CONCLUSION

The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

The ancestral type of the R-RAS protein has oncogenic potential

BACKGROUND

The R-RAS2 is a small GTPase highly similar to classical RAS proteins at the regulatory and signaling levels. The high evolutionary conservation of R-RAS2, its links to basic cellular processes and its role in cancer, make R-RAS2 an interesting research topic. To elucidate the evolutionary history of R-RAS proteins, we investigated and compared structural and functional properties of ancestral type R-RAS protein with human R-RAS2.

METHODS

Bioinformatics analysis were used to elucidate the evolution of R-RAS proteins. Intrinsic GTPase activity of purified human and sponge proteins was analyzed with GTPase-GloTM Assay kit. The cell model consisted of human breast cancer cell lines MCF-7 and MDA-MB-231 transiently transfected with EsuRRAS2-like or HsaRRAS2. Biological characterization of R-RAS2 proteins was performed by Western blot on whole cell lysates or cell adhesion protein isolates, immunofluorescence and confocal microscopy, MTT test, colony formation assay, wound healing and Boyden chamber migration assays.

RESULTS

We found that the single sponge R-RAS2-like gene/protein probably reflects the properties of the ancestral R-RAS protein that existed prior to duplications during the transition to Bilateria, and to Vertebrata. Biochemical characterization of sponge and human R-RAS2 showed that they have the same intrinsic GTPase activity and RNA binding properties. By testing cell proliferation, migration and colony forming efficiency in MDA-MB-231 human breast cancer cells, we showed that the ancestral type of the R-RAS protein, sponge R-RAS2-like, enhances their oncogenic potential, similar to human R-RAS2. In addition, sponge and human R-RAS2 were not found in focal adhesions, but both homologs play a role in their regulation by increasing talin1 and vinculin.

CONCLUSIONS

This study suggests that the ancestor of all animals possessed an R-RAS2-like protein with oncogenic properties similar to evolutionarily more recent versions of the protein, even before the appearance of true tissue and the origin of tumors. Therefore, we have unraveled the evolutionary history of R-RAS2 in metazoans and improved our knowledge of R-RAS2 properties, including its structure, regulation and function.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

Gene purging and the evolution of Neoave metabolism and longevity

Maintenance of the proteasome requires oxidative phosphorylation (ATP) and mitigation of oxidative damage, in an increasingly dysfunctional relationship with aging. SLC3A2 plays a role on both sides of this dichotomy as an adaptor to SLC7A5, a transporter of branched-chain amino acids (BCAA: Leu, Ile, Val), and to SLC7A11, a cystine importer supplying cysteine to the synthesis of the antioxidant glutathione. Endurance in mammalian muscle depends in part on oxidation of BCAA; however, elevated serum levels are associated with insulin resistance and shortened lifespans. Intriguingly, the evolution of modern birds (Neoaves) has entailed the purging of genes including SLC3A2, SLC7A5, -7, -8, -10, and SLC1A4, -5, largely removing BCAA exchangers and their interacting Na+/Gln symporters in pursuit of improved energetics. Additional gene purging included mitochondrial BCAA aminotransferase (BCAT2), pointing to reduced oxidation of BCAA and increased hepatic conversion to triglycerides and glucose. Fat deposits are anhydrous and highly reduced, maximizing the fuel/weight ratio for prolonged flight, but fat accumulation in muscle cells of aging humans contributes to inflammation and senescence. Duplications of the bidirectional α-ketoacid transporters SLC16A3, SLC16A7, the cystine transporters SLC7A9, SLC7A11, and N-glycan branching enzymes MGAT4B, MGAT4C in Neoaves suggests a shift to the transport of deaminated essential amino acid, and stronger mitigation of oxidative stress supported by the galectin lattice. We suggest that Alfred Lotka's theory of natural selection as a maximum power organizer (PNAS 8:151,1922) made an unusually large contribution to Neoave evolution. Further molecular analysis of Neoaves may reveal novel rewiring with applications for human health and longevity.

Noncompetitive immunoassay optimized for pharmacokinetic assessments of biologically active efruxifermin

Efruxifermin (EFX) is a homodimeric human IgG1 Fc-FGF21 fusion protein undergoing investigation for treatment of liver fibrosis due to nonalcoholic steatohepatitis (NASH), a prevalent and serious metabolic disease for which there is no approved treatment. Biological activity of FGF21 requires its intact C-terminus, which enables binding to its obligate co-receptor β-Klotho on the surface of target cells. This interaction is a prerequisite for FGF21 signal transduction through its canonical FGF receptors: FGFR1c, 2c, and 3c. Therefore, the C-terminus of each FGF21 polypeptide chain must be intact, with no proteolytic truncation, for EFX to exert its pharmacological activity in patients. A sensitive immunoassay for quantification of biologically active EFX in human serum was therefore needed to support pharmacokinetic assessments in patients with NASH. We present a validated noncompetitive electrochemiluminescent immunoassay (ECLIA) that employs a rat monoclonal antibody for specific capture of EFX via its intact C-terminus. Bound EFX is detected by a SULFO-TAG™-conjugated, affinity purified chicken anti-EFX antiserum. The ECLIA reported herein for quantification of EFX demonstrated suitable analytical performance, with a sensitivity (LLOQ) of 20.0 ng/mL, to support reliable pharmacokinetic assessments of EFX. The validated assay was used to quantify serum EFX concentrations in a phase 2a study of NASH patients (BALANCED) with either moderate-to-advanced fibrosis or compensated cirrhosis. The pharmacokinetic profile of EFX was dose-proportional and did not differ between patients with moderate-to-advanced fibrosis and those with compensated cirrhosis. This report presents the first example of a validated pharmacokinetic assay specific for a biologically active Fc-FGF21 fusion protein, as well as the first demonstration of use of a chicken antibody conjugate as a detection reagent specific for an FGF21 analog.

The conservation and diversity of the exons encoding the glycine and arginine rich domain of the fibrillarin gene in vertebrates, with special focus on reptiles and birds

The nucleolar rRNA 2'-O-methyltransferase fibrillarin (FBL) contains a highly conserved methyltransferase domain at the C-terminus and a diverse glycine arginine-rich (GAR) domain at the N-terminus in eukaryotes. We found that a nine-exon configuration of fbl and exon 2-3 encoded GAR domain are conserved and specific in vertebrates. All internal exons except exon 2 and 3 are of the same lengths in different vertebrate lineages. The lengths of exon 2 and 3 vary in different vertebrate species but the ones with longer exon 2 usually have shorter exon 3 complementarily, limiting lengths of the GAR domain within a certain range. In tetrapods except for reptiles, exon 2 appears to be longer than exon 3. We specifically analyzed different lineages of reptiles for their GAR sequences and exon lengths. The lengths of exon 2 in reptiles are around 80-130-nt shorter and the lengths of exon 3 in reptiles are around 50-90 nt longer than those in other tetrapods, all in the GAR-coding regions. An FSPR sequence is present at the beginning of the GAR domain encoded by exon 2 in all vertebrates, and a specific FXSP/G element (X can be K, R, Q, N, and H) exist in the middle of GAR with phenylalanine as the 3rd exon 3-encoded amino acid residue starting from jawfish. Snakes, turtles, and songbirds contain shorter exon 2 compared with lizards, indicating continuous deletions in exon 2 and insertions/duplications in exon 3 in these lineages. Specifically, we confirmed the presence the fbl gene in chicken and validated the RNA expression. Our analyses of the GAR-encoding exons of fbl in vertebrates and reptiles should provide the basis for further evolutionary analyses of more GAR domain encoding proteins.

Development of specific monoclonal antibodies for the detection of natural chicken tumor necrosis factor-alpha

Tumor necrosis factor alpha (TNF-α) is an important proinflammatory cytokine and the only known cytokine that can directly kill tumor cells. Unlike mammalian counterparts, chicken TNF-α (chTNF-α) gene has not been identified until very recently due to its high GC content (∼70%) and long GC fragments. The biological functions of this newly-identified cytokine and its detection methods remain to be further investigated. In this study, the extracellular domain of chTNF-α was cloned into prokaryotic vector after codon optimization and recombinant chTNF-α protein was expressed. Subsequently, using recombinant chTNF-ɑ as immunogen, rabbit polyclonal antibody (pAb) and eight clones of mouse anti-chTNF-ɑ monoclonal antibodies (mAbs) were produced, respectively. Both the pAb and mAbs specifically recognized recombinant chTNF-ɑ expressed in E.coli and transfected COS-7 cells. Further mapping the antigenic region showed that all the mAbs recognized a region of amino acid residues 195-285 of chTNF-ɑ. Furthermore, an antigen-capture enzyme-linked immunosorbent assay for the detection of chTNF-ɑ was established using one mAb and the pAb. This assay showed no cross-reactivity with irrelevant Trx-fused antigens and could detect natural chTNF-ɑ expressed by mitogen-activated chicken splenocytes in a dose-dependent manner, with a detection limit of 1 ng/mL. Collectively, our results indicated that the mAbs and pAb against chTNF-α are specific and could be used for the study of the biological functions of chTNF-ɑ and the detection of chTNF-ɑ.

New insights into the novel sequences of the chicken pan-genome by liquid chip

Increasing evidence indicates that the missing sequences and genes in the chicken reference genome are involved in many crucial biological pathways, including metabolism and immunity. The low detection rate of novel sequences by resequencing hindered the acquisition of these sequences and the exploration of the relationship between new genes and economic traits. To improve the capture ratio of novel sequences, a 48K liquid chip including 25K from the reference sequence and 23K from the novel sequence was designed. The assay was tested on a panel of 218 animals from 5 chicken breeds. The average capture ratio of the reference sequence was 99.55%, and the average sequencing depth of the target sites was approximately 187X, indicating a good performance and successful application of liquid chips in farm animals. For the target region in the novel sequence, the average capture ratio was 33.15% and the average sequencing depth of target sites was approximately 60X, both of which were higher than that of resequencing. However, the different capture ratios and capture regions among varieties and individuals proved the difficulty of capturing these regions with complex structures. After genotyping, GWAS showed variations in novel sequences potentially relevant to immune-related traits. For example, a SNP close to the differentiation of lymphocyte-related gene IGHV3-23-like was associated with the H/L ratio. These results suggest that targeted capture sequencing is a preferred method to capture these sequences with complex structures and genes potentially associated with immune-related traits.

False gene and chromosome losses in genome assemblies caused by GC content variation and repeats

BACKGROUND

Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements.

RESULTS

Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies.

CONCLUSIONS

Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

Annotations of novel antennae-expressed genes in male Glossina morsitans morsitans tsetse flies

Tsetse flies use antennal expressed genes to navigate their environment. While most canonical genes associated with chemoreception are annotated, potential gaps with important antennal genes are uncharacterized in Glossina morsitans morsitans. We generated antennae-specific transcriptomes from adult male G. m. morsitans flies fed/unfed on bloodmeal and/or exposed to an attractant (ε-nonalactone), a repellant (δ-nonalactone) or paraffin diluent. Using bioinformatics approach, we mapped raw reads onto G. m. morsitans gene-set from VectorBase and collected un-mapped reads (constituting the gaps in annotation). We de novo assembled these reads (un-mapped) into transcript and identified corresponding genes of the transcripts in G. m. morsitans gene-set and protein homologs in UniProt protein database to further annotate the gaps. We predicted potential protein-coding gene regions associated with these transcripts in G. m. morsitans genome, annotated/curated these genes and identified their putative annotated orthologs/homologs in Drosophila melanogaster, Musca domestica or Anopheles gambiae genomes. We finally evaluated differential expression of the novel genes in relation to odor exposures relative to no-odor control (unfed flies). About 45.21% of the sequenced reads had no corresponding transcripts within G. m. morsitans gene-set, corresponding to the gap in existing annotation of the tsetse fly genome. The total reads assembled into 72,428 unique transcripts, most (74.43%) of which had no corresponding genes in the UniProt database. We annotated/curated 592 genes from these transcripts, among which 202 were novel while 390 were improvements of existing genes in the G. m. morsitans genome. Among the novel genes, 94 had orthologs in D. melanogaster, M. domestica or An. gambiae while 88 had homologs in UniProt. These orthologs were putatively associated with oxidative regulation, protein synthesis, transcriptional and/or translational regulation, detoxification and metal ion binding, thus providing insight into their specific roles in antennal physiological processes in male G. m. morsitans. A novel gene (GMOY014237.R1396) was differentially expressed in response to the attractant. We thus established significant gaps in G. m. morsitans genome annotation and identified novel male antennae-expressed genes in the genome, among which > 53% (108) are potentially G. m. morsitans specific.

De novo assembly of 20 chicken genomes reveals the undetectable phenomenon for thousands of core genes on micro-chromosomes and sub-telomeric regions

The gene numbers and evolutionary rates of birds were assumed to be much lower than those of mammals, which is in sharp contrast to the huge species number and morphological diversity of birds. It is, therefore, necessary to construct a complete avian genome and analyze its evolution. We constructed a chicken pan-genome from 20 de novo assembled genomes with high sequencing depth, and identified 1,335 protein-coding genes and 3,011 long noncoding RNAs not found in GRCg6a. The majority of these novel genes were detected across most individuals of the examined transcriptomes but were seldomly measured in each of the DNA sequencing data regardless of Illumina or PacBio technology. Furthermore, different from previous pan-genome models, most of these novel genes were overrepresented on chromosomal subtelomeric regions and microchromosomes, surrounded by extremely high proportions of tandem repeats, which strongly blocks DNA sequencing. These hidden genes were proved to be shared by all chicken genomes, included many housekeeping genes, and enriched in immune pathways. Comparative genomics revealed the novel genes had 3-fold elevated substitution rates than known ones, updating the knowledge about evolutionary rates in birds. Our study provides a framework for constructing a better chicken genome, which will contribute toward the understanding of avian evolution and the improvement of poultry breeding.

The Discovery of Chicken Foxp3 Demands Redefinition of Avian Regulatory T Cells

Since the publication of the first chicken genome sequence, we have encountered genes playing key roles in mammalian immunology, but being seemingly absent in birds. One of those was, until recently, Foxp3, the master transcription factor of regulatory T cells in mammals. Therefore, avian regulatory T cell research is still poorly standardized. In this study we identify a chicken ortholog of Foxp3 We prove sequence homology with known mammalian and sauropsid sequences, but also reveal differences in major domains. Expression profiling shows an association of Foxp3 and CD25 expression levels in CD4⁺CD25⁺ peripheral T cells and identifies a CD4^-CD25⁺Foxp3^high subset of thymic lymphocytes that likely represents yet undescribed avian regulatory T precursor cells. We conclude that Foxp3 is existent in chickens and that it shares certain functional characteristics with its mammalian ortholog. Nevertheless, pathways for regulatory T cell development and Foxp3 function are likely to differ between mammals and birds. The identification and characterization of chicken Foxp3 will help to define avian regulatory T cells and to analyze their functional properties and thereby advance the field of avian immunology.

Three chromosome-level duck genome assemblies provide insights into genomic variation during domestication

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed "missing genes" in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Bird evolution by insulin resistance

Drift of oxygen concentrations in the atmosphere was one of the main drivers of the evolution of vertebrates. The drop in oxygen concentrations at the Permian-Triassic (PT) boundary may have been the biggest challenge to vertebrates. This hypoxic condition forced theropods to lose certain genes to maximize their efficiency of oxygen usage. Recent studies show that omentin and insulin-sensitive glucose transporter 4 (GLUT4) are missing in the bird genome. Since these gene products play essential roles in maintaining insulin sensitivity, this loss forced theropods to become insulin resistant. Insulin resistance may have been the key to allowing theropods to become hyperathletic under hypoxic conditions and to outcompete mammals during the Triassic period. A second challenge was the gradual increase in oxygen concentrations during the late Jurassic, Cretaceous, and Tertiary periods when reactive oxygen species (ROS) leakage from mitochondria became a problem. Since the simplest solution was the expansion of body size, some theropods became bigger to reduce ROS leakage per volume. Another solution was the development of a constitutively active countermeasure against ROS. A recent study shows that Neoaves have constitutively active nuclear factor erythroid 2-related factor 2 (NRF2) due to deletion of the C-terminal part of the KEAP1 protein, thus allowing Neoaves to express antioxidant enzymes to overcome ROS leakage.

Sequencing refractory regions in bird genomes are hotspots for accelerated protein evolution

Background

Approximately 1000 protein encoding genes common for vertebrates are still unannotated in avian genomes. Are these genes evolutionary lost or are they not yet found for technical reasons? Using genome landscapes as a tool to visualize large-scale regional effects of genome evolution, we reexamined this question.

Results

On basis of gene annotation in non-avian vertebrate genomes, we established a list of 15,135 common vertebrate genes. Of these, 1026 were not found in any of eight examined bird genomes. Visualizing regional genome effects by our sliding window approach showed that the majority of these "missing" genes can be clustered to 14 regions of the human reference genome. In these clusters, an additional 1517 genes (often gene fragments) were underrepresented in bird genomes. The clusters of “missing” genes coincided with regions of very high GC content, particularly in avian genomes, making them “hidden” because of incomplete sequencing. Moreover, proteins encoded by genes in these sequencing refractory regions showed signs of accelerated protein evolution. As a proof of principle for this idea we experimentally characterized the mRNA and protein products of four "hidden" bird genes that are crucial for energy homeostasis in skeletal muscle: ALDOA, ENO3, PYGM and SLC2A4.

Conclusions

A least part of the “missing” genes in bird genomes can be attributed to an artifact caused by the difficulty to sequence regions with extreme GC% (“hidden” genes). Biologically, these “hidden” genes are of interest as they encode proteins that evolve more rapidly than the genome wide average. Finally we show that four of these “hidden” genes encode key proteins for energy metabolism in flight muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01905-7.

Adipokines in metabolic and reproductive functions in birds: An overview of current knowns and unknowns

Adipose tissue is now recognized as an active endocrine organ, which synthesizes and secretes numerous peptides factors called adipokines. In mammals, they exert pleiotropic effects affecting energy metabolism but also fertility. In mammals, secretion of adipokines is altered in adipose tissue dysfunctions and may participate to obesity-associated disorders. Thus, adipokines are promising candidates both for novel pharmacological treatment strategies and as diagnostic tools. As compared to mammals, birds exhibit several unique physiological features, which make them an interesting model for comparative studies on endocrine control of metabolism and adiposity and reproductive functions. Some adipokines such as leptin and visfatin may have different roles in avian species as compared to mammals. In addition, some of them found in mammals such as CCL2 (chemokine ligand 2), resistin, omentin and FGF21 (Fibroblast Growth factor 21) have not yet been mapped to the chicken genome model and among its annotated gene models. This brief review aims to summarize data (structure, metabolic and reproductive roles and molecular mechanisms involved) related to main avian adipokines (leptin, adiponectin, visfatin, and chemerin) and we will briefly discuss the adipokines that are still lacking in avian species.

An Unbiased Molecular Approach Using 3'-UTRs Resolves the Avian Family-Level Tree of Life

Presumably, due to a rapid early diversification, major parts of the higher-level phylogeny of birds are still resolved controversially in different analyses or are considered unresolvable. To address this problem, we produced an avian tree of life, which includes molecular sequences of one or several species of ∼90% of the currently recognized family-level taxa (429 species, 379 genera) including all 106 family-level taxa of the nonpasserines and 115 of the passerines (Passeriformes). The unconstrained analyses of noncoding 3-prime untranslated region (3′-UTR) sequences and those of coding sequences yielded different trees. In contrast to the coding sequences, the 3′-UTR sequences resulted in a well-resolved and stable tree topology. The 3′-UTR contained, unexpectedly, transcription factor binding motifs that were specific for different higher-level taxa. In this tree, grebes and flamingos are the sister clade of all other Neoaves, which are subdivided into five major clades. All nonpasserine taxa were placed with robust statistical support including the long-time enigmatic hoatzin (Opisthocomiformes), which was found being the sister taxon of the Caprimulgiformes. The comparatively late radiation of family-level clades of the songbirds (oscine Passeriformes) contrasts with the attenuated diversification of nonpasseriform taxa since the early Miocene. This correlates with the evolution of vocal production learning, an important speciation factor, which is ancestral for songbirds and evolved convergent only in hummingbirds and parrots. As 3′-UTR-based phylotranscriptomics resolved the avian family-level tree of life, we suggest that this procedure will also resolve the all-species avian tree of life

Natural loss of function of ephrin-B3 shapes spinal flight circuitry in birds

Flight in birds evolved through patterning of the wings from forelimbs and transition from alternating gait to synchronous flapping. In mammals, the spinal midline guidance molecule ephrin-B3 instructs the wiring that enables limb alternation, and its deletion leads to synchronous hopping gait. Here, we show that the ephrin-B3 protein in birds lacks several motifs present in other vertebrates, diminishing its affinity for the EphA4 receptor. The avian ephrin-B3 gene lacks an enhancer that drives midline expression and is missing in galliforms. The morphology and wiring at brachial levels of the chicken embryonic spinal cord resemble those of ephrin-B3 null mice. Dorsal midline decussation, evident in the mutant mouse, is apparent at the chick brachial level and is prevented by expression of exogenous ephrin-B3 at the roof plate. Our findings support a role for loss of ephrin-B3 function in shaping the avian brachial spinal cord circuitry and facilitating synchronous wing flapping.

The Airway Pathobiome in Complex Respiratory Diseases: A Perspective in Domestic Animals

Respiratory infections in domestic animals are a major issue for veterinary and livestock industry. Pathogens in the respiratory tract share their habitat with a myriad of commensal microorganisms. Increasing evidence points towards a respiratory pathobiome concept, integrating the dysbiotic bacterial communities, the host and the environment in a new understanding of respiratory disease etiology. During the infection, the airway microbiota likely regulates and is regulated by pathogens through diverse mechanisms, thereby acting either as a gatekeeper that provides resistance to pathogen colonization or enhancing their prevalence and bacterial co-infectivity, which often results in disease exacerbation. Insight into the complex interplay taking place in the respiratory tract between the pathogens, microbiota, the host and its environment during infection in domestic animals is a research field in its infancy in which most studies are focused on infections from enteric pathogens and gut microbiota. However, its understanding may improve pathogen control and reduce the severity of microbial-related diseases, including those with zoonotic potential.

SLC2A12 of SLC2 Gene Family in Bird Provides Functional Compensation for the Loss of SLC2A4 Gene in Other Vertebrates

Avian genomes are small and lack some genes that are conserved in the genomes of most other vertebrates including nonavian sauropsids. One hypothesis stated that paralogs may provide biochemical or physiological compensation for certain gene losses; however, no functional evidence has been reported to date. By integrating evolutionary analysis, physiological genomics, and experimental gene interference, we clearly demonstrate functional compensation for gene loss. A large-scale phylogenetic analysis of over 1,400 SLC2 gene sequences identifies six new SLC2 genes from nonmammalian vertebrates and divides the SLC2 gene family into four classes. Vertebrates retain class III SLC2 genes but partially lack the more recent duplicates of classes I and II. Birds appear to have completely lost the SLC2A4 gene that encodes an important insulin-sensitive GLUT in mammals. We found strong evidence for positive selection, indicating that the N-termini of SLC2A4 and SLC2A12 have undergone diversifying selection in birds and mammals, and there is a significant correlation between SLC2A12 functionality and basal metabolic rates in endotherms. Physiological genomics have uncovered that SLC2A12 expression and allelic variants are associated with insulin sensitivity and blood glucose levels in wild birds. Functional tests have indicated that SLC2A12 abrogation causes hyperglycemia, insulin resistance, and high relative activity, thus increasing energy expenditures that resemble a diabetic phenotype. These analyses suggest that the SLC2A12 gene not only functionally compensates insulin response for SLC2A4 loss but also affects daily physical behavior and basal metabolic rate during bird evolution, highlighting that older genes retain a higher level of functional diversification.

Universal nomenclature for oxytocin-vasotocin ligand and receptor families

Oxytocin (OXT; hereafter OT) and arginine vasopressin or vasotocin (AVP or VT; hereafter VT) are neurotransmitter ligands that function through specific receptors to control diverse functions^1,2. Here we performed genomic analyses on 35 species that span all major vertebrate lineages, including newly generated high-contiguity assemblies from the Vertebrate Genomes Project^3,4. Our findings support the claim⁵ that OT (also known as OXT) and VT (also known as AVP) are adjacent paralogous genes that have resulted from a local duplication, which we infer was through DNA transposable elements near the origin of vertebrates and in which VT retained more of the parental sequence. We identified six major oxytocin-vasotocin receptors among vertebrates. We propose that all six of these receptors arose from a single receptor that was shared with the common ancestor of invertebrates, through a combination of whole-genome and large segmental duplications. We propose a universal nomenclature based on evolutionary relationships for the genes that encode these receptors, in which the genes are given the same orthologous names across vertebrates and paralogous names relative to each other. This nomenclature avoids confusion due to differential naming in the pre-genomic era and incomplete genome assemblies, furthers our understanding of the evolution of these genes, aids in the translation of findings across species and serves as a model for other gene families.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

Evidence for the loss of plasminogen receptor KT gene in chicken

The loss of conserved genes has the potential to alter phenotypes drastically. Screening of vertebrate genomes for lineage-specific gene loss events has identified numerous natural knockouts associated with specific phenotypes. We provide evidence for the loss of a multi-exonic plasminogen receptor KT (PLGRKT) protein-encoding gene located on the Z chromosome in chicken. Exons 1 and 2 are entirely missing; remnants of exon 3 and a mostly intact exon 4 are identified in an assembly gap-free region in chicken with conserved synteny across species and verified using transcriptome and genome sequencing. PLGRKT gene disrupting changes are present in representative species from all five galliform families. In contrast to this, the presence of an intact transcriptionally active PLGRKT gene in species such as mallard, swan goose, and Anolis lizard suggests that gene loss occurred in the galliform lineage sometime between 68 and 80 Mya. The presence of galliform specific chicken repeat 1 (CR1) insertion at the erstwhile exon 2 of PLGRKT gene suggests repeat insertion-mediated loss. However, at least nine other independent PLGRKT coding frame disrupting changes in other bird species are supported by genome sequencing and indicate a role for relaxed purifying selection before CR1 insertion. The recurrent loss of a conserved gene with a role in the regulation of macrophage migration, efferocytosis, and blood coagulation is intriguing. Hence, we propose potential candidate genes that might be compensating the function of PLGRKT based on the presence of a C-terminal lysine residue, transmembrane domains, and gene expression patterns.

Chromonomer: A Tool Set for Repairing and Enhancing Assembled Genomes Through Integration of Genetic Maps and Conserved Synteny

The pace of the sequencing and computational assembly of novel reference genomes is accelerating. Though DNA sequencing technologies and assembly software tools continue to improve, biological features of genomes such as repetitive sequence as well as molecular artifacts that often accompany sequencing library preparation can lead to fragmented or chimeric assemblies. If left uncorrected, defects like these trammel progress on understanding genome structure and function, or worse, positively mislead this research. Fortunately, integration of additional, independent streams of information, such as a marker-dense genetic map and conserved orthologous gene order from related taxa, can be used to scaffold together unlinked, disordered fragments and to restructure a reference genome where it is incorrectly joined. We present a tool set for automating these processes, one that additionally tracks any changes to the assembly and to the genetic map, and which allows the user to scrutinize these changes with the help of web-based, graphical visualizations. Chromonomer takes a user-defined reference genome, a map of genetic markers, and, optionally, conserved synteny information to construct an improved reference genome of chromosome models: a “chromonome”. We demonstrate Chromonomer’s performance on genome assemblies and genetic maps that have disparate characteristics and levels of quality.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Prioritizing sequence variants in conserved non-coding elements in the chicken genome using chCADD

The availability of genomes for many species has advanced our understanding of the non-protein-coding fraction of the genome. Comparative genomics has proven itself to be an invaluable approach for the systematic, genome-wide identification of conserved non-protein-coding elements (CNEs). However, for many non-mammalian model species, including chicken, our capability to interpret the functional importance of variants overlapping CNEs has been limited by current genomic annotations, which rely on a single information type (e.g. conservation). We here studied CNEs in chicken using a combination of population genomics and comparative genomics. To investigate the functional importance of variants found in CNEs we develop a ch(icken) Combined Annotation-Dependent Depletion (chCADD) model, a variant effect prediction tool first introduced for humans and later on for mouse and pig. We show that 73 Mb of the chicken genome has been conserved across more than 280 million years of vertebrate evolution. The vast majority of the conserved elements are in non-protein-coding regions, which display SNP densities and allele frequency distributions characteristic of genomic regions constrained by purifying selection. By annotating SNPs with the chCADD score we are able to pinpoint specific subregions of the CNEs to be of higher functional importance, as supported by SNPs found in these subregions are associated with known disease genes in humans, mice, and rats. Taken together, our findings indicate that CNEs harbor variants of functional significance that should be object of further investigation along with protein-coding mutations. We therefore anticipate chCADD to be of great use to the scientific community and breeding companies in future functional studies in chicken.

Genome size evolution: towards new model systems for old questions

Genome size (GS) variation is a fundamental biological characteristic; however, its evolutionary causes and consequences are the topic of ongoing debate. Whether GS is a neutral trait or one subject to selective pressures, and how strong these selective pressures are, may remain open questions. Fundamentally, the genomic sequences responsible for this variation directly impact the potential evolutionary outcomes and, equally, are the targets of different evolutionary pressures. For example, duplications and deletions of genic regions (large or small) can have immediate and drastic phenotypic effects, while an expansion or contraction of non-coding DNA is less likely to cause catastrophic phenotypic effects. However, in the long term, the accumulation or deletion of ncDNA is likely to have larger effects. Modern sequencing technologies are allowing for the dissection of these proximate causes, but a combination of these new technologies with more traditional evolutionary experiments and approaches could revolutionize this debate and potentially resolve many of these arguments. Here, I discuss an ambitious way forward for GS research, putting it in context of historical debates, theories and sometimes contradictory evidence, and highlighting the promise of combining new sequencing technologies and analytical developments with more traditional experimental evolution approaches.

New genome assembly of the barn owl (Tyto alba alba)

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single-Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

ZEBrA: Zebra finch Expression Brain Atlas-A resource for comparative molecular neuroanatomy and brain evolution studies

An in-depth understanding of the genetics and evolution of brain function and behavior requires a detailed mapping of gene expression in functional brain circuits across major vertebrate clades. Here we present the Zebra finch Expression Brain Atlas (ZEBrA; www.zebrafinchatlas.org, RRID: SCR_012988), a web-based resource that maps the expression of genes linked to a broad range of functions onto the brain of zebra finches. ZEBrA is a first of its kind gene expression brain atlas for a bird species and a first for any sauropsid. ZEBrA's >3,200 high-resolution digital images of in situ hybridized sections for ~650 genes (as of June 2019) are presented in alignment with an annotated histological atlas and can be browsed down to cellular resolution. An extensive relational database connects expression patterns to information about gene function, mouse expression patterns and phenotypes, and gene involvement in human diseases and communication disorders. By enabling brain-wide gene expression assessments in a bird, ZEBrA provides important substrates for comparative neuroanatomy and molecular brain evolution studies. ZEBrA also provides unique opportunities for linking genetic pathways to vocal learning and motor control circuits, as well as for novel insights into the molecular basis of sex steroids actions, brain dimorphisms, reproductive and social behaviors, sleep function, and adult neurogenesis, among many fundamental themes.

β2 adrenergic receptors and leptin interplay to decrease food intake in chicken

1. The present study was designed to examine the effects of intracerebroventricular (ICV) injection of different α and [Formula: see text] adrenergic receptor antagonists on leptin-induced hypophagia in broiler chickens.2. The study consisted of six experiments. In all experiments, chickens were deprived of feed for 3 h prior to the ICV injections and thereafter were returned immediately to the individual cages and cumulative feed intake, based on the percentage of body weight, was measured at 30, 60 and 120 min post-injection.3. In experiment 1, leptin (2.5, 5 or 10 µg) were injected in birds. In experiment 2, groups received either control solution, prazosin (10 nmol), leptin (10 µg) or a co-injection of prazosin (10 nmol) and leptin (10 µg). The other experiments were conducted as experiment 2, but instead of prazosine (10 nmol), yohimbine (13 nmol) was used in experiment 3, metoprolol (24 nmol) in experiment 4, ICI 118,551 (5 nmol) in experiment 5 and SR 59230R (5 nmol) in experiment 6 were injected either in a group or in combination with leptin (10 µg).4. The results of this study revealed a dose-dependent hypophagic effect of leptin and, in experiment 5, ICV co-injection of ICI118, 551 (5 nmol) and leptin (10 µg) significantly attenuated this effect (P˂0.5). These results suggest that the hypophagic effect of leptin is probably mediated by β2 adrenergic receptors in chickens.

Comparative Analyses of Vertebrate Gut Microbiomes Reveal Convergence between Birds and Bats

Diet and host phylogeny drive the taxonomic and functional contents of the gut microbiome in mammals, yet it is unknown whether these patterns hold across all vertebrate lineages. Here, we assessed gut microbiomes from ∼900 vertebrate species, including 315 mammals and 491 birds, assessing contributions of diet, phylogeny, and physiology to structuring gut microbiomes. In most nonflying mammals, strong correlations exist between microbial community similarity, host diet, and host phylogenetic distance up to the host order level. In birds, by contrast, gut microbiomes are only very weakly correlated to diet or host phylogeny. Furthermore, while most microbes resident in mammalian guts are present in only a restricted taxonomic range of hosts, most microbes recovered from birds show little evidence of host specificity. Notably, among the mammals, bats host especially bird-like gut microbiomes, with little evidence for correlation to host diet or phylogeny. This suggests that host-gut microbiome phylosymbiosis depends on factors convergently absent in birds and bats, potentially associated with physiological adaptations to flight. Our findings expose major variations in the behavior of these important symbioses in endothermic vertebrates and may signal fundamental evolutionary shifts in the cost/benefit framework of the gut microbiome.IMPORTANCE In this comprehensive survey of microbiomes of >900 species, including 315 mammals and 491 birds, we find a striking convergence of the microbiomes of birds and animals that fly. In nonflying mammals, diet and short-term evolutionary relatedness drive the microbiome, and many microbial species are specific to a particular kind of mammal, but flying mammals and birds break this pattern with many microbes shared across different species, with little correlation either with diet or with relatedness of the hosts. This finding suggests that adaptation to flight breaks long-held relationships between hosts and their microbes.

Evolution of the genes mediating phototransduction in rod and cone photoreceptors

This paper reviews current knowledge of the evolution of the multiple genes encoding proteins that mediate the process of phototransduction in rod and cone photoreceptors of vertebrates. The approach primarily involves molecular phylogenetic analysis of phototransduction protein sequences, combined with analysis of the syntenic arrangement of the genes. At least 35 of these phototransduction genes appear to reside on no more than five paralogons - paralogous regions that each arose from a common ancestral region. Furthermore, it appears that such paralogs arose through quadruplication during the two rounds of genome duplication (2R WGD) that occurred in a chordate ancestor prior to the vertebrate radiation, probably around 600 millions years ago. For several components of the phototransduction cascade, it is shown that distinct isoforms already existed prior to WGD, with the likely implication that separate classes of scotopic and photopic photoreceptor cells had already evolved by that stage. The subsequent quadruplication of the entire genome then permitted the refinement of multiple distinct protein isoforms in rods and cones. A unified picture of the likely pattern and approximate timing of all the important gene duplications is synthesised, and the implications for our understanding of the evolution of rod and cone phototransduction are presented.

Avian Leptin: Bird's-Eye View of the Evolution of Vertebrate Energy-Balance Control

Discovery of the satiety hormone leptin in 1994 and its characterization in mammals provided a key tool to deciphering the complex mechanism governing adipose tissue regulation of appetite and energy expenditure. Surprisingly, despite the perfectly logical notion of an energy-storing tissue announcing the amount of fat stores using leptin signaling, alternate mechanisms were chosen in bird evolution. This conclusion emerged based on the recent discovery and characterization of genuine avian leptin - after it had been assumed missing by some, and erroneously identified by others. Critical evaluation of the past and present indications of the role of leptin in Aves provides a new perspective on the evolution of energy-balance control in vertebrates; proposing a regulation strategy alternative to the adipostat mechanism.

Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies

Background

More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced.

Results

The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected.

Conclusions

High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.

Exploring the molecular basis of neuronal excitability in a vocal learner

BACKGROUND

Vocal learning, the ability to learn to produce vocalizations through imitation, relies on specialized brain circuitry known in songbirds as the song system. While the connectivity and various physiological properties of this system have been characterized, the molecular genetic basis of neuronal excitability in song nuclei remains understudied. We have focused our efforts on examining voltage-gated ion channels to gain insight into electrophysiological and functional features of vocal nuclei. A previous investigation of potassium channel genes in zebra finches (Taeniopygia guttata) revealed evolutionary modifications unique to songbirds, as well as transcriptional specializations in the song system [Lovell PV, Carleton JB, Mello CV. BMC Genomics 14:470 2013]. Here, we expand this approach to sodium, calcium, and chloride channels along with their modulatory subunits using comparative genomics and gene expression analysis encompassing microarrays and in situ hybridization.

RESULTS

We found 23 sodium, 38 calcium, and 33 chloride channel genes (HGNC-based classification) in the zebra finch genome, several of which were previously unannotated. We determined 15 genes are missing relative to mammals, including several genes (CLCAs, BEST2) linked to olfactory transduction. The majority of sodium and calcium but few chloride channels showed differential expression in the song system, among them SCN8A and CACNA1E in the direct motor pathway, and CACNG4 and RYR2 in the anterior forebrain pathway. In several cases, we noted a seemingly coordinated pattern across multiple nuclei (SCN1B, SCN3B, SCN4B, CACNB4) or sparse expression (SCN1A, CACNG5, CACNA1B).

CONCLUSION

The gene families examined are highly conserved between avian and mammalian lineages. Several cases of differential expression likely support high-frequency and burst firing in specific song nuclei, whereas cases of sparse patterns of expression may contribute to the unique electrophysiological signatures of distinct cell populations. These observations lay the groundwork for manipulations to determine how ion channels contribute to the neuronal excitability properties of vocal learning systems.

Measuring the impact of gene prediction on gene loss estimates in Eukaryotes by quantifying falsely inferred absences

In recent years it became clear that in eukaryotic genome evolution gene loss is prevalent over gene gain. However, the absence of genes in an annotated genome is not always equivalent to the loss of genes. Due to sequencing issues, or incorrect gene prediction, genes can be falsely inferred as absent. This implies that loss estimates are overestimated and, more generally, that falsely inferred absences impact genomic comparative studies. However, reliable estimates of how prevalent this issue is are lacking. Here we quantified the impact of gene prediction on gene loss estimates in eukaryotes by analysing 209 phylogenetically diverse eukaryotic organisms and comparing their predicted proteomes to that of their respective six-frame translated genomes. We observe that 4.61% of domains per species were falsely inferred to be absent for Pfam domains predicted to have been present in the last eukaryotic common ancestor. Between phylogenetically different categories this estimate varies substantially: for clade-specific loss (ancestral loss) we found 1.30% and for species-specific loss 16.88% to be falsely inferred as absent. For BUSCO 1-to-1 orthologous families, 18.30% were falsely inferred to be absent. Finally, we showed that falsely inferred absences indeed impact loss estimates, with the number of losses decreasing by 11.78%. Our work strengthens the increasing number of studies showing that gene loss is an important factor in eukaryotic genome evolution. However, while we demonstrate that on average inferring gene absences from predicted proteomes is reliable, caution is warranted when inferring species-specific absences.

Transcripts with systematic nucleotide deletion of 1-12 nucleotide in human mitochondrion suggest potential non-canonical transcription

Raw transcriptomic data contain numerous RNA reads whose homology with template DNA doesn't match canonical transcription. Transcriptome analyses usually ignore such noncanonical RNA reads. Here, analyses search for noncanonical mitochondrial RNAs systematically deleting 1 to 12 nucleotides after each transcribed nucleotide triplet, producing deletion-RNAs (delRNAs). We detected delRNAs in the human whole cell and purified mitochondrial transcriptomes, and in Genbank's human EST database corresponding to systematic deletions of 1 to 12 nucleotides after each transcribed trinucleotide. DelRNAs detected in both transcriptomes mapped along with 55.63% of the EST delRNAs. A bias exists for delRNAs covering identical mitogenomic regions in both transcriptomic and EST datasets. Among 227 delRNAs detected in these 3 datasets, 81.1% and 8.4% of delRNAs were mapped on mitochondrial coding and hypervariable region 2 of dloop. Del-transcription analyses of GenBank's EST database confirm observations from whole cell and purified mitochondrial transcriptomes, eliminating the possibility that detected delRNAs are false positives matches, cytosolic DNA/RNA nuclear contamination or sequencing artefacts. These detected delRNAs are enriched in frameshift-inducing homopolymers and are poor in frameshift-preventing circular code codons (a set of 20 codons which regulate reading frame detection, over- and underrepresented in coding and other frames of genes, respectively) suggesting a motif-based regulation of non-canonical transcription. These findings show that rare non-canonical transcripts exist. Such non canonical del-transcription does increases mitochondrial coding potential and non-coding regulation of intracellular mechanisms, and could explain the dark DNA conundrum.

Origin and evolution of the specialized forms of proteasomes involved in antigen presentation

Proteasomes are a multi-subunit protease complex that produces peptides bound by major histocompatibility complex (MHC) class I molecules. Phylogenetic studies indicate that two specialized forms of proteasomes, immunoproteasomes and thymoproteasomes, and the proteasome activator PA28αβ emerged in a common ancestor of jawed vertebrates which acquired adaptive immunity based on the MHC, T cell receptors, and B cell receptors ~ 500 million years ago. Comparative genomics studies now provide strong evidence that the genes coding for the immunoproteasome subunits emerged by genome-wide duplication. On the other hand, the gene encoding the thymoproteasome subunit β5t emerged by tandem duplication from the gene coding for the β5 subunit. Strikingly, birds lack immunoproteasomes, thymoproteasomes, and the proteasome activator PA28αβ, raising an interesting question of whether they have evolved any compensatory mechanisms.

Origin and diversification of the plasminogen activation system among chordates

Background

The plasminogen (PLG) activation system is composed by a series of serine proteases, inhibitors and several binding proteins, which together control the temporal and spatial generation of the active serine protease plasmin. As this proteolytic system plays a central role in human physiology and pathophysiology it has been extensively studied in mammals. The serine proteases of this system are believed to originate from an ancestral gene by gene duplications followed by domain gains and deletions. However, the identification of ancestral forms in primitive chordates supporting these theories remains elusive. In addition, evolutionary studies of the non-proteolytic members of this system are scarce.

Results

Our phylogenetic analyses place lamprey PLG at the root of the vertebrate PLG-group, while lamprey PLG-related growth factors represent the ancestral forms of the jawed-vertebrate orthologues. Furthermore, we find that the earliest putative orthologue of the PLG activator group is the hyaluronan binding protein 2 (HABP2) gene found in lampreys. The prime plasminogen activators (tissue- and urokinase-type plasminogen activator, tPA and uPA) first occur in cartilaginous fish and phylogenetic analyses confirm that all orthologues identified compose monophyletic groups to their mammalian counterparts. Cartilaginous fishes exhibit the most ancient vitronectin of all vertebrates, while plasminogen activator inhibitor 1 (PAI-1) appears for the first time in cartilaginous fishes and is conserved in the rest of jawed vertebrate clades. PAI-2 appears for the first time in the common ancestor of reptiles and mammals, and represents the latest appearing plasminogen activator inhibitor. Finally, we noted that the urokinase-type plasminogen activator receptor (uPAR)—and three-LU domain containing genes in general—occurred later in evolution and was first detectable after coelacanths.

Conclusions

This study identifies several primitive orthologues of the mammalian plasminogen activation system. These ancestral forms provide clues to the origin and diversification of this enzyme system. Further, the discovery of several members—hitherto unknown in mammals—provide new perspectives on the evolution of this important enzyme system.

Electronic supplementary material

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

Exploring the unmapped DNA and RNA reads in a songbird genome

BACKGROUND

A widely used approach in next-generation sequencing projects is the alignment of reads to a reference genome. Despite methodological and hardware improvements which have enhanced the efficiency and accuracy of alignments, a significant percentage of reads frequently remain unmapped. Usually, unmapped reads are discarded from the analysis process, but significant biological information and insights can be uncovered from these data. We explored the unmapped DNA (normal and bisulfite treated) and RNA sequence reads of the great tit (Parus major) reference genome individual. From the unmapped reads we generated de novo assemblies, after which the generated sequence contigs were aligned to the NCBI non-redundant nucleotide database using BLAST, identifying the closest known matching sequence.

RESULTS

Many of the aligned contigs showed sequence similarity to different bird species and genes that were absent in the great tit reference assembly. Furthermore, there were also contigs that represented known P. major pathogenic species. Most interesting were several species of blood parasites such as Plasmodium and Trypanosoma.

CONCLUSIONS

Our analyses revealed that meaningful biological information can be found when further exploring unmapped reads. For instance, it is possible to discover sequences that are either absent or misassembled in the reference genome, and sequences that indicate infection or sample contamination. In this study we also propose strategies to aid the capture and interpretation of this information from unmapped reads.

Revisiting avian 'missing' genes from de novo assembled transcripts

Background

Argument remains as to whether birds have lost genes compared with mammals and non-avian vertebrates during speciation. High quality-reference gene sets are necessary for precisely evaluating gene gain and loss. It is essential to explore new reference transcripts from large-scale de novo assembled transcriptomes to recover the potential hidden genes in avian genomes.

Results

We explored 196 high quality transcriptomic datasets from five bird species to reconstruct transcripts for the purpose of discovering potential hidden genes in the avian genomes. We constructed a relatively complete and high-quality bird transcript database (1,623,045 transcripts after quality control in five birds) from a large amount of avian transcriptomic data, and found most of the presumed missing genes (83.2%) could be recovered in at least one bird species. Most of these genes have been identified for the first time in birds. Our results demonstrate that 67.94% genes have GC content over 50%, while 2.91% genes are AT-rich (AT% > 60%). In our results, 239 (53.59%) genes had a tissue-specific expression index of more than 0.9 in chicken. The missing genes also have lower Ka/Ks values than average (genome-wide: Ka/Ks = 0.99; missing gene: Ka/Ks = 0.90; t-test = 1.25E-14). Among all presumed missing genes, there were 135 for which we did not find any meaningful orthologues in any of the 5 species studied.

Conclusion

Insufficient reference genome quality is the major reason for wrongly inferring missing genes in birds. Those presumably missing genes often have a very strong tissue-specific expression pattern. We show multi-tissue transcriptomic data from various species are necessary for inferring gene family evolution for species with only draft reference genomes.

Electronic supplementary material

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

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

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

Genome Sequence of Peacock Reveals the Peculiar Case of a Glittering Bird

The unique ornamental features and extreme sexual traits of Peacock have always intrigued scientists and naturalists for centuries. However, the genomic basis of these phenotypes are yet unknown. Here, we report the first genome sequence and comparative analysis of peacock with the high quality genomes of chicken, turkey, duck, flycatcher and zebra finch. Genes involved in early developmental pathways including TGF-β, BMP, and Wnt signaling, which have been shown to be involved in feather patterning, bone morphogenesis, and skeletal muscle development, revealed signs of adaptive evolution and provided useful clues on the phenotypes of peacock. Innate and adaptive immune genes involved in complement system and T-cell response also showed signs of adaptive evolution in peacock suggesting their possible role in building a robust immune system which is consistent with the predictions of the Hamilton–Zuk hypothesis. This study provides novel genomic and evolutionary insights into the molecular understanding toward the phenotypic evolution of Indian peacock.