Genome-wide candidate regions for selective sweeps revealed through massive parallel sequencing of DNA across ten turkey populations

Sequencing and characterization of complete mitogenome DNA of worldwide turkey (Meleagris gallopavo) populations

The history of turkey (Meleagris gallopavo) domestication can be traced back to the period between 700 and 200 BC in Mexico. This process involved multiple contributors and resulted in the development of modern local turkey breeds. This research investigates the complete mitochondrial diversity across a diverse range of local turkeys. Seventy-three turkeys were sampled from various populations, including autochthonous Italian breeds, an American breed (Narragansett), as well as wild turkeys from the USA and Mexico. The mitochondrial DNA (mtDNA) was employed as a powerful tool for biodiversity and breed phylogeny investigation. An analysis of the entire mtDNA was conducted to identify breed-specific unique traits, mitochondrial-specific characteristics, and the phylogenetic relationship among turkey populations. A total of 44 polymorphic sites were identified. Brianzolo and Narragansett birds were characterized as genetically homogeneous populations. Thirty-two different haplotypes were identified when our samples were compared with mtDNA D-loop of 96 online available turkeys from various geographical countries. H1 and H2, differing for one mutation, were the most abundant, comprising 132 of the 185 sequences. H1 included samples coming from every region, while H2 was predominantly characterized by Italian samples. USA and Mexican samples appear to be more variable in their mtDNA than the other populations.

Investigating inbreeding in the turkey (Meleagris gallopavo) genome

The detrimental effects of increased homozygosity due to inbreeding have prompted the development of methods to reduce inbreeding. The detection of runs of homozygosity (ROH), or contiguous stretches of homozygous marker genotypes, can be used to describe and quantify the level of inbreeding in an individual. The estimation of inbreeding coefficients can be calculated based on pedigree information, ROH, or the genomic relationship matrix. The aim of this study was to detect and describe ROH in the turkey genome and compare estimates of pedigree-based inbreeding coefficients (F_PED) with genomic-based inbreeding coefficients estimated from ROH (F_ROH) and the genomic relationship matrix (F_GRM). A total of 2,616,890 pedigree records were available. Of these records, 6,371 genotyped animals from three purebred turkey (Meleagris gallopavo) lines between 2013 and 2019 were available, and these were obtained using a dense single nucleotide polymorphism array (56,452 SNPs). The overall mean length of detected ROH was 2.87 ± 0.29 Mb with a mean number of 84.87 ± 8.79 ROH per animal. Short ROH with lengths of 1 to 2 Mb long were the most abundant throughout the genome. Mean ROH coverage differed greatly between chromosomes and lines. Considering inbreeding coefficient means across all lines, genomic derived inbreeding coefficients (F_ROH = 0.27; F_GRM = 0.32) were higher than coefficients estimated from pedigree records (F_PED = 0.14). Correlations between F_ROH and F_PED, F_ROH and F_GRM, and F_PED and F_GRM ranged between 0.19 to 0.31, 0.68 to 0.73, and 0.17 to 0.30, respectively. Additionally, correlations between F_ROH from different lengths and F_PED substantially increased with ROH length from -0.06 to 0.33. Results of the current research, including the distribution of ROH throughout the genome and ROH-derived inbreeding estimates, can provide a more comprehensive description of inbreeding in the turkey genome. This knowledge can be used to evaluate genetic diversity, a requirement for genetic improvement, and develop methods to minimize inbreeding in turkey breeding programs.

Genetic Diversity and Identification of Homozygosity-Rich Genomic Regions in Seven Italian Heritage Turkey (Meleagris gallopavo) Breeds

Italian autochthonous turkey breeds are an important reservoir of genetic biodiversity that should be maintained with an in vivo approach. The aim of this study, part of the TuBAvI national project on biodiversity, was to use run of homozygosity (ROH), together with others statistical approaches (e.g., Wright's F-statistics, principal component analysis, ADMIXTURE analysis), to investigate the genomic diversity in several heritage turkey breeds. We performed a genome-wide characterization of ROH-rich regions in seven autochthonous turkey breeds, i.e., Brianzolo (Brzl), Bronzato Comune Italiano (BrCI), Bronzato dei Colli Euganei (CoEu), Parma e Piacenza (PrPc), Nero d'Italia (NeIt), Ermellinato di Rovigo (ErRo) and Romagnolo (Roma). ROHs were detected based on a 650K SNP genotyping. ROH_islands were identified as homozygous ROH regions shared by at least 75% of birds (within breed). Annotation of genes was performed with DAVID. The admixture analyses revealed that six breeds are unique populations while the Roma breed consists in an admixture of founder populations. Effective population size estimated on genomic data shows a numeric contraction. ROH_islands harbour genes that may be interesting for target selection in commercial populations also. Among them the PTGS2 and PLA2G4A genes on chr10 were related to reproduction efficiency. This is the first study mapping genetic variation in autochthonous turkey populations. Breeds were genetically different among them, with the Roma breed proving to be a mixture of the other breeds. The ROH_islands identified harboured genes peculiar to the selection that occurred in heritage breeds. Finally, this study releases previously undisclosed information on existing genetic variation in the turkey species.

Copy Number Variants in Four Italian Turkey Breeds

Simple Summary

Hybrid Turkey selection is focusing on meat production traits characterized by high genetic heritability; the strong directional selection is well known to produce a constant loss in genetic diversity. Genetic characterization is one of the essential activities in the management of populations at risk of extinction. In addition, the genetic structure at the population level and the relationships between individuals are nowadays analysable at the genomic level. In this paper, the genome of 4 different Italian turkey breeds included in the Autochthonous Italian Poultry Breeds Register are analysed in order to obtain a genome-wide Copy Number Variant scan to ameliorate the existing knowledge of the genomic structure of Italian local turkey breeds. Differences have been described at genomic level for physiological, reproductive, and behavioral traits. The analyzed breeds are clearly distinguishable at the genomic level, and their relationships are clearly linked to their geographical origin and to the history of the rural structure of their developing regions. Genome information based on Copy Number Variant (CNV) detection has generated important information in this study concerning the uniqueness of the Italian local turkey breeds.

Abstract

Heritage breeds can be considered a genetic reservoir of genetic variability to be conserved and valorized considering their historical, cultural, and adaptive characteristics and possibly for their high potential in commercial hybrid genetic improvement by gene introgression. The aim of the present research is to investigate via Copy Number Variant (CNVs) the genomic makeup of 4 Italian autochthonous turkey breeds (Bronzato Comune—BrCI, 24; Ermellinato di Rovigo—ErRo, 24; Parma e Piacenza—PrPc, 25; Romagnolo—RoMa, 29). CNVs detection was performed using two different software and an interbreed CNVs comparison was carried out. A total of 1077 CNVs were identified in 102 turkeys, summarized into 519 CNV regions (CNVRs), which resulted after merging in 101 and 18 breed and shared regions. Biodiversity was analyzed using the effective information supplied by CNVs analysis, and BrCI and ErRo were characterized by a low mapped CNV number. Differences were described at a genomic level related to physiological, reproductive, and behavioral traits. The comparison with other three Italian turkey breeds (Brianzolo, Colle Euganei, and Nero Italiano) using a CNV data set available in the literature showed high clustering properties at the genomic level, and their relationships are strictly linked to the geographical origin and to the history of the rural structure of their native regions.

Integrating Selection Mapping With Genetic Mapping and Functional Genomics

Genomic scans for signatures of selection allow us to, in principle, detect variants and genes that underlie recent adaptations. By combining selection mapping with genetic mapping of traits known to be relevant to adaptation, we can simultaneously investigate whether genes and variants show signals of recent selection and whether they impact traits that have likely been selected. There are three ways to integrate selection mapping with genetic mapping or functional genomics: (1) To use genetic mapping data from other populations as a form of genome annotation. (2) To perform experimental evolution or artificial selection to be able to study selected variants when they segregate, either by performing genetic mapping before selection or by crossing the selected individuals to some reference population. (3) To perform a comparative study of related populations facing different selection regimes. This short review discusses these different ways of integrating selection mapping with genetic mapping and functional genomics, with examples of how each has been done.

Gene editing in birds takes flight

The application of gene editing (GE) technology to create precise changes to the genome of bird species will provide new and exciting opportunities for the biomedical, agricultural and biotechnology industries, as well as providing new approaches for producing research models. Recent advances in modifying both the somatic and germ cell lineages in chicken indicate that this species, and conceivably soon other avian species, has joined a growing number of model organisms in the gene editing revolution.