A multi-platform draft de novo genome assembly and comparative analysis for the Scarlet Macaw (Ara macao)

A Bird's-Eye View of Chromosomic Evolution in the Class Aves

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

Patterns of cytonuclear discordance and divergence between subspecies of the scarlet macaw (Ara macao) in Central America

The scarlet macaw, Ara macao, is a neotropical parrot that contains two described subspecies with broadly discrete geographical distributions. One subspecies, A. m. macao, is found from South America north into southwestern Costa Rica, while the second subspecies, A. m. cyanoptera, is found from eastern Costa Rica north into central Mexico. Our previous research using mitochondrial data to examine phylogeographical divergence across the collective range of these two subspecies concluded that they represent distinct evolutionary entities, with minimal contemporary hybridization between them. Here we further examine phylogenetic relationships and patterns of genetic variation between these two subspecies using a dataset of genetic markers derived from their nuclear genomes. Our analyses show clear nuclear divergence between A. m. macao and A. m. cyanoptera in Central America. Collectively however, samples from this region appear genetically more similar to one another than they do to the examined South American (Brazilian) A. m. macao sample. This observation contradicts our previous assessments based on mitochondrial DNA analyses that A. m. macao in Central and South America represent a single phylogeographical group that is evolutionarily distinct from Central American A. m. cyanoptera. Nonetheless, in agreement with our previous findings, ongoing genetic exchange between the two subspecies appears limited. Rather, our analyses indicate that incomplete lineage sorting is the best supported explanation for cytonuclear discordance within these parrots. High-altitude regions in Central America may act as a reproductive barrier, limiting contemporary hybridization between A. m. macao and A. m. cyanoptera. The phylogeographic complexities of scarlet macaw taxa in this region highlight the need for additional evolutionary examinations of these populations.

Genome survey sequencing and mining of genome-wide microsatellite markers in yellow-billed babbler (Turdoides affinis)

Turdoides affinis is a species of group dwelling old world passerine of family Leiothrichidae. Unavailability of genome-wide sequence and species-specific molecular markers have hindered comprehensive understanding of cooperative breeding behaviour in T. affinis. Therefore, we generated genome-wide microsatellite markers through whole genome short read sequencing of T. affinis. A total of 68.8 gigabytes of paired-end raw data were sequenced containing 195,067,054 reads. Total sequenced reads spanned a coverage of 17X with genome size of 1.18 Gb. A large number of microsatellite markers (265,297) were mined in the T. affinis genome using Krait, and 50 most informative markers were identified and validated further. In-silico PCR results validated 47 markers. Of these 47 markers, five were randomly selected and validated in-vitro in twelve individuals of T. affinis. Genotyping data on these five loci estimated observed heterozygosity (H₀) and expected heterozygosity (H_e) ratios between 0.333 - 0.833 and 0.851-0.906, respectively. Effective allele size ranged from 6.698 to 10.667, inbreeding coefficient of the population ranged from 0.080 to 0.631 and null allele frequency was calculated at 0.055 to 0.303. Polymorphic information content of all the five loci varied between 0.850 and 0.906. Probabilities of exclusion and identity across 5 loci was estimated to be 0.95 and 0.0036, respectively. All the loci showed significant adherence to Hardy-Weinberg equilibrium. The microsatellite markers reported in this study will facilitate future population genetics studies on T. affinis and other congeneric species.

Improved genome assembly provides new insights into the environmental adaptation of the American cockroach, Periplaneta americana

The synanthropic pest and a model organism for entomological research, American cockroach, Periplaneta americana (Linnaeus), can survive in unfavorable environments for humans. To investigate the genetic mechanisms of success in environmental adaptation of P. americana, we de novo reassembled its whole genome based on next-generation sequencing and PacBio sequencing. The final genome reassembly consisted of approximately 3.34 Gb with scaffold N50 of 465.51 Kb. The completeness (95.4%) of the complete genome was evaluated with single-copy orthologous genes using BUSCO. We identified 18,618 protein-coding genes, 16,443 (88.32%) of which were well supported by public protein databases. We identified 482.04 Mb (approximately 14.45%) repeat elements, 1,385,093 perfect microsatellites simple sequence repeats in P. americana genome, which was higher than other four Blattaria insects. Comparative genomics analysis revealed obvious expansion in the gene families associated with chemoreception (olfactory receptors, gustatory receptors, ionotropic glutamate receptors, chemosensory protein, and sensory neuron membrane protein), which provided the necessary information for functional characterization of the chemosensory receptors of P. americana, with potential for new or refined applications of semiochemicals-based control of this pest insect. Similarly, gene families (cytochrome P450s, carboxyl/choline esterases, and UDP-glycosyl-transferases) encoding receptors for bitter or toxic substances and detoxification enzymes were obviously expanded in P. americana, enabling its ability to detect and detoxify many toxins. Enrichment analysis of positively selected genes in P. americana revealed items associated with metabolic process and catalytic activity, which possibly contributed to the pesticide resistance of P. americana. We also analyzed the homologs to antimicrobial peptide genes reported in the Drosophila genome, and identified two attacins and seven defensins in P. americana. Our data and findings will substantially facilitate molecular studies in P. americana, including elucidation of detoxification mechanisms of xenobiotic, as well as development of new pest management strategies for the control of pests like P. americana.

Exploring genetic variation and population structure in a threatened species, Noturus placidus, with whole-genome sequence data

The Neosho madtom (Noturus placidus) is a small catfish, generally less than 3 inches in length, unique to the Neosho-Spring River system within the Arkansas River Basin. It was federally listed as threatened in 1990, largely due to habitat loss. For conservation efforts, we generated whole-genome sequence data from 10 Neosho madtom individuals originating from 3 geographically separated populations to evaluate genetic diversity and population structure. A Neosho madtom genome was de novo assembled, and genome size and content were assessed. Single nucleotide polymorphisms were assessed from de Bruijn graphs, and via reference alignment with both the channel catfish (Ictalurus punctatus) reference genome and Neosho madtom reference genome. Principal component analysis and structure analysis indicated weak population structure, suggesting fish from the 3 locations represent a single population. Using a novel method, genome-wide conservation and divergence between the Neosho madtom, channel catfish, and zebrafish (Danio rerio) was assessed by pairwise contig alignment, which demonstrated that genes important to embryonic development frequently had conserved sequences. This research in a threatened species with no previously published genomic resources provides novel genetic information to guide current and future conservation efforts and demonstrates that using whole-genome sequencing provides detailed information of population structure and demography using only a limited number of rare and valuable samples.

Positive selection in noncoding genomic regions of vocal learning birds is associated with genes implicated in vocal learning and speech functions in humans

Vocal learning, the ability to imitate sounds from conspecifics and the environment, is a key component of human spoken language and learned song in three independently evolved avian groups—oscine songbirds, parrots, and hummingbirds. Humans and each of these three bird clades exhibit specialized behavioral, neuroanatomical, and brain gene expression convergence related to vocal learning, speech, and song. To understand the evolutionary basis of vocal learning gene specializations and convergence, we searched for and identified accelerated genomic regions (ARs), a marker of positive selection, specific to vocal learning birds. We found avian vocal learner-specific ARs, and they were enriched in noncoding regions near genes with known speech functions or brain gene expression specializations in humans and vocal learning birds, including FOXP2, NEUROD6, ZEB2, and MEF2C, and near genes with major neurodevelopmental functions, including NR2F1, NRP2, and BCL11B. We also found enrichment near the SFARI class S genes associated with syndromic vocal communication forms of autism spectrum disorders. These findings reveal strong candidate noncoding regions near genes for the evolutionary adaptations that distinguish vocal learning species from their close vocal nonlearning relatives and provide further evidence of molecular convergence between birdsong and human spoken language.

Polly Wants a Genome: The Lack of Genetic Testing for Pet Parrot Species

Parrots are considered the third most popular pet species, after dogs and cats, in the United States of America. Popular birds include budgerigars, lovebirds and cockatiels and are known for their plumage and vocal learning abilities. Plumage colour variation remains the main driving force behind breeder selection. Despite the birds' popularity, only two molecular genetic tests-bird sexing and pathogen screening-are commercially available to breeders. For a limited number of species, parentage verification tests are available, but are mainly used in conservation and not for breeding purposes. No plumage colour genotyping test is available for any of the species. Due to the fact that there isn't any commercial plumage genotype screening or parentage verification tests available, breeders mate close relatives to ensure recessive colour alleles are passed to the next generation. This, in turn, leads to inbreeding depression and decreased fertility, lower hatchability and smaller clutch sizes, all important traits in commercial breeding systems. This review highlights the research carried out in the field of pet parrot genomics and points out the areas where future research can make a vital contribution to understanding how parrot breeding can be improved to breed healthy, genetically diverse birds.

New Perspectives on Avian Models for Studies of Basic Aging Processes

Avian models have the potential to elucidate basic cellular and molecular mechanisms underlying the slow aging rates and exceptional longevity typical of this group of vertebrates. To date, most studies of avian aging have focused on relatively few of the phenomena now thought to be intrinsic to the aging process, but primarily on responses to oxidative stress and telomere dynamics. But a variety of whole-animal and cell-based approaches to avian aging and stress resistance have been developed-especially the use of primary cell lines and isolated erythrocytes-which permit other processes to be investigated. In this review, we highlight newer studies using these approaches. We also discuss recent research on age-related changes in neural function in birds in the context of sensory changes relevant to homing and navigation, as well as the maintenance of song. More recently, with the advent of "-omic" methodologies, including whole-genome studies, new approaches have gained momentum for investigating the mechanistic basis of aging in birds. Overall, current research suggests that birds exhibit an enhanced resistance to the detrimental effects of oxidative damage and maintain higher than expected levels of cellular function as they age. There is also evidence that genetic signatures associated with cellular defenses, as well as metabolic and immune function, are enhanced in birds but data are still lacking relative to that available from more conventional model organisms. We are optimistic that continued development of avian models in geroscience, especially under controlled laboratory conditions, will provide novel insights into the exceptional longevity of this animal taxon.

Pre-Columbian transregional captive rearing of Amazonian parrots in the Atacama Desert

The feathers of tropical birds were one of the most significant symbols of economic, social, and sacred status in the pre-Columbian Americas. In the Andes, finely produced clothing and textiles containing multicolored feathers of tropical parrots materialized power, prestige, and distinction and were particularly prized by political and religious elites. Here we report 27 complete or partial remains of macaws and amazon parrots from five archaeological sites in the Atacama Desert of northern Chile to improve our understanding of their taxonomic identity, chronology, cultural context, and mechanisms of acquisition. We conducted a multiproxy archaeometric study that included zooarchaeological analysis, isotopic dietary reconstruction, accelerated mass spectrometry radiocarbon dating, and paleogenomic analysis. The results reveal that during the Late Intermediate Period (1100 to 1450 CE), Atacama oasis communities acquired scarlet macaws (Ara macao) and at least five additional translocated parrot species through vast exchange networks that extended more than 500 km toward the eastern Amazonian tropics. Carbon and nitrogen stable isotopes indicate that Atacama aviculturalists sustained these birds on diets rich in marine bird guano-fertilized maize-based foods. The captive rearing of these colorful, exotic, and charismatic birds served to unambiguously signal relational wealth in a context of emergent intercommunity competition.

The draft genome of the blood pheasant (Ithaginis cruentus): Phylogeny and high-altitude adaptation

The blood pheasant (Ithaginis cruentus), the only species in the genus Ithaginis, lives in an extremely inhospitable high-altitude environment, coping with hypoxia and ultraviolet (UV) radiation. To further investigate the phylogeny of Phasianidae species based on complete genomes and understand the molecular genetic mechanisms of the high-altitude adaptation of the blood pheasant, we de novo assembled and annotated the complete genome of the blood pheasant. The blood pheasant genome size is 1.04 Gb with scaffold N50 of 10.88 Mb. We identified 109.92 Mb (10.62%) repetitive elements, 279,037 perfect microsatellites, and 17,209 protein-coding genes. The phylogenetic tree of Phasianidae based on whole genomes revealed three highly supported major clades with the blood pheasant included in the "erectile clade." Comparative genomics analysis showed that many genes were positively selected in the blood pheasant, which was associated with response to hypoxia and/or UV radiation. More importantly, among these positively selected genes (PSGs) which were related to high-altitude adaptation, sixteen PSGs had blood pheasant-specific missense mutations. Our data and analysis lay solid foundation to the study of Phasianidae phylogeny and provided new insights into the potential adaptation mechanisms to the high altitude employed by the blood pheasant.

Chromosomal Evolution in the Phylogenetic Context: A Remarkable Karyotype Reorganization in Neotropical Parrot Myiopsitta monachus (Psittacidae)

Myiopsitta monachus is a small Neotropical parrot (Psittaciformes: Arini Tribe) from subtropical and temperate regions of South America. It has a diploid chromosome number 2n = 48, different from other members of the Arini Tribe that have usually 70 chromosomes. The species has the lowest 2n within the Arini Tribe. In this study, we combined comparative chromosome painting with probes generated from chromosomes of Gallus gallus and Leucopternis albicollis, and FISH with bacterial artificial chromosomes (BACs) selected from the genome library of G. gallus with the aim to shed light on the dynamics of genome reorganization in M. monachus in the phylogenetic context. The homology maps showed a great number of fissions in macrochromosomes, and many fusions between microchromosomes and fragments of macrochromosomes. Our phylogenetic analysis by Maximum Parsimony agree with molecular data, placing M. monachus in a basal position within the Arini Tribe, together with Amazona aestiva (short tailed species). In M. monachus many chromosome rearrangements were found to represent autopomorphic characters, indicating that after this species split as an independent branch, an intensive karyotype reorganization took place. In addition, our results show that M. monachus probes generated by flow cytometry provide novel cytogenetic tools for the detection of avian chromosome rearrangements, since this species presents breakpoints that have not been described in other species.

Phylogenetic Analysis and Karyotype Evolution in Two Species of Core Gruiformes: Aramides cajaneus and Psophia viridis

Gruiformes is a group with phylogenetic issues. Recent studies based on mitochondrial and genomic DNA have proposed the existence of a core Gruiformes, consisting of five families: Heliornithidae, Aramidae, Gruidae, Psophiidae and Rallidae. Karyotype studies on these species are still scarce, either by conventional staining or molecular cytogenetics. Due to this, this study aimed to analyze the karyotype of two species (Aramides cajaneus and Psophia viridis) belonging to families Rallidae and Psopiidae, respectively, by comparative chromosome painting. The results show that some chromosome rearrangements in this group have different origins, such as the association of GGA5/GGA7 in A. cajaneus, as well as the fission of GGA4p and association GGA6/GGA7, which place P. viridis close to Fulica atra and Gallinula chloropus. In addition, we conclude that the common ancestor of the core Gruiformes maintained the original syntenic groups found in the putative avian ancestral karyotype.

Evolutionary History, Genomic Adaptation to Toxic Diet, and Extinction of the Carolina Parakeet

As the only endemic neotropical parrot to have recently lived in the northern hemisphere, the Carolina parakeet (Conuropsis carolinensis) was an iconic North American bird. The last surviving specimen died in the Cincinnati Zoo in 1918 [1]. The cause of its extinction remains contentious: besides excessive mortality associated to habitat destruction and active hunting, their survival could have been negatively affected by its range having become increasingly patchy [2] or by the exposure to poultry pathogens [3, 4]. In addition, the Carolina parakeet showed a predilection for cockleburs, an herbaceous plant that contains a powerful toxin, carboxyatractyloside, or CAT [5], which did not seem to affect them but made the birds notoriously toxic to most predators [3]. To explore the demographic history of this bird, we generated the complete genomic sequence of a preserved specimen held in a private collection in Espinelves (Girona, Spain), as well as of a close extant relative, Aratinga solstitialis. We identified two non-synonymous genetic changes in two highly conserved proteins known to interact with CAT that could underlie a specific dietary adaptation to this toxin. Our genomic analyses did not reveal evidence of a dramatic past demographic decline in the Carolina parakeet; also, its genome did not exhibit the long runs of homozygosity that are signals of recent inbreeding and are typically found in endangered species. As such, our results suggest its extinction was an abrupt process and thus likely solely attributable to human causes.

A High-Quality Draft Genome Assembly of the Black-Necked Crane (Grus nigricollis) Based on Nanopore Sequencing

The black-necked crane (Grus nigricollis) which inhabits high-altitude areas has the largest body size of the world's 15 crane species, and is classified as threatened by the IUCN. To support further studies on population genetics and genomics, we present a high-quality genome assembly based on both Illumina and nanopore sequencing. In total, 54.59 Gb Illumina short reads and 116.5 Gb nanopore long reads were generated. The 1.23 Gb assembled genome has a high contig N50 of 17.89 Mb, and has a longest contig of 87.83 Mb. The completeness (97.7%) of the draft genome was evaluated with single-copy orthologous genes using BUSCO. We identified 17,789 genes and found that 8.11% of the genome is composed of repetitive elements. In total, 84 of the 2,272 one-to-one orthologous genes were under positive selection in the black-necked crane lineage. SNP-based inference indicated two bottlenecks in the recent demographic trajectories of the black-necked crane. The genome information will contribute to future study of crane evolutionary history and provide new insights into the potential adaptation mechanisms of the black-necked crane to its high-altitude habitat.

The Draft Genome of the Endangered Sichuan Partridge (Arborophila rufipectus) with Evolutionary Implications

The Sichuan partridge (Arborophila rufipectus, Phasianidae, Galliformes) is distributed in south-west China, and classified as endangered grade. To examine the evolution and genomic features of Sichuan partridge, we de novo assembled the Sichuan partridge reference genome. The final draft assembly consisted of approximately 1.09 Gb, and had a scaffold N50 of 4.57 Mb. About 1.94 million heterozygous single-nucleotide polymorphisms (SNPs) were detected, 17,519 protein-coding genes were predicted, and 9.29% of the genome was identified as repetitive elements. A total of 56 olfactory receptor (OR) genes were found in Sichuan partridge, and conserved motifs were detected. Comparisons between the Sichuan partridge genome and chicken genome revealed a conserved genome structure, and phylogenetic analysis demonstrated that Arborophila possessed a basal phylogenetic position within Phasianidae. Gene Ontology (GO) enrichment analysis of positively selected genes (PSGs) in Sichuan partridge showed over-represented GO functions related to environmental adaptation, such as energy metabolism and behavior. Pairwise sequentially Markovian coalescent analysis revealed the recent demographic trajectory for the Sichuan partridge. Our data and findings provide valuable genomic resources not only for studying the evolutionary adaptation, but also for facilitating the long-term conservation and genetic diversity for this endangered species.

Genome-wide analysis sheds light on the high-altitude adaptation of the buff-throated partridge (Tetraophasis szechenyii)

The buff-throated partridge (Tetraophasis szechenyii) is a hypoxia-tolerant bird living in an extremely inhospitable high-altitude environment, which has high ultraviolet (UV) radiation as well as a low oxygen supply when compared with low-altitude areas. To further understand the molecular genetic mechanisms of the high-altitude adaptation of the buff-throated partridges, we de novo assembled the complete genome of the buff-throated partridge. Comparative genomics revealed that positively selected hypoxia-related genes in the buff-throated partridge were distributed in the HIF-1 signaling pathway (map04066), response to hypoxia (GO:0001666), response to oxygen-containing compound (GO:1901700), ATP binding (GO:0005524), and angiogenesis (GO:0001525). Of these positively selected hypoxia-related genes, one positively selected gene (LONP1) had one buff-throated partridge-specific missense mutation which was classified as deleterious by PolyPhen-2. Moreover, positively selected genes in the buff-throated partridge were enriched in cellular response to DNA damage stimulus (corrected P value: 0.028006) and DNA repair (corrected P value: 0.044549), which was related to the increased exposure of the buff-throated partridge to UV radiation. Compared with other avian genomes, the buff-throated partridge showed expansion in genes associated with steroid hormone receptor activity and contractions in genes related to immune and olfactory perception. Furthermore, comparisons between the buff-throated partridge genome and red junglefowl genome revealed a conserved genome structure and provided strong evidence of the sibling relationship between Tetraophasis and Lophophorus. Our data and analysis contributed to the study of Phasianidae evolutionary history and provided new insights into the potential adaptation mechanisms to the high altitude employed by the buff-throated partridge.

In silico identification and characterization of novel microsatellite loci for the Blue-and-yellow Macaw Ara ararauna (Linnaeus, 1758) (Psittaciformes, Psittacidae)

The illegal trade is a major threat to many bird species, and parrots are common victims of this activity. Domestic and international pet markets are interested on different parrot species, such as the Blue-and-yellow Macaw (Ara ararauna). This South American macaw is not globally threatened, but is under protection from over-exploitation. This study aimed to identify and characterize novel microsatellite loci for population and parentage analysis of A. ararauna. Scaffold sequences of Ara macao available in the NCBI database were used for microsatellite searches using MsatCommander software. We tested a total of 28 loci, from which 25 were polymorphic, one was monomorphic, and two did not generated amplification products. For polymorphic loci, the mean number of alleles was 8.24 (4 – 15 alleles per locus), the observed heterozygosity ranged from 0.333 to 0.917, and the expected heterozygosity from 0.353 to 0.890. The paternity exclusion probability and identity probability were highly discriminatory. Thus, these novel microsatellite markers can be useful for population assignment and paternity tests, helping the authorities to manage macaws from the illegal trafficking and control commercial breeders.

Chromosome Painting in Neotropical Long- and Short-Tailed Parrots (Aves, Psittaciformes): Phylogeny and Proposal for a Putative Ancestral Karyotype for Tribe Arini

Most Neotropical Psittacidae have a diploid number of 2n = 70, and a dichotomy in chromosome patterns. Long-tailed species have biarmed macrochromosomes, while short-tailed species have telo/acrocentric macrochromosomes. However, the use of chromosome painting has demonstrated that karyotype evolution in Psittacidae includes a high number of inter/intrachromosomal rearrangements. To determine the phylogeny of long- and short-tailed species, and to propose a putative ancestral karyotype for this group, we constructed homology maps of Pyrrhura frontalis (PFR) and Amazona aestiva (AAE), belonging to the long- and short-tailed groups, respectively. Chromosomes were analyzed by conventional staining and fluorescent in situ hybridization using whole chromosome paints of Gallusgallus and Leucopternis albicollis. Conventional staining showed a karyotype with 2n = 70 in both species, with biarmed macrochromosomes in PFR and telo/acrocentric chromosomes in AAE. Comparison of the results with the putative avian ancestral karyotype (PAK) showed fusions in PFR of PAK1p/PAK4q (PFR1) and PAK6/PAK7 (PFR6) with a paracentric inversion in PFR6. However, in AAE, there was only the fusion between PAK6/7 (AAE7) with a paracentric inversion. Our results indicate that PFR retained a more basal karyotype than long-tailed species previously studied, and AAE a more basal karyotype for Neotropical Psittacidae analyzed so far.

Archaeogenomic evidence from the southwestern US points to a pre-Hispanic scarlet macaw breeding colony

Hundreds of scarlet macaw (Ara macao cyanoptera) skeletons have been recovered from archaeological contexts in the southwestern United States and northwestern Mexico (SW/NW). The location of these skeletons, >1,000 km outside their Neotropical endemic range, has suggested a far-reaching pre-Hispanic acquisition network. Clear evidence for scarlet macaw breeding within this network is only known from the settlement of Paquimé in NW dating between 1250 and 1450 CE. Although some scholars have speculated on the probable existence of earlier breeding centers in the SW/NW region, there has been no supporting evidence. In this study, we performed an ancient DNA analysis of scarlet macaws recovered from archaeological sites in Chaco Canyon and the contemporaneous Mimbres area of New Mexico. All samples were directly radiocarbon dated between 900 and 1200 CE. We reconstructed complete or near-complete mitochondrial genome sequences of 14 scarlet macaws from five different sites. We observed remarkably low genetic diversity in this sample, consistent with breeding of a small founder population translocated outside their natural range. Phylogeographic comparisons of our ancient DNA mitogenomes with mitochondrial sequences from macaws collected during the last 200 years from their endemic Neotropical range identified genetic affinity between the ancient macaws and a single rare haplogroup (Haplo6) observed only among wild macaws in Mexico and northern Guatemala. Our results suggest that people at an undiscovered pre-Hispanic settlement dating between 900 and 1200 CE managed a macaw breeding colony outside their endemic range and distributed these symbolically important birds through the SW.

Karyotype Evolution in Birds: From Conventional Staining to Chromosome Painting

In the last few decades, there have been great efforts to reconstruct the phylogeny of Neoaves based mainly on DNA sequencing. Despite the importance of karyotype data in phylogenetic studies, especially with the advent of fluorescence in situ hybridization (FISH) techniques using different types of probes, the use of chromosomal data to clarify phylogenetic proposals is still minimal. Additionally, comparative chromosome painting in birds is restricted to a few orders, while in mammals, for example, virtually all orders have already been analyzed using this method. Most reports are based on comparisons using Gallus gallus probes, and only a small number of species have been analyzed with more informative sets of probes, such as those from Leucopternis albicollis and Gyps fulvus, which show ancestral macrochromosomes rearranged in alternative patterns. Despite this, it is appropriate to review the available cytogenetic information and possible phylogenetic conclusions. In this report, the authors gather both classical and molecular cytogenetic data and describe some interesting and unique characteristics of karyotype evolution in birds.

Annotated Draft Genome Assemblies for the Northern Bobwhite (Colinus virginianus) and the Scaled Quail (Callipepla squamata) Reveal Disparate Estimates of Modern Genome Diversity and Historic Effective Population Size

Northern bobwhite (Colinus virginianus; hereafter bobwhite) and scaled quail (Callipepla squamata) populations have suffered precipitous declines across most of their US ranges. Illumina-based first- (v1.0) and second- (v2.0) generation draft genome assemblies for the scaled quail and the bobwhite produced N50 scaffold sizes of 1.035 and 2.042 Mb, thereby producing a 45-fold improvement in contiguity over the existing bobwhite assembly, and ≥90% of the assembled genomes were captured within 1313 and 8990 scaffolds, respectively. The scaled quail assembly (v1.0 = 1.045 Gb) was ∼20% smaller than the bobwhite (v2.0 = 1.254 Gb), which was supported by kmer-based estimates of genome size. Nevertheless, estimates of GC content (41.72%; 42.66%), genome-wide repetitive content (10.40%; 10.43%), and MAKER-predicted protein coding genes (17,131; 17,165) were similar for the scaled quail (v1.0) and bobwhite (v2.0) assemblies, respectively. BUSCO analyses utilizing 3023 single-copy orthologs revealed a high level of assembly completeness for the scaled quail (v1.0; 84.8%) and the bobwhite (v2.0; 82.5%), as verified by comparison with well-established avian genomes. We also detected 273 putative segmental duplications in the scaled quail genome (v1.0), and 711 in the bobwhite genome (v2.0), including some that were shared among both species. Autosomal variant prediction revealed ∼2.48 and 4.17 heterozygous variants per kilobase within the scaled quail (v1.0) and bobwhite (v2.0) genomes, respectively, and estimates of historic effective population size were uniformly higher for the bobwhite across all time points in a coalescent model. However, large-scale declines were predicted for both species beginning ∼15-20 KYA.

Ancient horizontal transfers of retrotransposons between birds and ancestors of human pathogenic nematodes

Parasite host switches may trigger disease emergence, but prehistoric host ranges are often unknowable. Lymphatic filariasis and loiasis are major human diseases caused by the insect-borne filarial nematodes Brugia, Wuchereria and Loa. Here we show that the genomes of these nematodes and seven tropical bird lineages exclusively share a novel retrotransposon, AviRTE, resulting from horizontal transfer (HT). AviRTE subfamilies exhibit 83–99% nucleotide identity between genomes, and their phylogenetic distribution, paleobiogeography and invasion times suggest that HTs involved filarial nematodes. The HTs between bird and nematode genomes took place in two pantropical waves, >25–22 million years ago (Myr ago) involving the Brugia/Wuchereria lineage and >20–17 Myr ago involving the Loa lineage. Contrary to the expectation from the mammal-dominated host range of filarial nematodes, we hypothesize that these major human pathogens may have independently evolved from bird endoparasites that formerly infected the global breadth of avian biodiversity.

Lymphatic filariasis and loiasis are diseases caused by insect-borne filarial nematodes. Here, Suh et al. identify a retrotransposon that is present in the genomes of these nematodes and seven tropical bird lineages, indicating two waves of horizontal gene transfer around 17–25 million years ago.

Comparative Cytogenetics of the Congo African Grey Parrot (Psittacus erithacus)

The Congo African grey parrot (Psittacus erithacus, PER) is an endemic species of Central Africa, valued for its intelligence and listed as vulnerable due to poaching and habitat destruction. Improved knowledge about the P. erithacus genome is needed to address key biological questions and conservation of this species. The P. erithacus genome was studied using conventional and molecular cytogenetic approaches including Zoo-FISH. P. erithacus has a 'typical' parrot karyotype with 2n = 62-64 and 8 pairs of macrochromosomes. A distinct feature was a sharp macro-microchromosome boundary. Telomeric sequences were present at all chromosome ends and interstitially in PER2q, the latter coinciding with a C-band. NORs mapped to 4 pairs of microchromosomes which is in contrast to a single NOR in ancestral type avian karyotypes. Zoo-FISH with chicken macrochromosomes GGA1-9 and Z revealed patterns of conserved synteny similar to many other avian groups, though neighboring synteny combinations of GGA6/7, 8/9, and 1/4 were distinctive only to parrots. Overall, P. erithacus shared more Zoo-FISH patterns with neotropical macaws than Australian species such as cockatiel and budgerigar. The observations suggest that Psittaciformes karyotypes have undergone more extensive evolutionary rearrangements compared to the majority of other avian genomes.

Complete mitochondrial genome of Blue-and-yellow Macaw (Ara ararauna): the species morphologically similar to Blue-throated Macaw (Ara glaucogularis)

Six genera of the Arini tribe form morphologically diverse group termed as Macaws. The presence of bare facial area distinguishes them from other members of the tribe. Genera and species of Macaw group differ in pattern of this bare skin as well as in body size. The genus Ara contains eight extant species. Two of them are classified as medium-sized Macaws and other six species are classified as large Macaws. Based on morphological similarities and differences, large Macaw species can be segregated into three pairs according to their plumage coloration. So far, representative mitochondrial genomes were sequenced only for Ara glaucogularis (blue and yellow coloration), Ara macao (predominantly red/scarlet) and Ara militaris (predominantly green) species. Ara ararauna is the second of two Ara species with predominant blue and yellow plumage. Hitherto, close phylogenetic relation of these two species was shown using only ND2 genes or incomplete combined COI/Cytb/ND2 sequences. Therefore, Ara ararauna mitogenome, which was sequenced in this study, will be indispensable to refine the phylogenetic relationships within Macaws group.

Chromosomal Diversity and Karyotype Evolution in South American Macaws (Psittaciformes, Psittacidae)

Most species of macaws, which represent the largest species of Neotropical Psittacidae, characterized by their long tails and exuberant colours, are endangered, mainly because of hunting, illegal trade and habitat destruction. Long tailed species seem to represent a monophyletic group within Psittacidae, supported by cytogenetic data. Hence, these species show karyotypes with predominance of biarmed macrochromosomes, in contrast to short tailed species, with a predominance of acro/telocentric macrochromosomes. Because of their similar karyotypes, it has been proposed that inversions and translocations may be the main types of rearrangements occurring during the evolution of this group. However, only one species of macaw, Ara macao, that has had its genome sequenced was analyzed by means of molecular cytogenetics. Hence, in order to verify the rearrangements, we analyzed the karyotype of two species of macaws, Ara chloropterus and Anodorhynchus hyacinthinus, using cross-species chromosome painting with two different sets of probes from chicken and white hawk. Both intra- and interchromosomal rearrangements were observed. Chicken probes revealed the occurrence of fusions, fissions and inversions in both species, while the probes from white hawk determined the correct breakpoints or chromosome segments involved in the rearrangements. Some of these rearrangements were common for both species of macaws (fission of GGA1 and fusions of GGA1p/GGA4q, GGA6/GGA7 and GGA8/GGA9), while the fissions of GGA 2 and 4p were found only in A. chloropterus. These results confirm that despite apparent chromosomal similarity, macaws have very diverse karyotypes, which differ from each other not only by inversions and translocations as postulated before, but also by fissions and fusions.

Complete mitochondrial genome of Military Macaw (Ara militaris): its comparison with mitogenomes of two other Ara species

The Military Macaw is one of the eight species of the genus Ara. The genus is one of six genera, which form morphologically diverse group termed as Macaws. Parrots of this group differ in body size on demand of the genus and species. Six of Ara species are classified as large Macaws. Based on morphological similarities and differences, these species can be segregated into three pairs according to their plumage coloration. Representative mitochondrial genomes were sequenced only for A. glaucogularis (blue and yellow coloration) and A. macao (predominantly red/scarlet). Ara militaris is one of two predominantly green species and full mitochondrial genome of considered species was sequenced in this study. It's comparison with A. glaucogularis and A. macao mitogenomes revealed higher degree of identity between militaris and macao sequences than between militaris and glaucogularis mtDNAs. Ara militaris mitogenome will be indispensable to refine the phylogenetic relationships within Macaw group.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

De novo genome assembly of the economically important weed horseweed using integrated data from multiple sequencing platforms

Horseweed (Conyza canadensis), a member of the Compositae (Asteraceae) family, was the first broadleaf weed to evolve resistance to glyphosate. Horseweed, one of the most problematic weeds in the world, is a true diploid (2n = 2x = 18), with the smallest genome of any known agricultural weed (335 Mb). Thus, it is an appropriate candidate to help us understand the genetic and genomic bases of weediness. We undertook a draft de novo genome assembly of horseweed by combining data from multiple sequencing platforms (454 GS-FLX, Illumina HiSeq 2000, and PacBio RS) using various libraries with different insertion sizes (approximately 350 bp, 600 bp, 3 kb, and 10 kb) of a Tennessee-accessed, glyphosate-resistant horseweed biotype. From 116.3 Gb (approximately 350× coverage) of data, the genome was assembled into 13,966 scaffolds with 50% of the assembly = 33,561 bp. The assembly covered 92.3% of the genome, including the complete chloroplast genome (approximately 153 kb) and a nearly complete mitochondrial genome (approximately 450 kb in 120 scaffolds). The nuclear genome is composed of 44,592 protein-coding genes. Genome resequencing of seven additional horseweed biotypes was performed. These sequence data were assembled and used to analyze genome variation. Simple sequence repeat and single-nucleotide polymorphisms were surveyed. Genomic patterns were detected that associated with glyphosate-resistant or -susceptible biotypes. The draft genome will be useful to better understand weediness and the evolution of herbicide resistance and to devise new management strategies. The genome will also be useful as another reference genome in the Compositae. To our knowledge, this article represents the first published draft genome of an agricultural weed.

High-coverage sequencing and annotated assemblies of the budgerigar genome

BACKGROUND

Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome.

FINDINGS

We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing.

CONCLUSIONS

Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Tracing the evolution of amniote chromosomes

A great deal of diversity in chromosome number and arrangement is observed across the amniote phylogeny. Understanding how this diversity is generated is important for determining the role of chromosomal rearrangements in generating phenotypic variation and speciation. Gaining this understanding is achieved by reconstructing the ancestral genome arrangement based on comparisons of genome organization of extant species. Ancestral karyotypes for several amniote lineages have been reconstructed, mainly from cross-species chromosome painting data. The availability of anchored whole genome sequences for amniote species has increased the evolutionary depth and confidence of ancestral reconstructions from those made solely from chromosome painting data. Nonetheless, there are still several key lineages where the appropriate data required for ancestral reconstructions is lacking. This review highlights the progress that has been made towards understanding the chromosomal changes that have occurred during amniote evolution and the reconstruction of ancestral karyotypes.

The genome sequence of a widespread apex predator, the golden eagle (Aquila chrysaetos)

Biologists routinely use molecular markers to identify conservation units, to quantify genetic connectivity, to estimate population sizes, and to identify targets of selection. Many imperiled eagle populations require such efforts and would benefit from enhanced genomic resources. We sequenced, assembled, and annotated the first eagle genome using DNA from a male golden eagle (Aquila chrysaetos) captured in western North America. We constructed genomic libraries that were sequenced using Illumina technology and assembled the high-quality data to a depth of ∼40x coverage. The genome assembly includes 2,552 scaffolds >10 Kb and 415 scaffolds >1.2 Mb. We annotated 16,571 genes that are involved in myriad biological processes, including such disparate traits as beak formation and color vision. We also identified repetitive regions spanning 92 Mb (∼6% of the assembly), including LINES, SINES, LTR-RTs and DNA transposons. The mitochondrial genome encompasses 17,332 bp and is ∼91% identical to the Mountain Hawk-Eagle (Nisaetus nipalensis). Finally, the data reveal that several anonymous microsatellites commonly used for population studies are embedded within protein-coding genes and thus may not have evolved in a neutral fashion. Because the genome sequence includes ∼800,000 novel polymorphisms, markers can now be chosen based on their proximity to functional genes involved in migration, carnivory, and other biological processes.

A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Wild populations of northern bobwhites (Colinus virginianus; hereafter bobwhite) have declined across nearly all of their U.S. range, and despite their importance as an experimental wildlife model for ecotoxicology studies, no bobwhite draft genome assembly currently exists. Herein, we present a bobwhite draft de novo genome assembly with annotation, comparative analyses including genome-wide analyses of divergence with the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes, and coalescent modeling to reconstruct the demographic history of the bobwhite for comparison to other birds currently in decline (i.e., scarlet macaw; Ara macao). More than 90% of the assembled bobwhite genome was captured within <40,000 final scaffolds (N50 = 45.4 Kb) despite evidence for approximately 3.22 heterozygous polymorphisms per Kb, and three annotation analyses produced evidence for >14,000 unique genes and proteins. Bobwhite analyses of divergence with the chicken and zebra finch genomes revealed many extremely conserved gene sequences, and evidence for lineage-specific divergence of noncoding regions. Coalescent models for reconstructing the demographic history of the bobwhite and the scarlet macaw provided evidence for population bottlenecks which were temporally coincident with human colonization of the New World, the late Pleistocene collapse of the megafauna, and the last glacial maximum. Demographic trends predicted for the bobwhite and the scarlet macaw also were concordant with how opposing natural selection strategies (i.e., skewness in the r-/K-selection continuum) would be expected to shape genome diversity and the effective population sizes in these species, which is directly relevant to future conservation efforts.

Advancing avian behavioral neuroendocrinology through genomics

Genome technologies are transforming all areas of biology, including the study of hormones, brain and behavior. Annotated reference genome assemblies are rapidly being produced for many avian species. Here we briefly review the basic concepts and tools used in genomics. We then consider how these are informing the study of avian behavioral neuroendocrinology, focusing in particular on lessons from the study of songbirds. We discuss the impact of having a complete "parts list" for an organism; the transformational potential of studying large sets of genes at once instead one gene at a time; the growing recognition that environmental and behavioral signals trigger massive shifts in gene expression in the brain; and the prospects for using comparative genomics to uncover the genetic roots of behavioral variation. Throughout, we identify promising new directions for bolstering the application of genomic information to further advance the study of avian brain and behavior.