Cargando…

Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here,...

Descripción completa

Detalles Bibliográficos
Autores principales: Gmel, Annik Imogen, Guichard, Matthieu, Dainat, Benjamin, Williams, Geoffrey Rhys, Eynard, Sonia, Vignal, Alain, Servin, Bertrand, Neuditschko, Markus
Formato: Online Artículo Texto
Lenguaje:English
Publicado: John Wiley and Sons Inc. 2023
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9843643/
https://www.ncbi.nlm.nih.gov/pubmed/36694553
http://dx.doi.org/10.1002/ece3.9723
_version_ 1784870454131949568
author Gmel, Annik Imogen
Guichard, Matthieu
Dainat, Benjamin
Williams, Geoffrey Rhys
Eynard, Sonia
Vignal, Alain
Servin, Bertrand
Neuditschko, Markus
author_facet Gmel, Annik Imogen
Guichard, Matthieu
Dainat, Benjamin
Williams, Geoffrey Rhys
Eynard, Sonia
Vignal, Alain
Servin, Bertrand
Neuditschko, Markus
author_sort Gmel, Annik Imogen
collection PubMed
description Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole‐genome queen genotypes, from a pool‐seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well‐defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F (ROH)) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F (ROH) (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F (ROH) (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F (ROH) of 2.33%. Within the breed‐specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh‐BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.
format Online
Article
Text
id pubmed-9843643
institution National Center for Biotechnology Information
language English
publishDate 2023
publisher John Wiley and Sons Inc.
record_format MEDLINE/PubMed
spelling pubmed-98436432023-01-23 Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data Gmel, Annik Imogen Guichard, Matthieu Dainat, Benjamin Williams, Geoffrey Rhys Eynard, Sonia Vignal, Alain Servin, Bertrand Neuditschko, Markus Ecol Evol Research Articles Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole‐genome queen genotypes, from a pool‐seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well‐defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F (ROH)) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F (ROH) (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F (ROH) (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F (ROH) of 2.33%. Within the breed‐specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh‐BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees. John Wiley and Sons Inc. 2023-01-17 /pmc/articles/PMC9843643/ /pubmed/36694553 http://dx.doi.org/10.1002/ece3.9723 Text en © 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. https://creativecommons.org/licenses/by/4.0/This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ (https://creativecommons.org/licenses/by/4.0/) License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
spellingShingle Research Articles
Gmel, Annik Imogen
Guichard, Matthieu
Dainat, Benjamin
Williams, Geoffrey Rhys
Eynard, Sonia
Vignal, Alain
Servin, Bertrand
Neuditschko, Markus
Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title_full Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title_fullStr Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title_full_unstemmed Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title_short Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole‐genome sequencing data
title_sort identification of runs of homozygosity in western honey bees (apis mellifera) using whole‐genome sequencing data
topic Research Articles
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9843643/
https://www.ncbi.nlm.nih.gov/pubmed/36694553
http://dx.doi.org/10.1002/ece3.9723
work_keys_str_mv AT gmelannikimogen identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT guichardmatthieu identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT dainatbenjamin identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT williamsgeoffreyrhys identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT eynardsonia identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT vignalalain identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT servinbertrand identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata
AT neuditschkomarkus identificationofrunsofhomozygosityinwesternhoneybeesapismelliferausingwholegenomesequencingdata