Here by the entire genome sequencing out of 55 honey bees and by design a top solution recombination chart in honey bee, i learned that crossovers is from the GC posts, nucleotide range, and you can gene thickness. I also verified the previous idea one genes indicated in the staff member brains enjoys oddly large CO rates. All of our studies contain the evaluate that diversity from personnel behavior, not resistant setting, try a drivers of the higher crossing-more price inside bees. We find no research the crossing-more rates is with a leading NCO price.
Measures and you can content
Five territories away from honeybees (Apis mellifera ligustica Twist) was in fact compiled out of a good bee farm within the Zhenjiang, Asia. For each colony consisted of one to queen, those drones, and you may hundreds of professionals. Bees out of about three colonies was indeed selected to own whole genome sequencing.
The DNA each and every personal is removed using phenol/chloroform/isoamyl liquor approach. To reduce the possibility of microbial contaminants, the brand new stomachs away from bees was basically removed prior to DNA extraction. On step three ?g off DNA from for every attempt were utilized getting entire genome resequencing as the remaining DNA try leftover for PCR and you will Sanger sequencing. Construction of your DNA libraries and you may Illumina sequencing were performed in the BGI-Shenzhen. Inside brief, paired-prevent sequencing libraries having enter measurements of five hundred bp have been developed per attempt according to the manufacturer’s advice. After that dos ? a hundred bp matched up-prevent reads was in fact generated for the IlluminaHiSEq 2000. The new queens had been sequenced at the around 67? visibility on average, drones in the as much as thirty five? publicity, and you can specialists during the everything 30? exposure (Table S1 in A lot more file dos). The fresh sequences was basically deposited regarding the GenBank database (accession no. SRP043350).
SNP getting in touch with and you may marker identification
Honeybee site genome is actually downloaded regarding NCBI . New sequencing checks out had been earliest mapped on to reference genome with bwa and realigned that have stampy . Following regional realignment around indels try did because of the Genome Studies Toolkit (GATK) , and you may variations was named by GATK UnifiedGenotyper.
Considering the lower reliability regarding getting in touch with indel variants, only recognized SNPs can be used just like the indicators. Basic, 920,528 in order to 960,246 hetSNPs had been called in the for every single queen (Table S2 for the Most file dos). Next, as much as 22% ones had been got rid of because those sites are hetSNPs for the at least one haploid drone (this might echo non-allelic series alignments for the reason that CNVs, sequencing error, otherwise low sequencing high quality). Similar size of the fresh new hetSNPs and additionally were found in people spunk sequencing . Ultimately, 671,690 to 740,763 reputable hetSNPs into the for each and every colony were utilized given that indicators to help you choose recombination situations https://datingranking.net/local-hookup/corpus-christi/ (Dining table S2 from inside the Most file 2).
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.