Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes.

BACKGROUND: The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions. RESULTS: We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest. CONCLUSIONS: The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

Genetic Mapping and Biochemical Basis of Yellow Feather Pigmentation in Budgerigars.

Parrot feathers contain red, orange, and yellow polyene pigments called psittacofulvins. Budgerigars are parrots that have been extensively bred for plumage traits during the last century, but the underlying genes are unknown. Here we use genome-wide association mapping and gene-expression analysis to map the Mendelian blue locus, which abolishes yellow pigmentation in the budgerigar. We find that the blue trait maps to a single amino acid substitution (R644W) in an uncharacterized polyketide synthase (MuPKS). When we expressed MuPKS heterologously in yeast, yellow pigments accumulated. Mass spectrometry confirmed that these yellow pigments match those found in feathers. The R644W substitution abolished MuPKS activity. Furthermore, gene-expression data from feathers of different bird species suggest that parrots acquired their colors through regulatory changes that drive high expression of MuPKS in feather epithelia. Our data also help formulate biochemical models that may explain natural color variation in parrots. VIDEO ABSTRACT.

A pioneering study on cytotoxicity in Australian parakeets (Melopsittacus undulates) exposed to tannery effluent.

Waste effluent from the tannery industry is a major source of environmental pollution. Considering that the bird intake of water contaminated with tannery effluent constitutes a potential genotoxic source, especially for birds inhabiting areas closest to tanning industries, the aim of this study is to assess the possible mutagenic effects that the intake may have on Melopsittacus undulatus (Australian parakeet). In order to do so, adult male and female M. undulatus were distributed in two experimental groups: control (drinking water) and TE (5%). After 60 days of exposure, the micronucleus test, as well as tests looking for other nuclear abnormalities in the peripheral blood of the birds were performed. The male and female birds exposed to the pollutant have presented the highest total number of nuclear abnormalities, as well as increased individual abnormalities such as nuclei with symmetrical constricted bi-lobed/bi-nucleated erythrocytes, indented nuclei and micro-lobed nuclei (top)/micro-nuclei (bottom). In addition, the exposure to TE has caused a nuclear variant increase rarely reported in the literature concerning poultry erythrocyte nuclei. The birds exposed to the pollutant have presented the highest frequency of displaced nuclei forming different rotation/displacement angles within the cells. Therefore, the current study confirmed the toxicological potential of TE and was pioneer in showing that male and female M. undulatus exposed to pollutant present the highest frequency of erythrocyte nuclear abnormalities, thus corroborating the initial hypothesis herein presented.

The complete mitochondrial genome sequence of the budgerigar, Melopsittacus undulatus.

Here, we describe the budgie's mitochondrial genome sequence, a resource that can facilitate this parrot's use as a model organism as well as for determining its phylogenetic relatedness to other parrots/Psittaciformes. The estimated total length of the sequence was 18,193 bp. In addition to the to the 13 protein and tRNA and rRNA coding regions, the sequence also includes a duplicated hypervariable region, a feature unique to only a few birds. The two hypervariable regions shared a sequence identity of about 86%.

Misassembly detection using paired-end sequence reads and optical mapping data.

MOTIVATION: A crucial problem in genome assembly is the discovery and correction of misassembly errors in draft genomes. We develop a method called misSEQuel that enhances the quality of draft genomes by identifying misassembly errors and their breakpoints using paired-end sequence reads and optical mapping data. Our method also fulfills the critical need for open source computational methods for analyzing optical mapping data. We apply our method to various assemblies of the loblolly pine, Francisella tularensis, rice and budgerigar genomes. We generated and used stimulated optical mapping data for loblolly pine and F.tularensis and used real optical mapping data for rice and budgerigar. RESULTS: Our results demonstrate that we detect more than 54% of extensively misassembled contigs and more than 60% of locally misassembled contigs in assemblies of F.tularensis and between 31% and 100% of extensively misassembled contigs and between 57% and 73% of locally misassembled contigs in assemblies of loblolly pine. Using the real optical mapping data, we correctly identified 75% of extensively misassembled contigs and 100% of locally misassembled contigs in rice, and 77% of extensively misassembled contigs and 80% of locally misassembled contigs in budgerigar. AVAILABILITY AND IMPLEMENTATION: misSEQuel can be used as a post-processing step in combination with any genome assembler and is freely available at http://www.cs.colostate.edu/seq/.

Mind the gap: upgrading genomes with Pacific Biosciences RS long-read sequencing technology.

Many genomes have been sequenced to high-quality draft status using Sanger capillary electrophoresis and/or newer short-read sequence data and whole genome assembly techniques. However, even the best draft genomes contain gaps and other imperfections due to limitations in the input data and the techniques used to build draft assemblies. Sequencing biases, repetitive genomic features, genomic polymorphism, and other complicating factors all come together to make some regions difficult or impossible to assemble. Traditionally, draft genomes were upgraded to "phase 3 finished" status using time-consuming and expensive Sanger-based manual finishing processes. For more facile assembly and automated finishing of draft genomes, we present here an automated approach to finishing using long-reads from the Pacific Biosciences RS (PacBio) platform. Our algorithm and associated software tool, PBJelly, (publicly available at https://sourceforge.net/projects/pb-jelly/) automates the finishing process using long sequence reads in a reference-guided assembly process. PBJelly also provides "lift-over" co-ordinate tables to easily port existing annotations to the upgraded assembly. Using PBJelly and long PacBio reads, we upgraded the draft genome sequences of a simulated Drosophila melanogaster, the version 2 draft Drosophila pseudoobscura, an assembly of the Assemblathon 2.0 budgerigar dataset, and a preliminary assembly of the Sooty mangabey. With 24× mapped coverage of PacBio long-reads, we addressed 99% of gaps and were able to close 69% and improve 12% of all gaps in D. pseudoobscura. With 4× mapped coverage of PacBio long-reads we saw reads address 63% of gaps in our budgerigar assembly, of which 32% were closed and 63% improved. With 6.8× mapped coverage of mangabey PacBio long-reads we addressed 97% of gaps and closed 66% of addressed gaps and improved 19%. The accuracy of gap closure was validated by comparison to Sanger sequencing on gaps from the original D. pseudoobscura draft assembly and shown to be dependent on initial reference quality.

Endogenous hepadnaviruses in the genome of the budgerigar (Melopsittacus undulatus) and the evolution of avian hepadnaviruses.

Endogenous hepadnaviruses (hepatitis B viruses [HBVs]) were recently discovered in the genomes of passerine birds. We mined six additional avian genomes and discovered multiple copies of endogenous HBVs in the budgerigar (order Psittaciformes), designated eBHBV. A phylogenetic analysis reveals that the endogenous hepadnaviruses are more diverse than their exogenous counterparts and that the endogenous and exogenous hepadnaviruses form distinct lineages even when sampled from the same avian order, indicative of multiple genomic integration events.

Comparative analysis of mineralocorticoid receptor expression among vocal learners (Bengalese finch and budgerigar) and non-vocal learners (quail and ring dove) has implications for the evolution of avian vocal learning.

Mineralocorticoid receptor is the receptor for corticosteroids such as corticosterone or aldosterone. Previously, we found that mineralocorticoid receptor was highly expressed in song nuclei of a songbird, Bengalese finch (Lonchura striata var. domestica). Here, to examine the relationship between mineralocorticoid receptor expression and avian vocal learning, we analyzed mineralocorticoid receptor expression in the developing brain of another vocal learner, budgerigar (Melopsittacus undulatus) and non-vocal learners, quail (Coturnix japonica) and ring dove (Streptopelia capicola). Mineralocorticoid receptor showed vocal control area-related expressions in budgerigars as Bengalese finches, whereas no such mineralocorticoid receptor expressions were seen in the telencephalon of non-vocal learners. Thus, these results suggest the possibility that mineralocorticoid receptor plays a role in vocal development of parrots as songbirds and that the acquisition of mineralocorticoid receptor expression is involved in the evolution of avian vocal learning.