Specific chicken repeat 1 (CR1) retrotransposon insertion suggests phylogenetic affinity of rockfowls (genus Picathartes) to crows and ravens (Corvidae)

While the monophyly of the order Passeriformes as well as its suborders suboscines (Tyranni) and oscines (Passeri) is well established, both on morphological and molecular grounds, lower phylogenetic relationships have been a continuous matter of debate, especially within oscines. This is particularly true for the rockfowls (genus Picathartes), which phylogenetic classification has been an ongoing puzzle. Sequence-based molecular studies failed in deriving unambiguously resolved and supported hypotheses. We present here a novel approach: use of retrotransposon insertions as phylogenetic markers in passerine birds. Chicken repeat 1 (CR1) is the most important non-LTR retrotransposon in birds. We present two truncated CR1 loci in passerine birds, not only found in representatives of Corvinae (jays, crows and allies), but also in the West-African Picathartes species which provide new evidence for a closer relationship of these species to Corvidae than has previously been thought. Additionally, we show that not only the absence/presence pattern of a CR1 insertion, but also the CR1 sequences themselves contain phylogenetic information.     

Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds

Using standard phylogenetic methods, it can be hard to resolve the order in which speciation events took place when new lineages evolved in the distant past and within a short time frame. As an example, phylogenies of galliform birds (including well-known species such as chicken, turkey, and quail) usually show low bootstrap support values at short internal branches, reflecting the rapid diversification of these birds in the Eocene. However, given the key role of chicken and related poultry species in agricultural, evolutionary, general biological and disease studies, it is important to know their internal relationships. Recently, insertion patterns of transposable elements such as long and short interspersed nuclear element markers have proved powerful in revealing branching orders of difficult phylogenies. Here we decipher the order of speciation events in a group of 27 galliform species based on insertion events of chicken repeat 1 (CR1) transposable elements. Forty-four CR1 marker loci were identified from the draft sequence of the chicken genome, and from turkey BAC clone sequence, and the presence or absence of markers across species was investigated via electrophoretic size separation of amplification products and subsequent confirmation by DNA sequencing. Thirty markers proved possible to type with electrophoresis of which 20 were phylogenetically informative. The distribution of these repeat elements supported a single homoplasy-free cladogram, which confirmed that megapodes, cracids, New World quail, and guinea fowl form outgroups to Phasianidae and that quails, pheasants, and partridges are each polyphyletic groups. Importantly, we show that chicken is an outgroup to turkey and quail, an observation which does not have significant support from previous DNA sequence- and DNA-DNA hybridization-based trees and has important implications for evolutionary studies based on sequence or karyotype data from galliforms. We discuss the potential and limitations of using a genome-based retrotransposon approach in resolving problematic phylogenies among birds.     

Biased distributions and decay of long interspersed nuclear elements in the chicken genome

The genomes of birds are much smaller than mammalian genomes, and transposable elements (TEs) make up only 10% of the chicken genome, compared with the 45% of the human genome. To study the mechanisms that constrain the copy numbers of TEs, and as a consequence the genome size of birds, we analyzed the distributions of LINEs (CR1's) and SINEs (MIRs) on the chicken autosomes and Z chromosome. We show that (1) CR1 repeats are longest on the Z chromosome and their length is negatively correlated with the local GC content; (2) the decay of CR1 elements is highly biased, and the 5'-ends of the insertions are lost much faster than their 3'-ends; (3) the GC distribution of CR1 repeats shows a bimodal pattern with repeats enriched in both AT-rich and GC-rich regions of the genome, but the CR1 families show large differences in their GC distribution; and (4) the few MIRs in the chicken are most abundant in regions with intermediate GC content. Our results indicate that the primary mechanism that removes repeats from the chicken genome is ectopic exchange and that the low abundance of repeats in avian genomes is likely to be the consequence of their high recombination rates.     

Marker utility of miniature inverted-repeat transposable elements for wheat biodiversity and evolution

Transposable elements (TEs) account for up to 80% of the wheat genome and are considered one of the main drivers of wheat genome evolution. However, the contribution of TEs to the divergence and evolution of wheat genomes is not fully understood. In this study, we have developed 55 miniature inverted-repeat transposable element (MITE) markers that are based on the presence/absence of an element, with over 60% of these 55 MITE insertions associated with wheat genes. We then applied these markers to assess genetic diversity among Triticum and Aegilops species, including diploid (AA, BB and DD genomes), tetraploid (BBAA genome) and hexaploid (BBAADD genome) species. While 18.2% of the MITE markers showed similar insertions in all species indicating that those are fossil insertions, 81.8% of the markers showed polymorphic insertions among species, subspecies, and accessions. Furthermore, a phylogenetic analysis based on MITE markers revealed that species were clustered based on genus, genome composition, and ploidy level, while 47.13% genetic divergence was observed between the two main clusters, diploids versus polyploids. In addition, we provide evidence for MITE dynamics in wild emmer populations. The use of MITEs as evolutionary markers might shed more light on the origin of the B-genome of polyploid wheat.     

Retroposon insertion patterns of neoavian birds: strong evidence for an extensive incomplete lineage sorting era

More than 150 Ma, the avian lineage separated from that of other dinosaurs and later diversified into the more than 10,000 species extant today. The early neoavian bird radiations most likely occurred in the late Cretaceous (more than 65 Ma) but left behind few if any molecular signals of their archaic evolutionary past. Retroposed elements, once established in an ancestral population, are highly valuable, virtually homoplasy-free markers of species evolution; after applying stringent orthology criteria, their phylogenetically informative presence/absence patterns are free of random noise and independent of evolutionary rate or nucleotide composition effects. We screened for early neoavian orthologous retroposon insertions and identified six markers with conflicting presence/absence patterns, whereas six additional retroposons established before or after the presumed major neoavian radiation show consistent phylogenetic patterns. The exceptionally frequent conflicting retroposon presence/absence patterns of neoavian orders are strong indicators of an extensive incomplete lineage sorting era, potentially induced by an early rapid successive speciation of ancestral Neoaves.     

Identification of CR1 retroposons in Arborophila rufipectus and their application to Phasianidae phylogeny

Chicken repeat 1 (CR1), a member of non‐LTR retroposon, is an important phylogenetic marker in avian systematics. In this study, we reported several characteristics of CR1 elements in a draft genome of Arborophila rufipectus (Sichuan partridge). According to the analyses of RepeatMasker, approximately 254 966 CR1 elements were identified in A. rufipectus, covering 6.7% of the genome. Subsequently, we selected eighteen novel CR1 elements by comparing the chicken genome, turkey genome and assembled A. rufipectus scaffolds. Here, a combined data set comprising of 22 CR1 loci, mitochondrial genomes and eight unlinked introns was analysed to infer the evolutionary relationships of twelve Phasianidae species. The applicability of CR1 sequences for inferring avian phylogeny relative to mtDNA and intron sequences was investigated as well. Our results elucidated the position of A. rufipectus in Phasianidae with robust supports that it presented a sister clade to Arborophila ardens/Arborophila brunneopectus, and implied that genus Arborophila was in a basal phylogenetic position within Phasianidae and a phylogenetic affinity between Meleagris gallopavo and Pucrasia macrolopha. Therefore, this work not only resolved some of the confounding relationships among Phasianidae, but also suggested CR1 sequences could provide powerful complementary data for phylogeny reconstruction.     

Tracking ancient polyploids: a retroposon insertion reveals an extinct diploid ancestor in the polyploid origin of belladonna

Polyploidy is a prominent process in plant evolution and adaptation, but molecular phylogenetic studies of polyploids based on DNA sequences have often been confounded by their complex gene and genome histories. We report here a retroposon insertion in the nuclear gene granule-bound starch synthase I (GBSSI or "waxy") that clearly reveals the ancient hybrid history of the medically important polyploid species belladonna (Atropa belladonna) and resolves the controversy over the taxonomic group to which it belongs, the tribe Hyoscyameae (Solanaceae). Our inferences based on the pattern of presence or absence of the retroposon insertion are corroborated by phylogenetic analyses of the GBSSI gene sequences. This case may suggest that retroposons are promising molecular markers to study polyploid evolution.     

Phylogeny of the macaques (Cercopithecidae: Macaca) based on Alu elements

Genus Macaca (Cercopithecidae: Papionini) is one of the most successful primate radiations. Despite previous studies on morphology and mitochondrial DNA analysis, a number of issues regarding the details of macaque evolution remain unsolved. Alu elements are a class of non-autonomous retroposons belonging to short interspersed elements that are specific to the primate lineage. Because retroposon insertions show very little homoplasy, and because the ancestral state (absence of the SINE) is known, Alu elements are useful genetic markers and have been utilized for analyzing primate phylogenentic relationships and human population genetic relationships. Using PCR display methodology, 298 new Alu insertions have been identified from ten species of macaques. Together with 60 loci reported previously, a total of 358 loci are used to infer the phylogenetic relationships of genus Macaca. With regard to earlier unresolved issues on the macaque evolution, the topology of our tree suggests that: 1) genus Macaca contains four monophyletic species groups; 2) within the Asian macaques, the silenus group diverged first, and members of the sinica and fascicularis groups share a common ancestor; 3) Macaca arctoides are classified in the sinica group. Our results provide a robust molecular phylogeny for genus Macaca with stronger statistical support than previous studies. The present study also illustrates that SINE-based approaches are a powerful tool in primate phylogenetic studies and can be used to successfully resolve evolutionary relationships between taxa at scales from the ordinal level to closely related species within one genus.     

Activity of Ancient RTE Retroposons during the Evolution of Cows, Spiral-Horned Antelopes, and Nilgais (Bovinae)

In the genome of Artiodactyla (cow, sheep, pigs, camels, and whales), a major retroposon group originated from a presumable horizontal transfer of BovB, a retrotransposon-like element retroposon, between 52 and 70 million years ago. Since then, BovB retroposons have proliferated and today occupy a quarter of the cow's genome sequence. The BovB-related short interspersed elements (SINEs) were used for resolving the phylogeny of Bovinae (cows, spiral-horned antelopes, and nilgais) and their relatives. In silico screening of 55,000 intronic retroposon insertions in the cow genome and experimental validation of 126 introns resulted in 29 informative retroposon markers for resolving bovine evolutionary relationships. A transposition-in-transposition analysis identifies three different phases of SINE activity and show how BovB elements have expanded in the cattle genome.     

Characterization of CR1 repeat random PCR markers for mapping the chicken genome

Polymerase chain reaction (PCR) primers complementary to portions of the chicken repetitive element CR1 have been used previously to generate useful markers on the chicken genome linkage map. To understand better the genetic basis for this technique and to convert CR1–PCR loci to markers useful in physical genome mapping, five polymorphic CR1–PCR-generated DNAs were cloned and partially sequenced. Inverse PCR was then employed to clone the corresponding region of the genomes of both the Jungle Fowl (JF) and White Leghorn (WL) parental DNA templates. Our results demonstrate that some of the CR1–PCR-generated DNAs arise from priming at an endogenous CR1 element, whereas others are due to chance complementarity between the CR1–PCR primer in use and random annealing sites in the genome, unrelated to a demonstrable CR1 element. In all five instances, it was possible to identify the sequence difference between the JF and WL parental DNAs that gave rise to the initial polymorphism and design allele-specific PCR primer sets that uniquely detect that polymorphism. In four of the five instances, the polymorphism was a one or two basepair sequence difference within the primer annealing site, but in the fifth case the responsible difference was outside, but very close to, the annealing site. In all instances the allele-specific PCR for the sequence polymorphism mapped identically with the corresponding CR1–PCR amplification polymorphism. We conclude that CR1–PCR provides an efficient and reliable mechanism for genome mapping in avians that can correlate linkage and physical mapping approaches.     

A comparison of synteny and gene order on the homologue of chicken chromosome 7 between two passerine species and between passerines and chicken

Wild populations of passerine birds have frequently been used in studies of ecology and evolution. However, the majority of species lack genetic linkage maps. The completion of a model avian genome sequence (the jungle fowl, Gallus gallus) has created an opportunity to advance genetic knowledge of other birds. Here we constructed genetic linkage maps of the homologue of chicken chromosome 7 for two passerine species, the zebra finch (Taeniopygia guttata) and the house sparrow (Passer domesticus). Linkage groups were constructed with a combination of SNP and microsatellite markers that were developed both in silico based on similarity to the chicken genome sequence, and de novo in the laboratory. Synteny was well conserved throughout the chromosome; however there were some gene rearrangements between chickens and passerines. This suggests that the transfer of map information from chicken to passerine and between different passerine species should not assume conserved gene order. The length of linkage maps of chromosome 7 were on average 60% that of chicken, consistent with the idea that passerines have a reduced recombination rate relative to chicken. Some evidence of heterochiasmy, i.e. a difference in the recombination rate between the sexes, was observed.     

CR1 retroposons provide a new insight into the phylogeny of Phasianidae species (Aves: Galliformes)

Chicken repeat 1 (CR1) elements, a class of retroposons belonging to non-long-terminal repeats, have been recognized as powerful tools for phylogenetic studies. Here we examine the phylogenetic relationships of 11 Phasianidae species based on CR1 retroposons. Together with 19 loci reported previously, a total of 99 CR1 loci were identified from chicken genome and turkey BAC clone sequences. 75 insertion events were used to address the branching order of 11 species in Phasianidae. The topology of our tree suggests that: 1) Gallus gallus possessed a basal phylogenetic position within Phasianidae and was related to Bambusicola thoracica (BSP=100%); 2) After the split of G. gallus and B. thoracica, Arborophila rufipectus diverged from Phasianidae (BSP=100%). Nine unambiguous insertion events supported a phylogenetic position of A. rufipectus different to previous mitochondrial data suggesting a hybrid origin or an ancient introgression of A. rufipectus; and 3) 22 CR1 insertion events strongly supported the eight phasianids under investigation sharing a common ancestor. Our study has revisited the phylogenetic position of G. gallus and A. rufipectus and provided a new insight into the phylogeny of Phasianidae birds. It showed that a CR1-based methodology has a great potential to be informative within Phasianidae in resolving relationships of closely related species whose radiation and speciation have occurred very recently.     

Phylogenetics of modern birds in the era of genomics

In the 14 years since the first higher-level bird phylogenies based on DNA sequence data, avian phylogenetics has witnessed the advent and maturation of the genomics era, the completion of the chicken genome and a suite of technologies that promise to add considerably to the agenda of avian phylogenetics. In this review, we summarize current approaches and data characteristics of recent higher-level bird studies and suggest a number of as yet untested molecular and analytical approaches for the unfolding tree of life for birds. A variety of comparative genomics strategies, including adoption of objective quality scores for sequence data, analysis of contiguous DNA sequences provided by large-insert genomic libraries, and the systematic use of retroposon insertions and other rare genomic changes all promise an integrated phylogenetics that is solidly grounded in genome evolution. The avian genome is an excellent testing ground for such approaches because of the more balanced representation of single-copy and repetitive DNA regions than in mammals. Although comparative genomics has a number of obvious uses in avian phylogenetics, its application to large numbers of taxa poses a number of methodological and infrastructural challenges, and can be greatly facilitated by a 'community genomics' approach in which the modest sequencing throughputs of single PI laboratories are pooled to produce larger, complementary datasets. Although the polymerase chain reaction era of avian phylogenetics is far from complete, the comparative genomics era-with its ability to vastly increase the number and type of molecular characters and to provide a genomic context for these characters-will usher in a host of new perspectives and opportunities for integrating genome evolution and avian phylogenetics.     

Microsatellite loci in Japanese quail and cross-species amplification in chicken and guinea fowl

In line with the Gifu University''s initiative to map the Japanese quail genome, a total of 100 Japanese quail microsatellite markers isolated in our laboratory were evaluated in a population of 20 unrelated quails randomly sampled from a colony of wild quail origin. Ninety-eight markers were polymorphic with an average of 3.7 alleles per locus and a mean heterozygosity of 0.423. To determine the utility of these markers for comparative genome mapping in Phasianidae, cross-species amplification of all the markers was tested with chicken and guinea fowl DNA. Amplification products similar in size to the orthologous loci in quail were observed in 42 loci in chicken and 20 loci in guinea fowl. Of the cross-reactive markers, 57.1% in chicken and 55.0% in guinea fowl were polymorphic when tested in 20 birds from their respective populations. Five of 15 markers that could cross-amplify Japanese quail, chicken, and guinea fowl DNA were polymorphic in all three species. Amplification of orthologous loci was confirmed by sequencing 10 loci each from chicken and guinea fowl and comparing with them the corresponding quail sequence. The microsatellite markers reported would serve as a useful resource base for genetic mapping in quail and comparative mapping in Phasianidae.     

Retropositional events consolidate the branching order among New World monkey genera

Due to contradicting relationships obtained from various morphological and genetic studies, phylogenetic relationships among New World monkey genera are highly disputed. In the present study, we analyzed the presence/absence pattern of 128 SINE integrations in all New World monkey genera. Among them, 70 were specific for only a single genus, whereas another 18 were present in all New World monkey genera. The 40 remaining insertions were informative to elucidate phylogenetic relationships among genera. Several of them confirmed the monophyly of the three families Cebidae, Atelidae and Pitheciidae as well as of the subfamily Callithrichinae. Further markers provided evidence for a sister grouping of Cebidae and Atelidae to the exclusion of Pitheciidae as well as for relationships among genera belonging to Callithrichinae and Atelidae. Although a close affiliation of Saimiri, Aotus and Cebus to Callithrichinae was shown, the relationships among the three genera remained unresolved due to three contradicting insertions.     

A genome-wide set of 106 microsatellite markers for the blue tit (Cyanistes caeruleus)

We have characterized a set of 106 microsatellite markers in 26-127 individual blue tits (Cyanistes caeruleus), and assigned their location on the zebra finch (Taeniopygia guttata) and on the chicken (Gallus gallus) genome on the basis of sequence homology. Thirty-one markers are newly designed from zebra finch EST (expressed sequence tags) sequences, 22 markers were developed by others from EST sequences using different methods and the remaining 53 loci were previously designed or modified passerine markers. The 106 microsatellite markers are distributed over 26 and 24 chromosomes in the zebra finch and in the chicken genome respectively and the number of alleles varies between 2 and 49. Eight loci deviate significantly from Hardy-Weinberg equilibrium and show a high frequency of null alleles, and three pairs of markers located in the same chromosome appear to be in linkage disequilibrium. With the exception of these few loci, the polymorphic microsatellite markers presented here provide a useful genome-wide resource for population and evolutionary genetic studies of the blue tit, in addition to their potential utility in other passerine birds.     

Incomplete Lineage Sorting and Hybridization Statistics for Large-Scale Retroposon Insertion Data

Ancient retroposon insertions can be used as virtually homoplasy-free markers to reconstruct the phylogenetic history of species. Inherited, orthologous insertions in related species offer reliable signals of a common origin of the given species. One prerequisite for such a phylogenetically informative insertion is that the inserted element was fixed in the ancestral population before speciation; if not, polymorphically inserted elements may lead to random distributions of presence/absence states during speciation and possibly to apparently conflicting reconstructions of their ancestry. Fortunately, such misleading fixed cases are relatively rare but nevertheless, need to be considered. Here, we present novel, comprehensive statistical models applicable for (1) analyzing any pattern of rare genomic changes, (2) testing and differentiating conflicting phylogenetic reconstructions based on rare genomic changes caused by incomplete lineage sorting or/and ancestral hybridization, and (3) differentiating between search strategies involving genome information from one or several lineages. When the new statistics are applied, in non-conflicting cases a minimum of three elements present in both of two species and absent in a third group are considered significant support (p<0.05) for the branching of the third from the other two, if all three of the given species are screened equally for genome or experimental data. Five elements are necessary for significant support (p<0.05) if a diagnostic locus derived from only one of three species is screened, and no conflicting markers are detected. Most potentially conflicting patterns can be evaluated for their significance and ancestral hybridization can be distinguished from incomplete lineage sorting by considering symmetric or asymmetric distribution of rare genomic changes among possible tree configurations. Additionally, we provide an R-application to make the new KKSC insertion significance test available for the scientific community at http://retrogenomics.uni-muenster.de:3838/KKSC_significance_test/.     

Retroposon insertions and the chronology of avian sex chromosome evolution

The vast majority of extant birds possess highly differentiated Z and W sex chromosomes. Nucleotide sequence data from gametologs (homologs on opposite sex chromosomes) suggest that this divergence occurred throughout early bird evolution via stepwise cessation of recombination between identical sex chromosomal regions. Here, we investigated avian sex chromosome differentiation from a novel perspective, using retroposon insertions and random insertions/deletions for the reconstruction of gametologous gene trees. Our data confirm that the CHD1Z/CHD1W genes differentiated in the ancestor of the neognaths, whereas the NIPBLZ/NIPBLW genes diverged in the neoavian ancestor and independently within Galloanserae. The divergence of the ATP5A1Z/ATP5A1W genes in galloanserans occurred independently in the chicken, the screamer, and the ancestor of duck-related birds. In Neoaves, this gene pair differentiated in each of the six sampled representatives, respectively. Additionally, three of our investigated loci can be utilized as universal, easy-to-use independent tools for molecular sexing of Neoaves or Neognathae.