PATTERNS OF DESCENT IN CLONAL LINEAGES AND THEIR MULTILOCUS FINGERPRINTS ARE RESOLVED WITH COMBINED GENE GENEALOGIES

Clonal lineages in the filamentous ascomycete (fungi) Sclerotinia sclerotiorum were determined by analysis of genealogies of four loci: the intergenic spacer of the nuclear ribosomal repeat (IGS; approximately 4 kb), the translation elongation factor (EF-1α; approximately 300 bp), an anonymous region (44.11; approximately 700 bp), and the calmodulin gene (CAL; approximately 400 bp). Three of the four loci are physically unlinked. The combined analysis of the four loci provided the best estimate of phylogeny, which is consistent with a pattern of some recombination among clonal lineages against a background of predominant clonality. Comparison of gene genealogies with a phylogeny inferred from DNA fingerprints and a combined phylogeny of the entire dataset identified convergent or parallel changes in fingerprints. Analysis of the entire data matrix allowed us to resolve patterns of descent among clonal lineages that could not be inferred from fingerprints alone and to discern recent episodes of divergence that were not detected in gene genealogies. Prerequisites for applying this approach to other systems are a haploid context for inferring multiple gene genealogies (such as the mitochondrial genome) that indicate limited recombination and another data matrix that identifies recently evolved genotypes.     

A microbial population-species interface: nested cladistic and coalescent inference with multilocus data

Using sequence data from seven nuclear loci in 385 isolates of the haploid, plant parasitic, ascomycete fungus, Sclerotinia, divergence times of populations and of species were distinguished. The evolutionary history of haplotypes on both population and species scales was reconstructed using a combination of parsimony, maximum likelihood and coalescent methods, implemented in a specific order. Analysis of site compatibility revealed recombination blocks from which alternative (marginal) networks were inferred, reducing uncertainty in the network due to recombination. Our own modifications of Templeton and co-workers' cladistic inference method and a coalescent approach detected the same phylogeographic processes. Assuming neutrality and a molecular clock, the boundary between divergent populations and species is an interval of time between coalescence (to a common ancestor) of populations and coalescence of species.     

Patterns of DNA Variability at X-Linked Loci in Mus Domesticus

Introns of four X-linked genes (Hprt, Plp, Glra2, and Amg) were sequenced to provide an estimate of nucleotide diversity at nuclear genes within the house mouse and to test the neutral prediction that the ratio of intraspecific polymorphism to interspecific divergence is the same for different loci. Hprt and Plp lie in a region of the X chromosome that experiences relatively low recombination rates, while Glra2 and Amg lie near the telomere of the X chromosome, a region that experiences higher recombination rates. A total of 6022 bases were sequenced in each of 10 Mus domesticus and one M. caroli. Average nucleotide diversity (π) for introns within M. domesticus was quite low (π = 0.078%). However, there was substantial variation in the level of heterozygosity among loci. The two telomeric loci, Glra2 and Amg, had higher ratios of polymorphism to divergence than the two loci experiencing lower recombination rates. These results are consistent with the hypothesis that heterozygosity is reduced in regions with lower rates of recombination, although sampling of additional genes is needed to establish whether there is a general correlation between heterozygosity and recombination rate as in Drosophila melanogaster.     

Clonal composition of the peach-potato aphid Myzus persicae (Homoptera: Aphididae) in France and Scotland: Comparative analysis with IGS fingerprinting and microsatellite markers

Fourteen colonies of the peach‐potato aphid, Myzus persicae, were taken either from French peach trees or weeds in 2001. Thirty five apomictic parthenogenetic lineages (APLs) were established. Ribosomal DNA intergenic spacer (IGS) fingerprinting was used to characterise these and 28 fingerprints were duly obtained. Those lineages with different fingerprints were considered different genotypes and those with the same fingerprint as the same. The genetic identity of APLs was further tested using four microsatellite loci. APLs that differed by IGS fingerprint had distinct microsatellite allele combinations and those that had the same IGS fingerprint had the same microsatellite allele combinations. The results confirmed that IGS types corresponded to different aphid genotypes. Independent APLs with identical IGS and microsatellite genotype were therefore considered different representatives of the same clone. APLs from M. persicae found on Scottish crops in 1995, 1996 and 2001, as well as a long‐term laboratory line were also examined by the same methods. Their IGS fingerprints were similar or identical suggesting that they all belonged to the same clone. Microsatellite markers also suggested that these lineages were derived from a single clone. Some field lineages exhibited slight modifications to their IGS fingerprints confirming that the IGS evolves more rapidly than these microsatellite alleles. Thus, IGS will continue to provide a useful marker for aphid fieldwork.     

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.     

Kinetic and segregational analysis of mitochondrial DNA recombination in yeast

A pair of yeast strains of opposite mating type was constructed to contain polymorphisms at three loci on the mitochondrial genome--the 21 S rRNA gene, var1, and cob--such that parental and recombinant forms of these genes could be easily detected by Southern blot analysis. These polymorphisms were used to measure in a single cross gene conversions at the 21 S rRNA and var1 loci and a reciprocal recombination at cob. For all three loci, recombination initiates at about the same time, 4 to 6 h after mixing cells, and increases with similar kinetics over a 24-h period. The segregation of parental and recombinant forms of these genes was then followed by pedigree analysis. The results, which show a high variance in the distribution of parental and recombinant forms of all three genes in cells derived from both the first bud and the mother zygote, are consistent with the segregation of a small number of mitochondrial DNA molecules from the zygote to diploid buds. Based on these results and previous experiments of this type, a limited "zone of mixing" of parental mitochondrial DNA molecules probably exists in the zygote. The extent of sampling from this zone, together with the intrinsic properties of the recombination events themselves, is likely to determine the observed pattern of recombination of mitochondrial DNA sequences at the population level.     

Recombination in RNA

The aphthovirus genome consists of a single molecule of single-stranded RNA that encodes all the virus-induced proteins. We isolated recombinant aphthoviruses from cells simultaneously infected with temperature-sensitive mutants of two different subtype strains. Analysis of the proteins induced by 16 independently generated recombinants revealed two types of protein pattern, which were consistent with single genetic crossovers on the 5' side and 3' side, respectively, of the central P34-coding region. Recombinants invariably inherited all four coat proteins from the same parent, and novel recombinant proteins were not observed. RNAase T1 fingerprints of virus RNA, prepared from representatives of each recombinant type, confirmed the approximate crossover sites that had been deduced from the inheritance of proteins. These fingerprints provide molecular evidence of recombination at the level of RNA and demonstrate the potential of RNA recombination for producing genetic diversity among picornaviruses.     

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated F_ST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated F_ST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of F_ST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased F_ST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated F_ST in geographically distant locations. We discuss possible explanations for this unexpected result.     

Indica and japonica crosses resulting in linkage block and recombination suppression on rice chromosome 12

Understanding linkage block size and molecular mechanisms of recombination suppression is important for plant breeding. Previously large linkage blocks ranging from 14 megabases to 27 megabases were observed around the rice blast resistance gene Pi-ta in rice cultivars and backcross progeny involving an indica and japonica cross. In the present study, the same linkage block was further examined in 456 random recombinant individuals of rice involving 5 crosses ranging from F(2) to F(10) generation, with and without Pi-ta containing genomic indica regions with both indica and japonica germplasm. Simple sequence repeat markers spanning the entire chromosome 12 were used to detect recombination break points and to delimit physical size of linkage blocks. Large linkage blocks ranging from 4.1 megabases to 10 megabases were predicted from recombinant individuals involving genomic regions of indica and japonica. However, a significantly reduced block from less than 800 kb to 2.1megabases was identified from crosses of indica with indica rice regardless of the existence of Pi-ta. These findings suggest that crosses of indica and japonica rice have significant recombination suppression near the centromere on chromosome 12.     

DNA variation in a conifer,Cryptomeria japonica (Cupressaceae sensu lato)

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.     

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F _ST) and sequence divergence (d _XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for speciation in the face of gene flow.     

Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune

The common split-gilled mushroom Schizophyllum commune is found throughout the world on woody substrates. This study addresses the dispersal and population structure of this fungal species by studying the phylogeny and evolutionary dynamics of ribosomal DNA (rDNA) spacer regions. Extensive sampling (n = 195) of sequences of the intergenic spacer region (IGS1) revealed a large number of unique haplotypes (n = 143). The phylogeny of these IGS1 sequences revealed strong geographic patterns and supported three evolutionarily distinct lineages within the global population. The same three geographic lineages were found in phylogenetic analysis of both other rDNA spacer regions (IGS2 and ITS). However, nested clade analysis of the IGS1 phylogeny suggested the population structure of S. commune has undergone recent changes, such as a long distance colonization of western North America from Europe as well as a recent range expansion in the Caribbean. Among all spacer regions, variation in length and nucleotide sequence was observed between but not within the tandem rDNA repeats (arrays). This pattern is consistent with strong within-array and weak among-array homogenizing forces. We present evidence for the suppression of recombination between rDNA arrays on homologous chromosomes that may account for this pattern of concerted evolution.     

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.     

Recombination and population mosaic of a multifunctional viral gene, adeno-associated virus cap

Homologous recombination is a dominant force in evolution and results in genetic mosaics. To detect evidence of recombination events and assess the biological significance of genetic mosaics, genome sequences for various viral populations of reasonably large size are now available in the GenBank. We studied a multi-functional viral gene, the adeno-associated virus (AAV) cap gene, which codes for three capsid proteins, VP1, VP2 and VP3. VP1-3 share a common C-terminal domain corresponding to VP3, which forms the viral core structure, while the VP1 unique N-terminal part contains an enzymatic domain with phospholipase A2 activity. Our recombinant detection program (RecI) revealed five novel recombination events, four of which have their cross-over points in the N-terminal, VP1 and VP2 unique region. Comparison of phylogenetic trees for different cap gene regions confirmed discordant phylogenies for the recombinant sequences. Furthermore, differences in the phylogenetic tree structures for the VP1 unique (VP1u) region and the rest of cap highlighted the mosaic nature of cap gene in the AAV population: two dominant forms of VP1u sequences were identified and these forms are linked to diverse sequences in the rest of cap gene. This observation together with the finding of frequent recombination in the VP1 and 2 unique regions suggests that this region is a recombination hot spot. Recombination events in this region preserve protein blocks of distinctive functions and contribute to convergence in VP1u and divergence of the rest of cap. Additionally the possible biological significance of two dominant VP1u forms is inferred.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

The relationship of nucleotide polymorphism, recombination rate and selection in wild tomato species

We analyzed the effects of mating system and recombination rate on single nucleotide polymorphisms using 14 single-copy nuclear loci from single populations of five species of wild tomatoes (Solanum section Lycopersicon). The taxa investigated comprise two self-compatible (SC) and three self-incompatible (SI) species. The observed reduction in nucleotide diversity in the SC populations compared to the SI populations is much stronger than expected under the neutral effects of the mating system on effective population size. Importantly, outgroup sequences available for 11 of the 14 loci yield strong positive correlations between silent nucleotide diversity and silent divergence, indicative of marked among-locus differences in mutation rates and/or selective constraints. Furthermore, using a physical estimate of local recombination rates, we find that silent nucleotide diversity (but not divergence) is positively correlated with recombination rate in two of the SI species. However, this correlation is not nearly as strong as in other well-characterized species (in particular, Drosophila). We propose that nucleotide diversity in Lycopersicon is dominated mainly by differences in neutral mutation rates and/or selective constraints among loci, demographic processes (such as population subdivision), and background selection. In addition, we hypothesize that the soil seed bank plays an important role in the maintenance of the large genetic diversity in the SI species (in particular L. peruvianum).     

Monitoring the mode and tempo of concerted evolution in the Drosophila melanogaster rDNA locus

Non-LTR retrotransposons R1 and R2 have persisted in rRNA gene loci (rDNA) since the origin of arthropods despite their continued elimination by the recombinational mechanisms of concerted evolution. This study evaluated the short-term evolutionary dynamics of the rDNA locus by measuring the divergence among replicate Drosophila melanogaster lines after 400 generations. The total number of rDNA units on the X chromosome of each line varied from 140 to 310, while the fraction of units inserted with R1 and R2 retrotransposons ranged from 37 to 65%. This level of variation is comparable to that found in natural population surveys. Variation in locus size and retrotransposon load was correlated with large changes in the number of uninserted and R1-inserted units, yet the numbers of R2-inserted units were relatively unchanged. Intergenic spacer (IGS) region length variants were also used to evaluate changes in the rDNA loci. All IGS length variants present in the lines showed significant increases and decreases of copy number. These studies, combined with previous data following specific R1 and R2 insertions in these lines, help to define the type and distribution, both within the locus and within the individual units, of recombinational events that give rise to the concerted evolution of the rDNA locus.     

Estimating divergence time and ancestral effective population size of Bornean and Sumatran orangutan subspecies using a coalescent hidden Markov model

Due to genetic variation in the ancestor of two populations or two species, the divergence time for DNA sequences from two populations is variable along the genome. Within genomic segments all bases will share the same divergence-because they share a most recent common ancestor-when no recombination event has occurred to split them apart. The size of these segments of constant divergence depends on the recombination rate, but also on the speciation time, the effective population size of the ancestral population, as well as demographic effects and selection. Thus, inference of these parameters may be possible if we can decode the divergence times along a genomic alignment. Here, we present a new hidden Markov model that infers the changing divergence (coalescence) times along the genome alignment using a coalescent framework, in order to estimate the speciation time, the recombination rate, and the ancestral effective population size. The model is efficient enough to allow inference on whole-genome data sets. We first investigate the power and consistency of the model with coalescent simulations and then apply it to the whole-genome sequences of the two orangutan sub-species, Bornean (P. p. pygmaeus) and Sumatran (P. p. abelii) orangutans from the Orangutan Genome Project. We estimate the speciation time between the two sub-species to be thousand years ago and the effective population size of the ancestral orangutan species to be , consistent with recent results based on smaller data sets. We also report a negative correlation between chromosome size and ancestral effective population size, which we interpret as a signature of recombination increasing the efficacy of selection.     

The Effect of Attachment Site Mutations on Strand Exchange in Bacteriophage λ Site-Specific Recombination

Recombination of phage λ attachment sites occurs by sequential exchange of the DNA strands at two specific locations. The first exchange produces a Holliday structure, and the second resolves it to recombinant products. Heterology for base substitution mutations in the region between the two strand exchange points (the overlap region) reduces recombination; some mutations inhibit the accumulation of Holliday structures, others inhibit their resolution to recombinant products. To see if heterology also alters the location of the strand exchange points, we determined the segregation pattern of three single and one multiple base pair substitution mutations of the overlap region in crosses with wild type sites. The mutations are known to differ in the severity of their recombination defect and in the stage of strand exchange they affect. The three single mutations behaved similarly: each segregated into both products of recombination,, and the two products of a single crossover were frequently nonreciprocal in the overlap region. In contrast, the multiple mutation preferentially segregated into one of the two recombinant products, and the two products of a single crossover appeared to be fully reciprocal. The simplest explanation of the segregation pattern of the single mutations is that strand exchanges occur at the normal locations to produce recombinants with mismatched base pairs that are frequently repaired. The segregation pattern of the multiple mutation is consistent with the view that both strand exchanges usually occur to one side of the mutant site. We suggest that the segregation pattern of a particular mutation is determined by which stage of strand exchange it inhibits and by the severity of the inhibition.     

Speciation genomics and a role for the Z chromosome in the early stages of divergence between Mexican ducks and mallards

Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and Mexican duck (A. [p.] diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (Φ_ST = 0.088) were 5.2 times more differentiated than autosomal DNA (Φ_ST = 0.017) but comparable to mtDNA (Φ_ST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (Φ_ST) and absolute (d_XY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (Φ_ST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.