Every five perceiver scored very specimens, specifically those off around the latest hybrid area
Morphological evaluations
We obtained mature sample peels handled in the STRIBC (N =67), the latest Western Art gallery off Natural Record (AMNH, Letter =25), and the United states National Art gallery bird range (USNM, Letter = 10) from Costa Rica, Panama, and you can Colombia for every single off three phenotypic traits: i) dorsal plumage, ii) ventral plumage, and you may iii) face shield ornamentation. Following the , MJM, Sel-, or more to 3 other perceiver (come across Additional file step 1: Table S1) scored for every bird by themselves for the a size off 0 – step 1, where 0 portrayed a good phenotypically absolute Jacana spinosa, and 1 illustrated an effective phenotypically pure Jacana jacana (based on reviews with peels off outside of the get in touch with zone), and we averaged the newest ratings for each attribute as phenotypic hybrid list for every single specimen . Specimens with phenotypic crossbreed indices less than 0.10 were categorized as “pure” spinosa, when you’re a collection from 0.ninety or maybe more is obtained because “pure” jacana. A score between 0.10 and you can 0.25 try classified given that “spinosa-like”, whereas a directory anywhere between 0.75 and you can 0.ninety are categorized because “jacana-like”. Specimens that have indices anywhere between 0.twenty-five and 0.75 have been categorized because “hybrids”.
Genetic data
I sequenced a great 652 legs pair fragment of your own mitochondrial cytochrome oxidase We gene (COI: the fresh new DNA barcode fragment) off 59 Panamanian jacanas, together with one or two specimens off eastern Honduras sent to our study because of the J. Klicka while the Marjorie Barrick Art gallery. I chose which gene when planning on taking advantageous asset of established sequence study toward Genbank getting jacanas away from varied geographical countries, that happen to be together with put into the newest dataset (Letter = 10). Research methods follow standard actions which have been explained somewhere else . We tried however, were not able to produce series research from toe mat extractions throughout the Smith and you will Emlen series of specimens.
This new COI investigation were utilized to help you guess new hereditary distinctiveness of both jacana types and infer enough time of speciation ranging from J. spinosa and J. jacana, and assess the knowledge that geography restrictions gene circulate inside the jacanas (we.elizabeth. isolation-by-distance). We generated average-joining community inside PopART . I put Super so you can determine the common distance utilising the HKY make of molecular development so you’re able to contrast divergence inside jacanas to other avian species where genetic length could have been computed with the HKY design .
We tested towards the effect of geographic point on genetic variations (i.e., separation from the point) through an excellent Mantel try followed regarding system Alleles in space , hence compares hereditary and you will geographic ranges among anybody in place of populations. For J. spinosa and you can J. jacana mtDNA succession sets, we synchronised uncorrected pairwise sequence divergence that have pairwise geographic point to have every you can easily sets. Importance of new relationship are influenced by evaluating the brand new noticed correlation so you’re able to haphazard correlations considering ten,100000 permutations.
In order Finn mer to match the newest mtDNA-built analyses off hereditary introgression across the jacana hybrid zone, we sequenced several autosomal introns (gapdh 11: , and you may 10551: ) getting STRIBC specimens accumulated across Panama. Centered on an initial evaluation out-of sequence variation out-of a few products away from for each geographical endpoint (outside the hybrid zone) in our sampling design, i recognized a single nucleotide polymorphism for each and every intron: gapdh step 3- anti-sense status 136 (cytosine-thymine: C-T change), and you may 10551 status 153 (adenine-guanine: A-Grams transition). Positions is in line with aimed sequences of each and every intron on J. spinosa try MJM8238 (Genbank accessions: gapdh: KM891734, 10551: KM891735).
Cline figure analyses
We estimated the shape of the phenotypic and genotypic clines using ClineFit , which applies the cline fitting formulas developed by Szymura and Barton . The shape of a cline is determined by three formulas. The first formula (1) defines a symmetric S-shaped curve around the center of the hybrid zone, and relies on two parameters: c – for the center of the curve, and w – which is 1/maximum slope of the cline, and represents the cline width. Two additional formulae define the cline tails (2,3) and add four additional parameters: z L and z R , which represent the distance from c to a vertical asymptote for the exponential decay of allele frequencies to the left and right sides of the hybrid zone, while ? L and ? R define the rates of this exponential decay relative to the shape of the central cline (1).