Search

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., (© 2023. The Author(s).)

A species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole-genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population-specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process., (© 2023 John Wiley & Sons Ltd.)

Introduction: Highly pathogenic avian influenza (HPAI) viruses, such as H5N1, continue to pose a serious threat to animal agriculture, wildlife and to public health. Controlling and mitigating this disease in domestic birds requires a better understanding of what makes some species highly susceptible (such as turkey and chicken) while others are highly resistant (such as pigeon and goose). Susceptibility to H5N1 varies both with species and strain; for example, species that are tolerant of most H5N1 strains, such as crows and ducks, have shown high mortality to emerging strains in recent years. Therefore, in this study we aimed to examine and compare the response of these six species, to low pathogenic avian influenza (H9N2) and two strains of H5N1 with differing virulence (clade 2.2 and clade 2.3.2.1) to determine how susceptible and tolerant species respond to HPAI challenge., Methods: Birds were challenged in infection trials and samples (brain, ileum and lung) were collected at three time points post infection. The transcriptomic response of birds was examined using a comparative approach, revealing several important discoveries., Results: We found that susceptible birds had high viral loads and strong neuro-inflammatory response in the brain, which may explain the neurological symptoms and high mortality rates exhibited following H5N1 infection. We discovered differential regulation of genes associated with nerve function in the lung and ileum, with stronger differential regulation in resistant species. This has intriguing implications for the transmission of the virus to the central nervous system (CNS) and may also indicate neuro-immune involvement at the mucosal surfaces. Additionally, we identified delayed timing of the immune response in ducks and crows following infection with the more deadly H5N1 strain, which may account for the higher mortality in these species caused by this strain. Lastly, we identified candidate genes with potential roles in susceptibility/resistance which provide excellent targets for future research., Discussion: This study has helped elucidate the responses underlying susceptibility to H5N1 influenza in avian species, which will be critical in developing sustainable strategies for future control of HPAI in domestic poultry., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Morris, Mishra, Raut, Gaunt, Borowska, Kuo, Wang, Vijayakumar, Chingtham, Dutta, Baillie, Digard, Vervelde, Burt and Smith.)

Background: The Australian black swan (Cygnus atratus) is an iconic species with contrasting plumage to that of the closely related northern hemisphere white swans. The relative geographic isolation of the black swan may have resulted in a limited immune repertoire and increased susceptibility to infectious diseases, notably infectious diseases from which Australia has been largely shielded. Unlike mallard ducks and the mute swan (Cygnus olor), the black swan is extremely sensitive to highly pathogenic avian influenza. Understanding this susceptibility has been impaired by the absence of any available swan genome and transcriptome information., Results: Here, we generate the first chromosome-length black and mute swan genomes annotated with transcriptome data, all using long-read based pipelines generated for vertebrate species. We use these genomes and transcriptomes to show that unlike other wild waterfowl, black swans lack an expanded immune gene repertoire, lack a key viral pattern-recognition receptor in endothelial cells and mount a poorly controlled inflammatory response to highly pathogenic avian influenza. We also implicate genetic differences in SLC45A2 gene in the iconic plumage of the black swan., Conclusion: Together, these data suggest that the immune system of the black swan is such that should any avian viral infection become established in its native habitat, the black swan would be in a significant peril., (© 2023. The Author(s).)

Marek's Disease (MD) has a significant impact on both the global poultry economy and animal welfare. The disease pathology can include neurological damage and tumour formation. Sexual dimorphism in immunity and known higher susceptibility of females to MD makes the chicken Z chromosome (GGZ) a particularly attractive target to study the chicken MD response. Previously, we used a Hy-Line F 6 population from a full-sib advanced intercross line to map MD QTL regions (QTLRs) on all chicken autosomes. Here, we mapped MD QTLRs on GGZ in the previously utilized F 6 population with individual genotypes and phenotypes, and in eight elite commercial egg production lines with daughter-tested sires and selective DNA pooling (SDP). Four MD QTLRs were found from each analysis. Some of these QTLRs overlap regions from previous reports. All QTLRs were tested by individuals from the same eight lines used in the SDP and genotyped with markers located within and around the QTLRs. All QTLRs were confirmed. The results exemplify the complexity of MD resistance in chickens and the complex distribution of p -values and Linkage Disequilibrium (LD) pattern and their effect on localization of the causative elements. Considering the fragments and interdigitated LD blocks while using LD to aid localization of causative elements, one must look beyond the non-significant markers, for possible distant markers and blocks in high LD with the significant block. The QTLRs found here may explain at least part of the gender differences in MD tolerance, and provide targets for mitigating the effects of MD.

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long-read transcriptome of an Australian-sourced European starling (S. vulgaris vAU), and a second, North American, short-read genome assembly (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterization. S. vulgaris vAU combined 10× genomics linked-reads, low-coverage Nanopore sequencing, and PacBio Iso-Seq full-length transcript scaffolding to generate a 1050 Mb assembly on 6222 scaffolds (7.6 Mb scaffold N50, 94.6% busco completeness). Further scaffolding against the high-quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Species-specific transcript mapping and gene annotation revealed good gene-level assembly and high functional completeness. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso-Seq transcript data and complementary homology-based annotation of sequential assembly steps (assessed using a new tool, saaga) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counterintuitive behaviour in traditional busco metrics, and present buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base-calling quality. This work expands our knowledge of avian genomes and the available toolkit for assessing and improving genome quality. The new genomic resources presented will facilitate further global genomic and transcriptomic analysis on this ecologically important species., (© 2022 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

OASs play critical roles in immune response against virus infection by polymerizing ATP into 2-5As, which initiate the classical OAS/RNase L pathway and induce degradation of viral RNA. OAS members are functionally diverged in four known innate immune pathways (OAS/RNase L, OASL/IRF7, OASL/RIG-I, and OASL/cGAS), but how they functionally diverged is unclear. Here, we focus on evolutionary patterns and explore the link between evolutionary processes and functional divergence of Tetrapod OAS1. We show that Palaeognathae and Primate OAS1 genes are conserved in genomic and protein structures but differ in function. The former (i.e., ostrich) efficiently synthesized long 2-5A and activated RNase L, while the latter (i.e., human) synthesized short 2-5A and did not activate RNase L. We predicted and verified that two in-frame indels and one positively selected site in the active site pocket contributed to the functional divergence of Palaeognathae and Primate OAS1. Moreover, we discovered and validated that an in-frame indel in the C-terminus of Palaeognathae OAS1 affected the binding affinity of dsRNA and enzymatic activity, and contributed to the functional divergence of Palaeognathae OAS1 proteins. Our findings unravel the molecular mechanism for functional divergence and give insights into the emergence of novel functions in Tetrapod OAS1., (© 2021. Science China Press and Springer-Verlag GmbH Germany, part of Springer Nature.)

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., (© 2021. The Author(s).)

Avian Leukosis Virus subgroup E (ALVE) integrations are endogenous retroviral elements found in the chicken genome. The presence of ALVE has been reported to have negative impacts on multiple traits, including egg production and body weight. The recent development of rapid, inexpensive and specific ALVE detection methods has facilitated their characterization in elite commercial egg production lines across multiple generations. The presence of 20 ALVE was examined in 8 elite lines, from 3 different breeds. Seventeen of these ALVE (85%) were informative and found to be segregating in at least one of the lines. To test for an association between specific ALVE inserts and traits, a large genotype by phenotype study was undertaken. Genotypes were obtained for 500 to 1500 males per line, and the phenotypes used were sire-daughter averages. Phenotype data were analyzed by line with a linear model that included the effects of generation, ALVE genotype and their interaction. If genotype effect was significant, the number of ALVE copies was fitted as a regression to estimate additive ALVE gene substitution effect. Significant associations between the presence of specific ALVE inserts and 18 commercially relevant performance and egg quality traits, including egg production, egg weight and albumen height, were observed. When an ALVE was segregating in more than one line, these associations did not always have the same impact (negative, positive or none) in each line. It is hypothesized that the presence of ALVE in the chicken genome may influence production traits by 3 mechanisms: viral protein production may modulate the immune system and impact overall production performance (virus effect); insertional mutagenesis caused by viral integration may cause direct gene alterations or affect gene regulation (gene effect); or the integration site may be within or adjacent to a quantitative trait region which impacts a performance trait (linkage disequilibrium, marker effect)., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1-5 Gbp) and idiosyncratic genome features., Results: Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species., Conclusions: Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity 1-4 . Sparse taxon sampling has previously been proposed to confound phylogenetic inference 5 , and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Background: The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise., Results: We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models., Conclusions: Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

Marek's disease (MD) represents a significant global economic and animal welfare issue. Marek's disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify QTL regions (QTLR) influencing resistance to MDV, including an F 6 population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, miRNAs, lncRNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

The annual photoperiod cycle provides the critical environmental cue synchronizing rhythms of life in seasonal habitats. In 1936, Bünning proposed a circadian-based coincidence timer for photoperiodic synchronization in plants. Formal studies support the universality of this so-called coincidence timer, but we lack understanding of the mechanisms involved. Here we show in mammals that long photoperiods induce the circadian transcription factor BMAL2, in the pars tuberalis of the pituitary, and triggers summer biology through the eyes absent/thyrotrophin (EYA3/TSH) pathway. Conversely, long-duration melatonin signals on short photoperiods induce circadian repressors including DEC1, suppressing BMAL2 and the EYA3/TSH pathway, triggering winter biology. These actions are associated with progressive genome-wide changes in chromatin state, elaborating the effect of the circadian coincidence timer. Hence, circadian clock-pituitary epigenetic pathway interactions form the basis of the mammalian coincidence timer mechanism. Our results constitute a blueprint for circadian-based seasonal timekeeping in vertebrates.

Background: Endogenous retroviruses (ERVs) are the remnants of retroviral infections which can elicit prolonged genomic and immunological stress on their host organism. In chickens, endogenous Avian Leukosis Virus subgroup E (ALVE) expression has been associated with reductions in muscle growth rate and egg production, as well as providing the potential for novel recombinant viruses. However, ALVEs can remain in commercial stock due to their incomplete identification and association with desirable traits, such as ALVE21 and slow feathering. The availability of whole genome sequencing (WGS) data facilitates high-throughput identification and characterisation of these retroviral remnants., Results: We have developed obsERVer, a new bioinformatic ERV identification pipeline which can identify ALVEs in WGS data without further sequencing. With this pipeline, 20 ALVEs were identified across eight elite layer lines from Hy-Line International, including four novel integrations and characterisation of a fast feathered phenotypic revertant that still contained ALVE21. These bioinformatically detected sites were subsequently validated using new high-throughput KASP assays, which showed that obsERVer was highly precise and exhibited a 0% false discovery rate. A further fifty-seven diverse chicken WGS datasets were analysed for their ALVE content, identifying a total of 322 integration sites, over 80% of which were novel. Like exogenous ALV, ALVEs show site preference for proximity to protein-coding genes, but also exhibit signs of selection against deleterious integrations within genes., Conclusions: obsERVer is a highly precise and broadly applicable pipeline for identifying retroviral integrations in WGS data. ALVE identification in commercial layers has aided development of high-throughput diagnostic assays which will aid ALVE management, with the aim to eventually eradicate ALVEs from high performance lines. Analysis of non-commercial chicken datasets with obsERVer has revealed broad ALVE diversity and facilitates the study of the biological effects of these ERVs in wild and domesticated populations., Competing Interests: Competing interestsThe authors declare they have no competing interests., (© The Author(s) 2020.)

Background: Dinoflagellates are taxonomically diverse and ecologically important phytoplankton that are ubiquitously present in marine and freshwater environments. Mostly photosynthetic, dinoflagellates provide the basis of aquatic primary production; most taxa are free-living, while some can form symbiotic and parasitic associations with other organisms. However, knowledge of the molecular mechanisms that underpin the adaptation of these organisms to diverse ecological niches is limited by the scarce availability of genomic data, partly due to their large genome sizes estimated up to 250 Gbp. Currently available dinoflagellate genome data are restricted to Symbiodiniaceae (particularly symbionts of reef-building corals) and parasitic lineages, from taxa that have smaller genome size ranges, while genomic information from more diverse free-living species is still lacking., Results: Here, we present two draft diploid genome assemblies of the free-living dinoflagellate Polarella glacialis, isolated from the Arctic and Antarctica. We found that about 68% of the genomes are composed of repetitive sequence, with long terminal repeats likely contributing to intra-species structural divergence and distinct genome sizes (3.0 and 2.7 Gbp). For each genome, guided using full-length transcriptome data, we predicted > 50,000 high-quality protein-coding genes, of which ~40% are in unidirectional gene clusters and ~25% comprise single exons. Multi-genome comparison unveiled genes specific to P. glacialis and a common, putatively bacterial origin of ice-binding domains in cold-adapted dinoflagellates., Conclusions: Our results elucidate how selection acts within the context of a complex genome structure to facilitate local adaptation. Because most dinoflagellate genes are constitutively expressed, Polarella glacialis has enhanced transcriptional responses via unidirectional, tandem duplication of single-exon genes that encode functions critical to survival in cold, low-light polar environments. These genomes provide a foundational reference for future research on dinoflagellate evolution.

Background: The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research., Results: We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1., Conclusions: We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

In macrolecithal species, cryopreservation of the oocyte and zygote is not possible due to the large size and quantity of lipid deposited within the egg. For birds, this signifies that cryopreserving and regenerating a species from frozen cellular material are currently technically unfeasible. Diploid primordial germ cells (PGCs) are a potential means to freeze down the entire genome and reconstitute an avian species from frozen material. Here, we examine the use of genetically engineered (GE) sterile female layer chicken as surrogate hosts for the transplantation of cryopreserved avian PGCs from rare heritage breeds of chicken. We first amplified PGC numbers in culture before cryopreservation and subsequent transplantation into host GE embryos. We found that all hatched offspring from the chimera GE hens were derived from the donor rare heritage breed broiler PGCs, and using cryopreserved semen, we were able to produce pure offspring. Measurement of the mutation rate of PGCs in culture revealed that 2.7 × 10 -10 de novo single-nucleotide variants (SNVs) were generated per cell division, which is comparable with other stem cell lineages. We also found that endogenous avian leukosis virus (ALV) retroviral insertions were not mobilized during in vitro propagation. Taken together, these results show that mutation rates are no higher than normal stem cells, essential if we are to conserve avian breeds. Thus, GE sterile avian surrogate hosts provide a viable platform to conserve and regenerate avian species using cryopreserved PGCs., Competing Interests: The authors declare no competing interest.

Missing in Metastasis (MIM), or Metastasis Suppressor 1 (MTSS1), is a highly conserved protein, which links the plasma membrane to the actin cytoskeleton. MIM has been implicated in various cancers, however, its modes of action remain largely enigmatic. Here, we performed an extensive in silico characterisation of MIM to gain better understanding of its function. We detected previously unappreciated functional motifs including adaptor protein (AP) complex interaction site and a C-helix, pointing to a role in endocytosis and regulation of actin dynamics, respectively. We also identified new functional regions, characterised with phosphorylation sites or distinct hydrophilic properties. Strong negative selection during evolution, yielding high conservation of MIM, has been combined with positive selection at key sites. Interestingly, our analysis of intra-molecular co-evolution revealed potential regulatory hotspots that coincided with reduced potentially pathogenic polymorphisms. We explored databases for the mutations and expression levels of MIM in cancer. Experimentally, we focused on chronic lymphocytic leukaemia (CLL), where MIM showed high overall expression, however, downregulation on poor prognosis samples. Finally, we propose strong conservation of MTSS1 also on the transcriptional level and predict novel transcriptional regulators. Our data highlight important targets for future studies on the role of MIM in different tissues and cancers.

Toll-like receptors (TLRs) are key sensor molecules in vertebrates triggering initial phases of immune responses to pathogens. The avian TLR family typically consists of ten receptors, each adapted to distinct ligands. To understand the complex evolutionary history of each avian TLR, we analyzed all members of the TLR family in the whole genome assemblies and target sequence data of 63 bird species covering all major avian clades. Our results indicate that gene duplication events most probably occurred in TLR1 before synapsids diversified from sauropsids. Unlike mammals, ssRNA-recognizing TLR7 has duplicated independently in several avian taxa, while flagellin-sensing TLR5 has pseudogenized multiple times in bird phylogeny. Our analysis revealed stronger positive, diversifying selection acting in TLR5 and the three-domain TLRs (TLR10 [TLR1A], TLR1 [TLR1B], TLR2A, TLR2B, TLR4) that face the extracellular space and bind complex ligands than in single-domain TLR15 and endosomal TLRs (TLR3, TLR7, TLR21). In total, 84 out of 306 positively selected sites were predicted to harbor substitutions dramatically changing the amino acid physicochemical properties. Furthermore, 105 positively selected sites were located in the known functionally relevant TLR regions. We found evidence for convergent evolution acting between birds and mammals at 54 of these sites. Our comparative study provides a comprehensive insight into the evolution of avian TLR genetic variability. Besides describing the history of avian TLR gene gain and gene loss, we also identified candidate positions in the receptors that have been likely shaped by direct molecular host-pathogen coevolutionary interactions and most probably play key functional roles in birds., (© The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Host cells develop the OAS/RNase L [2'-5'-oligoadenylate synthetase (OAS)/ribonuclease L] system to degrade cellular and viral RNA, and/or the OASL/RIG-I (2'-5'-OAS like/retinoic acid inducible protein I) system to enhance RIG-I-mediated IFN induction, thus providing the first line of defense against viral infection. The 2'-5'-OAS-like (OASL) protein may activate the OAS/RNase L system using its typical OAS-like domain (OLD) or mimic the K63-linked pUb to enhance antiviral activity of the OASL/RIG-I system using its two tandem ubiquitin-like domains (UBLs). We first describe that divergent avian (duck and ostrich) OASL inhibit the replication of a broad range of RNA viruses by activating and magnifying the OAS/RNase L pathway in a UBL-dependent manner. This is in sharp contrast to mammalian enzymatic OASL, which activates and magnifies the OAS/RNase L pathway in a UBL-independent manner, similar to 2'-5'-oligoadenylate synthetase 1 (OAS1). We further show that both avian and mammalian OASL can reversibly exchange to activate and magnify the OAS/RNase L and OASL/RIG-I system by introducing only three key residues, suggesting that ancient OASL possess 2-5A [p x 5'A(2'p5'A) n ; x = 1-3; n ≥ 2] activity and has functionally switched to the OASL/RIG-I pathway recently. Our findings indicate the molecular mechanisms involved in the switching of avian and mammalian OASL molecules to activate and enhance the OAS/RNase L and OASL/RIG-I pathways in response to infection by RNA viruses.

Background: Japanese quail (Coturnix japonica), a recently domesticated poultry species, is important not only as an agricultural product, but also as a model bird species for genetic research. However, most of the biological questions concerning genomics, phylogenetics, and genetics of some important economic traits have not been answered. It is thus necessary to complete a high-quality genome sequence as well as a series of comparative genomics, evolution, and functional studies., Results: Here, we present a quail genome assembly spanning 1.04 Gb with 86.63% of sequences anchored to 30 chromosomes (28 autosomes and 2 sex chromosomes Z/W). Our genomic data have resolved the long-term debate of phylogeny among Perdicinae (Japanese quail), Meleagridinae (turkey), and Phasianinae (chicken). Comparative genomics and functional genomic data found that four candidate genes involved in early maturation had experienced positive selection, and one of them encodes follicle stimulating hormone beta (FSHβ), which is correlated with different FSHβ levels in quail and chicken. We re-sequenced 31 quails (10 wild, 11 egg-type, and 10 meat-type) and identified 18 and 26 candidate selective sweep regions in the egg-type and meat-type lines, respectively. That only one of them is shared between egg-type and meat-type lines suggests that they were subject to an independent selection. We also detected a haplotype on chromosome Z, which was closely linked with maroon/yellow plumage in quail using population resequencing and a genome-wide association study. This haplotype block will be useful for quail breeding programs., Conclusions: This study provided a high-quality quail reference genome, identified quail-specific genes, and resolved quail phylogeny. We have identified genes related to quail early maturation and a marker for plumage color, which is significant for quail breeding. These results will facilitate biological discovery in quails and help us elucidate the evolutionary processes within the Phasianidae family.

Chicken anaemia virus (CAV) is a lymphotropic virus that causes anaemia and immunosuppression in chickens. Previously, we proposed that CAV evades host antiviral responses in vivo by disrupting T-cell signalling, but the precise cellular targets and modes of action remain elusive. In this study, we examined gene expression in Marek's disease virus-transformed chicken T-cell line MSB-1 after infection with CAV using both a custom 5K immune-focused microarray and quantitative real-time PCR at 24, 48 and 72 h post-infection. The data demonstrate an intricate equilibrium between CAV and the host gene expression, displaying subtle but significant modulation of transcripts involved in the T-cell, inflammation and NF-κB signalling cascades. CAV efficiently blocked the induction of type-I interferons and interferon-stimulated genes at 72 h. The cell expression pattern implies that CAV subverts host antiviral responses and that the transformed environment of MSB-1 cells offers an opportunistic advantage for virus growth.

As one of the most economically important species and a unique model organism for biological and medical research, the chicken represents the first non-mammalian amniotic species to have its genome sequenced; and so far, the chicken reference genome represents the best assembled and annotated avian genome. Since the release of the first draft genome sequence, the chicken genome assembly has improved greatly in coverage, contiguity and accuracy owing to the continuous efforts made by the chicken genomics community to generate extensive new data using novel sequencing technologies. Transcriptome sequencing, especially the recent effort to characterise full-length transcripts in chicken tissues, has provided key insights into the complexity of structure and function of the chicken genome. In this article, we review the progress in chicken genome assembly and annotation, and recent advances in comparative genomics in birds. Limitations of current data and plans of research are also discussed.

Correction: After the publication of this work [1] an error was noticed in one of the author surnames. The author name Leif Anderson should be spelt as Leif Andersson.

Increased thyrotrophin-stimulating hormone β (TSHβ) expression in the pars tuberalis is assumed to be an early step in the neuroendocrine mechanism transducing photoperiodic information. The present study aimed to determine the relationship between long-photoperiod (LP) and diurnal TSHβ gene expression in the juvenile chicken by comparing LP-photostimulated birds with groups kept on a short photoperiod (SP) for 1 or 12 days. TSHβ expression increased by 3- and 23-fold after 1 and 12 days of LP-photostimulation both during the day and at night. Under both SP and LP conditions, TSHβ expression was between 3- and 14-fold higher at night than in the day, suggesting that TSHβ expression cycles in a diurnal pattern irrespective of photoperiod. The ratio of DIO2/3 was decreased on LPs, consequent to changes in DIO3 expression, although there was no evidence of any diurnal effect on DIO2 or DIO3 expression. Plasma prolactin concentrations revealed both an effect of LPs and time-of-day. Thus, TSHβ expression changes in a dynamic fashion both diurnally and in response to photoperiod., (© 2017 The Authors. Journal of Neuroendocrinology published by John Wiley & Sons Ltd on behalf of British Society for Neuroendocrinology.)

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.

Background: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences., Results: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences., Conclusions: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

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., (Copyright © 2017 Warren et al.)

The analysis of linkage disequilibrium (LD) underpins the development of effective genotyping technologies, trait mapping and understanding of biological mechanisms such as those driving recombination and the impact of selection. We apply the Malécot-Morton model of LD to create additive LD maps that describe the high-resolution LD landscape of commercial chickens. We investigated LD in chickens (Gallus gallus) at the highest resolution to date for broiler, white egg and brown egg layer commercial lines. There is minimal concordance between breeds of fine-scale LD patterns (correlation coefficient <0.21), and even between discrete broiler lines. Regions of LD breakdown, which may align with recombination hot spots, are enriched near CpG islands and transcription start sites (P<2.2 × 10 -16 ), consistent with recent evidence described in finches, but concordance in hot spot locations between commercial breeds is only marginally greater than random. As in other birds, functional elements in the chicken genome are associated with recombination but, unlike evidence from other bird species, the LD landscape is not stable in the populations studied. The development of optimal genotyping panels for genome-led selection programmes will depend on careful analysis of the LD structure of each line of interest. Further study is required to fully elucidate the mechanisms underlying highly divergent LD patterns found in commercial chickens.

Background: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact., Results: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor., Conclusions: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size.

Organizers are regions of the embryo that can both induce new fates and impart pattern on other regions. So far, surprisingly few organizers have been discovered, considering the number of patterned tissue types generated during development. This may be because their discovery has relied on transplantation and ablation experiments. Here we describe a new approach, using chick embryos, to discover organizers based on a common gene expression signature, and use it to uncover the anterior intestinal portal (AIP) endoderm as a putative heart organizer. We show that the AIP can induce cardiac identity from non-cardiac mesoderm and that it can pattern this by specifying ventricular and suppressing atrial regional identity. We also uncover some of the signals responsible. The method holds promise as a tool to discover other novel organizers acting during development.

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., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Rapid growth in broilers is associated with susceptibility to metabolic disorders such as pulmonary hypertension syndrome (ascites) and sudden death. This study describes a genome search for QTL associated with relative weight of cardio respiratory and metabolically important organs (heart, lungs, liver and gizzard), and hematocrit value in a Brazilian broiler-layer cross. QTL with similar or different effects across sexes were investigated. At 42 days of age after fasted for 6 h, the F2 chickens were weighed and slaughtered. Weights and percentages of the weight relative to BW42 of gizzard, heart, lungs, liver and hematocrit were used in the QTL search. Parental, F1 and F2 individuals were genotyped with 128 genetic markers (127 microsatellites and 1 SNP) covering 22 linkage groups. QTL mapping analyses were carried out using mixed models. A total of 11 genome-wide significant QTL and five suggestive linkages were mapped. Thus, genome-wide significant QTL with similar effects across sexes were mapped to GGA2, 4 and 14 for heart weight, and to GGA2, 8 and 12 for gizzard %. Additionally, five genome-wide significant QTL with different effects across sexes were mapped to GGA 8, 19 and 26 for heart weight; GGA26 for heart % and GGA3 for hematocrit value. Five QTL were detected in chromosomal regions where QTL for similar traits were previously mapped in other F2 chicken populations. Seven novel genome-wide significant QTL are reported here, and 21 positional candidate genes in QTL regions were identified.

Avian pathogens are responsible for major costs to society, both in terms of huge economic losses to the poultry industry and their implications for human health. The health and welfare of millions of birds is under continued threat from many infectious diseases, some of which are increasing in virulence and thus becoming harder to control, such as Marek's disease virus and avian influenza viruses. The current era in animal genomics has seen huge developments in both technologies and resources, which means that researchers have never been in a better position to investigate the genetics of disease resistance and determine the underlying genes/mutations which make birds susceptible or resistant to infection. Avian genomics has reached a point where the biological mechanisms of infectious diseases can be investigated and understood in poultry and other avian species. Knowledge of genes conferring disease resistance can be used in selective breeding programmes or to develop vaccines which help to control the effects of these pathogens, which have such a major impact on birds and humans alike.

Persistent free-running circannual (approximately year-long) rhythms have evolved in animals to regulate hormone cycles, drive metabolic rhythms (including hibernation), and time annual reproduction. Recent studies have defined the photoperiodic input to this rhythm, wherein melatonin acts on thyrotroph cells of the pituitary pars tuberalis (PT), leading to seasonal changes in the control of thyroid hormone metabolism in the hypothalamus. However, seasonal rhythms persist in constant conditions in many species in the absence of a changing photoperiod signal, leading to the generation of circannual cycles. It is not known which cells, tissues, and pathways generate these remarkable long-term rhythmic processes. We show that individual PT thyrotrophs can be in one of two binary states reflecting either a long (EYA3(+)) or short (CHGA(+)) photoperiod, with the relative proportion in each state defining the phase of the circannual cycle. We also show that a morphogenic cycle driven by the PT leads to extensive re-modeling of the PT and hypothalamus over the circannual cycle. We propose that the PT may employ a recapitulated developmental pathway to drive changes in morphology of tissues and cells. Our data are consistent with the hypothesis that the circannual timer may reside within the PT thyrotroph and is encoded by a binary switch timing mechanism, which may regulate the generation of circannual neuroendocrine rhythms, leading to dynamic re-modeling of the hypothalamic interface. In summary, the PT-ventral hypothalamus now appears to be a prime structure involved in long-term rhythm generation., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: EAV-HP is an ancient retrovirus pre-dating Gallus speciation, which continues to circulate in modern chicken populations, and led to the emergence of avian leukosis virus subgroup J causing significant economic losses to the poultry industry. We mapped EAV-HP integration sites in Ethiopian village chickens, a Silkie, Taiwan Country chicken, red junglefowl Gallus gallus and several inbred experimental lines using whole-genome sequence data., Results: An average of 75.22 ± 9.52 integration sites per bird were identified, which collectively group into 279 intervals of which 5 % are common to 90 % of the genomes analysed and are suggestive of pre-domestication integration events. More than a third of intervals are specific to individual genomes, supporting active circulation of EAV-HP in modern chickens. Interval density is correlated with chromosome length (P < 2.31(-6)), and 27 % of intervals are located within 5 kb of a transcript. Functional annotation clustering of genes reveals enrichment for immune-related functions (P < 0.05)., Conclusions: Our results illustrate a non-random distribution of EAV-HP in the genome, emphasising the importance it may have played in the adaptation of the species, and provide a platform from which to extend investigations on the co-evolutionary significance of endogenous retroviral genera with their hosts.

Background: Infectious Bronchitis is a highly contagious respiratory disease which causes tracheal lesions and also affects the reproductive tract and is responsible for large economic losses to the poultry industry every year. This is due to both mortality (either directly provoked by IBV itself or due to subsequent bacterial infection) and lost egg production. The virus is difficult to control by vaccination, so new methods to curb the impact of the disease need to be sought. Here, we seek to identify genes conferring resistance to this coronavirus, which could help in selective breeding programs to rear chickens which do not succumb to the effects of this disease., Methods: Whole genome gene expression microarrays were used to analyse the gene expression differences, which occur upon infection of birds with Infectious Bronchitis Virus (IBV). Tracheal tissue was examined from control and infected birds at 2, 3 and 4 days post-infection in birds known to be either susceptible or resistant to the virus. The host innate immune response was evaluated over these 3 days and differences between the susceptible and resistant lines examined., Results: Genes and biological pathways involved in the early host response to IBV infection were determined andgene expression differences between susceptible and resistant birds were identified. Potential candidate genes for resistance to IBV are highlighted., Conclusions: The early host response to IBV is analysed and potential candidate genes for disease resistance are identified. These putative resistance genes can be used as targets for future genetic and functional studies to prove a causative link with resistance to IBV.

Background: β-defensins and cathelicidins are two families of cationic antimicrobial peptides (AMPs) with a broad range of antimicrobial activities that are key components of the innate immune system. Due to their important roles in host defense against rapidly evolving pathogens, the two gene families provide an ideal system for studying adaptive gene evolution. In this study we performed phylogenetic and selection analyses on β-defensins and cathelicidins from 53 avian species representing 32 orders to examine the evolutionary dynamics of these peptides in birds., Results and Conclusions: Avian β-defensins are found in a gene cluster consisting of 13 subfamiles. Nine of these are conserved as one to one orthologs in all birds, while the others (AvBD1, AvBD3, AvBD7 and AvBD14) are more subject to gene duplication or pseudogenisation events in specific avian lineages. Avian cathelicidins are found in a gene cluster consisting of three subfamilies with species-specific duplications and gene loss. Evidence suggested that the propiece and mature peptide domains of avian cathelicidins are possibly co-evolving in such a way that the cationicity of the mature peptide is partially neutralised by the negative charge of the propiece prior to peptide secretion (further evidence obtained by repeating the analyses on primate cathelicidins). Negative selection (overall mean dN < dS) was detected in most of the gene domains examined, conserving certain amino acid residues that may be functionally crucial for the avian β-defensins and cathelicidins, while episodic positive selection was also involved in driving the diversification of specific codon sites of certain AMPs in avian evolutionary history. These findings have greatly improved our understanding of the molecular evolution of avian AMPs and will be useful to understand their role in the avian innate immune response. Additionally, the large dataset of β-defensin and cathelicidin peptides may also provide a valuable resource for translational research and development of novel antimicrobial agents in the future.