Here because of the whole genome sequencing of 55 honey bees and also by creating a high resolution recombination map when you look at the honey bee, i unearthed that crossovers is for the GC articles, nucleotide range, and gene density. I including verified the former suggestion that family genes shown during the staff member brains have strangely high CO pricing. The investigation contain the examine that variation out of staff member choices, although not immune means, is a motorist of the highest crossing-over rate within the bees. We find no research your crossing-more price was accompanied by a leading NCO rate.

Procedures and you may content

Five colonies off honeybees (Apis mellifera ligustica Twist) was indeed built-up regarding a bee ranch when you look at the Zhenjiang, China. Each nest contained one queen, dozens of drones, and you can a huge selection of experts. Bees out of three territories was basically picked having whole genome sequencing.

Brand new DNA of each personal are extracted using phenol/chloroform/isoamyl alcohol means. To minimize the possibility of microbial toxic contamination, the newest stomachs off bees was indeed removed ahead of DNA extraction. On the step three ?g from DNA out-of for each try were utilized for entire genome resequencing as the leftover DNA was remaining for PCR and Sanger sequencing. Build of your DNA libraries and Illumina sequencing were did in the BGI-Shenzhen. Inside short-term, paired-avoid sequencing libraries which have insert size of 500 bp was indeed created each test depending on the manufacturer’s rules. Then 2 ? 100 bp matched-avoid reads was produced toward IlluminaHiSEq 2000. New queens had been sequenced on everything 67? publicity normally, drones at as much as thirty five? coverage, and gurus from the just as much as 29? coverage (Table S1 within the A lot more file dos). The new sequences was indeed transferred on GenBank database (accession no. SRP043350).

SNP getting in touch with and marker personality

Honeybee site genome are installed of NCBI . Brand new sequencing reads was very first mapped on to site genome with bwa right after which realigned which have stampy . Then regional realignment doing indels is did of the Genome Study Toolkit (GATK) , and you may variations had been named by GATK UnifiedGenotyper.

Considering the down precision out of contacting indel variants, just recognized SNPs are utilized as the markers. Earliest, 920,528 so you’re able to 960,246 hetSNPs was basically titled inside for every single queen (Dining table S2 into the Additional document dos). Following, approximately twenty two% of those was basically eliminated due to the fact that web sites are hetSNPs in the one haploid drone (this might reflect non-allelic succession alignments because of CNVs, sequencing error, otherwise reduced sequencing quality). Similar dimensions of the newest hetSNPs and have been noticed in human spunk sequencing . Finally, 671,690 to 740,763 credible hetSNPs inside each nest were used because indicators to position recombination events (Table S2 from San Francisco free hookup website inside the A lot more document dos).

Haploid phasing

For each colony, the identified markers were used for haploid phasing. The linkage of every two adjacent markers was inferred to determine the two chromosome haplotypes of the queen by comparing the SNP linkage information across all drones from the same colony. Detailed methods were described in Lu’s study . In brief, for each pair of adjacent hetSNPs, for example A/G and C/T, there could be two types of link in the queen ‘A-C, G-T’ or ‘A-T, G-C’. Assuming recombination events are low probability, if more ‘A-C, G-T’ drones are found than ‘A-T, G-C’ drones, then ‘A-C, G-T’ is assumed to be the correct link in the queen and vice versa. The two haplotypes can be clearly discriminated between >99% of ple). For linkage of the <1% markers, as shown in Additional file 1: Figure S2B, between markers at ‘LG1:20555174' and ‘LG1:20555456' , there are 14 ‘A-A or G-G' type drones against 1 ‘A-G or G-A' type drone, so ‘A-A, G-G' is assumed to be the correct link in queen and a recombination event is identified at this site in sample I-9.