Avian Genome Research Articles

Current Approaches to and the Application of Intracytoplasmic Sperm Injection (ICSI) for Avian Genome Editing.

Poultry are one of the most valuable resources for human society. They are also recognized as a powerful experimental animal for basic research on embryogenesis. Demands for the supply of low-allergen eggs and bioreactors have increased with the development of programmable genome editing technology. The CRISPR/Cas9 system has recently been used to produce transgenic animals and various animals in the agricultural industry and has also been successfully adopted for the modification of chicken and quail genomes. In this review, we describe the successful establishment of genome-edited lines combined with germline chimera production systems mediated by primordial germ cells and by viral infection in poultry. The avian intracytoplasmic sperm injection (ICSI) system that we previously established and recent advances in ICSI for genome editing are also summarized.

Comprehensive Comparative Analysis Sheds Light on the Patterns of Microsatellite Distribution across Birds Based on the Chromosome-Level Genomes.

Microsatellites (SSRs) are widely distributed in the genomes of organisms and are an important genetic basis for genome evolution and phenotypic adaptation. Although the distribution patterns of microsatellites have been investigated in many phylogenetic lineages, they remain unclear within the morphologically and physiologically diverse avian clades. Here, based on high-quality chromosome-level genomes, we examined the microsatellite distribution patterns for 53 birds from 16 orders. The results demonstrated that each type of SSR had the same ratio between taxa. For example, the frequency of imperfect SSRs (I-SSRs) was 69.90-84.61%, while perfect SSRs (P-SSRs) were 14.86-28.13% and compound SSRs (C-SSRs) were 0.39-2.24%. Mononucleotide SSRs were dominant for perfect SSRs (32.66-76.48%) in most bird species (98.11%), and A(n) was the most abundant repeat motifs of P-SSRs in all birds (5.42-68.22%). Our study further confirmed that the abundance and diversity of microsatellites were less effected by evolutionary history but its length. The number of P-SSRs decreased with increasing repeat times, and longer P-SSRs motifs had a higher variability coefficient of the repeat copy number and lower diversity, indicating that longer motifs tended to have more stable preferences in avian genomes. We also found that P-SSRs were mainly distributed at the gene ends, and the functional annotation for these genes demonstrated that they were related to signal transduction and cellular process. In conclusion, our research provided avian SSR distribution patterns, which will help to explore the genetic basis for phenotypic diversity in birds.

In silico investigation of uncoupling protein function in avian genomes.

The uncoupling proteins (UCPs) are involved in lipid metabolism and belong to a family of mitochondrial anionic transporters. In poultry, only one UCP homologue has been identified and experimentally shown to be associated with growth, feed conversion ratio, and abdominal fat according to its predominant expression in bird muscles. In endotherm birds, cell metabolic efficiency can be tuned by the rate of mitochondrial coupling. Thus, avUCP may be a key contributor to controlling metabolic rate during particular environmental changes. This study aimed to perform a set of in-silico investigations primarily focused on the structural, biological, and biomimetic functions of avUCP. Thereby, using in silico genome analyses among 8 avian species (chicken, turkey, swallow, manakin, sparrow, wagtail, pigeon, and mallard) and a series of bioinformatic approaches, we provide phylogenetic inference and comparative genomics of avUCPs and investigate whether sequence variation can alter coding sequence characteristics, the protein structure, and its biological features. Complementarily, a combination of literature mining and prediction approaches was also applied to predict the gene networks of avUCP to identify genes, pathways, and biological crosstalk associated with avUCP function. The results showed the evolutionary alteration of UCP proteins in different avian species. Uncoupling proteins in avian species are highly conserved trans membrane proteins as seen by sequence alignment, physio-chemical parameters, and predicted protein structures. Taken together, avUCP has the potential to be considered a functional marker for the identification of cell metabolic state, thermogenesis, and oxidative stress caused by cold, heat, fasting, transfer, and other chemical stimuli stresses in birds. It can also be deduced that avUCP, in migrant or domestic birds, may increase heat stress resistance by reducing fatty acid transport/b-oxidation and thermoregulation alongside antioxidant defense mechanisms. The predicted gene network for avUCP highlighted a cluster of 21 genes involved in response to stress and 28 genes related to lipid metabolism and the proton buffering system. Finally, among 11 enriched pathways, crosstalk of 5 signaling pathways including MAPK, adipocytokine, mTOR, insulin, ErbB, and GnRH was predicted, indicating a possible combination of positive or negative feedback among pathways to regulate avUCP functions. Genetic selection for fast-growing commercial poultry has unintentionally increased susceptibility to many kinds of oxidative stress, and so avUCP could be considered as a potential candidate gene for balancing energy expenditure and reactive oxygen species production, especially in breeding programs. In conclusion, avUCP can be introduced as a pleiotropic gene that requires the contribution of regulatory genes, hormones, pathways, and genetic crosstalk to allow its finely-tuned function.

Dynamic Evolution of Avian RNA Virus Sensors: Repeated Loss of RIG-I and RIPLET

Retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated protein 5 (MDA5) are key RNA virus sensors belonging to the RIG-I-like receptor (RLR) family. The activation of the RLR inflammasome leads to the establishment of antiviral state, mainly through interferon-mediated signaling. The evolutionary dynamics of RLRs has been studied mainly in mammals, where rare cases of RLR gene losses were described. By in silico screening of avian genomes, we previously described two independent disruptions of MDA5 in two bird orders. Here, we extend this analysis to approximately 150 avian genomes and report 16 independent evolutionary events of RIG-I inactivation. Interestingly, in almost all cases, these inactivations are coupled with genetic disruptions of RIPLET/RNF135, an ubiquitin ligase RIG-I regulator. Complete absence of any detectable RIG-I sequences is unique to several galliform species, including the domestic chicken (Gallus gallus). We further aimed to determine compensatory evolution of MDA5 in RIG-I-deficient species. While we were unable to show any specific global pattern of adaptive evolution in RIG-I-deficient species, in galliforms, the analyses of positive selection and surface charge distribution support the hypothesis of some compensatory evolution in MDA5 after RIG-I loss. This work highlights the dynamic nature of evolution in bird RNA virus sensors.

Full-length fruit transcriptomes of southern highbush (Vaccinium sp.) and rabbiteye (V. virgatum Ait.) blueberry

BackgroundBlueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries.ResultsWe generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the ‘Draper’ genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in ‘Suziblue’ and ‘Powderblue’ respectively by combining the reads mapped to Northern Highbush blueberry ‘Draper’ genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both ‘Suziblue’ and ‘Powderblue’ with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in ‘Suziblue’ and ‘Powderblue’, respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant.Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in ‘Powderblue’. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics.ConclusionsA full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality–regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.

Avian surfactant protein (SP)-A2 first arose in an early tetrapod before the divergence of amphibians and gradually lost the collagen domain

The air-liquid interface of the mammalian lung is lined with pulmonary surfactants, a mixture of specific proteins and lipids that serve a dual purpose-enabling air-breathing and protection against pathogens. In mammals, surfactant proteins A (SP-A) and D (SP -D) are involved in innate defence of the lung. Birds seem to lack the SP-D gene, but possess SP-A2, an additional SP-A-like gene. Here we investigated the evolution of the SP-A and SP-D genes using computational gene prediction, homology, simulation modelling and phylogeny with published avian and other vertebrate genomes. PCR was used to confirm the identity and expression of SP-A analogues in various tissue homogenates of zebra finch and turkey. In silico analysis confirmed the absence of SP-D-like genes in all 47 published avian genomes. Zebra finch and turkey SP-A1 and SP-A2 sequences, confirmed by PCR of lung homogenates, were compared with sequenced and in silico predicted vertebrate homologs to construct a phylogenetic tree. The collagen domain of avian SP-A1, especially that of zebra finch, was dramatically shorter than that of mammalian SP-A. Amphibian and reptilian genomes also contain avian-like SP-A2 protein sequences with a collagen domain. NCBI Gnomon-predicted avian and alligator SP-A2 proteins all lacked the collagen domain completely. Both avian SP-A1 and SP-A2 sequences form separate clades, which are most closely related to their closest relatives, the alligators. The C-terminal carbohydrate recognition domain (CRD) of zebra finch SP-A1 was structurally almost identical to that of rat SP-A. In fact, the CRD of SP-A is highly conserved among all the vertebrates. Birds retained a truncated version of mammalian type SP-A1 as well as a non-collagenous C-type lectin, designated SP-A2, while losing the large collagenous SP-D lectin, reflecting their evolutionary trajectory towards a unidirectional respiratory system. In the context of zoonotic infections, how these evolutionary changes affect avian pulmonary surface protection is not clear.

Environmental niche and flight intensity are associated with molecular evolutionary rates in a large avian radiation

BackgroundMetabolic activity and environmental energy are two of the most studied putative drivers of molecular evolutionary rates. Their extensive study, however, has resulted in mixed results and has rarely included the exploration of interactions among various factors impacting molecular evolutionary rates across large clades. Taking the diverse avian family Furnariidae as a case study, we examined the association between several estimates of molecular evolutionary rates with proxies of metabolic demands imposed by flight (wing loading and wing shape) and proxies of environmental energy across the geographic ranges of species (temperature and UV radiation).ResultsWe found weak evidence of a positive effect of environmental and morphological variables on mitochondrial substitution rates. Additionally, we found that temperature and UV radiation interact to explain molecular rates at nucleotide sites affected by selection and population size (non-synonymous substitutions), contrary to the expectation of their impact on sites associated with mutation rates (synonymous substitutions). We also found a negative interaction between wing shape (as described by the hand-wing index) and body mass explaining mitochondrial molecular rates, suggesting molecular signatures of positive selection or reduced population sizes in small-bodied species with greater flight activity.ConclusionsOur results suggest that the demands of flight and environmental energy pose multiple evolutionary pressures on the genome either by driving mutation rates or via their association with natural selection or population size. Data from whole genomes and detailed physiology across taxa will bring a more complete picture of the impact of metabolism, population size, and the environment on avian genome evolution.

Transcript- and annotation-guided genome assembly of the European starling.

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.6Mb 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.

Differential Viral Genome Diversity of Healthy and RSS-Affected Broiler Flocks.

The intestinal virus community contributes to health and disease. Runting and stunting syndrome (RSS) is associated with enteric viruses and leads to economic losses in the poultry industry. However, many viruses that potentially cause this syndrome have also been identified in healthy animals. To determine the difference in the virome of healthy and diseased broilers, samples from 11 healthy and 17 affected broiler flocks were collected at two time points and analyzed by Next-Generation Sequencing. Virus genomes of Parvoviridae, Astroviridae, Picornaviridae, Caliciviridae, Reoviridae, Adenoviridae, Coronaviridae, and Smacoviridae were identified at various days of poultry production. De novo sequence analysis revealed 288 full or partial avian virus genomes, of which 97 belonged to the novel genus Chaphamaparvovirus. This study expands the knowledge of the diversity of enteric viruses in healthy and RSS-affected broiler flocks and questions the association of some viruses with the diseases.

A high-quality assembly reveals genomic characteristics, phylogenetic status, and causal genes for leucism plumage of Indian peafowl.

BackgroundThe dazzling phenotypic characteristics of male Indian peafowl (Pavo cristatus) are attractive both to the female of the species and to humans. However, little is known about the evolution of the phenotype and phylogeny of these birds at the whole-genome level. So far, there are no reports regarding the genetic mechanism of the formation of leucism plumage in this variant of Indian peafowl.ResultsA draft genome of Indian peafowl was assembled, with a genome size of 1.05 Gb (the sequencing depth is 362×), and contig and scaffold N50 were up to 6.2 and 11.4 Mb, respectively. Compared with other birds, Indian peafowl showed changes in terms of metabolism, immunity, and skeletal and feather development, which provided a novel insight into the phenotypic evolution of peafowl, such as the large body size and feather morphologies. Moreover, we determined that the phylogeny of Indian peafowl was more closely linked to turkey than chicken. Specifically, we first identified that PMEL was a potential causal gene leading to the formation of the leucism plumage variant in Indian peafowl.ConclusionsThis study provides an Indian peafowl genome of high quality, as well as a novel understanding of phenotypic evolution and phylogeny of Indian peafowl. These results provide a valuable reference for the study of avian genome evolution. Furthermore, the discovery of the genetic mechanism for the development of leucism plumage is both a breakthrough in the exploration of peafowl plumage and also offers clues and directions for further investigations of the avian plumage coloration and artificial breeding in peafowl.

Nicotinic acetylcholine receptors in a songbird brain.

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that mediate fast synaptic transmission and cell signaling, which contribute to learning, memory, and the execution of motor skills. Birdsong is a complex learned motor skill in songbirds. Although the existence of 15 nAChR subunits has been predicted in the avian genome, their expression patterns and potential contributions to song learning and production have not been comprehensively investigated. Here, we cloned all the 15 nAChR subunits (ChrnA1-10, B2-4, D, and G) from the zebra finch brain and investigated the mRNA expression patterns in the neural pathways responsible for the learning and production of birdsong during a critical period of song learning. Although there were no detectable hybridization signals for ChrnA1, A6, A9, and A10, the other 11 nAChR subunits were uniquely expressed in one or more major subdivisions in the song nuclei of the songbird brain. Of these 11 subunits, ChrnA3-5, A7, and B2 were differentially regulated in the song nuclei compared with the surrounding anatomically related regions. ChrnA5 was upregulated during the critical period of song learning in the lateral magnocellular nucleus of the anterior nidopallium. Furthermore, single-cell RNA sequencing revealed ChrnA7 and B2 to be the major subunits expressed in neurons of the vocal motor nuclei HVC and robust nucleus of the arcopallium, indicating the potential existence of ChrnA7-homomeric and ChrnB2-heteromeric nAChRs in limited cell populations. These results suggest that relatively limited types of nAChR subunits provide functional contributions to song learning and production in songbirds.

De Novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric 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.

Evolution of RNA sensing receptors in birds.

Birds are important hosts for many RNA viruses, including influenza A virus, Newcastle disease virus, West Nile virus and coronaviruses. Innate defense against RNA viruses in birds involves detection of viral RNA by pattern recognition receptors. Several receptors of different classes are involved, such as endosomal toll-like receptors and cytoplasmic retinoic acid–inducible gene I-like receptors, and their downstream adaptor proteins. The function of these receptors and their antagonism by viruses is well established in mammals; however, this has received less attention in birds. These receptors have been characterized in a few bird species, and the completion of avian genomes will permit study of their evolution. For each receptor, functional work has established ligand specificity and activation by viral infection. Engagement of adaptors, regulation by modulators and the supramolecular organization of proteins required for activation are incompletely understood in both mammals and birds. These receptors bind conserved nucleic acid agonists such as single- or double-stranded RNA and generally show purifying selection, particularly the ligand binding regions. However, in birds, these receptors and adaptors differ between species, and between individuals, suggesting that they are under selection for diversification over time. Avian receptors and signalling pathways, like their mammalian counterparts, are targets for antagonism by a variety of viruses, intent on escape from innate immune responses.

Организация некодирующих элементов в геномах птиц

The review focuses on non-coding elements in bird genomes. Like most living organisms, birdscarry a significant number of non-coding sequences in their genome. Technically, these sequences cannot beclassified as genes. They include DNA regions with structural functions and high level of conservatismas well as a considerable number of sequences, including repetitive fractions, whose functions are not yetknown. Apart from a general overview of tandem repeats in avian genomes, the paper provides a special focuson differences in copy number of individual repeats. Interestingly, this number directly correlates witha reduction in genome size in birds. Besides, the review provides a detailed description of how repetitivesequences are involved in structural and functional organisation of centromeric and telomeric regions witha special emphasis given to the organisation of pericentromeric and subtelomeric heterochromatin. It alsodescribes what makes these regions different from centromeric and telomeric sequences per se. The data onthe distribution of repeats in bird genomes and their specific organisation are compared with the level of theirconservatism and their role in maintaining structural and functional organisation of the genome.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity.

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Repeated MDA5 Gene Loss in Birds: An Evolutionary Perspective

Two key cytosolic receptors belonging to the retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) family sense the viral RNA-derived danger signals: RIG-I and melanoma differentiation-associated protein 5 (MDA5). Their activation establishes an antiviral state by downstream signaling that ultimately activates interferon-stimulated genes (ISGs). While in rare cases RIG-I gene loss has been detected in mammalian and avian species, most notably in the chicken, MDA5 pseudogenization has only been detected once in mammals. We have screened over a hundred publicly available avian genome sequences and describe an independent disruption of MDA5 in two unrelated avian lineages, the storks (Ciconiiformes) and the rallids (Gruiformes). The results of our RELAX analysis confirmed the absence of negative selection in the MDA5 pseudogene. In contrast to our prediction, we have shown, using multiple dN/dS-based approaches, that the MDA5 loss does not appear to have resulted in any compensatory evolution in the RIG-I gene, which may partially share its ligand-binding specificity. Together, our results indicate that the MDA5 pseudogenization may have important functional effects on immune responsiveness in these two avian clades.

Genetically engineered birds; pre-CRISPR and CRISPR era†.

Generating biopharmaceuticals in genetically engineered bioreactors continues to reign supreme. Hence, genetically engineered birds have attracted considerable attention from the biopharmaceutical industry. Fairly recent genome engineering methods have made genome manipulation an easy and affordable task. In this review, we first provide a broad overview of the approaches and main impediments ahead of generating efficient and reliable genetically engineered birds, and various factors that affect the fate of a transgene. This section provides an essential background for the rest of the review, in which we discuss and compare different genome manipulation methods in the pre-clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR era in the field of avian genome engineering.

Metagenomic sequencing determines complete infectious bronchitis virus (avian Gammacoronavirus) vaccine strain genomes and associated viromes in chicken clinical samples.

Infectious bronchitis virus (IBV, genus Gammacoronavirus) causes an economically important and highly contagious disease in chicken. Random primed RNA sequencing was applied to two IBV positive clinical samples and one in ovo-passaged virus. The virome of a cloacal swab pool was dominated by IBV (82% of viral reads) allowing de novo assembly of a GI-13 lineage complete genome with 99.95% nucleotide identity to vaccine strain 793B. In addition, substantial read counts (16% of viral reads) allowed the assembly of a near-complete chicken astrovirus genome, while lower read counts identified the presence of chicken calicivirus and avian leucosis virus. Viral reads in a respiratory/intestinal tissue pool were distributed between IBV (22.53%), Sicinivirus (Picornaviridae, 24%), and avian leucosis virus (37.04%). A complete IBV genome with 99.95% nucleotide identity to vaccine strain H120 (lineage GI-1), as well as a near-complete avian leucosis virus genome and a partial Sicinivirus genome were assembled from the tissue sample data. Lower read counts identified chicken calicivirus, Avibirnavirus (infectious bursal disease virus, assembling to 98.85% of segment A and 69.66% of segment B closely related to D3976/1 from Germany, 2017) and avian orthoreovirus, while three avian orthoavulavirus 1 reads confirmed prior real-time RT-PCR result. IBV sequence variation analysis identified both fixed and minor frequency variations in the tissue sample compared to its in ovo-passaged virus. Metagenomic methods allow the determination of complete coronavirus genomes from clinical chicken samples while providing additional insights in RNA virus sequence diversity and coinfecting viruses potentially contributing to pathogenicity.Supplementary InformationThe online version contains supplementary material available at 10.1007/s11262-021-01872-7.

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

BackgroundApproximately 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.ResultsOn 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.ConclusionsA 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.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.