Right here from the entire genome sequencing away from 55 honey bees and by design a high quality recombination map in honey-bee, i learned that crossovers was associated with GC stuff, nucleotide diversity, and gene thickness. I also verified the former suggestion you to definitely family genes expressed for the staff brains possess unusually large CO prices. The research hold the take a look at you to definitely diversity out-of employee choices, yet not protected means, try a motorist of the large crossing-more rate within the bees. We find no proof the crossing-more than rate is actually followed closely by a premier NCO rates.

Procedures and content

Four territories off honeybees (Apis mellifera ligustica Spin) was indeed amassed away from a great bee farm when you look at the Zhenjiang, China. For each and every colony consisted of you to definitely king, all those drones, and numerous pros. Bees of about three territories was indeed chosen to possess whole genome sequencing.

The brand new DNA of each personal are extracted having fun with phenol/chloroform/isoamyl alcohol approach. To attenuate the possibility of bacterial contaminants, the new abdomens out-of bees was indeed removed before DNA removal. In the 3 ?g out-of DNA from per test were used to possess whole genome resequencing as remaining DNA are remaining to own PCR and you may Sanger sequencing. Construction of your own DNA libraries and you will Illumina sequencing was performed within BGI-Shenzhen. Inside the short-term, paired-end sequencing libraries having type sized five-hundred bp was in fact created for each decide to try according to the maker’s instructions. Following dos ? one hundred bp coordinated-stop checks out were produced with the IlluminaHiSEq 2000. The newest queens had been sequenced at as much as 67? publicity normally, drones from the approximately thirty-five? publicity, and you can workers within just as much as 29? visibility (Table S1 during the Extra document 2). Brand new sequences had been placed from the GenBank database (accession no. SRP043350).

SNP getting in touch with and marker identity

Honeybee resource genome try installed of NCBI . The brand new sequencing checks out were very first mapped onto reference genome that have bwa and realigned having stampy . Following local realignment as much as indels try performed of the Genome Studies Toolkit (GATK) , and you can variations was in fact entitled of the GATK UnifiedGenotyper.

Due to the all the way down reliability regarding contacting indel variations, only understood SNPs can be used while the markers. Very first, 920,528 to 960,246 hetSNPs was entitled during the for every single king (Table S2 from inside the More file 2). After that, whenever 22% of these have been got rid of because the internet sites are hetSNPs when you look at the one or more haploid drone (this may mirror non-allelic sequence alignments as a result of CNVs, sequencing error, or reduced sequencing quality). Comparable size of the brand new hetSNPs and have been present in person jizz sequencing . Finally, 671,690 in order to 740,763 reliable hetSNPs when you look at the for each and every colony were utilized once the markers so you’re able to discover recombination events (Dining table S2 inside the Extra file 2).

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.