Phenotypic and genetic characterization of Africanized Apis mellifera colonies with natural tolerance to Varroa destructor and contrasting defensive behavior

Africanized Apis mellifera colonies with promising characteristics for beekeeping have been detected in northern Argentina (subtropical climate) and are considered of interest for breeding programs. Integral evaluation of this feral material revealed high colony strength and resistance/tolerance to brood diseases. However, these Africanized honeybees (AHB) also showed variable negative behavioral traits for beekeeping, such as defensiveness, tendency to swarm and avoidance behavior. We developed a protocol for the selection of AHB stocks based on defensive behavior and characterized contrasting colonies for this trait using NGS technologies. For this purpose, population and behavioral parameters were surveyed throughout a beekeeping season in nine daughter colonies obtained from a mother colony (A1 mitochondrial haplotype) with valuable characteristics (tolerance to the mite Varroa destructor , high colony strength and low defensiveness). A Defensive Behavior Index was developed and tested in the colonies under study. Mother and two daughter colonies displaying contrasting defensive behavior were analyzed by ddRADseq. High-quality DNA samples were obtained from 16 workers of each colony. Six pooled samples, including two replicates of each of the three colonies, were processed. A total of 12,971 SNPs were detected against the reference genome of A. mellifera , 142 of which showed significant differences between colonies. We detected SNPs in coding regions, lncRNA, miRNA, rRNA, tRNA, among others. From the original data set, we also identified 647 SNPs located in protein-coding regions, 128 of which are related to 21 genes previously associated with defensive behavior, such as dop3 and dopR2 , CaMKII and ADAR , obp9 and obp10 , and members of the 5-HT family. We discuss the obtained results by considering the influence of polyandry and paternal lineages on the defensive behavior in AHB and provide baseline information to use this innovative molecular approach, ddRADseq, to assist in the selection and evaluation of honey bee stocks showing low defensive behavior for commercial uses.