Search

Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (F ST : 0.011-0.043) and negligible inbreeding coefficients (F IS : -0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92-6.85 Mb, with short ROH segments (1-5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20-40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features., (© 2024 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (F ST ), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics., (© 2024 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Celtic-Iberian pig breeds were majority in Spain and Portugal until the first half of the 20th century. In the 1990s, they were nearly extinct as a result of the introduction of foreign improved pig breeds. Despite its historical importance, the genetic background of the Celtic-Iberian pig strain is poorly documented. In this study, we have identified genomic regions that might contain signatures of selection peculiar of the Celtic-Iberian genetic lineage. A total of 153 DNA samples of Celtic-Iberian pigs (Spanish Gochu Asturcelta and Portuguese Bísara breeds), Iberian pigs (Spanish Iberian and Portuguese Alentejano breeds), Cinta Senese pig, Korean local pig and Cosmopolitan pig (Hampshire, Landrace and Large White individuals) were analysed. A pairwise-comparison approach was applied: the Gochu Asturcelta and the Bísara samples as test populations and the five other pig populations as reference populations. Three different statistics (XP-EHH, F ST and ΔDAF) were computed on each comparison. Strict criteria were used to identify selection sweeps in order to reduce the noise brought on by the Gochu Asturcelta and Bísara breeds' severe population bottlenecks. Within test population, SNPs used to construct potential candidate genomic areas under selection were only considered if they were identified in four of ten two-by-two pairwise comparisons and in at least two of three statistics. Genomic regions under selection constructed within test population were subsequently overlapped to construct candidate regions under selection putatively unique to the Celtic-Iberian pig strain. These genomic regions were finally used for enrichment analyses. A total of 39 candidate regions, mainly located on SSC5 and SSC9 and covering 3130.5 kb, were identified and could be considered representative of the ancient genomic background of the Celtic-Iberian strain. Enrichment analysis allowed to identify a total of seven candidate genes (NOL12, LGALS1, PDXP, SH3BP1, GGA1, WIF1, and LYPD6). Other studies reported that the WIF1 gene is associated with ear size, one of the characteristic traits of the Celtic-Iberian pig strain. The function of the other candidate genes could be related to reproduction, adaptation and immunity traits, indirectly fitting with the rusticity of a non-improved pig strain traditionally exploited under semi-extensive conditions., (© 2023 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The analysis of the genomic landscape of inbreeding using runs of homozygosity (ROH) patterns is becoming an interesting tool to partially understand phenotypic differences among individuals. In this study, we analysed genome-wide ROH patterns in two groups of Florida goats. We first determined the inbreeding levels of each individual by calculating ROH-based inbreeding coefficients (F ROH ). Then, the individuals were divided into two groups based on F ROH : high inbreeding (HI, F ROH  >0.1) and low inbreeding (LI, F ROH  <0.03). Finally, we performed an extensive in-depth analysis of ROH distribution in each group separately. We found a higher abundance of short ROH in LI, whereas long ROH was more frequent in HI. Furthermore, ROH abundance was not evenly distributed among chromosomes within groups, with some chromosomes showing larger numbers of ROH, like CHI6, CHI7 and CHI27. A different landscape was observed in recent inbreeding (ROH >8 Mb), with significant increases in CHI6, CHI11 and CHI28. Determination of genomic regions with significantly increased ROH (ROH islands-ROHi) showed 13 ROHi related to whole inbreeding and five ROHi associated with recent inbreeding analysis. Within these genomic regions, 123 and 101 genes were identified in HI and LI, respectively, including 10 and seven candidate genes previously related to production, fertility and heat resistance in goats and livestock species., (© 2023 John Wiley & Sons Ltd.)

In Vietnam, local chicken breeds account for over 70% of the national poultry population. Although these breeds are abundant, their productivity is low and their use is threatened by the extensive importation of foreign productive breeds. In this context, conservation programmes targeting several emblematic breeds have been established. The goal of these programmes was to characterize endangered breeds and maintain a pool of characteristic birds for preserving their genetic heritage. To contribute to these programmes, we comprehensively characterized four Vietnamese local chicken breeds (Dong Tao, Ho, Mia and Mong) at the genomic level using high-density single-nucleotide polymorphism (SNP) genotyping. Despite originating in geographically close areas, Dong Tao and Ho were evidently different from each other as well as from Mong and Mia, which shared a more recent common ancestor. The genomic inbreeding coefficient revealed high homozygosity amongst the four breeds (10%-20%). The observation of clear differentiation at the genomic level supported the presence of distinct breeds; nonetheless, the occurrence of crossbred birds in a presumably purebred sample demonstrated the need to apply genomic tools to unambiguously assign the birds to the correct breed. Moreover, the occurrence of substantial inbreeding and the presence of subgroups in certain breeds warranted attention to create future nuclei for use in the conservation of these local breeds., (© 2022 Wiley-VCH GmbH.)

In the past 50 years, there has been a steep increase in the demand for poultry products, met by increasing production along with genetic selection for improved growth, efficiency, health and reproduction. The selection tends to reduce the number and type of genetic resources contributing to the majority of production. The University of Alberta maintains 10 heritage chicken lines (Brown Leghorn (BL), Light Sussex (LS), New Hampshire (NH), Saskatchewan Barred Rock (SaskBR), Shaver Barred Rock (ShaverBR), Shaver Rhode Island Red (RIR), White Leghorn (WL) and three commercial crosses called Alberta Meat Control strains 1957 (AMC-1957), 1978 sire line (AMC-1978-20S) and 1978 dam line (AMC-1978-30D), that played a large role in the evolution of the poultry industry in Canada. Since these lines have not been subjected to the same intensive selection pressures as commercial counterparts, they may contain unique genetic variants lost in commercial lines. Thus, for conservation management of these lines, the first step is to assess their genetic diversity. 71 male samples from across 10 lines were analysed using whole-genome sequencing and patterns of genetic diversity and relatedness among these lines were explored. AMC-1978-30D showed the highest genetic diversity as reflected in observed and expected heterozygosity (0.327 and 0.250), percentage of polymorphic markers (~ 65%) and average recent inbreeding coefficient (-0.039), followed by AMC-1978-20S and AMC-1957. BL showed the lowest genetic diversity as reflected in observed and expected heterozygosity (0.130 and 0.116), percentage of polymorphic markers (~31%) and average recent inbreeding coefficient (0.577), followed by LS, WL and NH. Our findings highlight the need for special attention for the populations of BL, WL, LS and NH, with the largest levels of inbreeding. Our results can be used to develop a breeding strategy to optimize and conserve the genetic variation present in heritage lines., (© 2022 Wiley-VCH GmbH.)

Advances in long-read sequencing technology and genome assembly provide an opportunity to improve the pig genome and reveal the full range of structural variations (SVs) between local Chinese and European pigs. To date, little is known about the genomes of some unique Chinese indigenous breeds, such as the Ningxiang pig. Here, we report the sequencing and assembly of a highly contiguous Ningxiang pig genome (NX) via an integration of PacBio single-molecule real-time sequencing, Illumina next-generation sequencing, BioNano optical mapping and Hi-C (chromosome conformation capture) approaches. The assembled genome comprises 2.44 Gb with a contig N50 of 26.1 Mb and 418 contigs in total. These contigs are organized into 121 scaffolds with a scaffold N50 of 139.0 Mb. More than 99.1% of the assembled sequence could be localized to 19 pseudochromosomes and is annotated with 20,914 protein-coding genes and 34.04% repetitive sequences. Comparisons between the NX and European Duroc assemblies revealed many SVs in genes involved in the immune system, nervous system, lipid metabolism and environmental adaptation. The genetic variants include 47 Chinese domestic pig-specific SVs and the associated 74 genes may contribute to the differences in domestic traits compared to European pigs. Moreover, single nucleotide polymorphisms (SNPs) identified from whole genome resequencing data of 73 Chinese pigs, representing 17 geographically isolated breeds, showed their specific genetic variations, population structure and evolutionary patterns. Finally, we explore transcriptional regulation in the first intron of the MYL4 gene, as the genomic SV (281-bp deletion) in Ningxiang pig promotes its subcutaneous fat compared to European pig breeds. This work identifies a set of Asian-specific SVs and SNPs, which will be important resources for modern pig breeding and genetic conservation., (© 2021 John Wiley & Sons Ltd.)

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species., (© 2021 John Wiley & Sons Ltd.)

There are wide genomic and phenotypic differences between Asian and European pig breeds, yet the current reference genome is the European Duroc pig genome. A high-quality pig genome is lacking for genetic analysis of agricultural traits in Asian pigs. Here, using a hybrid approach, a high-quality reference genome (MSCAAS v1) for the Asian Meishan breed is assembled with a contig N50 size of 48.05 Mb. MSCAAS v1 outperforms the Duroc genome as a reference genome for Asian breeds. Genomic comparison reveals 49,103 structural variations (SVs) between Meishan and Duroc, 4.02% of which are Asian-specific SVs (AP-SVs). Notably, a 30-Mb hotspot for AP-SVs on chromosome X enriched for genes associated with Asian-pig-specific phenotypes is present in Asian domestic pig breeds, but absent in Asian wild boars, suggesting that Asian domestic breeds share a common ancestor. Interbreed transcriptomics reveals transcriptional suppression roles of AP-SVs in multiple tissues. Finally, transcriptional regulation in the intron of IGF2R is reported, as genomic SV (274-bp deletion) in Tibetan pig limits its growth compared to domestic pig breeds. In summary, this study provides insights regarding the genetic changes underlying pig domestication and presents a benchmark-setting resource for the utilization of agricultural valuable loci in Asian pigs., (© 2021 John Wiley & Sons Ltd.)

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep., (© 2024 Wiley‐VCH GmbH. Published by John Wiley & Sons Ltd.)

Biological information regarding markers and gene association may be used to attribute different weights for single nucleotide polymorphism (SNP) in genome-wide selection. Therefore, we aimed to evaluate the predictive ability and the bias of genomic prediction using models that allow SNP weighting in the genomic relationship matrix (G) building, with and without incorporating biological information to obtain the weights. Firstly, we performed a genome-wide association studies (GWAS) in data set containing single- (SL) or a multi-line (ML) pig population for androstenone, skatole and indole levels. Secondly, 1%, 2%, 5%, 10%, 30% and 50% of the markers explaining the highest proportions of the genetic variance for each trait were selected to build gene networks through the association weight matrix (AWM) approach. The number of edges in the network was computed and used to derive weights for G (AWM-WssGBLUP). The single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used as standard scenarios. All scenarios presented predictive abilities different from zero; however, the great overlap in their confidences interval suggests no differences among scenarios. Most of scenarios of based on AWM provide overestimations for skatole in both SL and ML populations. On the other hand, the skatole and indole prediction were no biased in the ssGBLUP (S1) in both SL and ML populations. Most of scenarios based on AWM provide no biased predictions for indole in both SL and ML populations. In summary, using biological information through AWM matrix and gene networks to derive weights for genomic prediction resulted in no increase in predictive ability for boar taint compounds. In addition, this approach increased the number of analyses steps. Thus, we can conclude that ssGBLUP is most appropriate for the analysis of boar taint compounds in comparison with the weighted strategies used in the present work., (© 2020 Wiley-VCH GmbH.)

The Pacific oyster (Crassostrea gigas) is a widely distributed marine bivalve of great ecological and economic importance. In this study, we provide a high-quality chromosome-level genome assembled using Pacific Bioscience long reads and Hi-C-based and linkage-map-based scaffolding technologies and a high-resolution variation map constructed using large-scale resequencing analysis. The 586.8 Mb genome consists of 10 pseudochromosome sequences ranging from 38.6 to 78.9 Mb, containing 301 contigs with an N50 size of 3.1 Mb. A total of 30,078 protein-coding genes were predicted, of which 22,757 (75.7%) were high-reliability annotations supported by a homologous match to a curated protein in the SWISS-PROT database or transcript expression. Although a medium level of repeat components (57.2%) was detected, the genomic content of the segmental duplications reached 26.2%, which is the highest among the reported genomes. By whole genome resequencing analysis of 495 Pacific oysters, a comprehensive variation map was built, comprised of 4.78 million single nucleotide polymorphisms, 0.60 million short insertions and deletions, and 49,333 copy number variation regions. The structural variations can lead to an average interindividual genomic divergence of 0.21, indicating their crucial role in shaping the Pacific oyster genome diversity. The large amount of mosaic distributed repeat elements, small variations, and copy number variations indicate that the Pacific oyster is a diploid organism with an extremely high genomic complexity at the intra- and interindividual level. The genome and variation maps can improve our understanding of oyster genome diversity and enrich the resources for oyster molecular evolution, comparative genomics, and genetic research., (© 2021 John Wiley & Sons Ltd.)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses., (© 2019 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The house sparrow is an important model species for studying physiological, ecological and evolutionary processes in wild populations. Here, we present a medium density, genome wide linkage map for house sparrow (Passer domesticus) that has aided the assembly of the house sparrow reference genome, and that will provide an important resource for ongoing mapping of genes controlling important traits in the ecology and evolution of this species. Using a custom house sparrow 10 K iSelect Illumina SNP chip we have assigned 6,498 SNPs to 29 autosomal linkage groups, based on a mean of 430 informative meioses per SNP. The map was constructed by combining the information from linkage with that of the physical position of SNPs within scaffold sequences in an iterative process. Averaged between the sexes; the linkage map had a total length of 2,004 cM, with a longer map for females (2,240 cM) than males (1,801 cM). Additionally, recombination rates also varied along the chromosomes. Comparison of the linkage map to the reference genomes of zebra finch, collared flycatcher and chicken, showed a chromosome fusion of the two avian chromosomes 8 and 4A in house sparrow. Lastly, information from the linkage map was utilized to conduct analysis of linkage disequilibrium (LD) in eight populations with different effective population sizes (N e ) in order to quantify the background level LD. Together, these results aid the design of future association studies, facilitate the development of new genomic tools and support the body of research that describes the evolution of the avian genome., (© 2020 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (F PED ) information, on ROH identified using the SW method (F ROHSW ) and on ROH identified using the CR method (F ROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between F PED and F ROH were high, particularly for F ROHSW . Moreover, correlations between F ROHSW and F PED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of F ROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with F PED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness., (© 2023 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

Anthropogenic hybridization, or higher and non-natural rates of gene flow directly and indirectly induced by human activities, is considered a significant threat to biodiversity. The primary concern for conservation is the potential for genomic extinction and loss of adaptiveness for native species due to the extensive introgression of non-native genes. To alleviate or reverse trends for such scenarios requires the direct integration of genomic data within a model framework for effective management. Towards this end, we developed the simRestore R program as a decision-making tool that integrates ecological and genomic information to simulate ancestry outcomes from optimized conservation strategies. In short, the program optimizes supplementation and removal strategies across generations until a set native genetic threshold is reached within the studied population. Importantly, in addition to helping with initial decision-making, simulations can be updated with the outcomes of ongoing efforts, allowing for the adaptive management of populations. After demonstrating functionality, we apply and optimize among actionable management strategies for the endangered Hawaiian duck for which the current primary threat is genetic extinction through ongoing anthropogenic hybridization with feral mallards. Simulations demonstrate that supplemental and removal efforts can be strategically tailored to move the genetic ancestry of Hawaii's hybrid populations towards Hawaiian duck without the need to completely start over. Further, we discuss ecological parameter sensitivity, including which factors are most important to ensure genetic outcomes (i.e. number of offspring). Finally, to facilitate use, the program is also available online as a Shiny Web application., (© 2023 John Wiley & Sons Ltd.)

The field of molecular biology is advancing fast with new powerful technologies, sequencing methods and analysis software being developed constantly. Commonly used tools originally developed for research on humans and model species are now regularly used in ecological and evolutionary research. There is also a growing interest in the causes and consequences of epigenetic variation in natural populations. Studying ecological epigenetics is currently challenging, especially for vertebrate systems, because of the required technical expertise, complications with analyses and interpretation, and limitations in acquiring sufficiently high sample sizes. Importantly, neglecting the limitations of the experimental setup, technology and analyses may affect the reliability and reproducibility, and the extent to which unbiased conclusions can be drawn from these studies. Here, we provide a practical guide for researchers aiming to study DNA methylation variation in wild vertebrates. We review the technical aspects of epigenetic research, concentrating on DNA methylation using bisulfite sequencing, discuss the limitations and possible pitfalls, and how to overcome them through rigid and reproducible data analysis. This review provides a solid foundation for the proper design of epigenetic studies, a clear roadmap on the best practices for correct data analysis and a realistic view on the limitations for studying ecological epigenetics in vertebrates. This review will help researchers studying the ecological and evolutionary implications of epigenetic variation in wild populations., (© 2022 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

A high-quality reference genome can be a valuable resource for threatened species by providing a foundation to assess their evolutionary potential to adapt to future pressures such as environmental change. We assembled the genome of a female hihi (Notiomysits cincta), a threatened passerine bird endemic to Aotearoa New Zealand. The assembled genome is 1.06 Gb, and is of high quality and highly contiguous, with a contig N50 of 7.0 Mb, estimated QV of 44 and a BUSCO completeness of 96.8%. A male assembly of comparable quality was generated in parallel. A population linkage map was used to scaffold the autosomal contigs into chromosomes. Female and male sequence coverage and comparative genomics analyses were used to identify Z-, and W-linked contigs. In total, 94.6% of the assembly length was assigned to putative nuclear chromosome scaffolds. Native DNA methylation was highly correlated between sexes, with the W chromosome contigs more highly methylated than autosomal chromosomes and Z contigs. 43 differentially methylated regions were identified, and these may represent interesting candidates for the establishment or maintenance of sex differences. By generating a high-quality reference assembly of the heterogametic sex, we have created a resource that enables characterization of genome-wide diversity and facilitates the investigation of female-specific evolutionary processes. The reference genomes will form the basis for fine-scale assessment of the impacts of low genetic diversity and inbreeding on the adaptive potential of the species and will therefore enable tailored and informed conservation management of this threatened taonga (treasured) species., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Swedish Red (SR) and Swedish Holstein (SH) are the dominating commercial dairy cattle breeds in Sweden. Both breeds have undergone substantial changes during the last half century due to intensive selection for breeding goal traits, but also resulting from increased international exchange of breeding animals and genetic drift. The aim of this study was to learn more about changes in genomic diversity and inbreeding in these two breeds over time. Therefore, semen samples from old bulls were genotyped using the 150K Genomic Profiler SNP array and combined with 50K SNP array genotype data, obtained for more recent bulls from the Nordic Cattle Genetic Evaluation. Different measures of level of homozygosity, genomic inbreeding, relatedness and changes in allele frequency were estimated for bulls born during different time periods from the 1950s until 2020. In total, more than 33,000 SNPs for 9737 SR and 5041 SH bulls were included in the analysis using PLINK v1.9. The results showed higher average homozygosity for SR than for SH bulls up to around 2000, but the difference was very small after that. The average inbreeding coefficients based on deviation from expected homozygosity as well as on runs of homozygosity decreased until the early 1980s in both breeds, whereafter they started to increase again for SH, but stayed more stable for SR. From the 1990s onwards, SH displayed higher average inbreeding coefficients than SR. In the last studied birth year group (2015-2020), the mean inbreeding coefficient based on runs of homozygosity was 5.9% for SH and 3.7% for SR. A principal component analysis showed a pattern of genetic relationships related to the birth year period of the bulls, illustrating the gradual change of the genetic material within each breed. The change in allele frequency over time was generally larger for SH than for SR. The results show that the inbreeding level was higher half a century ago than at present, and the inbreeding levels were lower than in some other studied populations. Still, the increase seen for inbreeding coefficients and homozygosity, especially in SH during recent years, should be considered in future breeding strategies., (© 2023 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

Disentangling the interaction between the genetic basis and environmental context underlying phenotypic variation is critical for understanding organismal evolution. Environmental change, such as increased rates of urbanization, can induce shifts in phenotypic plasticity with some individuals adapting to city life while others are displaced. A key trait that can facilitate adaptation is the degree at which animals respond to stressors. This stress response, which includes elevation of baseline circulating concentrations of glucocorticoids, has a heritable component and exhibits intra- and inter-individual variation. However, the mechanisms behind this variability and whether they might be responsible for adaptation to different environments are not known. Variation in DNA methylation can be a potential mechanism that mediates environmental effects on the stress response, as early-life stressors increase glucocorticoid concentrations and change adult phenotype. We used an inter- and intra-environmental cross-foster experiment to analyse the contribution of DNA methylation to early-life phenotypic variation. We found that at hatching, urban house wren (Troglodytes aedon) offspring had higher methylation frequencies compared with their rural counterparts. We also observed age-related patterns in offspring methylation, indicating the developmental effects of the rearing environment on methylation. At fledgling, differential methylation analyses showed that cellular respiration genes were differentially methylated in broods of different origins and behavioural and metabolism genes were differentially methylated in broods of different rearing environments. Lastly, hyper-methylation of a single gene (CNTNAP2) is associated with decreased glucocorticoid levels and the rearing environment. These differential methylation patterns linked to a specific physiological phenotype suggest that DNA methylation may be a mechanism by which individuals adjust to novel environments during their lifespan. Characterizing genetic and environmental influences on methylation is critical for understanding the role of epigenetic mechanisms in evolutionary adaptation., (© 2022 The Authors. Journal of Animal Ecology © 2022 British Ecological Society.)

This study aimed at investigating the influence of genetic and non-genetic factors on immune traits to inform on possibilities of genetic improvement of disease resistance traits in local chicken of Kenya. Immune traits such as natural and specific antibodies are considered suitable indicators of an individual's health status and consequently, used as indicator traits of disease resistance. In this study, natural antibodies binding to Keyhole Limpet Hemocyanin (KLH-NAbs) was used to measure general disease resistance. Specific antibodies binding to Newcastle disease virus (NDV-IgG) post vaccination was used to measure specific disease resistance. Titers of KLH-NAbs isotypes (KLH-IgM, KLH-IgG and KLH-IgA) and NDV-IgG were measured in 1,540 chickens of different ages ranging from 12 to 56 weeks. A general linear model was fitted to determine the effect of sex, generation, population type, phylogenetic cluster, line, genotype and age on the antibody traits. A multivariate animal mixed model was fitted to estimate heritability and genetic correlations among the antibody traits. The model constituted of non-genetic factors found to have a significant influence on the antibody traits as fixed effects, and animal and residual effects as random variables. Overall mean (±SE) concentration levels for KLH-IgM, KLH-IgG, KLH-IgA and NDV-IgG were 10.33 ± 0.04, 9.08 ± 0.02, 6.00 ± 0.02 and 10.12 ± 0.03, respectively. Sex, generation and age (linear covariate) significantly (p < 0.05) influenced variation across all the antibody traits. Genotype effects (p < 0.05) were present in all antibody traits, apart from KLH-IgA. Interaction between generation and line was significant (p < 0.05) in KLH-IgM and NDV-IgG while nesting phylogenetic cluster within population significantly (p < 0.05) influenced all antibody traits, apart from KLH-IgA. Heritability estimates for KLH-IgM, KLH-IgG, KLH-IgA and NDV-IgG were 0.28 ± 0.08, 0.14 ± 0.06, 0.07 ± 0.04 and 0.31 ± 0.06, respectively. There were positive genetic correlations (0.40-0.61) among the KLH-NAbs while negative genetic correlations (-0.26 to -0.98) were observed between the KLH-NAbs and NDV-IgG. Results from this study indicate that non-genetic effects due to biological and environmental factors influence natural and specific antibodies and should be accounted for to reduce bias and improve accuracy when evaluating the traits. Subsequently, the moderate heritability estimates in KLH-IgM and NDV-IgG suggest selection possibilities for genetic improvement of general and specific immunity, respectively, and consequently disease resistance. However, the negative correlations between KLH-NAbs and NDV-IgG indicate the need to consider a suitable approach that can optimally combine both traits in a multiple trait selection strategies., (© 2022 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

This study aimed to integrate GWAS and structural variants to propose possible molecular biomarkers related to gastrointestinal nematode resistance traits in Santa Inês sheep. The phenotypic records FAMACHA, haematocrit, white blood cell count, red blood cell count, haemoglobin, platelets and egg counts per gram of faeces were collected from 700 naturally infected animals, belonging to four Brazilian flocks. A total of 576 animals were genotyped using the Ovine SNP12k BeadChip and were imputed using a reference population with Ovine SNP50 BeadChip. The GWAS approaches were based on SNPs, haplotypes, CNVs and ROH. The overlapping between the significant genomic regions detected from all approaches was investigated, and the results were integrated using a network analysis. Genes related to the immune system were found, such as ABCB1, IL6, WNT5A and IRF5. Genomic regions containing candidate genes and metabolic pathways involved in immune responses, inflammatory processes and immune cells affecting parasite resistance traits were identified. The genomic regions, biological processes and candidate genes uncovered could lead to biomarkers for selecting more resilient sheep and improving herd welfare and productivity. The results obtained are the start point to identify molecular biomarkers related to indicator traits of gastrointestinal nematode resistance in Santa Inês sheep., (© 2022 Wiley-VCH GmbH.)

Sample mix-ups occur when samples have accidentally been duplicated, mislabelled or swapped. When samples are subsequently genotyped or sequenced, this can lead to individual IDs being incorrectly linked to genetic data, resulting in incorrect or biased research results, or reduced power to detect true biological patterns. We surveyed the community and found that almost 80% of responding researchers have encountered sample mix-ups. However, many recent studies in the field of molecular ecology do not appear to systematically report individual assignment checks as part of their publications. Although checks may be done, lack of consistent reporting means that it is difficult to assess whether sample mix-ups have occurred or been detected. Here, we present an easy-to-follow sample verification framework that can utilise existing metadata, including species, population structure, sex and pedigree information. We demonstrate its application to a data set representing individuals of a threatened Aotearoa New Zealand bird species, the hihi, genotyped on a 50K SNP array. We detected numerous incorrect genotype-ID associations when comparing observed and genetic sex or comparing to relationships in a verified microsatellite pedigree. The framework proposed here helped to confirm 488 individuals (39%), correct another 20 bird-genotype links, and detect hundreds of incorrect sample IDs, emphasizing the value of routinely checking genetic and genomic data sets for their accuracy. We therefore promote the implementation and reporting of this simple yet effective sample verification framework as a standardized quality control step for studies in the field of molecular ecology., (© 2021 John Wiley & Sons Ltd.)

Copy number variations (CNV) are an important source of genetic variation. CNV has been increasingly studied and frequently associated with diseases and productive traits in livestock animals. However, CNV-based genome-wide association studies (GWAS) in Santa Inês sheep, one of the principal sheep breeds in Brazil, have not yet been reported. Thus, the aim of this study was to investigate the association between CNV and growth, efficiency and carcass traits in sheep. The Illumina OvineSNP50 BeadChip array was used to detect CNV in 491 Santa Inês individuals. Then, CNV-based GWAS was performed with a linear mixed model approach considering a genomic relationship matrix, for ten traits: (1) growth: body weight at three (W3) and six (W6) months of age; (2) efficiency: residual feed intake (RFI) and feed efficiency (FE) and (3) carcass: external carcass length (ECL), leg length (LL), carcass yield (CY), commercial cuts weight (CCW), loin eye area (LEA) and subcutaneous fat thickness (SFT). We identified 1,167 autosomal CNV in 438 sheep, with 294 non-redundant CNV, ranging from 21.8 to 861.9 kb, merged into 216 distinct copy number variation regions (CNVRs). One significant CNV segment (p FDR -value<0.05) in OAR3 was associated with CY, while another significant CNV in OAR6 was associated with RFI. Additionally, another 5 CNV segments were considered relevant for investigation in the future studies. The significant segments overlapped 4 QTLs and spanned 8 genes, including the SPAST,TGFA and ADGRL3 genes, involved in cell differentiation and energy metabolism. Therefore, the results of the present study increase knowledge about CNV in sheep, their possible impacts on productive traits, and provide information for future investigations, being especially useful for those interested in structural variations in the sheep genome., (© 2022 Wiley-VCH GmbH.)

Quantifying the level of linkage disequilibrium (LD), non-random association of alleles at two or more loci, is important to determine the number of markers needed for genomic selection. The aims of this study were to evaluate the extent of LD in Dual-Purpose Belgian Blue (DPBB) and to compare the level of LD in DPBB with that of Walloon Holstein. Data of 28,427 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), of 639 DPBB and 398 Holstein bulls were used. The level of LD between pairwise SNPs separated by up to 10 Mb was evaluated, separately for each breed, using the squared correlation of the alleles at two loci. The analysis of molecular variance showed that the percentage of variation within populations (85.48%) was higher than between populations (14.52%). However, permutation tests showed a significant genetic differentiation between the two studied populations (p < .01). The average LD found between adjacent SNP pairs in DPBB (0.16 (SD = 0.22)) was generally lower than in Holstein (0.23 (SD = 0.27)). The proportion of SNPs in useful LD (r 2  > 0.30) within a genomic distance of ≤0.10 Mb between SNPs was 18.58% and 28.23% in DPBB and Holstein bulls, respectively. In both breeds, the effective population size decreased over generations; however, the decline was greater in DPBB than that in Holstein. Based on results, it can be concluded that at least 68,000 SNPs are needed for implementing genomic selection in DPBB cattle with enough accuracy., (© 2021 Wiley-VCH GmbH.)

With the rapid growth of the number of sequenced ancient genomes, there has been increasing interest in using this new information to study past and present adaptation. Such an additional temporal component has the promise of providing improved power for the estimation of natural selection. Over the last decade, statistical approaches for the detection and quantification of natural selection from ancient DNA (aDNA) data have been developed. However, most of the existing methods do not allow us to estimate the timing of natural selection along with its strength, which is key to understanding the evolution and persistence of organismal diversity. Additionally, most methods ignore the fact that natural populations are almost always structured, which can result in an overestimation of the effect of natural selection. To address these issues, we introduce a novel Bayesian framework for the inference of natural selection and gene migration from aDNA data with Markov chain Monte Carlo techniques, co-estimating both timing and strength of natural selection and gene migration. Such an advance enables us to infer drivers of natural selection and gene migration by correlating genetic evolution with potential causes such as the changes in the ecological context in which an organism has evolved. The performance of our procedure is evaluated through extensive simulations, with its utility shown with an application to ancient chicken samples., (© 2021 John Wiley & Sons Ltd.)

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, F ST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. F ST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using F ST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both F ST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation., (© 2021 Wiley-VCH GmbH.)

Qingyuan partridge chicken is an important indigenous chicken in China. In its breeding schemes, chickens are usually selected at the age of 105-day-old for five traits, including body weight (BW), shank length (SL), shank girth (SG), comb height (CH) and feather maturity (FM). At present, genetic parameters of the aforementioned traits have still not been studied in Qingyuan partridge chickens. The objectives of this study were the following: (1) to investigate whether the optimal statistical models of these traits need to consider maternal genetic and permanent environmental effects in late-feathering Qingyuan partridge hens, and (2) to estimate genetic parameters for these traits based on the optimal models. The numbers of records for BW, SL, SG, CH and FM were 13,721, 13,671, 13,670, 13,669 and 13,672, respectively. Variance components were estimated using average information-restricted maximum likelihood method, and the optimal model was determined based on Bayesian information criterion. More specifically, the optimal model for BW considered maternal genetic and permanent environmental effects in addition to direct additive genetic effect; SL, SG and FM considered direct and maternal genetic effects; and CH considered direct and maternal genetic effects, and the covariance between them. The direct heritabilities of these traits estimated using the optimal models were 0.21 ± 0.04, 0.30 ± 0.05, 0.40 ± 0.05, 0.59 ± 0.09 and 0.09 ± 0.04, respectively; the maternal heritabilities were 0.01 ± 0.04, 0.05 ± 0.05, 0.04 ± 0.05, 0.09 ± 0.09 and 0.03 ± 0.04, respectively. Maternal genetic effect evidently played an important part in FM and maternal heritability accounted for 30 per cent of total heritability. Furthermore, the direct and maternal genetic effects for CH were estimated to be negatively and moderately correlated (-0.51 ± 0.11). For all traits, neglecting existent maternal effects biased the estimation of direct heritability. Therefore, to implement optimum breeding strategies for improvement of these traits in Qingyuan partridge hens, maternal effects should be taken into consideration., (© 2020 Wiley-VCH GmbH.)

Mud crabs, found throughout the Indo-Pacific region, are coastal species that are important fisheries resources in many tropical and subtropical Asian countries. Here, we present a chromosome-level genome assembly of a mud crab (Scylla paramamosain). The genome is 1.55 Gb (contig N50 191 kb) in length and encodes 17,821 proteins. The heterozygosity of the assembled genome was estimated to be 0.47%. Effective population size analysis suggested that an initial large population size of this species was maintained until 200 thousand years ago. The contraction of cuticle protein and opsin genes compared with Litopenaeus vannamei is assumed to be correlated with shell hardness and light perception ability, respectively. Furthermore, the analysis of three chemoreceptor gene families, the odorant receptor (OR), gustatory receptor (GR) and ionotropic receptor (IR) families, suggested that the mud crab has no OR genes and shows a contraction of GR genes and expansion of IR genes. The numbers of the three gene families were similar to those in three other decapods but different from those in two nondecapods and insects. In addition, IRs were more diversified in decapods than in nondecapod crustaceans, and most of the expanded IRs in the mud crab genome were clustered with the antennal IR clades. These findings suggested that IRs might exhibit more diverse functions in decapods than in nondecapods, which may compensate for the smaller number of GR genes. Decoding the S. paramamosain genome not only provides insight into the genetic changes underpinning ecological traits but also provides valuable information for improving the breeding and aquaculture of this species., (© 2021 John Wiley & Sons Ltd.)

The common pheasant Phasianus colchicus, belonging to the order Galliformes and family Phasianidae, is the most widespread species. Despite a long history of captivity, the domestication of this bird is still at a preliminary stage. Recently, the demand for accelerating its transformation to poultry for meat and egg production has been increasing. In this study, we assembled high quality, chromosome scale genome of the common pheasant by using PacBio long reads, next-generation short reads, and Hi-C technology. The primary assembly has contig N50 size of 1.33 Mb and scaffold N50 size of 59.46 Mb, with a total size of 0.99 Gb, resolving most macrochromosomes into single scaffolds. A total of 23,058 genes and 10.71 Mb interspersed repeats were identified, constituting 30.31% and 10.71% of the common pheasant genome, respectively. Our phylogenetic analysis revealed that the common pheasant shared common ancestors with turkey about 24.7-34.5 million years ago (Ma). Rapidly evolved gene families, as well as branch-specific positively selected genes, indicate that calcium-related genes are potentially related to the adaptive and evolutionary change of the common pheasant. Interestingly, we found that the common pheasant has a unique major histocompatibility complex B locus (MHC-B) structure: three major inversions occurred in the sequence compared with chicken MHC-B. Furthermore, we detected signals of selection in five breeds of domestic common pheasant, several of which are production-oriented., (© 2020 John Wiley & Sons Ltd.)

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes., (© 2020 John Wiley & Sons Ltd.)

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (F ST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock., (© 2020 Wiley-VCH GmbH.)

Runs of homozygosity (ROH) are continuous segments of the genome that arose as a result of inbreeding, resulting in the inheritance of identical haplotypes from both parents who shared a common ancestor. In the present study, we performed a detailed characterization and comparison of ROH in four pig breeds, including intensively selected Polish Landrace as well as native unselected animals of Puławska and two Złotnicka breeds (White and Spotted). We used a medium-density PorcineSNP60 BeadChip assay (Illumina) and cgaTOH software to detect ROH covering a minimum of 30 adjacent SNPs and maintaining a size over 1 Mb. By analysing ROH distribution and frequency across the genome, we also identified genomic regions with high ROH frequency (so-called "ROH hotspots"). The obtained results showed that the analysed conserved breeds were characterized by a higher ROH span and higher ROH-based inbreeding coefficients (F ROH ), which likely result from past population bottlenecks, increasing the overall inbreeding level within these populations. The analysis of ROH distribution across the genomes revealed the presence of both shared and breed-specific ROH hotspots. These hotspots, presumably representing genome regions under selection, overlapped with a variety of genes associated with processes connected with immune system functioning, reproduction, glucose homeostasis and metabolism. The genome regions with ROH hotspots overlapping in all analysed populations, located on SSC4 (51.9-55.9 Mb) and 13 (92.6-97.8 Mb), covered thirty-one different genes, including MMP16, SLC7A13, ATP6V0D2, CNGB3, WWiP1, RiMDN1 and CPNE3. These genes are primarily associated with biological regulation and metabolism, processes that could be responsible for the variety of the selected production and functional features., (© 2020 Wiley-VCH GmbH.)

Equine insect bite hypersensitivity (IBH) is a pruritic skin allergy caused primarily by biting midges, Culicoides spp. IBH susceptibility has polygenic inheritance and occurs at high frequencies in several horse breeds worldwide, causing increased costs and reduced welfare of affected horses. The aim of this study was to identify and validate single nucleotide polymorphisms (SNPs) associated with equine IBH susceptibility. After quality control, 33,523 SNPs were included in a Bayesian genome-wide association study on 177 affected and 178 unaffected Icelandic horses. We report associated regions in E. caballus (ECA) 1, 3, 15 and 18, overlapping with known IBH QTLs in horses, and novel regions containing several genes, together explaining 11.46% of the total genetic variance. For validation, three SNPs on ECA 1 and ECA X (explaining the largest percentage of genetic variance) within 1-mb genomic windows for IBH were genotyped in an independent population of 280 Exmoor ponies. The associated genomic region (152-153 mb) on ECA 1 was confirmed in Exmoor ponies and contains the AQR gene involved in splicing processes and a long non-coding RNA. This study confirms the polygenic nature of IBH susceptibility and suggests a role of transcriptional regulatory mechanisms (e.g., alternative splicing) for IBH predisposition in these horse breeds., (© 2019 Blackwell Verlag GmbH.)

Most traits in animal breeding, including feed efficiency traits in pigs, are affected by many genes with small effect and have a moderately high heritability between 0.1 and 0.5, which enables efficient selection. Since the microbiota composition in the gastrointestinal tract is also partly heritable and was shown to have a substantial effect on feed efficiency, the host genes affect the phenotype not only directly by altering metabolic pathways, but also indirectly by changing the microbiota composition. The effect m i of the microbiota composition on the breeding value g i of an animal i is the conditional expectation of its breeding value, given the vector φ i with microbiota frequencies, that is m i = E g i | φ i . The breeding value g i of an animal can therefore be decomposed into a heritable contribution m i that arises from an altered microbiota composition and a heritable contribution p i that arises from altered metabolic pathways within the animal, so g i = m i + p i . Instead of selecting for breeding value g ^ i , an index comprising the two components m ^ i and p ^ i with appropriate weights, that is I i = λ 1 m ^ i + λ 2 p ^ i , can be used. The present study shows how this breeding strategy can be applied in pig genomic selection breeding scheme for two feed efficiency traits and daily gain., (© 2019 Blackwell Verlag GmbH.)

Recovering the native genetic background of a breed and increasing the founder genome equivalent (FGE) that is contributed by a breed to the gene pool of the species can increase its value for conservation. The suitability of several strategies was compared, whereby a hypothetical multi-breed population, the core set, was used to approximate the genetic diversity of the species. Twenty-five generations of management were simulated based on genotypes of German Angler cattle. The scenarios were compared when the kinship reached 0.10. The native contribution (NC) increased in a population with 400 births per generation from 0.317 to 0.706, whereas 1,000 births enabled to reach 0.894. This scenario maximized the NC, constrained the native kinship, and the kinship of the core set so that its genetic diversity could not decrease. It increased the proportions of mainstream breeds because their genes were removed from the target breed. A substantial increase of the FGE was achieved in some other scenarios, which arose from reduced genetic overlap and from increased diversity within the breed. The latter factor is especially important for breeds with high contributions to the core set., (© 2019 Blackwell Verlag GmbH.)

The plasma very low-density lipoprotein (VLDL) concentration is an effective blood biochemical indicator that could be used to select lean chicken lines. In the current study, we used Genome-wide association study (GWAS) method to detect SNPs with significant effects on plasma VLDL concentration. As a result, 38 SNPs significantly associated with plasma VLDL concentration were identified using at least one of the three mixed linear model (MLM) packages, including GRAMMAR, EMMAX and GEMMA. Nearly, all these SNPs with significant effects on plasma VLDL concentration (except Gga&#95;rs16160897) have significantly different allele frequencies between lean and fat lines. The 1-Mb regions surrounding these 38 SNPs were extracted, and twelve important regions were obtained after combining the overlaps. A total of 122 genes in these twelve important regions were detected. Among these genes, LRRK2, ABCD2, TLR4, E2F1, SUGP1, NCAN, KLF2 and RAB8A were identified as important genes for plasma VLDL concentration based on their basic functions. The results of this study may supply useful information to select lean chicken lines., (© 2019 Blackwell Verlag GmbH.)