The brand new four chromosomes found in the LD check and additionally displayed inversion-typical models from the PCA (Fig
Concept component analyses
1b, d, f, h, idea component loadings: A lot more document step one: Rates S2–S5). The three autosomal inversions got two fundamental homozygote haplotype clusters (for the heterozygous someone in-between) in addition to sex chromosome put into three main homozygote haplotype clusters (to your heterozygous anyone in between). New groups was in fact well-defined to your autosomes however, on the chromosome TguZ at least popular haplotype (haplotype C during the Fig. 1h) seemed to create specific recombination with every of the two almost every other haplotypes, deciding to make the groups far more diffuse. not, both low mediocre heterozygosity in this for every group out of homozygotes opposed in order to heterozygotes (Dining table 2) and you can average-joining channels (having fun with Community v4.six.1.step one that have basic setup ) into phased SNP data on inversion breakpoint (playing with Beagle v3.step three.dos ; Most document step 1: Shape S6) then support the interpretation that the LD regions depict inversion polymorphisms. It has to also be noted you to chromosomes Tgu5 and you can TguZ got come before https://hookupdaddy.net/ receive cytogenetically to take pericentric inversions additionally the breakpoints match accurately towards the LD part limitations [forty-five, 48–50].
In the newest analyses we do not understand with certainty and this arrangement are ancestral and now we therefore name him or her predicated on its allele frequency (An effective = major haplotype, B = small haplotype, C = minimum common haplotype toward chromosome TguZ; Fig. 1b, d, f, h; Dining table 2). The major alleles of all four inversion polymorphisms presented remarkably similar wavelengths varying between 0.53 and 0.sixty (Desk dos). With the chromosome TguZ, minimum of common allele (haplotype C) are unusual (regularity 0.074; Table 2). Every inversion polymorphisms had been for the Sturdy–Weinberg equilibrium (HWE; Desk 2) there is no LD between the two, which means they segregate alone (Additional document dos: Table S1).
Pooled heterozygosity and you will lesser allele counts during the inversion breakpoints
We calculated pooled heterozygosity (ZHp) in 50-kb non-overlapping sliding windows along each chromosome (Fig. 2a). Low values of ZHp are indicative of regions with a high degree of fixation, for example, due to positive selection ; high values of ZHp are expected, for example, in regions of local population structure (like inversions) or under balancing selection . We found pronounced peaks in ZHp at the presumed breakpoints of the inversions on chromosomes Tgu5, Tgu11, and Tgu13, whereas ZHp dropped to almost genome-wide average values in the interior of the inversions. Chromosome Tgu11 had only one such peak, suggesting that the proximal breakpoint is missing in the current genome assembly. Diversity (SNPs per site in a 50-kb window; Additional file 1: Figure S7) was slightly reduced at the presumed breakpoints of every inversion compared to the inversion’s interior (mean SNPs per site ± standard deviation at breakpoints versus interior of 0.017 ± 0.005 versus 0.020 ± 0.005 for Tgu5, 0.0057 ± 0.0036 versus 0.018 ± 0.004 for Tgu11, and 0.016 ± 0.006 versus 0.022 ± 0.004 for Tgu13; 0.021 ± 0.007 collinear autosomal genome-wide average SNPs per site). On chromosome TguZ, the entire inversion interior had high ZHp values, which only dropped to the genome-wide average outside the inverted region. Further, diversity on TguZ was markedly reduced all along the inverted region, including the presumed breakpoints, and increased to the genome-wide average only outside the inversion (0.0021 ± 0.0015 versus 0.022 ± 0.009, respectively).
a Pooled heterozygosity (ZHp) in 50-kb windows along each chromosome in the zebra finch genome. b–e For the highlighted areas in a, which are the presumed inversion breakpoints on the autosomes and the entire inversion interior on the sex chromosome, the minor allele count frequency (MAC) spectra are shown for chromosome Tgu5 with a local maximum at 0.34–0.36 and a frequency of the minor (B) haplotype in the sample of 0.35 (b), Tgu11 with a local maximum at 0.48–0.50 and a frequency of minor (B) haplotype in the sample of 0.47 (c), Tgu13 with a local maximum at 0.48–0.50 and a frequency of minor (B) haplotype in the sample of 0.50 (d), and TguZ with two local maxima at 0.28–0.30 and 0.42–0.44 and a frequency of the B haplotype in the sample of 0.30 and frequency of the major (A) haplotype in the sample of 0.63 (e). f For comparison, the MAC of all remaining SNPs peaks at an allele frequency of around 0.1 because SNPs with a lower frequency were not unambiguously called
