- Open Access
Reconstruction of avian ancestral karyotypes reveals differences in the evolutionary history of macro- and microchromosomes
- Joana Damas†1,
- Jaebum Kim†2,
- Marta Farré1,
- Darren K Griffin†3 and
- Denis M Larkin†1Email authorView ORCID ID profile
© The Author(s). 2018
- Received: 12 January 2018
- Accepted: 18 September 2018
- Published: 5 October 2018
Reconstruction of ancestral karyotypes is critical for our understanding of genome evolution, allowing for the identification of the gross changes that shaped extant genomes. The identification of such changes and their time of occurrence can shed light on the biology of each species, clade and their evolutionary history. However, this is impeded by both the fragmented nature of the majority of genome assemblies and the limitations of the available software to work with them. These limitations are particularly apparent in birds, with only 10 chromosome-level assemblies reported thus far. Algorithmic approaches applied to fragmented genome assemblies can nonetheless help define patterns of chromosomal change in defined taxonomic groups.
Here, we make use of the DESCHRAMBLER algorithm to perform the first large-scale study of ancestral chromosome structure and evolution in birds. This algorithm allows us to reconstruct the overall genome structure of 14 key nodes of avian evolution from the Avian ancestor to the ancestor of the Estrildidae, Thraupidae and Fringillidae families.
Analysis of these reconstructions provides important insights into the variability of rearrangement rates during avian evolution and allows the detection of patterns related to the chromosome distribution of evolutionary breakpoint regions. Moreover, the inclusion of microchromosomes in our reconstructions allows us to provide novel insights into the evolution of these avian chromosomes, specifically.
- Ancestral karyotypes
- Chromosome evolution
- Evolutionary breakpoint regions
- Homologous synteny blocks
Reconstructions of ancestral chromosome structures, utilising traditional cytogenetics or genetic map comparisons, offered the first insights into the evolutionary events that shaped extant animal gross genome organisation and the mechanisms that drive chromosome evolution. Indeed, cytogenetics (‘zoo-FISH’) information proved useful to describe basic ancestral patterns, e.g. in placental mammals  and birds  but has limited sensitivity. Differences in rearrangement rates and chromosome structures revealed from these reconstructions, nonetheless, pointed to variation in levels of genome reshuffling and prevailing types of chromosome changes in the evolution of distinct and related phylogenetic clades . Notable examples include widespread interchromosomal rearrangements in mammals, illustrated by a high variation in chromosome numbers, from six in Indian muntjac (Muntiacus muntjac ) to 102 in Vizcacha rat (Tympanoctomys barrerae ). Contrastingly, in birds, interchromosomal rearrangements are very rare with around two thirds of all species studied so far having similar chromosome numbers and karyotypic patterns . These observations raise questions regarding possible differences in the mechanisms that drive the evolution of these lineages.
The low resolution of cytogenetic methodologies inevitably leads to undetected intrachromosomal rearrangements and limited usefulness of karyotypic comparisons for the reconstruction of older ancestors (e.g. eutherian, avian and amniote ancestors). These restrictions can be overcome by the inclusion of higher-resolution genome sequence maps, which (a) expand the evolutionary depth of ancestral karyotype reconstructions and (b) increase the resolution and often accuracy at which genome rearrangements are identified. The first attempts to reconstruct ancestral karyotypes from sequence data were performed for mammals [6, 7]. These reconstructions allowed the identification of novel genome rearrangements and the detection of varied rates of change in different animal lineages. They also revealed that evolutionary breakpoint regions (EBRs) are often reused in evolution, that EBRs locate in gene-dense regions and that lineage-specific EBRs are usually associated with the location of segmental duplications in mammals [6–8]. These findings demonstrated the importance of genome sequence comparisons for the detection of the overall pattern of chromosomal events that shaped extant genomes and stimulated the development of several algorithms to perform the reconstruction of ancestral karyotypes based on genome sequence data. Most algorithms used for the inference of ancestral karyotypes, such as InferCARs  or ANGES , are optimised for chromosome-level genome assemblies (i.e. one scaffold per chromosome) and their suitability to deal with sub-chromosome (e.g. scaffold-level) genome assemblies is very limited . While there have been many newly sequenced genomes released in the last few years, only a small number of them were assembled to chromosomes and were of use for ancestral karyotype reconstruction . To overcome these limitations, Kim and colleagues developed DESCHRAMBLER . This algorithm, in contrast to those aforementioned, is optimised to generate reconstructed ancestral chromosome fragments (RACFs) using information from both chromosome- and scaffold-level genome assemblies . DESCHRAMBLER has been previously applied to the genomes of 19 mammalian species (12 chromosome-level and seven scaffold-level) to reconstruct the chromosome structure of seven eutherian ancestors. DESCHRAMBLER detected a significant number of intrachromosomal changes that cytogenetic studies could not identify. Moreover, it showed that, amongst the 10 orders studied, the primates had the largest number of chromosomes structurally identical to the Eutherian ancestor (orangutan exhibiting the largest) with chimps displaying more structural changes than humans .
Limited availability of chromosome-level assemblies for Aves has hitherto restricted the study of chromosome evolution in this class. Indeed, to date, the reconstruction of avian ancestral chromosome structures has been based on molecular cytogenetic comparisons  with sequence-based reconstructions limited to several macrochromosomes only and based on just six (four chromosome- and two scaffold-based) assemblies . Thus, the extent to which individual chromosomes (particularly microchromosomes) have remained unchanged or rearranged intrachromosomally remains unknown for most avian lineages. Moreover, the sequence features affecting the stability and dynamics of avian chromosomes during evolution remained under-explored. Nonetheless, the current availability of ~ 60 avian genomes, of which ~ 30 have low enough assembly fragmentation to be suitable for DESCHRAMBLER ancestral genome reconstructions, has the potential to rectify this problem.
In this study, we report the use of DESCHRAMBLER  for the first large-scale study of ancestral chromosome structure and evolution in birds. The large number of avian genome assemblies included in this study and extended sampling of the avian phylogenetic tree allowed the reconstruction of the likely overall genome structure of 14 nodes in avian phylogeny, from the Avian ancestor to the ancestor of Passeriformes (the most species-rich avian order) through to the zebra finch. The analysis of these reconstructions provided detailed insights into the evolutionary history of the majority of Avian ancestor chromosomes, revealed differences in rates and evolutionary times of occurrence of structural changes in micro- and macrochromosomes and identified genomic features that are associated with remarkable evolutionary stability of several avian chromosomes.
Reconstructed ancestral chromosome fragments
Statistics of the reconstructed ancestors (100 Kbp resolution)
Total length RACFs (Kbp)
Longest RACF (Kbp)
Shortest RACF (Kbp)
Longest SF (Kbp)
Telluraves & Aequornithia & Gruae
Telluraves & Aequornithia
Passeroidea & Paroidea
Estrildidae & Thraupidae & Fringillidae
Ancestors’ chromosome structure reconstructions generated at 100 Kbp resolution for syntenic fragments (SFs) (see reference genome selection and SF resolution selection criteria in Additional file 1: Supplementary Information) resulted in RACFs for the Avian ancestor and going through Neognathae, Neoavian, landbirds and Passeriformes ancestors to Estrildidae & Thraupidae & Fringillidae ancestor (Fig. 1, Table 1). The number of RACFs ranged from 46 (Eufalconimorphae ancestor) to 89 (Passerea ancestor) covering 78–99% of the reference (zebra finch) genome (Table 1). A lower fragmentation (no. RACFs ≤ 56) was observed for those ancestors for which both sides of the speciation node contained chromosome- or scaffold-level assemblies with N50 > 9 Mbp (ETF, PAO, PPA, PAR, EUF, TEL, NEA and NEO ancestors; Table 1) except for Psittacopasserae (PSI) reconstruction which contained 65 RACFs.
Avian ancestors’ chromosomes
Evolutionary history of the Avian ancestor chromosomes
The Neognathae ancestor karyotype was similar to the Avian ancestor one, with five chromosomes (AVI2, 7, 13, 14 and Z) affected by rearrangements (Fig. 3). The largest number of rearrangements was found in chromosome Z. The Neoavian ancestor had two more chromosomes rearranged (AVI4 and 12) and additional inversions in chromosome Z (Fig. 3). The Telluraves & Aequornithia & Gruae to Psittacopasserae lineage was characterised by a low number of changes detected only in AVI1, 5, 26 and 28. In contrast, the Passeriformes ancestor had multiple additional chromosomes affected (AVI4A, 9, 11, 18 and 20; Fig. 3). This elevated rate was maintained in the Passeri ancestor affecting chromosomes AVI3, 8, 15, 21, 23, 24 and 27 (Fig. 3). As a result, only three Avian ancestor chromosomes (AVI17, 19 and 22) were found intact in all reconstructed ancestors implying their maintenance for ~ 92 MY of avian evolution (Fig. 3).
In fact, our data suggest that small chromosomes (AVI15–28) tend to be affected by rearrangements at faster rates than larger chromosomes (AVI1–14 and Z) but in more recent ancestors (p value < 0.001; Fig. 4). This difference in the evolutionary history of smaller and larger chromosomes is detected independently of the tree topology as it was observed for every phylogenetic tree used to generate DESCHRAMBLER reconstructions (Additional file 1: Figure S2 and S3).
Rates of chromosome rearrangements in avian genome evolution
To estimate and compare rates of chromosome rearrangement during avian evolution in the lineage leading to zebra finch, we calculated the number of evolutionary breakpoint regions (EBRs) for each branch of the avian phylogenetic tree (Fig. 1). To minimise the influence of potential errors introduced while merging RACFs into ancestor chromosomes, we only counted EBRs located within the original RACFs. We detected a total of 201 EBRs occurring during the ~ 100 MY of avian evolution from the Avian ancestor to zebra finch. The average rearrangement rate was estimated as 2.01 EBRs/MY. We observed that different avian lineages evolved at different rates. The Neoavian to Passerea, Telluraves to Eufalconimorphae, Passeriformes to Passeri, Estrildidae & Thraupidae & Fringillidae to zebra finch branches had rearrangement rates significantly higher than the average (> 3.5 EBRs/MY; FDR-corrected p value < 0.006; Fig. 1; Additional file 1: Table S9). The opposite trend was observed for Neognathae to Neoavian, Passerea to Telluraves & Aequornithia, Passeri to Passeroidea & Paroidea branches, which had genome rearrangement rates significantly lower than the average (< 1 EBR/MY; FDR-corrected p value < 0.03; Fig. 1; Additional file 1: Table S9). By definition, the rates of rearrangements were dependent on the branch lengths of the phylogenetic trees used in the reconstruction and therefore need to be considered with caution if the TENT tree branch lengths are doubted (Additional file 1: Table S10 and S11). However, a significantly increased rate of rearrangements after the Neoavian ancestor was supported by all the phylogenies. A low rate before the Telluraves ancestor and a high rate after the Passeriformes ancestor were supported by both the TENT  and Prum  topologies.
The genome rearrangements in man and mouse (GRIMM) webserver  allowed us to detect types and number of chromosomal rearrangements that have occurred between the Avian, the least fragmented Eufalconimorphae ancestor (no. RACFs = 46), the Passeriformes ancestor, the zebra finch and the chicken. As mentioned previously, only one interchromosomal rearrangement was observed, which corresponds to the fission of the AVI1 to form Eufalconimorphae ancestor chromosomes (EUF) 1 and 1A. The remaining rearrangements were chromosomal inversions. Consistent with the counts of EBRs, we observed an increased rate of chromosomal inversions (number of inversions per MY) for the ancestors phylogenetically closer to the zebra finch. From the Avian to the Eufalconimorphae ancestor, the rate of inversions was 0.77 inversions/MY; from Eufalconimorphae to Passeriformes, it increased to 1.64 inversions/MY; and from Passeriformes to zebra finch, we observed the highest rate of 2.58 inversions/MY. The number of inversions detected in our reconstruction in the five largest Avian ancestor chromosomes (AVI1–5; N = 59) and their zebra finch homeologues was consistent with the number reported by Romanov and colleagues (N = 54) . In contrast, the number of inversions detected between AVI1-5 and chicken chromosomes was twofold higher using our reconstructions (N = 53) than that reported by Romanov and colleagues (N = 22) .
EBR distribution in Avian ancestor chromosomes
It has been proposed that rearrangements in avian microchromosomes are rare and these chromosomes represent highly conserved blocks of synteny . To test if this hypothesis holds for the lineage leading to zebra finch, we estimated distributions of EBRs in the Avian ancestor chromosomes detected from all reconstructed ancestors and the zebra finch genome.
EBR distribution and fraction within genes, CNEs and TEs for each Avian ancestor chromosome
Chromosome fraction within
EBRs per Mbp
Average EBR distance (Mbp)
We also tested differences in EBR density by averaging the distance between EBRs and between the last/first EBR end/start of the chromosomes. We observed that AVI26, 27 and 28 had average distances between EBRs significantly lower than the genome-wide average (FDR-corrected p value = 4.16E−05; Table 2), in agreement with their higher than average EBR density per Mbp. We also observed that the chromosomes with a lower EBR density per Mbp had a higher than average distance between EBRs (FDR corrected p value = 4.59E−06; Table 2).
EBR distribution and association with DNA sequence features
Several DNA sequence features were previously found associated with positions of EBRs implying that EBRs are “hotspots” of chromosome evolution, gene birth and death and changes in gene regulation [7, 9, 21, 22]. A distinct feature of avian genomes is a strong negative association of EBRs with DNA conserved non-coding elements (CNEs) with EBRs being found in CNE-sparse regions . Most previous studies of EBRs, however, have focused on pairwise comparisons of extant genome assemblies. These studies potentially suffer from misidentification of EBRs some of which could be assembly artefacts. To investigate what DNA features would be associated with the distribution of EBRs in reconstructed Avian ancestor chromosomes, we tested the association between avian CNEs, zebra finch transposable elements (TEs) and chicken gene content (due to a more complete annotation than exists for zebra finch) with EBR densities on Avian ancestor chromosomes.
Are AVI17, 19 and 22 ancestral homologous synteny blocks?
Fraction of bases within genes, TEs and CNEs in Avian ancestor intact chromosomes, msHSBs > 1.5 Mbp and the rest of the genome
Average fraction of bases within
msHSB > 1.5 Mbp
Rest of the genome
Gene ontology enrichment analysis of AVI17, 19 and 22
Using a combination of chromosome- and scaffold-level genome assemblies, we reconstructed, for the first time, the most likely chromosome structure of avian ancestors at 14 key phylogenetic nodes in the Passeriformes lineage leading to the zebra finch. None of these ancestors was reconstructed previously to the level of detail and coverage of the reference genome presented in this work. The reconstructed karyotypes allowed tracing differences in patterns of structural evolution between larger and smaller avian chromosomes for ~ 100 MY of evolution, as well as the identification of structurally stable microchromosomes and features of avian genomes that are associated with the conservation of chromosomes in evolution.
Contrary to previous studies, in which the proposed ancestral Avian karyotype included only macrochromosomes [2, 12], all the 14 reconstructions presented herein include microchromosomes orthologous to zebra finch chromosomes 11 to 28 (except 16 and 25). Our results demonstrate that smaller chromosomes (AVI15–27) become affected by rearrangements later in evolution than larger chromosomes, which started changing structurally as early as in the Neognathae ancestor (89 MYA). The reasons behind this phenomenon remain unclear; however, the data suggest that selection against the fixation of rearrangements could be responsible for the maintenance of the structural integrity of these chromosomes up to the Psittacopasserae (55 MYA) or even Passeriformes (32 MYA) ancestors. Due to their higher recombination rates (twice as high as in macrochromosomes and five times higher than in mammalian chromosomes), the requirement of at least one chiasma in each chromosome pair for correct segregation during meiosis [25, 26], and a positive association between the recombination rates and rearrangements in birds , microchromosomes are expected to have (and most of them demonstrate in extant genomes) higher rearrangement rates than larger chromosomes. Previous studies confirmed that microchromosomes were present as individual chromosomes as early as in the tetrapods’ ancestor [18, 25] and were not parts of larger chromosomes, therefore, could not escape the need for a chiasma in meiosis. The likely reason, therefore, why microchromosomes stayed intact for ~ 100 MY years of avian evolution could be the purifying selection against fixation of rearrangements affecting these chromosomes. Another possibility is that transposable elements (TEs) or other repetitive sequences could have been expanded in the Psittacopasserae lineage starting about 47 MYA providing more oportunities for non-allelic recombination in microchromosomes. This option implies that in the earlier ancestror genomes (> 47 MYA) repetitive elements were distributed differently between smaller and larger chromosomes (the later started rearranging in the Neognanthae ancestor). Our results demonstrate that zebra finch TE density translated to the Avian ancestor chromosomes indeed correlates with EBRs distribution along these chromosomes in the lineage leading to zebra finch. Therefore, it is possible that densities of TEs and CNEs acted in opposite directions: TEs increased oppportunities for aberrations in germ cells, while the presence of CNEs could have caused selection to disregard these changes. Further support to this hypothesis is provided by the fact that evolutionary stable AVI17, 19 and 22 (maintained intact in all reconstructed ancestors and majority of extant genomes) combined possess a significantly higher fraction of both TEs and genes than the rest of the genome suggesting that under no selection scenario, they should have accumulated internal rearrangements at the same rate as other chromosomes. Our data imply that the selection pressure has decreased about 47 MYA with 50% of microchromosomes rearranged in the Psittacopasserae ancestor and additional ones in the Passeri ancestor coinciding with the burst of diversification of Passeriformes  and reported TE expansion in passerine ancestors . It is tempting to speculate that changes in these chromosomes could have an important role in the generation of phenotypic diversity of extant Passeriformes by forming new regulatory networks (e.g. due to changes of positions of CNEs in chromosomes). An analogous burst of chromosomal rearrangements in mammals was observed in gibbons (also accompanied with activity of TEs), but it involved interchromosomal changes in a single species  while in birds the karyotypes were maintained in the majority of clades.
Multiple studies demonstrated that evolutionary breakpoint regions are found in chromosomal intervals enriched for TEs and genes, but sparse in CNEs that often play gene-regulatory or structural roles in the cell [7, 9, 21–23]. The fact that the density of descendant ancestral genome and zebra finch EBRs was inversely correlated with the density of avian CNEs on the Avian ancestor chromosomes, suggests that a higher CNE fraction found in some chromosomes could account for the evolutionary stability of these chromosomes, likely by preventing fixation of germ cell chromosomal aberrations. This finding is supported by our earlier discoveries that CNEs located near genes with development-related roles in avian and other reptilian phenotypes are likely to contribute to the formation of avian multispecies homologous synteny blocks, because they contained novel regulatory elements (e.g. transcription factor binding sites ) and that avian EBRs are found in CNE-sparse genome intervals . Our results, however, suggest that this mechanism could be extended to complete chromosomes, rather than individual homologous synteny blocks or EBRs within chromosomes. The difference in rearrangement rates between larger and smaller chromosomes, the fact that all three evolutionary stable Avian ancestor chromosomes found in this study were microchromosomes, and their sequence feature content, suggest that these chromosomes might behave as large homologous synteny blocks and their further breakage/rearrangements would likely have significant biological effects. Also, the fact that the TE density in Avian ancestor chromosomes correlates with their EBR density supports the role of TEs in the genome evolution of birds by promoting DNA breakage and/or joining [31, 32]. The location in these evolutionary stable chromosomes of many genes essential for the correct development of embryos as indicated by the GO analysis is supporting this hypothesis as well, as changes in the organisation of these genomic regions could disturb gene regulation and have deleterious functional effects, leading to their removal by purifying selection. Rearrangements found in these three chromosomes in a small number of extant genomes are therefore surprising suggesting that either these changes might had important biological implications or are assembly artefacts to be fixed in future genome assembly improvement experiments.
High variability of genome-wide chromosomal rearrangement rates between avian phylogenetic clades reported here agrees with the data reported by Zhang and colleagues  with discrepancies found for the Passeriformes to Passeri and the Avian to Neognathae branches. This might be related to the higher number of chromosome-level assemblies and overall higher continuity of genomes present in our dataset, which facilitated the detection of rearrangements unidentified in the previous work or to a large number of EBRs that did not pass Zhang and colleagues’ conservative filtering threshold . We also cannot exclude the possibility that the merging of reconstructed fragments to chromosomes could make the Avian ancestor chromosomes, in some cases, more structurally similar to its descendant genomes than they should be.
The number of inversions detected between the five largest Avian ancestor chromosomes (1 to 5) and their zebra finch homeologues is highly consistent between the current analysis and that of our previous studies . However, the opposite is observed for the chicken where our current reconstruction allows the detection of twice as many inversions compared to what we found previously . This inconsistency might be caused by an underrepresentation of avian clades in our earlier work, leading to a bias of our reconstructions to Galloanserae genome structures.
Despite the utility of the predicted ancestor genome structures to better understand avian chromosome evolution, they are not free of limitations. Due to the use of only one reference genome to define syntenic fragments, it is possible that some ancestral sequences that are not present in the reference genome, zebra finch in our case, were omitted from the reconstructions. Moreover, the predominance of scaffold-level assemblies in the descendant species results in fragmentation of predicted ancestral chromosomes. Indeed, we observed a lower number of RACFs for those ancestors to which both sides of the speciation node contained chromosome-level assemblies or scaffold-level assemblies with N50 > 9 Mbp, which reinforces the importance of having high continuity genome assemblies to facilitate the study of chromosome evolution. Another limitation is that the phylogenetic relationship of some avian clades is not well resolved as suggested by the disparities between avian phylogenetic trees proposed in different studies [13, 19, 34]. Our reconstructions, however, imply that the TENT tree topology from Jarvis and colleagues  was the most consistent with the underlying genomic data resulting in the least fragmented ancestral genome reconstructions compared to the other two tested trees. High agreement in RACF structures for well-established ancestral nodes shared between different trees suggests that the DESCHRAMBLER reconstructions are highly stable. Exclusion or inclusion of a small number of additional species had a small effect on the reconstruction of ancestors from ‘sister’ nodes. On the other hand, differences in estimations of branch lengths and evolutionary time of nodes between trees influenced the rates of rearrangements on individual phylogenetic nodes but had no effect on the patterns of evolutionary history of individual Avian ancestor chromosomes.
In this work, the reconstruction of the Avian ancestor and 13 additional avian ancestral karyotypes offered valuable novel insights into the history and patterns of chromosome rearrangements in the avian lineage. Our finding of significant differences in evolutionary histories of micro- and macrochromosomes in the ancestral Passeriformes lineage suggests that purifying selection was likely to act supporting the structural integrity of the majority of microchromosomes for ~ 50 MY of avian evolution. Why this pattern has changed ~ 47 MY ago with only three microchromosomes maintained intact in the most recent ancestor is still unknown. It could however be hypothesised that this event relates to the high density and diversity of retrotransposons observed in extant Passerine birds which activity was previously proposed to be linked to genome instability and species diversification in birds [29, 35].
The next questions to ask are as follows: (i) if this phenomenon is limited to Passeriformes or if it is a general signature of avian evolution and (ii) if it can be related to active speciation within Passeriformes? Moreover, the ancestral karyotype reconstructions presented herein provide an excellent resource for tracing structural changes in all avian lineages and to study their influence on the biology of the extant ~ 10,000 species.
To work towards answering these questions we need more high-quality (ideally chromosome-level) assemblies for all avian species. New sequencing and mapping technologies make this feasible. However, even the best assemblies are not free from errors. Reconstruction of ancestral genomes might assist the identification of potential problematic assembly regions and the production of more accurate genome assemblies. These, in turn, would increase the accuracy of the subsequent ancestral genome reconstructions, until one would eventually have access to a comprehensive catalogue of the events that shaped extant avian genomes and reveal their implications on biology of the species.
Avian and outgroup genome assemblies
The chicken (Gallus gallus; ICGSC Gallus_gallus 4.0 ), zebra finch (Taeniopygia guttata; WUGSC 3.2.4 ) and turkey (Meleagris gallopavo; TGC Turkey_2.01 ) chromosome assemblies were downloaded from the UCSC Genome Browser . The collared flycatcher (Ficedula albicollis; FicAlb1.5 ), peregrine falcon (Falco peregrinus) and rock pigeon (Columba livia)  chromosome assemblies were downloaded from NCBI. The Pekin duck (Anas platyrhynchos; BGI_duck_1.0 ) assembly was downloaded from NCBI and upgraded to chromosome level using radiation hybrid map obtained from Dr. Thomas Faraut (INRA). The hooded crow (Corvus cornix; Hooded_crow_genome ), canary (Serinus canaria; SCA1 ), Tibetan ground tit (Pseudopodoces humilis; PseHum1.0 ), golden eagle (Aquila chrysaetos; Aquila_chrysaetos-1.0.2 ) and bald eagle (Haliaeetus leucocephalus; Haliaeetus_leucocephalus-4.0) scaffold assemblies were obtained from NCBI. All remaining scaffold-based assemblies were downloaded from the GigaScience Database [46, 47]. Chromosome assemblies of outgroup genomes: anole lizard (Anolis carolinensis; AnoCar2.0 ) and opossum (Monodelphis domestica; MonDom5 ), the scaffolds assemblies of the Chinese alligator (Alligator sinensis; ASM45574v1 ) and the painted turtle (Chrysemys picta; Chrysemys_picta_bellii-3.0.1 ) were obtained from NCBI. General assembly statistics for each genome used are presented in Additional file 1: Table S1. Divergence times and topologies were obtained from the total evidence nucleotide tree reported by Jarvis and colleagues (2014) . Clade nomenclature was based on Jarvis and colleagues , Suh and colleagues , Yuri and colleagues  and taxonomy from Flux webpage .
We selected zebra finch as the reference genome for the reconstruction of avian ancestors as one of the best avian genome assemblies currently available, and because DESCHRAMBLER  requires the reference to be a descendant species for all reconstructed ancestors. Moreover, zebra finch is a representative of Passeriformes, the avian clade with the highest number of extant species, and which species also exhibit a high phenotypic diversity. Pairwise alignments using zebra finch chromosome assembly as the reference and all other genomes as targets were generated with LastZ (version 1.02.00 ) using the following parameters: C = 0 E = 30 H = 2000 K = 3000 L = 2200 O = 400. The pairwise alignments were converted into the UCSC “chain” and “net” alignment formats with axtChain (parameters: -minScore = 1000 -verbose = 0 -linearGap = loose for anole lizard and opossum, and -minScore = 1000 -verbose = 0 -linearGap = medium for all other species) followed by chainAntiRepeat, chainSort, chainPreNet, chainNet and netSyntenic, all with default parameters .
Reconstructed ancestral chromosome fragments
First, DESCHRAMBLER  was used to reconstruct RACFs of the Neognathae ancestor with a subset of species, as indicated in Additional file 1: Table S1. This experiment was performed at 100, 300 and 500 Kbp SFs resolution. After the selection of the best SF resolution for avian ancestral chromosomes reconstruction (100 Kbp; see resolution selection criteria in Additional file 1: Supplementary Information), DESCHRAMBLER was run with the full set of species to generate RACFs for all ancestors leading to zebra finch lineage, starting with the Avian ancestor.
After the resolution of the analysis was established, we performed test DESCHRAMBLER reconstructions using three phylogenetic trees. As the topology at the base of Neoaves is not completely resolved, we chose (1) a TENT tree from Jarvis and colleagues as a tree with the largest number of loci used , (2) a tree from Reddy and colleagues with a large taxon sampling and including both coding and non-cording loci in the analysis  and (3) a tree from Prum and colleagues due to a large taxon sampling . These trees disagree in placement of some species on the phylogenetic nodes. DESCHRAMBLER produced the least fragmented reconstructions using the Jarvis and colleagues (2014) tree. RACFs reconstructed using different trees were compared and found to be highly consistent (see Additional file 1: Supplementary Information) with all trees resulting in a similar number of EBRs. With the least fragmented reconstructions, the TENT tree results were chosen to be presented in the main text of the paper with other reconstructions being used to support the results or to indicate differences.
Detection of EBRs and chromosome rearrangements
We detected EBRs relative to the Avian ancestor in all other ancestors’ RACFs and the zebra finch using a previously published methodology . Breakpoint rates (EBRs/MY) for each branch leading to zebra finch were calculated dividing the number of detected EBRs by the length of the branch (in MY as in Jarvis and colleagues (2014) ). Differences in breakpoint rates compared to the average of all branches were tested as previously described . Differences in cumulative fractions of smaller and larger Avian ancestor chromosomes rearranged in reconstructed ancestors were tested by analysis of covariance in R (version 3.4.2 ).
We used the genome rearrangements in man and mouse (GRIMM) webserver  to predict the minimum number and the type of chromosomal rearrangements distinguishing the Avian ancestor chromosomes structure from those of the Eufalconimorphae and Passeriformes ancestors, the zebra finch and chicken genomes.
Avian ancestors’ chromosomes
The number of RACFs reconstructed by DESCHRAMBLER was higher than the number of Avian ancestor chromosomes previously proposed based on FISH experiments . This fragmentation is mostly due to the predominance of scaffold-level assemblies for the descendant species, resulting in a reduction of adjacency support. To reduce the fragmentation of the reconstructed avian ancestors’ genomes, we ordered RACFs by connecting RACFs which adjacency was supported by outgroup genomes or other, phylogenetically close and less fragmented (at the same position) ancestors. Specifically, for the Avian ancestor, we first merged those RACFs which adjacencies were supported (spanned) by an outgroup chromosome or scaffold. For the remaining Avian RACFs adjacencies with no support from outgroup genomes and the other reconstructed ancestors, we merged RACFs which adjacency was supported by other ancestor RACF, assuming that no rearrangement occurred between the target and the descendant ancestor in between RACFs. For each RACF adjacency, we used the support from the spanning RACF belonging to ancestors successively more distant. That is, we used first the support from the closest ancestor (e.g. Neognathae for the Avian ancestor) and successively more recent ancestors on the avian phylogenetic tree.
Fraction of rearranged avian ancestor chromosomes
For each Avian ancestor chromosomes, we calculated its fraction involved in intrachromosomal rearrangements on the other ancestors and extant descendant and outgroup species. We first established the ancestral state chromosome orientation by detecting which orientation would imply the least number of rearrangements. Then, we calculated its fraction involved in rearrangements by dividing the non-ancestral orientation by the cumulative length of the blocks mapped into that chromosome.
The fraction of the Avian ancestor chromosomes affected by interchromosomal rearrangements was calculated by dividing the cumulative length of the blocks of each target chromosome by the total length of the target blocks mapped into the Avian ancestor chromosome. The represented fraction corresponds to the lowest obtained value.
EBR rates and DNA sequence feature associations on Avian ancestor chromosomes and HSBs
We measured EBR density and distribution for the Avian ancestor chromosomes using the number of EBRs identified between the Avian ancestor and zebra finch. These measurements were obtained as the number of EBRs per Mbp and the average distance between EBRs. Differences between chromosomes for each of the analysed features were tested as previously reported .
Avian conserved non-coding elements (CNEs) were obtained from Farré and colleagues (2016) . Chicken gene (version of 27/04/2014) and zebra finch repetitive sequence (version of 08/05/2014) annotations were downloaded from the UCSC genome browser . We calculated the density of each of these features (CNEs, genes and TEs) for each Avian ancestor chromosome. The association between each sequence feature and chromosome-specific EBR density and distribution was tested using the Pearson’s correlation coefficient.
Avian ancestor chromosomes were divided into 1 Kbp non-overlapping windows. All intervals were assigned to either msHSBs, detected across the reconstructed ancestral chromosome and extant species genome assemblies, intact chromosomes (AVI17, 19 and 22) or intervals found in the rest of the genome. The CNE, TE and gene density were calculated for each window type using bedtools (version 2.20-1 ). Differences between each of these two sets for each analysed feature were tested as previously reported .
Gene ontology enrichment analysis
The basic version of gene ontology (GO) annotations (version 8 April 2017) was downloaded from the GO Consortium website . Sequence coordinates and Ensembl identifiers for chicken genes were obtained from Ensembl Biomart (version 74 ). All chicken genes located in regions included in the Avian ancestor chromosomes were used as the background list. To evaluate gene functional enrichment in the Avian ancestor chromosomes that were maintained intact during avian evolution, we assigned genes from the background list to these chromosomes. We used the GO::TermFinder Perl module  to detect GO terms overrepresented in our gene sets. We considered as significantly enriched the terms with p value < 0.05 and false discovery rate (FDR) < 5%.
This work was supported by the Biotechnology and Biological Sciences Research Council (grants BB/K008226/1 and BB/J010170/1 to DML and BB/K008161/1 to DKG)
Availability of data and materials
The Avian ancestor reconstruction datasets and alignment of ingroup and outgroup genomes against the Avian ancestor are publicly available from the Evolution Highway comparative chromosome browser: http://eh-demo.ncsa.uiuc.edu/birds [7, 17]. All other data generated or analysed during this study are included in this published article and its additional files.
This study utilised genome assemblies for 31 species which are publicly available from NCBI assembly database , with the following GenBank accession numbers:
Taeniopygia guttata , GCA_000151805; Geospiza fortis , GCA_000277835.1; Serinus canaria , GCA_000534875.1; Pseudopodoces humilis , GCA_000331425.1; Corvus brachyrhynchos , GCA_000691975.1; Corvus cornix , GCA_000738735.2; Ficedula albicollis , CA_000247815.2; Manacus vitellinus , GCA_000692015.2; Melopsittacus undulatus , GCA_000238935.1; Falco peregrinus , GCA_001887755.1; Aquila chrysaetos , GCA_000766835.1; Haliaeetus leucocephalus , GCA_000737465.1; Picoides pubescens , GCA_000699005.1; Pygoscelis adeliae , GCA_000699105.1; Aptenodytes forsteri , GCA_000699145.1; Nipponia nippon , GCA_000708225.1; Egretta garzetta , GCA_000687185.1; Opisthocomus hoatzin , GCA_000692075.1; Charadrius vociferus , GCA_000708025.2; Cuculus canorus , GCA_000709325.1; Chaetura pelagica , GCA_000747805.1; Calypte anna , GCA_000699085.1; Columba livia , GCA_001887795.1; Gallus gallus , GCA_000002315.2; Meleagris gallopavo , GCA_000002315.2; Anas platyrhynchos , GCA_000355885.1; Struthio camelus , GCA_000698965.1; Alligator sinensis , GCA_000455745.1; Chrysemys picta , GCA_000241765.2; Anolis carolinensis , GCA_000090745.2; and Monodelphis domestica , GCF_000002295.2.
JD prepared DESCHARMBER input data, performed the data analysis and drafted the manuscript. JK generated the DESCHARMBLER reconstructions and revised the manuscript. MF assisted with data analysis and revised the manuscript. DML and DKG conceived the study. DKG revised the manuscript. DML supervised the project, assisted with data analysis and revised the manuscript. All authors approved the final version of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Richard F, Lombard M, Dutrillaux B. Reconstruction of the ancestral karyotype of eutherian mammals. Chromosom Res. 2003;11:605–18.View ArticleGoogle Scholar
- Griffin DK, Robertson LB, Tempest HG, Skinner BM. The evolution of the avian genome as revealed by comparative molecular cytogenetics. Cytogenet Genome Res. 2007;117:64–77.View ArticleGoogle Scholar
- Ruiz-Herrera A, Farré M, Robinson TJ. Molecular cytogenetic and genomic insights into chromosomal evolution. Heredity. 2012;108:28–36.View ArticleGoogle Scholar
- Wurster DH, Benirschke K. Indian Muntjac, Muntiacus muntiak: a deer with a low diploid chromosome number. Science. 1970;168:1364–6.View ArticleGoogle Scholar
- Gallardo MH, González CA, Cebrián I. Molecular cytogenetics and allotetraploidy in the red vizcacha rat, Tympanoctomys barrerae (Rodentia, Octodontidae). Genomics. 2006;88:214–21.View ArticleGoogle Scholar
- Bourque G, Zdobnov EM, Bork P, Pevzner PA, Tesler G. Comparative architectures of mammalian and chicken genomes reveal highly variable rates of genomic rearrangements across different lineages. Genome Res. 2005;15:98–110.View ArticleGoogle Scholar
- Murphy WJ, Larkin DM, Everts-van der Wind A, Bourque G, Tesler G, Auvil L, et al. Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps. Science. 2005;309:613–7.View ArticleGoogle Scholar
- Skinner BM, Griffin DK. Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints. Heredity. 2012;108:37–41.View ArticleGoogle Scholar
- Ma J, Zhang L, Suh BB, Raney BJ, Burhans RC, Kent WJ, et al. Reconstructing contiguous regions of an ancestral genome. Genome Res. 2006;16:1557–65.View ArticleGoogle Scholar
- Jones BR, Rajaraman A, Tannier E, Chauve C. ANGES: reconstructing ANcestral GEnomeS maps. Bioinformatics. 2012;28:2388–90.View ArticleGoogle Scholar
- Kim J, Farré M, Auvil L, Capitanu B, Larkin DM, Ma J, et al. Reconstruction and evolutionary history of eutherian chromosomes. Proc Natl Acad Sci. 2017;114:E5379–88.View ArticleGoogle Scholar
- Romanov MN, Farré M, Lithgow PE, Fowler KE, Skinner BM, O’Connor R, et al. Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor. BMC Genomics. 2014;15:1–18.View ArticleGoogle Scholar
- Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, et al. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346:1320–31.View ArticleGoogle Scholar
- Ishishita S, Tsuruta Y, Uno Y, Nakamura A, Nishida C, Griffin DK, et al. Chromosome size-correlated and chromosome size-uncorrelated homogenization of centromeric repetitive sequences in New World quails. Chromosom Res. 2014;22:15–34.View ArticleGoogle Scholar
- de Oliveira EH, Habermann FA, Lacerda O, Sbalqueiro IJ, Wienberg J, Muller S. Chromosome reshuffling in birds of prey: the karyotype of the world's largest eagle (Harpy eagle, Harpia harpyja) compared to that of the chicken (Gallus gallus). Chromosoma. 2005;114:338–43.View ArticleGoogle Scholar
- Nanda I, Karl E, Griffin DK, Schartl M, Schmid M. Chromosome repatterning in three representative parrots (Psittaciformes) inferred from comparative chromosome painting. Cytogenet Genome Res. 2007;117:43–53.View ArticleGoogle Scholar
- Damas J, Kim J, Farré M, Griffin DK, Larkin DM. Avian ancestor comparative data. Evolution Highway Comparative Chromosome Browser. http://eh-demo.ncsa.uiuc.edu/birds; 2018.
- Uno Y, Nishida C, Tarui H, Ishishita S, Takagi C, Nishimura O, et al. Inference of the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes from comparative gene mapping. PLoS One. 2012;7:e53027.View ArticleGoogle Scholar
- Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, et al. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature. 2015;526:569–73.View ArticleGoogle Scholar
- Tesler G. GRIMM: genome rearrangements web server. Bioinformatics. 2002;18:492–3.View ArticleGoogle Scholar
- Farré M, Narayan J, Slavov GT, Damas J, Auvil L, Li C, et al. Novel insights into chromosome evolution in birds, archosaurs, and reptiles. Genome Biol Evol. 2016;8:2442–51.View ArticleGoogle Scholar
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491:393–8.Google Scholar
- Damas J, O'Connor R, Farré M, Lenis VPE, Martell HJ, Mandawala A, et al. Upgrading short-read animal genome assemblies to chromosome level using comparative genomics and a universal probe set. Genome Res. 2017;27:875–84.View ArticleGoogle Scholar
- Larkin DM, Pape G, Donthu R, Auvil L, Welge M, Lewin HA. Breakpoint regions and homologous synteny blocks in chromosomes have different evolutionary histories. Genome Res. 2009;19:770–7.View ArticleGoogle Scholar
- Burt DW. Origin and evolution of avian microchromosomes. Cytogenet Genome Res. 2002;96:97–112.View ArticleGoogle Scholar
- Rodionov AV, Chel'sheva LA, Solovei IV, Miakoshina Iu A. Chiasma distribution in the lampbrush chromosomes of the chicken Gallus gallus domesticus: hot spots of recombination and their possible role in proper dysjunction of homologous chromosomes at the first meiotic division. Genetika. 1992;28:151–60.PubMedGoogle Scholar
- Völker M, Backström N, Skinner BM, Langley EJ, Bunzey SK, Ellegren H, et al. Copy number variation, chromosome rearrangement, and their association with recombination during avian evolution. Genome Res. 2010;20:503–11.View ArticleGoogle Scholar
- Barker FK, Cibois A, Schikler P, Feinstein J, Cracraft J. Phylogeny and diversification of the largest avian radiation. Proc Natl Acad Sci U S A. 2004;101:11040–5.View ArticleGoogle Scholar
- Suh A, Bachg S, Donnellan S, Joseph L, Brosius J, Kriegs JO, et al. De-novo emergence of SINE retroposons during the early evolution of passerine birds. Mob DNA. 2017;8:21.View ArticleGoogle Scholar
- Carbone L, Harris RA, Gnerre S, Veeramah KR, Lorente-Galdos B, Huddleston J, et al. Gibbon genome and the fast karyotype evolution of small apes. Nature. 2014;513:195–201.View ArticleGoogle Scholar
- Gregory TR. The evolution of the genome. Burlington: Elsevier Academic; 2005.Google Scholar
- Wessler SR. Transposable elements and the evolution of eukaryotic genomes. Proc Natl Acad Sci. 2006;103:17600–1.View ArticleGoogle Scholar
- Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, et al. Comparative genomics reveals insights into avian genome evolution and adaptation. Science. 2014;346:1311–20.View ArticleGoogle Scholar
- Reddy S, Kimball RT, Pandey A, Hosner PA, Braun MJ, Hackett SJ, et al. Why do phylogenomic data sets yield conflicting trees? Data type influences the avian tree of life more than taxon sampling. Syst Biol. 2017;66:857–79.View ArticleGoogle Scholar
- Kapusta A, Suh A. Evolution of bird genomes—a transposon’s-eye view. Ann N Y Acad Sci. 2017;1389:164–85.View ArticleGoogle Scholar
- Consortium ICGS. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004;432:695–716.View ArticleGoogle Scholar
- Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, et al. The genome of a songbird. Nature. 2010;464:757–62.View ArticleGoogle Scholar
- Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Blomberg Le A, et al. Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis. PLoS Biol. 2010;8:e1000475.View ArticleGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12:996–1006.View ArticleGoogle Scholar
- Ellegren H, Smeds L, Burri R, Olason PI, Backstrom N, Kawakami T, et al. The genomic landscape of species divergence in Ficedula flycatchers. Nature. 2012;491:756–60.View ArticleGoogle Scholar
- Huang YH, Li YR, Burt DW, Chen HL, Zhang Y, Qian WB, et al. The duck genome and transcriptome provide insight into an avian influenza virus reservoir species. Nature Genet. 2013;45:776.View ArticleGoogle Scholar
- Poelstra JW, Vijay N, Bossu CM, Lantz H, Ryll B, Muller I, et al. The genomic landscape underlying phenotypic integrity in the face of gene flow in crows. Science. 2014;344:1410–4.View ArticleGoogle Scholar
- Frankl-Vilches C, Kuhl H, Werber M, Klages S, Kerick M, Bakker A, et al. Using the canary genome to decipher the evolution of hormone-sensitive gene regulation in seasonal singing birds. Genome Biol. 2015;16:19.View ArticleGoogle Scholar
- Cai Q, Qian X, Lang Y, Luo Y, Xu J, Pan S, et al. Genome sequence of ground tit Pseudopodoces humilis and its adaptation to high altitude. Genome Biol. 2013;14:R29.View ArticleGoogle Scholar
- Doyle JM, Katzner TE, Bloom PH, Ji Y, Wijayawardena BK, DeWoody JA. The genome sequence of a widespread apex predator, the golden eagle (Aquila chrysaetos). PLoS One. 2014;9:e95599.View ArticleGoogle Scholar
- Sneddon TP, Li P, Edmunds SC. GigaDB: announcing the GigaScience database. GigaScience. 2012;1:11.View ArticleGoogle Scholar
- Zhang G, Li B, Li C, Gilbert MTP, Jarvis ED, Wang J. Comparative genomic data of the Avian Phylogenomics Project. GigaScience. 2014;3:26.View ArticleGoogle Scholar
- Alfoldi J, Di Palma F, Grabherr M, Williams C, Kong L, Mauceli E, et al. The genome of the green anole lizard and a comparative analysis with birds and mammals. Nature. 2011;477:587–91.View ArticleGoogle Scholar
- Mikkelsen TS, Wakefield MJ, Aken B, Amemiya CT, Chang JL, Duke S, et al. Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences. Nature. 2007;447:167–77.View ArticleGoogle Scholar
- Wan QH, Pan SK, Hu L, Zhu Y, Xu PW, Xia JQ, et al. Genome analysis and signature discovery for diving and sensory properties of the endangered Chinese alligator. Cell Res. 2013;23:1091–105.View ArticleGoogle Scholar
- Shaffer HB, Minx P, Warren DE, Shedlock AM, Thomson RC, Valenzuela N, et al. The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage. Genome Biol. 2013;14:R28.View ArticleGoogle Scholar
- Suh A, Paus M, Kiefmann M, Churakov G, Franke FA, Brosius J, et al. Mesozoic retroposons reveal parrots as the closest living relatives of passerine birds. Nat Commun. 2011;2:443.View ArticleGoogle Scholar
- Yuri T, Kimball RT, Harshman J, Bowie RC, Braun MJ, Chojnowski JL, et al. Parsimony and model-based analyses of indels in avian nuclear genes reveal congruent and incongruent phylogenetic signals. Biology (Basel). 2013;2:419–44.Google Scholar
- Taxonomy in Flux http://www.jboyd.net/Taxo/List.html. Accessed Nov 2016.
- Harris RS. Improved pairwise alignment of genomic DNA. The Pennsylvania State University; 2007.Google Scholar
- Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D. Evolution’s cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci. 2003;100:11484–9.View ArticleGoogle Scholar
- CoreTeam R. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2012. http://www.R-project.org. 2015Google Scholar
- Rosenbloom KR, Armstrong J, Barber GP, Casper J, Clawson H, Diekhans M, et al. The UCSC Genome Browser database: 2015 update. Nucleic Acids Res. 2015;43:D670–81.View ArticleGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.View ArticleGoogle Scholar
- The Gene Ontology Consortium. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2015;43:D1049–56.View ArticleGoogle Scholar
- Kinsella RJ, Kahari A, Haider S, Zamora J, Proctor G, Spudich G, et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database. 2011;2011:bar030.View ArticleGoogle Scholar
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder—open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20:3710–5.View ArticleGoogle Scholar
- Kitts PA, Church DM, Thibaud-Nissen F, Choi J, Hem V, Sapojnikov V, et al. Assembly: a resource for assembled genomes at NCBI. Nucleic Acids Res. 2016;44:D73–80.View ArticleGoogle Scholar