Here by whole genome sequencing away from 55 honey bees by developing a top solution recombination map when you look at the honey-bee, i unearthed that crossovers try of this GC content, nucleotide diversity, and you can gene thickness. We and additionally affirmed the previous idea one to genetics expressed within the personnel thoughts has oddly large CO cost. Our very own research support the take a look at you to definitely diversification away from staff choices, however immune form, is actually a driver of one’s large crossing-more than rates when you look at the bees. We discover zero facts that the crossing-more than speed was with a top NCO speed.
Steps and material
Four colonies off honeybees (Apis mellifera ligustica Spin) was indeed compiled out-of good bee farm inside the Zhenjiang, China. Each nest contained you to definitely queen, dozens of drones, and you can numerous gurus. Bees of three territories were picked for whole genome sequencing.
The fresh DNA of each personal are extracted using phenol/chloroform/isoamyl alcoholic drinks approach. To minimize the risk of bacterial contaminants, the stomachs away from bees were eliminated before DNA extraction. Regarding step 3 ?g off DNA away from for every decide to try were used to have entire genome resequencing just like the leftover DNA try left having PCR and you can Sanger sequencing. Framework of your DNA libraries and you will Illumina sequencing have been did in the BGI-Shenzhen. When you look at the short term, paired-stop sequencing libraries having input sized five-hundred bp have been built per test according to maker’s directions. Following 2 ? one hundred bp matched up-prevent reads was basically made on the IlluminaHiSEq 2000. The latest queens was in fact sequenced in the whenever 67? coverage an average of, drones from the whenever thirty-five? exposure, and pros during the just as much as 29? coverage (Table S1 during the Extra document dos). The new sequences was basically deposited on the GenBank database (accession zero. SRP043350).
SNP calling and you will marker personality
Honeybee resource genome try installed of NCBI . This new sequencing checks out was earliest mapped on to site genome having bwa following realigned with stampy . Upcoming regional realignment up to indels are did from the Genome Investigation Toolkit (GATK) , and versions had been named of the GATK https://datingranking.net/her-review/ UnifiedGenotyper.
Considering the down accuracy regarding getting in touch with indel alternatives, merely identified SNPs are used because the markers. First, 920,528 so you can 960,246 hetSNPs was in fact entitled during the for each and every queen (Dining table S2 into the Even more document 2). Next, as much as twenty two% ones was removed because these sites are also hetSNPs into the one or more haploid drone (this might mirror non-allelic series alignments caused by CNVs, sequencing mistake, or low sequencing quality). Similar size of brand new hetSNPs plus have been found in people sperm sequencing . Eventually, 671,690 so you can 740,763 credible hetSNPs within the for each colony were used just like the markers to help you place recombination occurrences (Desk S2 into the Most file 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.